非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

高超声速复杂气动问题数值方法研究进展

高超声速流场具有复杂流动特征,其中真实气体效应、磁流体干扰效应和力热结构耦合效应等对气动力分析产生了重要影响。将流体力学研究扩展到分子动力学、电磁流体力学以及流固耦合等交叉学科领域,这给数值模拟方法带来了巨大挑战。针对高超声速气动力/热分析的热点问题,重点关注高温效应与低密度流动效应、磁流体干扰效应和力热结构耦合效应等,结合算例分析了相应的数值求解技术;在气动热方面主要比较了3类求解方法(纯工程方法、纯数值方法和基于Prandtl边界层理论的方法),并给出了相应算例;对于气动力/热/结构耦合问题,从耦合模型及耦合计算方法两方面开展了分析。最后指出了高超声速复杂气动问题数值求解技术未来需重点关注的几个方面。     

高超声速风洞气动布局设计

在分析国内外高超声速风洞发展现状的基础上,根据南京航空航天大学高超声速风洞(Nanjing Universityof Aeronautics & Astronautics Hypersonic Wind Tunnel,NHW)总体技术指标和要求,对该风洞气动布局设计方案和备部件的气动设计进行了研究.风洞气动布局设计点为马赫数5和8、设计总压为1 Mpa、总温685 K;风洞驱动方式采用高压下吹-真空吸气式方案,运行时间大于10 s、高压气源容积为32 m3、真空容积为650 m3;风洞加热方式采用金属板蓄热式加热器方案;风洞试验马赫数获取方式采用φ0.5 m口径的马赫数5,6,7和8的型面喷管方案.     

基于基因算法与博奕论的翼型跨音速Euler方程气动优化

将基因算法与博奕论中的 Nash平衡相结合 ,构建了多目标基因优化算法—— Nash基因算法 (NashGAs) ,并对初始翼型为 NACA0 0 1 2二维翼型进行给定跨音速流动下的形状增升优化。计算中应用 Bézier曲线对翼型形状进行参数化 ,避免了非流线型的不合理形状产生 ;采用有限元非结构网格 AUSM+通量分裂格式Euler方程数值解进行个体适应度值评估 ;应用动态网格技术调整计算网络 ,节约了 CPU机时 ;最后给出了优化结果     

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格，该方法的特点是在生成单元的同时引入结点，与其他预布结点的方法相比，其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构，有效地提高了程序运行效果，最后给出了Ｍ６机翼的数值算例。     

The Full Flowpath Analysis of a Hypersonic Vehicle

对一种类X43-A吸气式高超飞行器全流道开展了M7一级的三维数值仿真研究,深入探索了进气道处于起动状态和不起动状态时全流道的冷流流场结构和和气动力特性,且部分结果与实验数据进行了对比。研究结果表明:(1)前体横截面上存在显著的展向压强梯度,使得经过预压缩的气流偏离了进气道进口,但同时也减少了进入内通道的边界层气流,提高了进口流场的品质;(2)后体喷流股的膨胀过程受到了周围外流的显著干扰,因而沿流动方向其截面形状不断发生变化,如在喷口附近为近似矩形,而在后体末端附近则演化为近似三角形;(3)当进气道处于不起动状态时,其全流道流动结构发生了显著变化,进气道的外部压缩波系往复振荡,尾喷管出口的喷流股也在不停的膨胀和收缩,具有极强的非定常特征;(4)当进气道处于不起动状态时,全机的气动力特性呈周期性变化,升阻比的变化幅度较大,在最大、最小值分别可达起动状态的2倍和1/4倍。另外,升力、阻力系数的变化曲线之间存在一定的相位差;(5)与实验结果的对照表明,所采用的数值仿真方法具有较高的精度。     

计算跨音速三维机翼非定常载荷的核函数法

 本文给出了在跨音速流动中,计算三维机翼非定常载荷分布的核函数法和计算结果。此方法是基于线化升力面理论,利用在亚、超音速情况下,联系载荷分布与下洗的积分方程的核函数,采用局部线化法来处理跨音速的非定常绕流问题。在计算中引进了线分布的冲波偶极子,并嵌入适当的正冲波边界条件。 我们从统一的超、亚音速核函数表达式出发,对亚音速区数值积分作了新的简便的处理,对于亚、超音速区相互的诱导下洗给出了自己的数值积分过程。我们算了三个例题,和文献上的数值符合较好。