一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

冲压发动机进气道压力振荡过程的数值研究

针对来流马赫数3.5,0°迎角飞行条件下冲压发动机的简化模型,数值模拟了进气道整流罩开启过程的非定常流动。计算采用AUSMPW有限体积格式,基于结构与非结构的混合网格技术,求解了非定常N-S方程。计算结果表明,在进气道整流罩打开以后,超声速气流会在燃烧室通道内形成振幅衰减较慢的低频剧烈振荡。在同样计算条件下采用Euler方程进行对比性研究,数值实验发现对这类大雷诺数、强非定常性的问题,Euler方程和N-S方程得到的结果相差不大。     

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

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

自适应直角切割网格民机增升装置绕流数值模拟

提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明所述方法的正确性.     

基于伴随方程方法的非结构网格自适应技术及应用

对基于伴随方程的网格自适应技术的基本原理进行了研究,构造并实现了基于流场伴随方程的非结构网格自适应探测器,利用局部加密网格的办法,建立了基于Euler方程的无粘伴随自适应能力。对于粘性流动问题,提出首先应用无粘计算和伴随自适应技术从基准网格自动生成具有较好网格分辨率的非结构自适应网格,然后采用层推进法自动生成粘性计算网格的求解策略。采用RAE2822二维翼型和ONERA-M6三维机翼对所建立的无粘伴随自适应方法进行了验证,并应用所提出的方法和策略,对第二届国际涡流试验项目(VFE-2)中的65°尖前缘三角翼大攻角涡流场进行了数值模拟。对比计算结果表明,所建立的非结构网格伴随自适应方法是可行的和有效的。     

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

三维Euler方程非结构自适应网格投影和光顺技术研究

研究了三维Euler方程非结构自适应网格的投影和光顺技术.以ONERA M6机翼为计算模型,非结构网格为计算网格,守恒变量的范数为自适应加密准则,对采用一般自适应计算策略得到的网格和数值结果进行了分析,并发展了一种考虑网格投影和光顺的自适应计算策略.数值试验结果表明:通过对自适应所生成网格的投影与光顺,实现了自适应计算过程中的几何保真性,获得了相对光滑的压力分布曲线.     

目标函数误差估算及网格自适应处理

发展了一种基于目标函数误差估算的网格自适应准则,进而通过网格自适应处理提高目标函数计算的效率和准确性。首先描述了目标函数的误差估算及修正的方法,该方法通过伴随方程将原方程的残值误差与目标函数联系起来。然后,基于误差估算建立网格自适应准则,以减少目标函数修正后的剩余误差,提高目标函数计算结果的准确性,并将方法进一步发展至多目标问题。最后将该准则应用于Euler方程控制的NACA0012翼型无粘可压流动的网格自适应模拟。数值计算成功地捕获了与升力、阻力和力矩等目标函数相关的特征流动区域,计算出符合指定精度要求的目标函数值,验证了本文所发展的方法。     

基于Euler方程的三维自适应笛卡尔网格应用研究

本文采用非结构自适应笛卡尔网格和有限体积法求解了三维Euler方程。在网格生成过程中，采用了基于模型几何数据、模型表面曲率和流场特征为基础的网格自适应生成及加密方法，采用模型特征恢复技术保证了计算模型与原始外形的一致。在计算中，将以中心差分为基础的Jameson有限体积法在三维笛卡尔网格中进行了推广。本文中对导弹、歼击机等复杂问题进行了数值实验，计算结果与风洞实验结果符合良好，并在工程实际中成功应用。     

模拟存在运动激波流场的无网格自适应算法

将无网格自适应算法发展用于求解存在运动激波的非定常流动问题。为了降低初始布点要求,同时提高对运动激波的分辨率,提出了一种结合点云加密和点云粗化的动态点云技术。采用该技术,在非定常计算的每一个时间步,首先通过基于压强梯度的流场探针识别出激波所在位置,然后对激波附近点云进行加密,从而提高对激波的分辨率,同时对远离激波的点云进行粗化,尽可能减少冗余节点,提高计算效率。采用所提出的无网格自适应算法求解了N-S方程,首先对一维运动激波问题进行了数值模拟,通过将数值模拟结果与精确解比较,验证了所提算法求解存在运动激波的非定常流动的可行性,在此基础上,将方法推广用于二维非定常流动问题的求解,分别对瞬时激波碰撞圆柱以及双弧翼型激波振荡流动进行了数值模拟,计算结果表明,方法可以显著提高对激波的分辨率,同时相对于直接采用较密布点的非自适应算法,自适应算法在计算效率上具有优势。     

基于粘性笛卡尔网格的N-S方程计算

本文采用基于叉树数据结构的粘性笛卡尔网格求解N-S方程。为了使网格的分布更为合理,采用基于物体几何外形的方法对网格进行自适应加密。在物面处采用光顺、投影、分割加密的方法保证了物面处网格的正交性与贴体性,从而生成高质量粘性笛卡尔网格。将有限体积法应用于粘性笛卡尔网格求解N-S方程,采用了SA紊流模型。算例计算结果与实验数据符合良好。     

二维舵面绕流的N-S方程数值计算研究

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

翼型扑翼运动流场的数值模拟

