双模态冲压发动机燃烧室流场数值模拟研究

在三维、粘性、湍流及有化学反应的Navier-Stokes方程基础上，通过有限化学反应速率／涡扩散模型模化湍流燃烧，对以H2为燃料的双模态冲压发动机燃烧室流场进行了研究，分析了空燃比、燃料入射角、飞行马赫数对燃烧室工作模态的影响，并分析了燃烧室隔离段的作用。     

某型导弹尾喷流形状的数值模拟

从三维薄层近似N—S方程出发，采用EN0差分格式，对某型导弹层流场进行了数值模拟。计算模拟了飞行马赫数Mα0．8—1．6，海拔高度H＝2km—15km，燃烧室总压为20MPa和40MPa条件下导弹的层流场。计算结果表明与飞行速度相比，海拔高度和燃烧室总压对层流宽度有显著影响。通过对流场结构的分析与测量，建立了层流宽度与上述参数的拟合公式，为工程估算这一类型导弹层流影响边界提供参考。     

动静圆柱坐标系下N-S方程的分区算法及子域干涉处理研究

耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程,统一了动子域和静子域的控制方程;在子域计算中,采用高分辨率的三阶ENN格式及LU-SGS隐式解法;在子域干涉处理中,研究并对比了动静子域干涉面的重合处理与重叠处理;最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气机叶栅内三维紊流流场的数值模拟,验证了该分区算法的有效性,计算结果与实验数据吻合较好。     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数滴流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采和了自适应网格技术。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法，建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题，旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法，通过添加冗余约束，避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块，对基于运动嵌套网格 的空间流场进行求解，湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块，基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解，验证了网格变形模块和CFD模块的有效性，然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上，计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷，并与飞行测试数据进行了对比。计算结果表明，文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析，提高了旋翼气弹载荷的预测精度。     

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.     

A coupling VWM/CFD/CSD method for rotor airload prediction

A coupling fluid-structure method with a combination of viscous wake model(VWM),computational fluid dynamics(CFD) and comprehensive structural dynamics(CSD) modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method' is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。