A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

一种基于谱方法的三维热声不稳定模型

研究流体的流动参数,包括热源热释放的周向不均匀分布对热声不稳定现象的影响,将欧拉方程进行线化处理之后,对轴向和径向方向进行基于切比雪夫多项式的谱离散,周向方向则进行基于傅里叶多项式的谱离散,建立了一种三维热声不稳定数值模型.通过对圆柱管道共振频率以及声场的求解,进行了一系列算例的验证,证实了该方法的可靠性,并且通过热释...     

基于六自由度方程的飞机结冰问题仿真

基于六自由度方程对飞机的结冰飞行动力学问题进行了仿真研究。首先针对飞机遭遇结冰时的结冰飞行过程进行仿真,开发了一种结冰过程中飞机结冰严重程度及气动参数的计算方法。通过在各仿真时刻对飞机气动参数进行结冰影响修正以实现对飞机结冰飞行过程的模拟,同时分析了结冰过程对飞行特性的影响。进而对飞机结冰后带冰飞行时不同结冰严重程度下的操纵面阶跃响应进行了仿真,并对仿真结果进行了分析。     

一种求解S形进气道三维粘性流场的多重网格算法

提出了一种基于多重网格技术,可快速求解Ｓ形进气道粘性内流场的实用工程算法：从Ｂｒａｎｄｔ的多重网格理论出发,采用隐式Ｂｅａｍ－Ｗａｒｍｉｎｇ时间推进格式,结合“锯齿式”迭代技术交替在粗、细两种网格上对雷诺平均三维非定常Ｎ－Ｓ方程进行求解。用此算法数值模拟了某型号Ｓ形进气道的内流场,与试验结果对比显示：计算得到的流场符合一般Ｓ形管道流动特性,所捕获的二次流图谱清晰合理,进气道出口截面处总压恢复系数与试验值较为接近。与普通单层网格计算方法相比,计算结果精度高,收敛速度快。     

双欧法在克服伞-弹系统欧拉方程奇异性中的应用

文章针对伞弹系统欧拉方程奇异性 ,利用正、反欧拉方程精华区、奇异区呈倒挂的原理 ,建立了一套与飞机不同的反欧拉系统。将这套双欧系统应用在克服伞 -弹系统欧拉方程奇异性中 ,取得了满意的结果。同时对双欧法与四元数法计算结果进行了对比验证 ,从而更加清楚地看到了四元数法的方法缺陷和不足 ,为进一步能将双欧法应用于伞 -弹系统六自由度弹道实时仿真进行了有益的尝试     

用N-S方程有限体积法计算弹体绕流/底喷流

用单元中心有限体积法在物理空间曲线网格离散积分形式N-S方程，用坐标变换计算粘性通量项中的各偏导数值，用Runge-Kutta三步显式格式进行时间推进求解。采用了人工粘性、当地时间步长及隐式残值光顺。计算了旋成体大迎角涡流和带底部喷流的绕流场。计算中发现，各向异性隐式残值光顺较之各向同性隐式残值光顺能更快地加速迭代过程的收敛。计算的物面压力分布与实验基本符合，计算的旋成体大迎角涡流图谱合理。计算结果说明底部喷流对旋成体绕流及空气动力载荷有显著影响。     

高速风洞纵、横向支架干扰数值修正研究

本文提出了一种计算高速风洞模型支架形式对飞行器模型纵、横向气动力干扰量的数值计算方法。文中从跨音速全位势积分方程出发，编制了适用于计算支架系统对飞行器模型亚、跨、超音速纵、横向气动力干扰量的计算程序。通过对各种飞机模型纵、横向实验结果的修正计算表明，该计算程序使用方便，计算速度快，在中、小攻角和小侧滑角范围内计算结果与实验值有较好的一致性。可以做为选择模型支架形式和几何外形参数的理论工具，同时也为高速风洞纵、横向实验数据的支架干扰修正提供了可行的技术方法。该计算程序还具有飞行器模型几何外形和表面等压线以及压力分布云图显示功能。     

含裂纹固体推进剂装药的内流场数值模拟

利用“N－S方程＋粒子仿真”方法来计算含裂纹翼柱型装药内流场。计算结果表明：环形裂纹内的流动对主流的影响大于纵向裂纹，环形裂纹向上游位置流动对主流的影响较下游位置的大，并可保持到喷管喉部，翼槽附近环形裂纹内的流动对主流的影响相对较弱。     

High-order implicit discontinuous Galerkin schemes for unsteady compressible Navier–Stokes equations

Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.     

模型燃烧室紊流燃烧的大涡模拟

采用了两种不同的亚网格尺度燃烧模型对带 V型稳定器的模型燃烧室紊流化学反应流动进行了大涡模拟 ,用 k-ε方程亚网格尺度模型确定亚网格紊流粘性 ,为了考虑热辐射对燃烧室壁温和气流温度的影响 ,运用热流法辐射模型估算热辐射通量 ,用 SIMPLE算法和混合差分求解大涡模拟各守恒方程 ,通过对两种不同亚网格尺度燃烧模型数值模拟结果与实验值的比较表明 ,两种燃烧模型都与实验值较吻合 ,但 G方程小火焰模型要比亚网格 EBU燃烧模型符合得更好些