基于N-S方程的直升机旋翼翼型反设计方法

为了在实际旋翼翼型设计中考虑粘性影响,采用N-S方程为主控方程,建立了旋翼翼型流场求解方法。使用以Poisson方程为控制方程的网格生成程序生成围绕旋翼翼型的N-S粘性贴体网格,该网格生成程序能在反设计迭代过程中自动更新过渡翼型的粘性网格。在翼型流场计算及网格生成基础上,采用了MGM方程作为翼型反设计方程,建立了一套能够考虑粘性影响的直升机旋翼翼型的反设计方法。应用该方法,分别对二维NACA系列翼型、OA系列翼型、超临界翼型进行了反设计分析,获得了满足要求的二维旋翼翼型,并与目标值(压力分布、升力系数、阻力系数、力矩系数)吻合良好。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

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

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

复杂空气动力外形网格生成技术

通过对用于网格生成的Ｌａｐｌａｃｅ方程进行抛物化差分离散,发展了一种快速抛物型方程贴体及与边界正交网格生成方法,另外,采用源项修正技术求解用于网格生成的Ｐｏｉｓｓｏｎ方程,发展了一种可实现对网格线与边界间距及正交性进行双重控制的椭圆型方程方法,运用这种两种基本方法以及这两种方法的混合,分别生成了真实歼击机全机外形,带十字叉翼导弹以及双发,下单翼运输机外形的计算网格。并已成功地秀于跨声速得杂流场的娄     

液氢/液氧火箭发动机燃烧过程二维数值模拟分析

本文应用数值模拟方法对液氢/液氧火箭发动机燃烧过程进行了计算,给出了发动机燃烧室内燃气参数分布。计算中通过随机喷雾模型确定液氧滴进入燃烧室的初始条件,应用高压蒸发模型计算液氧蒸发速率及Arrhenius公式计算氢氧气相化学反应速率。贴体网格由TTM方法生成,数值方法采用一步隐式预测、两步显式校正的PISO算法,控制方程在非交错网格下离散并且对流项采用了二阶迎风格式。      

涡轮发动机复杂区域的结构化网格生成技术

通过对二维泊松方程增加伪时间项，采用时间推进法求解生成二维网格，利用拉格朗日插值生成三维网格，给出了一种在涡轮发动机复杂区域中生成结构化网格的计算方法。该方法很好地保证了网格的质量和正交性。大量的计算算例表明该网格的生成方法和区域分解算法相结合可以对复杂区域流场的流动细节进行很好的模拟。     

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

全机复杂流场网格生成及N-S方程数值模拟

全机复杂流场网格生成对计算结果起着举足轻重的作用。本文首先通过合理的拓扑结构,分块对接的方法对复杂流场进行了网格划分,然后采用N-S方程有限体积法及分区对接法对翼—身—尾全机流场进行数值模拟。计算结果验证了网格生成的有效性和计算结果的合理性。     

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

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

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。     

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

基于CFD/CSD方法的跨声速静气动弹性数值模拟应用研究

针对大型飞机跨声速静气动弹性问题,基于CFD/CSD流固耦合方法开展了数值模拟应用研究。具体方法是:采用RANS方程与静气动弹性平衡方程作为计算控制方程,CFD采用多块对接结构化网格进行分区并行计算,并利用多重网格技术加速收敛;使用RBF结合TFI方法进行网格变形;通过TPS插值实现多场数据交换。基于上述方法开发了计算程序,并与典型风洞试验结果开展了对比研究,验证了方法和程序的有效性。最后通过数值模拟结果,分析了跨声速时静气动弹性对典型大型飞机机翼的几何变形、表面压力及气动性能的影响特性。     

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 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

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

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

通过自动生成复杂构型的非结构混合网格。采用格心格式的有限体积法以适应不同的网格单元，应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟，将部分计算铰链力矩结果与风洞试验结果进行了对比分析。说明了计算方法对粘性流动良好的适应性,分析了迎角、马赫数、雷诺数对舵面铰链力矩的影响以及计算方法的可靠性和计算结果的合理性。     

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

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

基于隐式算法的悬停旋翼粘性绕流高效CFD分析方法

为显著提高旋翼粘性绕流CFD模拟效率,建立了一套高效的基于隐式LU-SGS算法和OpenMP并行策略的旋翼悬停流场求解方法。首先,求解Poisson方程生成桨叶剖面翼型的贴体正交网格,并通过剖面网格插值、翻折方法生成桨叶C-O型贴体网格;在此基础上,采用基于"扰动衍射"挖洞方法与"Inverse map"相结合的洞边界划定与贡献单元搜索方法,解决了嵌套网格技术中的相关瓶颈问题。然后,以非惯性系下耦合S-A湍流模型的RANS方程为流场主控方程,对流通量采用三阶Roe-MUSCL格式进行离散,时间推进采用高效的隐式LU-SGS方法,同时采用基于数据共享的OpenMP并行策略加速流场求解。最后,运用所建立的方法分别对不同旋翼翼型和悬停状态"Caradonna-Tung"以及UH-60A旋翼的流场及气动特性进行了计算,并给出了根据涡核位置加密网格来提高桨尖涡捕捉精度的方法,同时将计算结果与试验值进行了对比,验证了该方法在旋翼CFD流场模拟中的高效性和高精度特征。     

运输机组合体绕流非结构混合网格N-S方程数值模拟

本文利用通用网格生成软件生成运输类飞机翼身组合体（DLR-F6）非结构混合网格，在边界附面层区域内，生成三棱柱形网格，在其他流动区域采用四面体网格进行填充。粘性流场求解器采用基于N-S方程的CFX5计算流体软件求解，通过计算结果与风洞试验所得到的数据对比分析，验证了网格生成方法及网格数据的正确性和实用性。本文研究结果也为运输类飞机高速风洞试验提供了对比数据。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。