基于旋转坐标系的高阶间断有限元方法非定常湍流数值模拟  

Numerical Simulations of Turbulent Flow Based on Rotating Reference Frame Using High-Order Discontinuous Galerkin Method

作　　者：张涛[1] 吕宏强[1] 秦望龙 陈正武[2] ZHANG Tao;Lü Hongqiang;QIN Wanglong;CHEN Zhengwu(College of Aerospace Engineering,Nanjing University of Aeronautics&Astronautics,Nanjing,210016,China;China Aerodynamics Research and Development Center,Mianyang,621000,China)

机构地区：[1]南京航空航天大学航空学院,南京210016 [2]中国空气动力研究与发展中心,绵阳621000

出　　处：《南京航空航天大学学报》2019年第4期474-485,共12页Journal of Nanjing University of Aeronautics & Astronautics

基　　金：国家高技术研究发展计划(“八六三”计划)(2015AA015303)资助项目;国家自然科学基金(11272152)资助项目;航空基金(20152752033)资助项目;气动噪声控制重点实验室开放课题资助项目

摘　　要：基于高阶间断有限元方法(Discontinuous Galerkin method,DGM),对旋转非惯性系下耦合了修正的一方程S?A模型的RANS方程进行了离散求解。为了在稀疏网格上获得更贴近真实的物面形状,使用了多层高阶弯曲网格方法对物面进行拟合。非定常时间推进采用了隐式双时间步方法,每个时间步产生的线性系统采用预处理的方法,即广义最小残差方法(Generalized minimal residual method,GMRES)来求解。计算了旋转圆柱绕流以及经典翼型振荡算例的升力和力矩迟滞曲线,与实验结果以及前人的计算结果对比验证了本文方法的正确性和有效性。The RANS equations coupled with modified S?A model in a rotating reference frame are solved using the high?order discontinuous Galerkin method(DGM).In order to obtain more accurate wall boundaries on relatively coarse grids,a multilevel high?order curved mesh method is employed.An implicit dual time stepping method is adopted for the unsteady simulation,while generalized minimal residual method(GMRES)is deployed to solve the large linear system generated in each time step.The lift and moment hysteresis curves of rotating cylinder flow and classical airfoil oscillation are calculated.The correctness and validity of the proposed method are verified by comparing with the experimental results and previous calculation results.

关 键 词：间断有限元方法 弯曲网格 S-A模型 旋转坐标系 隐式双时间步 

分 类 号：V211.3[航空宇航科学与技术—航空宇航推进理论与工程]