Navier—Stokes方程粘性项的有限体积离散法研究

Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程粘性项的有限体积离散法研究封建湖西北工业大学应用数学系［摘要］：三维Ｎ－Ｓ方程是描述真实气体流动的非线性偏微分方程组，如何快速准确地计算其粘性项至关重要。这不仅因为粘性项的计算量非常巨大，而且粘性项对急剧变化的流动梯度的...     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

湍流脉动作用理论在轴对称流动计算中的应用

本文在湍流压力脉动远距离作用理论及其平面问题计算基础上,将其应用于轴对称流动的计算,并与轴对称湍流射流的实验结果进行了比较,表明轴对称流动计算比用平面计算更与实验一致。     

高超声速流动转捩的数值研究

基于线性稳定性的结果，以类似于处理湍流脉动的方式来模化流场脉动，用以研究高超声速的流动转捩，包括转捩起始位置及转捩区长度．研究中采用Faver平均Navier—Stokes方程k-ε湍流模型，分别以Roe二阶流通昔差分分裂格式和Chakravarthy三阶迎风TVD格式进行离散，数值模拟了椭圆锥的高超声速绕流流场。分析比较了不同的湍流转捩机制，并研究了雷诺数的影响规律。结果表明：Bypass转捩机制是传统的高湍流度风洞中模型流动转捩的主控因素。     

一种姿态控制发动机扩散段流场的数值模拟

采用CFD方法对一类火箭弹姿态控制发动机扩散段流场进行了数值模拟，并对几种典型流场进行了分析。计算中的控制方程使用雷诺平均Navier-Stokes方程，湍流模式选用Realizable k-g端流模式，近壁区用二层分区模型处理。通过与试验结果的比较，证明了在较大的主气源压力范围内计算结果合理，其推力计算的精度可满足火箭弹姿控发动机工程计算的要求。     

套室内气流的三维数值模拟

借助流场计算软件PHOENICS，对单元住宅室内的气流在夏季工况条件下进行三维模拟计算。考虑了外界气温和太阳辐射等因素对室内气流的影响，得到室内气流的速度矢量和温度等参数的分布。依据供暖、通风和空气调节在夏季工况下的标准，对室内的气流温度、气流速度做出评判。通过算例验证了数值模拟的可信性。     

较低湍流度范围湍流度对风洞实验结果的影响

变湍流变试验结果表明，湍流度对诸如洞壁边界层、平板边界层，各类翼型边界层和尾流的时均特性，湍流特性，转捩情况以及对翼型和旋成体压力分布等气动特性有显著影响。这种影响有一定规律，量值较大，不能忽视。     

受限共轴射流的湍流混合过程数值模拟

本文以下Ｓｔｉｒｌｉｎｇ发动机引射管设计为研究对象，对亚声速受限共轴射流中的两种气体混合过程，用数值模拟方法求解了轴对称雷诺平均Ｎ－Ｓ方程组。提出的一种改进的显式有限体积法，包括基于网格角点的具有高分辨的中心／迎风ＱＵＩＣＫ算法，四阶人工耗散、局部时间步长和隐式余数光滑技术的应用。为加速收敛，采用一种类似于总焓阻尼的方法，引入质量流量阻尼和气流中化学组份质量分数阻尼，以快速获得足够近似的解，适合工     

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

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

共轴异向旋转筒间不可压流的临界主不稳定性

基于长隙比无穷大和Tollmien-Schilchting 行谐波假设，通过导出坐标系下不可压一阶改型边疆特征值问题，实现了异向旋转筒流动临界主不稳定性的线性分析，在稳定性方程离散中引入了一类高精度对称紧致差分格式，在特征值搜索中发展了一个高效率复矩阵双重反迭代算法。双临界点分析揭示了几何参数、、流动参数、流动参数、扰动参数相互间的变化规律。理论结果同实验测量数据基本一致。