TY - JOUR T1 - Solution Remapping Method with Lower Bound Preservation for Navier-Stokes Equations in Aerodynamic Shape Optimization AU - Zhang , Bin AU - Yuan , Weixiong AU - Wang , Kun AU - Wang , Jufang AU - Liu , Tiegang JO - Communications in Computational Physics VL - 5 SP - 1381 EP - 1408 PY - 2023 DA - 2023/06 SN - 33 DO - http://doi.org/10.4208/cicp.OA-2023-0020 UR - https://global-sci.org/intro/article_detail/cicp/21765.html KW - Aerodynamic shape optimization, solution remapping technique, direct discontinuous Galerkin method, lower bound preservation, Navier-Stokes equations. AB -

It is found that the solution remapping technique proposed in [Numer. Math. Theor. Meth. Appl., 2020, 13(4)] and [J. Sci. Comput., 2021, 87(3): 1-26] does not work out for the Navier-Stokes equations with a high Reynolds number. The shape deformations usually reach several boundary layer mesh sizes for viscous flow, which far exceed one-layer mesh that the original method can tolerate. The direct application to Navier-Stokes equations can result in the unphysical pressures in remapped solutions, even though the conservative variables are within the reasonable range. In this work, a new solution remapping technique with lower bound preservation is proposed to construct initial values for the new shapes, and the global minimum density and pressure of the current shape which serve as lower bounds of the corresponding variables are used to constrain the remapped solutions. The solution distribution provided by the present method is proven to be acceptable as an initial value for the new shape. Several numerical experiments show that the present technique can substantially accelerate the flow convergence for large deformation problems with 70%-80% CPU time reduction in the viscous airfoil drag minimization.