AUSM+方法在黏性内流模拟中的应用

简要阐述了三维雷诺平均Navier-Stokes(RANS)方程求解程序中使用的对流迎风分裂格式改进形式(AUSM+)以及Spalart-Allmaras(S-A)一方程湍流模型,然后对Sajben收缩扩散通道内跨声速流动和高速透平静叶排内部流动进行了数值模拟,并将数值结果与实验结果进行了比较,结果表明自主开发的三维R...     

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

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

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

扩压器-水泵内非稳定流动研究

水泵内部流动实质上是复杂的三维非稳定流动 ,它对水泵性能及结构振动有重要影响。本文介绍了一种求解这种复杂内流动的数值方法。三维雷诺数平均的纳维斯托克斯方程 ( 3-DReynolds-averaged Navier-Stokes,RANS)以及标准 k-ε的方程用于描述水泵内非定常紊流流场。系统特性方程与水泵的 CFD模型相结合以考虑流体在管道中的加速 ;采用任意滑移网格界面模拟叶轮和静止部件之间的相互干涉 ;将整个叶轮作为分析对象 ,以考虑失速情况下流动的周向非对称。这些技术的结合包括了水泵内非稳定流动的物理实质。一台实验数据比较齐全的离心式 -扩压器水泵被用于验证所提出的数值方法     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

进口气流角对三维振荡涡轮叶片非定常流动影响的数值模拟研究

进一步发展了非定常雷诺平均N—S方程求解程序,并基于影响系数法,在三种不同振型下就进口气流角对三维涡轮振荡叶片绕流的影响进行了数值模拟研究。数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;进口气流角对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,且在不同模态下,进口气流角对振荡叶片非定常流动的影响规律表现不同。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

三维振荡叶栅内部非定常流动数值模拟研究

进一步发展了非定常雷诺平均N-S方程求解程序,并基于影响系数法对三种不同振型下的三维涡轮振荡叶片绕流问题进行了数值模拟研究.数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;振型对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,其中叶片运动导致叶栅各通道中流量的再分配以及出口气流角的变化等对叶片的负荷分布会产生较为显著的影响.     

粘性流场三维NS方程有限体积法数值解

提出了有限体积法中的隐式ＮＮＤ格式数值求解三维全ＮＳ方程的薄层三维ＮＳ方程的方法，对绕带座舱的再入飞行器高超声速粘性流动进行数值模拟。     

无尾布局支撑干扰数值模拟

采用TRIP2．0．SOLVER软件，开展了高、低速来流条件下，无尾飞翼布局支撑干扰的数值模拟技术研究。高速来流状态下支撑干扰的数值模拟，以尾支撑作为主支撑，腹支撑作为辅助支撑，完全模拟了实验中采用的两步法；低速来流状态下的支撑干扰的数值模拟，采用了有无尾支撑的方式得到支撑干扰量。介绍了支撑干扰的典型高、低速数值模拟结果，并与高速实验结果做了初步的对比。     

多块多重网格法及其跨声速转子内流并行模拟

为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。     

二维埋入式进气道的数值分析

采用并改进了适用于分离流动的 Ｊｏｈｎｓｏｎ Ｋｉｎｇ 紊流模型,结合 Ｌ ＵＳ Ｇ Ｓ近似隐式格式和基于 Ｍ Ｕ Ｓ Ｃ Ｌ非线性插值的 Ｒｏｅ 格式编写了雷诺平均的二维 Ｎ-Ｓ方程的求解程序。然后借助于该流场计算程序完成了二维埋入式进气道的数值研究,分析了几何设计参数对进气道性能的影响。     

定常吹/吸气控制凸包分离的数值研究

针对凸包流动分离主动控制开展了参数化数值研究。凸包模型为NASA设计用于流动分离主动控制的内外流通用模型。重点研究运行参数变化时定常吹/吸气影响。首先采用试验数据标定了所开发程序,然后进行了一系列参数化数值模拟。结果显示:定常吸气作用明显不同于定常吹气。对于所给定几何模型,不管吸气流量增加到多大,吸气都无法消除分离,而吹气可以做到。详细展示流动结构并讨论了流动机理。      

涡流发生器对Bump进气道性能影响数值研究

以一种Bump进气道为研究对象,通过在S弯扩压段入口处布置涡流发生器来控制流动分离,减小出口总压畸变。采用CFD数值计算软件对Bump进气道在设计点(Ma=2.0)与非设计点(Ma=1.8,0.8)工况下内、外流场进行计算,分析不同涡流发生器方案的效果。计算结果表明:在设计点工况下,安装涡流发生器能够抑制流动分离,改善进气道流场品质,减小出口总压畸变;在一些非设计点工况下会增大Bump进气道出口总压畸变;Bump进气道总压损失有所增大,不同叶片间距的涡流发生器对总压损失的影响相当。      

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

采用CFD方法对某型外压式超声速进气道流场进行了二维数值计算,叙述了网格生成、边界条件确定方法,控制方程求解收敛技巧以及计算结果的分析.数值计算结果与理论相符较好.     

一种改进的轴流涡轮叶尖对泄漏流影响的数值研究

 结合基于压力修正的采用雷诺应力湍流模型加壁面函数的三维计算流体力学程序,通过改进涡轮叶尖表面结构以期减小叶尖间隙泄漏损失,即在传统叶尖边缘镶嵌肋条的基础上倾斜压力面肋条,对这种改进叶尖的某一轴流涡轮转子的泄漏流场进行了数值研究,并详细分析了不同肋条高度和肋条深度对泄漏流场的影响,最后计算了涡轮效率。结果表明：通过改进的涡轮叶尖对泄漏流动进行被动控制,压力面斜肋条可以阻碍泄漏流动,同时叶尖区域回流区也减小泄漏流动,从而涡轮效率提高;肋条高度对涡轮效率有较大影响,且有个最佳肋条高度值,在最佳肋条高度下涡轮效率提高0.215%;肋条深度对流场影响不大,但小肋条深度可以稍微提高涡轮效率。     

跨声速轴流压气机特性全三维粘性流动分析

应用Steger-Ｗarming通量分裂技术,将守恒型方程中的流通向量分裂成两部分,应用 NND格式对无粘项进行数值离散,粘性项采用中心差分格式,采用LU-SGS隐式推进迭代建立了跨声速轴流压气机性能及内流场三维N-S方程的高分辨率和高效率的数值方法。数值模拟给出NASA ROTOR37和NASA ROTOR67性能及内流场流动结构。并进一步基于混合平面法建立了多级轴流压气机性能计算模型,详细分析了某一级跨声速轴流压气机性能和内流特性,计算与实验结果比较表明了本模型与方法是可行的和可靠的。     

基于螺纹锥形搅拌头的搅拌摩擦焊三维流场仿真研究

本文利用CFD软件包FLUENT建立了带螺纹表面的锥形搅拌头搅拌摩擦焊三维流动模型。通过数值模拟得到了焊缝沿厚度方向的水平切面速度矢量图,截面的速度矢量图以及用于观察材料流动行为的流动线图。同时分析了搅拌头表面螺纹对塑化金属流动行为的影响及焊缝塑化金属的流动规律。模拟结果对于搅拌摩擦焊过程中焊缝处塑化金属流动以及焊缝形成机理研究具有重要的指导意义。