Solution Remapping Method with Lower Bound Preservation for Navier-Stokes Equations in Aerodynamic Shape Optimization  

作　　者：Bin Zhang Weixiong Yuan Kun Wang Jufang Wang Tiegang Liu 

机构地区：[1]LMIB and School of Mathematical Sciences,Beihang University,Beijing 100191,China [2]Computational Aerodynamic Research Institute,China Aerodynamics Research and Development Center,Mianyang 621000,China

出　　处：《Communications in Computational Physics》2023年第5期1381-1408,共28页计算物理通讯（英文）

基　　金：This project is supported by the National Natural Science Foundation of China(No.12001031).

摘　　要：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 problemswith 70%-80%CPU time reduction in the viscous airfoil drag minimization.

关 键 词：Aerodynamic shape optimization solution remapping technique direct discontinuous Galerkin method lower bound preservation Navier-Stokes equations 

分 类 号：O175[理学—数学]