多重网格法求解离心压气机内部流场

本文提出了一种新的多重网格法,并将之与显式时间推进法、有限体积差分格式、粘性体积力法配合使用求解离心压气机内部流场。数值计算表明,作者多重网格法的计算结果与试验结果吻合比较好,同时该多重网格法的收敛速度比不用多重网格法加快了将近两倍。     

分离流动对矢量喷管性能的影响

采用有限体积法计算了轴对称矢量喷管三维粘性流场，通过改进人工粘性及其边界条件提高了计算的收敛性和数值模拟的精度。利用数值模拟研究了矢量喷管的分离流动对矢量推进性能的影响，通过同实验结果对比，验证了数值模拟的精度。研究表明，分离流动对矢量喷管的性能参数尤其是有效矢量角具有较大影响。     

微射流作动器三维外流场数值模拟

采用有限体积法求解Faver平均N—s方程，对三维、粘性、非定常、不可压微作动器外流场进行数值模拟。计算结果揭示了微射流产生、发展与耗散的过程，和实验中观察到的现象比较吻合。旋涡对沿展向的不稳定性是涡环破裂、旋涡对的相干性消失的主要原因。二次涡的发现为以后研究微射流旋涡对的细微结构打下良好的基础。     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

Simulation of 3D parachute fluid–structure interaction based on nonlinear finite element method and preconditioning finite volume method

A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual(GMRES) method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor(HHT) time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC(Harten–Lax–van Leer contact) scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel(LUSGS) approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.     

旋翼桨叶绕流Euler方程数值模拟

本研究应用格心格式有限体积法求解欧拉方程,模拟绕旋翼浆叶的流动。旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角。选择O—H型式生成弦向和展向网格。通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验。     

贫油预混预蒸发燃烧室回火现象的大涡模拟

为研究贫油预混预蒸发(LPP)燃烧室的回火现象及形成机理,基于大涡模拟(LES)及部分预混燃烧模型,对Schneider试验进行数值计算,分别将冷态及燃烧温度场与试验结果进行对比,验证了模型的准确性。将该模型应用于LPP燃烧室,并利用详细的化学机理对航空煤油火焰传播速度进行修正,成功重现了回火现象的发生并向主燃级上游传递的动态过程,通过流场分析发现局部形成的不规则涡结构在回火现象中扮演着重要角色。进一步从涡量方程的角度对涡结构进行分析,发现燃烧引起的斜压及体积膨胀是诱发回火现象的主要原因。     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

圆柱尾迹影响平板边界层转捩的实验研究

利用氢气泡和染色液流动显示的方法,在水槽中研究了圆柱尾迹对平板边界层转捩的影响。重点讨论圆柱与平板的流向距离和法向距离以及圆柱直径等几何参数对边界层转捩的不同作用,发现平板边界层转捩开始点的位置由尾迹与边界层作用的初始位置决定,这一位置由圆柱与平板的相对几何位置来确定。转捩区长度与边界层内部流体性质、扰动强度等有关,同时也受到圆柱、平板相对几何位置的影响。对转捩区近壁二次涡的产生和演化进行分析,发现尾迹对平板边界层的扰动是通过在近壁面诱导二次涡并影响二次涡的后续运动来体现,二次涡及其演化是边界层转捩位置大为提前的重要原因。