Topology optimization of flow domains using the lattice Boltzmann method

We consider the optimal design of two- (2D) and three-dimensional (3D) flow domains using the lattice Boltzmann method (LBM) as an approximation of Navier-Stokes (NS) flows. The problem is solved by a topology optimization approach varying the effective porosity of a fictitious material. The boundaries of the flow domain are represented by potentially discontinuous material distributions. NS flows are traditionally approximated by finite element and finite volume methods. These schemes, while well established as high-fidelity simulation tools using body-fitted meshes, are effected in their accuracy and robustness when regular meshes with zero-velocity constraints along the surface and in the interior of obstacles are used, as is common in topology optimization. Therefore, we study the potential of the LBM for approximating low Mach number incompressible viscous flows for topology optimization. In the LBM the geometry of flow domains is defined in a discontinuous manner, similar to the approach used in material-based topology optimization. In addition, this non-traditional discretization method features parallel scalability and allows for high-resolution, regular fluid meshes. In this paper, we show how the variation of the porosity can be used in conjunction with the LBM for the optimal design of fluid domains, making the LBM an interesting alternative to NS solvers for topology optimization problems. The potential of our topology optimization approach will be illustrated by 2D and 3D numerical examples.