非正交曲线坐标湍流旋流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标系统，用坐标变换的方法处理锥形燃烧室内线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。紊流模型采用工程上常用的ｋ－ε模型。还运用贴坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提出了理论依据。     

低速高湍流度90°弯管流动数值模拟

分析了曲率对弯曲管道流动的影响 ,给出了曲率修正的双层 k-ε湍流模型 ,并数值模拟了非均匀来流低速高湍流度 90°弯管内流动。文中比较了曲率修正双层 k- ε湍流模型的计算结果与实验结果 ,以及双层 k- ε湍流模型和逆压力梯度修正双层 k- ε湍流模型的计算结果。通过比较 ,发现经过曲率修正后的双层 k- ε湍流模型更好地模拟了非均匀来流低速高湍流度 90°弯管内流动。同时比较了壁面参数的插值方法 ,发现在弯曲管道流动模拟中对壁面参数进行一阶插值比零阶插值能更好地模拟真实流场     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

环形燃烧室性能计算

在任意曲线坐标系下对包括扩压器和火焰筒在内的环形燃烧室三维两相反应流进行了数值模拟。燃烧室性能采用多维经验分析法估算所用的数学模型有 κ-ε,紊流模型 ,EBU -Arrhenius紊流燃烧室模型 ,六通量热辐射模型以及颗粒群轨道模型。在非交错网格体系下 ,气相采用 SIMPLE算法和液相采用 PSIC算法求解 ,通过两种工况和三种扩压器进口速度分布计算表明 ,所得的流场各气流参数分布和燃烧室性能较合理 ,本文所建立的程序可为燃烧室优化设计和研制提供有用数据。     

来流马赫数对座舱气动加热影响的数值模拟

采用数值手段研究了来流马赫数对飞机座舱附近区域气动加热的影响,计算的来流马赫数为0.8~2.0,基于翼展的雷诺数为107。通过分析飞机蒙皮表面的温度分布、表面摩擦因数、壁面法向速度梯度和热流量等物理量,给出了飞机头部座舱附近区域的气动热分布情况。计算结果表明,随着来流马赫数的增加,壁面法向速度梯度增大,飞机对称面上的平均温度、表面摩擦因数和热流量随之增大。此外,针对某型号飞机的计算结果表明,座舱附近区域上表面中心线上的平均温度与来流马赫数、自由来流温度之间呈一定的函数关系。     

燃烧室冷态流场的数值研究

在任意曲线坐标系下对环形燃烧室三维冷态流场进行数值模拟。采用偏微分方程法和区域法生成三维贴体网格,紊流模型采用标准k-ε模型和重整化群k-(εRNGk-ε)模型,在非交错网格坐标系中采用SIMPLE算法求解。数值分析两种紊流模型对火焰筒内气流流动的影响,并与实验结果相比较。结果表明,与标准k-ε模型相比,RNGk-ε模型更适用于模拟紊流流动。     

双流混合边界的正交贴体网格生成

通过数值实验表明，用泊松方程变换形式生成的贴体坐标网格，其边界点的β值与源项呈连续单调的函数关系。根据这种关系，可以用方程迭代求根方法中的增值寻根法和对分法，使网格边界正交化。由薄壁结构相隔的两种流体混合边界正交网格生成的实例表明，这种方法对于生成适应于红外抑制器全流场计算的网格是很有效的。     

基于格子Boltzmann方法的可压缩翼型绕流模拟

采用能够恢复可压缩Navier-Stokes方程的耦合双分布(DDF)格子Boltzmann方法(LBM)模拟了贴体网格下的NACA0012翼型绕流.首先在Ma＝0.5,α＝0.0,Re＝5 000状态下通过模拟低雷诺数流动检验了该方法;然后分别在Ma＝0.5,0.85和1.2,α＝-0.05,1.0和0.0状态下进行了模拟,并在Ma＝0.85时通过加密网格得到了更好的结果.为了避免转捩问题采用了无粘边界条件进行模拟,但得到的压力系数与实验结果吻合良好,证明了该方法的可靠性并显示了其模拟可压缩的潜力.该方法在翼型绕流上的应用为湍流流动的数值模拟提供了基础.     

高速列车绕流的数值模拟

给出了用低速N－S方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Train3D，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

三维管内激波/湍流附面层干扰流场的数值模拟

给出了可压缩性及曲率修正两方程湍流模型，用它对跨音速三维喷管管内激波／湍流附面层干扰流场进行了数值模拟；给出了不同纵向截面上的激波结构、等马赫线图和四周固壁上的摩擦力线谱。将计算所得的摩擦力线谱、激波结构与实验结果进行了比较。分析了横截面上的流动结构，计算与实验吻合较好，进一步证实可压缩性及曲率修正两方程湍流模型能够较好地模拟激波／湍流附面层干扰流场，满足工程设计和分析要求。     

