可压缩后台阶流的数值模拟

基于可压缩Navier-Stokes方程,采用Beam-Warming矢通量分裂格式以及vanLeer的矢通量分裂方法,分别通过Johnson-King及Baldwin-Lomax湍流模型对亚音速和超音速的后台阶流进行了计算,计算结果与实验符合较好。     

双燃式冲压发动机中超燃燃烧室冷态流场数值模拟

采用Ｍａｃ Ｃｏｒｍ ａｃｋ 矢通分裂格式, 求解雷诺平均的Ｎ－Ｓ方程组, 采用Ｂａｌｄｗ ｉｎ－Ｌｏ－ｍ ａｘ 代数湍流模型和混合长度模型, 模拟超燃燃烧室内两股平行超声速来流的冷态流场, 取得了与试验相一致的结果, 计算表明回流区和射流流动参数对超燃燃烧室流场内波系的变化有着较大的影响。     

碳氢燃料超声速燃烧的化学动力学研究 (Ⅰ)甲烷超声速燃烧的简化化学动力学模型

采用雷诺平均Ｎ－Ｓ方程对碳氢燃料甲烷（ＣＨ４） 在超声速流动中混合燃烧过程进行了数值研究, 提出了简化的甲烷燃烧的反应机理及反应动力学模型, 建立了Ｂａｌｄｗ ｉｎ－Ｌｏｍ ａｘ 代数涡粘性湍流模型, 采用了高精度的ＴＶＤ隐式格式和化学源项隐式解法, 计算结果和氢气／空气超燃的计算结果进行了比较, 表明所采用的反应机理及计算模型合理, 算法可行, 就算例结果而言,碳氢化合物作为一种超燃燃料有其独特的优势, 值得进一步深入研究     

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

球心非对称的钝圆头体近区绕流场气动特性分析研究

本文采用数值模拟对球心非对称的钝圆头体附近流场气动特性进行分析研究。数值求解三维可压缩NS方程，采用NND格式，以修正的Baldwin-Lomax代数模型对不同飞行状态进行数值模拟，分析了球心非对称钝圆头体对流场的影响，与部分试验比较结果满意。     

喷管化学反应粘性流场的数值模拟

用时间相关法的显格式求解氢氧发动机喷管的粘性化学反应流动,采用６种组分,８个基元反应有限速率的化学反应模型,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流动参数和产物的质量分数在喷管流场中的分布。计算结果表明,采用显格式方法求解喷管化学反应粘性流动是可行的。     

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

SST湍流模型在高超声速绕流中的改进

为模拟高超声速湍流问题,对剪切应力输运(SST)湍流模型系数进行了修正。数值格式采用改进的总变差递减(TVD)格式,并对湍流模型的负值强制项进行了隐式处理。在此基础上计算了绕平板以及具有分离、再附、激波/边界层干扰等复杂流动结构的压缩拐角的高超声速流动。计算结果与试验数据及半经验公式的对比表明:SST湍流模型引入的雷诺剪切应力与湍动能之比为常数(Bradshaw数)在高超声速绕流中并不成立。Bradshaw数修正后的SST湍流模型与原模型相比,所计算的壁面压力、摩擦阻力和壁面热流分布更接近试验结果。     

Unsteady pressure fluctuations in the boattail region of a heat shield

Unsteady, turbulent, compressible, axisymmetric, Reynolds-averaged Navier–Stokes equations are solved for the flow past a bulbous payload shroud for a freestream Mach number of 0.95 and Reynolds number of 16.40×10⁶. A time-dependent numerical simulation is carried out using a finite-volume discretization technique in conjunction with a multistage Runge–Kutta time-stepping scheme. The closure of the system of equations is achieved using the Baldwin–Lomax turbulence model. Comparisons have been made with experimental results such as the schlieren picture and surface pressure distribution. A good agreement is found between them. A separated flow zone on the boattail region is observed and is found to be highly unsteady. Standard deviation, higher-order moments, histograms, and spectrum of surface pressure levels are analysed in the separated region of the boattail.     

三维机翼大迎角低速绕流及其涡控制的数值模拟