基于非定常N-S方程,利用动网格方法和预处理方法对翼型在低速度下扑翼运动进行数值模拟。采用Jameson中心格式的有限体积法、双时间推进法求解非定常N-S方程;紊流模型采用基于SA紊流模型的DES模型。本文计算模拟翼型在低雷诺数下扑翼运动过程,计算所得推力系数结果与实验参考结果吻合较好,对扑翼飞行器设计有一定参考价值。     

高超声速连续/稀薄流自适应界面推进重叠网格方法及应用研究

根据高超声速飞行器在过渡流域飞行时绕流流场具有连续流与稀薄流并存的特点,提出了一种适用于混合流场自适应界面推进的多尺度重叠网格方法。用两套独立生成不同尺度的结构网格和非结构网格覆盖整个流场。根据连续/稀薄流界面失效函数的截止值在结构网格内进行自适应"挖洞"。在挖洞后的稀薄流区利用非结构网格DSMC方法进行计算;在连续流区的结构网格内采用基于N-S方程的CFD算法。两种算法进行耦合计算和求解,在连续/稀薄流分界面进行流场信息的插值、交换和处理。根据流场实时耦合计算结果,利用界面失效函数自适应改变两套网格的重叠部分。最后,通过对超声速圆柱和高超声速带扩张角圆管绕流的数值模拟,验证了方法的可行性和有效性。     

三维直角叉树切割网格Euler方程自适应算法

基于叉树数据结构，实现了一种用于三维直角叉树切割网格的自适应算法，包括对几何外,形以及对流场计算的自适应网格加密技术。在初始网格的生成过程中，根据相邻网格的物面法向向量间的差值，进行针对外形的自适应网格加密；在流场计算中，根据相邻网格间选定物理量梯度的变化，进行针对流场的自适应网格加密。详细地描述了三维直角叉树切割网格的生成过程，以及以任意网格的切割细分算法。在自适应过程中，分别采用了八叉树和全叉树的数据结构，八叉树是基本的数据结构，而全叉树的采用，使网格具有了各向异性的特征，从而大大的减少了自适应网格的数量。采用中心有限体积法，求解Euler方程，并运用上术方法，完成了对外形和流场的逢适应网格加密算法，获得了较好的数值计算结果，证明了自适应算法的正确性，体现了直角叉树切割网格自适应技术的有效性和实用性。     

螺旋桨流场Euler/NS数值解对比分析

采用多块点对点网格生成技术,生成螺旋桨流场空间计算网格。针对不同前进比,分别求解Euler和NS方程对螺旋桨流场进行数值模拟,并与实验结果进行了对比分析。研究发现无粘和有粘的计算结果均接近于实验结果,因此基于Euler方程的流场数值模拟更适用于工程设计。     

基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

基于旋转爆震三维流场结构分析的计算模型对比研究

为了研究不同计算模型对三维旋转爆震数值模拟的影响,分别基于Euler方程、N-S方程、RANS方法和IDDES方法并结合滑移、无滑移壁面边界条件,耦合氢气/空气的有限化学反应速率模型（7组分8基元反应）,采用高分辨率的五阶有限差分格式WENO-PPM5离散对流项,对圆环筒型燃烧室内的三维旋转爆震波进行数值模拟,得到了不同模型下旋转爆震波的流场结构、传播特性和推力性能。对比了不同计算模型对流场结构的影响。当采用滑移壁面边界条件时,Euler方程和N-S方程的计算结果较为一致。当使用无滑移壁面边界条件时,边界层的存在会导致可燃混气与燃烧产物之间接触面上的爆燃燃烧区域沿壁面向上游渗透,增大爆燃区域的范围;且爆震波锋面非受限侧变窄、严重变形,不同计算方法计算的变形程度有所不同。     

自适应直角坐标网格方法在对流换热数值模拟中的应用

复杂几何边界是对流换热计算面临的主要难点之一。自适应直角坐标网格方法可以比较方便的处理复杂几何边界,但其应用目前主要集中于可压无粘流的Euler方程、可压粘性流的Navier-Stokes方程的外部绕流计算。本文应用自适应直角坐标网格方法计算了两个不可压内流对流换热算例,并把计算结果同商用CFD软件STAR-CD的计算结果进行了对比。自适应直角坐标网格方法在本文的两个算例中得到了同STAR-CD一致的结果。这说明采用该方法数值模拟对流换热问题是适用的。      

一种适合于旋翼涡流场计算的非结构自适应嵌套网格方法

针对旋翼涡流场的特点,提出了一种新的非结构自适应运动嵌套网格方法。在该运动嵌套网格方法中,贴体网格采用非结构网格,背景网格使用可自适应的笛卡尔网格,以便发挥非结构网格的形状描述能力和笛卡尔网格的空间自适应优势,方便网格生成过程并提高计算精度,同时克服在物体表面附近自适应笛卡尔网格生成工作量大和过程烦琐的不足。在此基础上,使用求解N-S主控方程的方法,分别对旋翼的悬停和前飞流场进行了数值模拟,计算结果表明本文的非结构自适应运动嵌套网格方法对提高旋翼涡流场的细节捕捉能力和计算精度是很有效的。