用非交错网格求解曲线坐标系下的流动控制方程

本文采用非交错网格法,求解任意曲线坐标系下的流动控制方程。文中阐述了压力振荡的根本原因和相应的消除措施,导出了压力修正方程中控制体边界上曲线速度分量的修正表达式。在这些表达式中,出现了一个表征一个网格间距和两个网格间距压力梯度差分值的附加差值项。当出现压力振荡时,该差值很大,可以有效地消除振荡,而当压力场没有振荡时,这一差值又很小。作为计算方法可靠性的检验,本文分别计算了二维直通道和渐缩形混合管内的气流流动问题,结果是令人满意的。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

多重网格在突扩燃烧室流场计算中的应用

多重网格法是一种具有快速收敛性特点的新计算技术,该法在求解偏微分方程时用一系列逐步加密或减疏的网格去离散求解区域,不同粗细网格可以消除不同波长的误差,从而加速了收敛。本文将此法应用于计算突扩燃烧室流场。SIMPLE算法和紊流k-ε模型被采用来预估旋流器后气流速度u,v,w和紊流参数k,ε分布以及回流区尺寸。计算结果表明,在达到同样计算精度条件下,与单层网格相比,此法可以加快收敛速度,缩短机时,减少计算费用,从而使数值计算在工程应用方面得到进一步推广。     

弯曲管道内任意曲率正交曲线坐标系下的流体流动基本方程

本文基于任意正交曲线坐标系下的流体流动基本方程的一般表达式,针对内流流体力学数值分析中经常遇到的具有任意曲率S弯管情况,通过几何分析,建立了适用于矩形、圆形及椭圆形变曲率S弯管道的曲线坐标系,进而导出相应坐标系下的三维可压流体流动的基本方程。     

大涡模拟在加力燃烧室冷态流场中的应用

对带V型稳定器模型加力燃烧室二维冷态流动进行了大涡模拟。采用了代数和k方程两种亚网格尺度模型，并用SIMPLE算法和混合差分格式求解离散方程，用壁面函数处理固壁边界条件，得到了模型加力燃烧室二维瞬态流场的模拟结果。通过对不同亚网格尺度模型和不同入口速度下的瞬态流场的比较，揭示了不同亚网格尺度模型和不同入口速度分布对带V型稳定器燃烧室二维冷态瞬态流场的影响。     

串列叶栅粘性紊流流场有限分析数值研究

本文将近年发展起来的有限分析方法(Finite Analytic Method)用于曲线坐标系上紊流N-S方程的数值计算。计算了来流雷诺数为2.0×10~5,气流攻角分别为0°,10°,-10°三种情况的二维粘性紊流串列叶栅流场。文中用k—ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。数值计算结果与实验结果相比较,吻合程度令人满意。     

竖直上升通道内微气泡湍流流动的数值模拟

为了提高气泡浓度分布预测精度,考虑流体和气泡之间的双向耦合作用,对竖直上升通道内含有微小气泡的湍流流动进行了数值模拟研究。采用直接数值模拟方法模拟流体流动,并用拉格朗日法跟踪气泡受流体和重力作用所产生的运动轨迹。气泡受力包括相间阻力、虚质量力、升力、壁面升力、压力梯度和重力等。数值模拟结果表明,通过加入壁面升力和采用双向耦合方法,可以较准确地预测气泡浓度分布和揭示气泡流动对流体湍流结构的影响。     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型,采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式,对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟,探讨了内外喷流及喷管内的流动特征,并与实验数据进行了比较,计算结果与实验结果吻合较好,验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性,从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

驻涡燃烧室驻涡区三维冷态流动特性数值研究

驻涡燃烧室驻涡区内的流动是影响驻涡燃烧室性能的一个关键因素。首先将数值模拟结果进行与试验测量结果进行对比,确定了最佳的湍流模型,并在此基础上通过数值方法深入研究了驻涡燃烧室驻涡区冷态流动特性。结果发现:通过在驻涡区前壁进气缝中设置一定的矩形挡片可以在驻涡区一定范围内诱发以反向旋转的涡对形式存在的流向涡。对比研究了挡片阻塞比BR(挡片面积与前壁开缝面积之比)分别为0,0.2,0.4时驻涡区内的流动结构,分析了流向涡产生的原因。定量结果表明,当BR=0.2,0.4时,(1)流向涡涡量大小比BR=0时提高了将近100%,BR=0.2对应的流向涡涡量比BR=0.4对应的流向涡涡量略大;(2)在流向涡混合层内,流向涡涡量沿轴向呈先增大后减小的趋势;(3)燃烧室总压损失约比BR=0时大1%。