Topology optimization for shallow-flow channel design using depth-integrated model

A topology optimization method is proposed for the design of shallow-flow channels based on quasi-three-dimensional flow models of laminar and turbulent flows. The models for laminar flow and turbulent flow are derived from the Navier– Stokes equations and the Reynolds-Averaged Navier–Stokes (RANS) equations, respectively, by integrating along the direction of channel thickness. The thickness is employed as the design variable in the topology optimization. The design variables are updated using a time-dependent diffusion equation with a design sensitivity which is calculated by a discrete adjoint approach. Numerical examples for minimizing dissipation energy or variance of flow velocity magnitude using the topology optimization demonstrates that the proposed method is capable of finding optimal solutions that satisfy the KKT conditions. In the former example, the design domain was clearly divided into domains where the thickness was either near the upper limit or near the lower limit. However, in the latter example, the thickness was at an intermediate level in almost the whole the design domain. The distribution of the thickness varied depending on the Reynolds number in both examples.

[1]  Yafei Jia,et al.  Numerical Model for Channel Flow and Morphological Change Studies , 1999 .

[2]  C. Othmer A continuous adjoint formulation for the computation of topological and surface sensitivities of ducted flows , 2008 .

[3]  A. Gersborg-Hansen Topology optimization of 3D Stokes flow problems , 2006 .

[4]  Atsushi Kawamoto,et al.  Drag minimization and lift maximization in laminar flows via topology optimization employing simple objective function expressions based on body force integration , 2012 .

[5]  Charles G. Speziale,et al.  ANALYTICAL METHODS FOR THE DEVELOPMENT OF REYNOLDS-STRESS CLOSURES IN TURBULENCE , 1990 .

[6]  Ken-ichi Abe,et al.  A new turbulence model for predicting fluid flow and heat transfer in separating and reattaching flows—I. Flow field calculations , 1995 .

[7]  Joel H. Ferziger,et al.  Zonal modeling of turbulent flows - Philosophy and accomplishments , 1990 .

[8]  L. H. Olesen,et al.  A high‐level programming‐language implementation of topology optimization applied to steady‐state Navier–Stokes flow , 2004, physics/0410086.

[9]  J. J. Leendertse,et al.  Aspects of a computational model for long-period water-wave propagation , 1967 .

[10]  İ. Tosun,et al.  Critical Reynolds number for Newtonian flow in rectangular ducts , 1988 .

[11]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[12]  C. H. Marchi,et al.  The lid-driven square cavity flow: numerical solution with a 1024 x 1024 grid , 2009 .

[13]  J. Petersson,et al.  Topology optimization of fluids in Stokes flow , 2003 .

[14]  Mitsuru Kitamura,et al.  Porous composite with negative thermal expansion obtained by photopolymer additive manufacturing , 2015, 1504.07724.

[15]  Keh-Chia Yeh,et al.  Bend-Flow Simulation Using 2D Depth-Averaged Model , 1999 .

[16]  Ping Zhang,et al.  Topology optimization of unsteady incompressible Navier-Stokes flows , 2011, J. Comput. Phys..

[17]  D. I. Papadimitriou,et al.  CONSTRAINED TOPOLOGY OPTIMIZATION FOR LAMINAR AND TURBULENT FLOWS , INCLUDING HEAT TRANSFER , 2011 .

[18]  Shintaro Yamasaki,et al.  Topology optimization by a time‐dependent diffusion equation , 2013 .

[19]  W. Rodi,et al.  Predictions of Heat and Mass Transfer in Open Channels , 1978 .

[20]  O. Sigmund,et al.  Topology optimization of channel flow problems , 2005 .

[21]  R. B. Dean Reynolds Number Dependence of Skin Friction and Other Bulk Flow Variables in Two-Dimensional Rectangular Duct Flow , 1978 .

[22]  M. Bendsøe,et al.  Generating optimal topologies in structural design using a homogenization method , 1988 .

[23]  Xianbao Duan,et al.  Shape-topology optimization for Navier-Stokes problem using variational level set method , 2008 .

[24]  N. Kasagi,et al.  Turbulence measurement in a separated and reattaching flow over a backward-facing step with the aid of three-dimensional particle tracking velocimetry , 1993 .