Topology optimization of fluid-to-fluid heat exchangers
The optimization of heat exchangers is a recent challenge, drawing more and more attention from industrial designers. These are multiphysics devices used to cool down or to heat up fluids by conveying them in the vicinity of another refrigerating or heating gas, liquid, or solid. The implementation of the constraint that both fluid must not mix during the optimization process is known to be quite challenging.
On this page, we show a few pictures of 2D and 3D test cases obtained with the methodology presented in the publication
Abstract:
We present a topology optimization approach for the design of
fluid-to-fluid heat exchangers which rests on an explicit meshed
discretization of the phases at stake, at every iteration of the optimization
process. The considered physical situations involve a weak coupling between the
Navier--Stokes equations for the velocity and the pressure in the fluid, and the
convection--diffusion equation for the temperature field. The proposed framework
combines several recent techniques from the field of shape and topology
optimization, and notably a level-set based mesh evolution algorithm for tracking
shapes and their deformations, an efficient optimization algorithm for
constrained shape optimization problems, and a numerical method to handle a wide
variety of geometric constraints such as thickness constraints and non-penetration
constraints. Our strategy is applied to the optimization of various types
of heat exchangers. At first, we consider a simplified 2D cross-flow model where
the optimized boundary is the section of the hot fluid phase flowing in the
transverse direction, which is naturally composed of multiple holes. A minimum
thickness constraint is imposed on the cross-section so as to account for
manufacturing and maximum pressure drop constraints. In a second part, we
optimize the design of 2D and 3D heat exchangers composed of two types of fluid
channels (hot and cold), which are separated by a solid body. A non-mixing
constraint between the fluid components containing the hot and cold phases is
enforced by prescribing a minimum distance between them. Numerical results are
presented on a variety of test cases, demonstrating the efficiency of our
approach in generating new, realistic, and unconventional heat exchanger designs.
@article{FEPPON2021113638,
title = "Body-fitted topology optimization of 2D and 3D fluid-to-fluid heat exchangers",
journal = "Computer Methods in Applied Mechanics and Engineering",
volume = "376",
pages = "113638",
year = "2021",
issn = "0045-7825",
doi = "https://doi.org/10.1016/j.cma.2020.113638",
url = "http://www.sciencedirect.com/science/article/pii/S0045782520308239",
author = "F. Feppon and G. Allaire and C. Dapogny and P. Jolivet",
}
A more pedagogical introduction to the method is also proposed in the following lecture notes:
[CHX] Feppon, F. Shape and topology optimization applied to Compact Heat Exchangers (2021). Submitted. HAL preprint hal-03207863.
(abstract)(bibtex)
Abstract:
The purpose of these notes is to offer a comprehensive introduction to topology
optimization for automated generation of complex heat exchanger designs, based on
the methodof Hadamard whereby the design variable is the shape of the fluid-solid
interfaces and is updated iteratively until convergence to a nearly optimal
design. The material is intended to be an introductory
exposure to our recent work (Feppon et al.(2021)) and PhD
thesis (Feppon, 2019).
@unpublished{feppon:hal-03207863,
TITLE = {{Shape and topology optimization applied to Compact Heat Exchangers}},
AUTHOR = {Feppon, Florian},
URL = {https://hal.archives-ouvertes.fr/hal-03207863},
NOTE = {working paper or preprint},
YEAR = {2021},
MONTH = Apr,
PDF = {https://hal.archives-ouvertes.fr/hal-03207863/file/vki_template.pdf},
HAL_ID = {hal-03207863},
HAL_VERSION = {v1},
}
The goal of the optimal design problem is to find the shape of two channels conveying respectively hot and cold fluid phase that maximize the heat exchanged between both components. Furthermore, we seek to limitate the static pressure drop constraint on each of the two channel, and we impose that both phases do not interpenetrate. The latter constraint is enforced conveniently in the framework of Hadamard by prescribing a minimum distance between the two fluid subdomains.
The topology optimization is achieved with a level-set based mesh evolution method, see this page and the following reference for more information.
Abstract:
An efficient framework is described
for the shape and topology optimization of realistic three-dimensional,
weakly-coupled fluid-thermal-mechanical systems. At the theoretical
level, the proposed methodology relies on the boundary variation of
Hadamard for describing the sensitivity of functions with respect to the
domain. From the numerical point of view, three key ingredients are
used:
(i) a level set based mesh evolution method allowing to describe large
deformations of the shape while maintaining an adapted, high-quality mesh of
the latter at every stage of the optimization process;
(ii) an efficient constrained optimization algorithm which is very well
adapted to the infinite-dimensional shape optimization context;
(iii) efficient preconditioning techniques for the solution of large finite
element systems in a reasonable computational time.
The performance of our strategy is illustrated with two examples of coupled
physics: respectively fluid--structure interaction and convective heat
transfer. Before that, we perform three other test cases, involving a
single physics (structural, thermal and aerodynamic design), for comparison
purposes and for assessing our various tools: in particular, they prove the
ability of the mesh evolution technique to capture very thin bodies or
shells in 3D.
@article{FEPPON2020109574,
title = "Topology optimization of thermal fluid–structure systems using body-fitted meshes and parallel computing",
journal = "Journal of Computational Physics",
volume = "417",
pages = "109574",
year = "2020",
issn = "0021-9991",
doi = "https://doi.org/10.1016/j.jcp.2020.109574",
url = "http://www.sciencedirect.com/science/article/pii/S002199912030348X",
author = "F. Feppon and G. Allaire and C. Dapogny and P. Jolivet",
keywords = "Shape and topology optimization, Fluid–structure interaction, Convective heat transfer, Aerodynamic design, Mesh adaptation, Distributed computing",
}
Some 2D optimization histories of heat exchanger design optimization with a non-mixing constraint
2D topology optimization of two non-mixing fluid channels for a Navier-Stokes flow (co-current heat exchange, fluid in white color)
2D topology optimization of two non-mixing fluid channels for a Navier-Stokes flow in a different setting (counter-current heat exchange, fluid in white color)
Some 3D optimization histories of heat exchanger design optimization with a non-mixing constraint
3D topology optimization of two non-mixing fluid channels (stokes flow) with very small pressure drop allowed
3D topology optimization of two non-mixing fluid channels (stokes flow) with larger pressure drops allowed (cut of the cubic domain)
3D topology optimization of two non-mixing fluid channels for a Navier-Stokes flow.
3D topology optimization of two non-mixing fluid channels for a Navier-Stokes flow with a smaller prescribed solid wall thickness.