混压式超声速进气道流场数值模拟

采用有限体积法对积分形式的Navier-Stokes方程组进行空间离散,所获得的半离散常微分方程采用多步Runge-Kutta方法推进求解。为增大稳定性允许的CFL数和加快收敛速度,采用了当地时间步长、残差光滑技术和多重网格法等策略。数值模拟了二维混压式进气道,给出了该进气道的内外流场结构和表面压力分布。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

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

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

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

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

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

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.      

风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0),采用刚性运动网格技术和动态拼接网格技术,开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究,影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程,采用有限体积法离散控制方程,离散方程的时间方向采用“双时间步”方法求解,空间方向无粘项离散采用 MUSCL-Roe 格式,湍流模型采用 SA 和 SST 湍流模型,并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构,进而影响吸力面的压力分布,而对压力面的流动结构和压力分布基本没有影响;采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

CFD技术在带动力飞机气动设计中的应用

应用数值模拟的方法,研究了飞机带动力发动机的流场模拟,基于此计算分析了其性能特性.采用基于雷诺时均的N-S方程,动力方程中的雷诺应力通过SST二方程湍流模型和自动壁面处理计算;同时,采用守恒隐式有限体积法.离散方程使用有限高精度对流程序类似于Barth和Jesperson进行计算.质量流动通过Rhie和Chow算法压力-速度耦合进行计算.     

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

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

涡流发生器数值计算方法研究

利用RANS方程和SA湍流模型解算器,采用多重网格法和预处理技术,对叶片式涡流发生器进行了数值模拟研究,初步探讨了涡流发生器的安装方式、剖面形状、几何尺度等因素在机翼分离流动控制中的影响规律和设计原则,并初步研究了安装涡流发生器对超临界机翼气动性能的影响。     

用于粘性可压缩流动数值计算的SIMPLE方法

对SIMPLE方法加以推广,使之适用于可压缩粘性流动的数值计算。采用非交错网格技术,推出三维任意曲线坐标系下可压缩形式SIMPLE方法的计算方程。对亚音速、跨音速及超音速等5个流场进行了数值计算并与有关文献及实验数据进行比较。      

求解可压缩流动的同位网格SIMPLE方法研究

在Rhie-Chou动量插值的基础上,推导了同位网格可压缩SIMPLE算法.经过无粘流超音速凸包算例和激波/湍流边界层干扰算例计算发现,如果对流项采用高阶有界HLPA格式,密度插值采用一阶迎风和中心差分的混合格式,这种算法能够很好地模拟凸包超音流的流动现象,在采用了新型GAO-YONG湍流模型后也能够较好地模拟激波/湍流边界层干扰.      

非结构同位网格SIMPLE类算法收敛性能比较

以非结构网格的SIMPLE算法为基础,将算法扩展为SIMPLEC算法.利用30°角的斜方腔流动计算成果,分析了非结构同位网格的SIMPLE/SIMPLEC算法的收敛性能;比较了因网格的非正交而引入的非正交项的取舍对该算法收敛性能的影响;并采用显式校正步法对SIMPLEC算法进行了显式校正.比较表明,在非结构同位网格SIMPLEC算法中可忽略非正交项,但有必要对压力作亚松弛.显式校正步法可显著地加速在非结构网格上求解N-S方程的收敛性能,而且在不同的松弛因子组合下,均有较好的收敛速率.     

突扩后突缩的三维流场计算

根据有关固体火箭发动机装药后端到喷管间的空间特征,对三维突扩后突缩内部流场计算进行了探索和研究。计算中成功地应用了SIMPLE算法解决复杂的三维流场计算问题,并在初场构成、源项离散及松弛系数的选取等方面都进行了调试改进。通过建立计算模型、划分网格、对方程及边界条件离散化和非线性处理,以及线性方程组求解等,计算得到了符合物理规律的数值解,指出了空间涡及二次涡的存在。     

可压缩流动SIMPLE算法压力修正方程边界条件处理探析

可压缩粘性流体Nivier-Stokes方程SIMPLE算法求解中压力修正方程起着关键的作用.提出从连续性方程本质出发确定压力修正方程边界条件的方法.通过对双圆弧(凸包)通道内流场的计算验证了方法的正确性.     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

转、静盘腔径向出口流场与换热研究

采用SIMPLE算法对具有中心入射流的开式转静盘腔中的流动与换热进行了计算研究。计算了流场中 3D速度分量、温度场 ,得出了转静盘系层流流动的规律 ,分析了流场中出现旋涡的变化规律、出口截面的速度变化情况以及扩散项和对流项对温度场的影响     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

液体发动机燃烧室流场模拟的并行SIMPLE算法

用串行和并行ＳＩＭＰＬＥ算法对液体火箭发动机燃烧室内的复杂流动进行了数值模拟。控制方程组用欧拉坐标系下的Ｎ－Ｓ方程组描述。在并行虚拟机环境下研究了交错网格系统中参数的传输和接受模式, 从计算结果可以看出并行计算的效率较高, 是解决此类复杂流动大规模数值模拟问题的有效手段之一。