通过数值求解三维非定常可压缩雷诺平均Navier-Stokes方程，研究了三维机翼大迎角低速绕流及其带表面周期性吹吸气的流动控制，并探讨了非定常质量引射作为外激发手段对机翼气动力特性的影响，计算表明，在机翼前缘附近的小区域内施加周期性吹吸气进行流动控制，这类外激发手段对三维机翼的大迎角分离流动的控制是可行的，通过选择适当的控制参数，尤其是合适的周期性激振频率，可以有效地改善气动力特性。     

超音速进气道建模方法研究

对某可调混压进气道在不同攻角、不同马赫数、不同斜板角度下进行了大量的数值计算。给出了设计状态下进气道内外流场特征;分析了攻角变化对进气道流场的影响;以数值仿真结果为基础,利用B样条理论建立了反映攻角、马赫数及可调斜板角度变化的超声速进气道数学模型。根据此数学模型,分析了攻角和进口气流马赫数对进气道性能的影响,同时给出了斜板对进气道性能的影响。      

时间推进法求解叶轮机械内不可压流动

发展了一种用时间推进方法求解叶轮机械内部不可压流动的计算程序。在这个程序中, 使用分布体力方法简化湍流粘性的求解过程, 根据Ｂａｌｄｗ ｉｎ－Ｌｏａｍ ｘ 两层代数湍流模型确定粘性应力的大小。程序中还使用局部时间步长和多重网格数值技术加速收敛。火箭燃料主泵前诱导轮内部流场的实例计算表明, 采用这一计算程序能够正确预测叶轮机械内部不可压流场。     

高空台扩压器的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发,利用时间相关法,采用ＭａｃＣｏｒｍａｃｋ两步显格式,结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,数值模拟了高空台扩压器的内部流场,并计算了扩压器的性能。结果表明,采用求解Ｎ－Ｓ方程的方法计算高空台扩压器性能是一种有效的方法。     

混压式进气道与弹体一体化流场数值模拟

基于Favre平均的N－S方程和B／L代数湍流模型,采用Jameson格式和矩阵人工粘性,对X布局的混压式超声速进气道与弹体的一体化流场进行了数值模拟,研究了解决了数值计算中的附面层抽吸及混压式进气 启动问题,得出了进气道内外流场的马赫数分布和速度分布,讨论了绕弹体非均匀来流条件下,进气道的位置,附面层抽吸及攻角对进气道性能的影响和进气道外型形状对弹性气动力的影响,结果表明,对混压式进气道必须进行附面层抽吸,进气道的位置对进气道的性能有很大的影响,进气道的外型型面对弹体的气动性能有一定的影响。     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

有侧壁干扰的机翼半模型风洞实验的N-S方程数值模拟

应用三维可压、雷诺平均Navier-Stokes方程数值模拟了机翼半模实验风洞侧壁干扰和三维机翼半模与安装侧壁结合部流场。计算采用中心有限体积多步Runge-Kutta时间步长格式,将Baldwin-Lomax两层代数紊流模型扩展应用于三维拐角区流动。用本文方法计算了3个算例 ,并与国外的有关实验数据进行比较,计算和实验吻合良好。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

跨声速和超声速流中激波/边界层干扰数值模拟

对SST湍流模型中的Bradshaw常数a1进行了修正,并对跨声速和超声速流中激波/边界层干扰进行了数值模拟研究,空间离散采用二阶精度差值的低耗散通量分裂格式(LDFSS),时间离散采用对称高斯-赛德尔(SGS)算法。结果表明:在跨声速流动中,计算得到的壁面压力分布、分离区长度和速度剖面都与实验值吻合较好,而且很好地模拟了典型的λ激波结构;在超声速流动中,修正后模型的计算精度较原始模型有了较大改善,计算得到壁面压力分布和分离点的位置都和实验值吻合较好。     

叶轮机转子内流动的三维跨声速Navier-Stokes方程组数值计算及试验验证

应用有限体积法及Baldwin-Lomax代数湍流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛.采用隐式算法对对流通量进行求解,而耗散通量采用显式算法.采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算耗散通量.关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解.应用该算法模拟了转子Rotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等.在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验有较大的差别.