基于直角叉树切割网格的复杂外形N-S方程数值模拟

对于粘性绕流的数值模拟，在直角叉树切割网格的基础上，结合三角形非结构网格和矩形结构化网格，利用其各自的优势和特点，提出了一种生成混合染交网格的思路和方法，在物面附近生成了适合粘性流计算的大长宽比矩形网格，地远场分布直角叉树切割网格，快速离散计算空间。对于复杂的多体问题，采用三角形网格来连接 各体网格，并运用网格合并的方法，保证了各网格之间的光滑过渡过连接，提高了网格的质量，针对复杂多段翼型绕流问题，在上述网格基础上，采用B－L代数湍流模型和中心有限体积法，完成了Navier-Stokes方程（N－S方程）的数值模拟，计算结果表明，网格生成和流场计算都是正确的。     

多段翼型结冰数值模拟研究

文章采用计算流体力学方法对多段翼型进行了结冰数值模拟研究。应用多块网格技术生成多段翼型网格。引入水滴流场的概念,通过求解空气-水滴两相流场的欧拉方法,获得翼型表面水滴收集特性。流场采用有限体积法及LU-SGS隐式时间推进数值求解,运用流场分区求解算法,完成了多段翼型绕流流场及水滴流场的数值求解。建立了冰层的法向增长模型,对冰形进行了预测。     

基于两相流欧拉方法的翼型结冰数值模拟

文章通过计算流体力学的方法对翼型结冰进行了预测。针对拉格朗日法求解水滴收集特性计算工作量太大的缺点,引入水滴流场的概念,通过求解空气-水滴两相流场的欧拉方法,获得翼型表面水滴收集特性。结冰计算过程中网格采用变形网格技术得到,既保证了网格质量又减少了网格生成的工作量,采用有限体积法及LU-SGS隐式时间推进数值求解N-S及水滴运动方程。建立了对应于霜冰和光冰的冰增长模型。通过对NACA0012翼型霜冰和光冰的计算验证了该方法是有效的。     

非结构混合网格下的三维舵面流动计算研究

运用非结构混合网格上N S方程的求解技术对三维舵面的粘性绕流进行了计算研究。在物面粘性作用区,运用改进的推进层方法生成三棱柱、金字塔和四面体型混合网格;在其它流动区域采用阵面推进方法生成四面体网格。采用基于重构中心格式的格心有限体积法和SPALART-A LLM ARAS一方程湍流模型对N S方程进行空间离散,运用通量线化假设和最大特征值方法进行雅可比矩阵分裂来实现隐式高斯-赛德尔时间迭代。以某带舵面偏转的机翼模型为数值算例,系统地计算和分析了迎角、舵偏角及马赫数对舵面气动特性和舵面流动分离的影响,并将计算得到的舵面铰链力矩与实验结果进行了对比验证。数值算例表明:该方法可以为三维舵面气动力和铰链力矩的计算与分析提供有效的数值工具。     

翼型跨音速定常及非定常粘性绕流的数值模拟

发展了一种计算翼型跨音速定常及非定常粘性绕流的数值方法。采用ＬＵ分解隐式三阶迎风格式求解二维非定常可压ＮＳ方程组，紊流模型采用ＢａｌｄｗｉｎＬｏｍａｘ双层代数模型。计算网格采用相对翼型固定的贴体运动网格。数值模拟了ＮＡＣＡ００１２翼型的几种跨音速定常及非定常粘性绕流状态，计算结果同已有的实验数据吻合较好     

基于影响函数控制的结构动网格生成方法

在静态结构网格的基础上,根据流场边界的运动提出了一种新的结构网格变形技术.通过影响函数的作用,将运动边界的信息向周围的网格传递,流场中的网格节点通过影响函数获取移动变形信息,构建出新的网格.相对于弹性网格变形方法,影响函数控制方法具有更强的网格变形能力和更好的网格质量.该方法能够满足生成贴体的粘性动网格,适用于Navier-Stokes方程.耦合Navier-Stokes方程与强迫运动方程,数值模拟了NACA0012翼型做强迫俯仰振荡(AGARD-CT1)的粘性动态流场.计算结果与试验值的对比,表明了该方法的有效性.     

基于N-S方程的翼型结冰数值模拟

为了研究结冰对翼型气动性能的影响，采用计算流体力学的方法对NACA0012翼型前缘的结冰进行了预测。在流场计算中，对二维可压粘流的时均N-S方程进行离散求解；采用四步龙格一库塔法求解水滴轨迹运动方程；建立了对应于霜冰和光冰的冰增长模型并假定冰沿着与当地翼型表面法向一致的方向增长来预测结冰的形状。分析了网格数和冰层时间步长对最终冰形的影响和结冰对翼型气动性能的影响，最后将计算的结果与文献中的实验值和计算值进行比较，结果一致。     

结冰对翼型和多段翼型绕流及气动特性影响研究

对翼型结冰问题，在结冰条件下，利用数值方法分析研究翼型和多段翼型绕流流场及气动特性的影响。对于翼型和多段翼型，采用分区多块网格技术，结合法向外推的代数方法和求解椭圆型方程方法，生成高质量计算网格。采用中心有限体积法，以及Runge-Kutta显式时间推进格式，结合B—L代数湍流模型，运用流场分区求解算法，完成了绕流流场的N-S方程数值模拟，进一步针对3种不同的冰型：钝头体、双角体和尖头体冰型，分析不同形状的冰型对翼型和多段翼型绕流流场及气动特性的影响。     

翼型突体三维绕流场的数值模拟计算

采用雷诺平均ＮＳ方程为基本方程,应用二层κ－ε湍流模型封闭雷诺平均ＮＳ方程,并用有限体积法离散方程,对翼型突体的三维绕流场进行了数值模拟计算。不同于通常的网格划分,采用了贴体坐标与空度法相结合的方法对翼型突体的几何形状及其对流场的影响进行描述,计算结果与实验结果进行了比较,结果表明,所用方法,可以对翼型突体三维绕流场给出合理的描述。     

Navier-Stokes方程预处理方法及其对翼型绕流数值模拟的应用

在与定常方程相容的情况下,通过对N av ier-Stokes方程的时间导数项实施Cho i及Turkel矩阵预处理,发展了适用于低、亚、跨声速粘性流动的有效的预处理方法。对预处理后的控制方程采用了Jam eson格式进行有限体积空间离散和Runge-K utta显示时间推进求解,以及采用了FA S多重网格方法加速收敛。对RAE 2822、GAW-1等不同类型翼型进行了低速和跨声速流动的数值模拟。算例表明:Cho i预处理方法和Turkel预处理方法均有效改善了时间推进方法对低马赫数流动计算的收敛性,而且提高了计算精度;预处理方法同样适用于跨声速非线性流动的计算,且表现出了一定的加速收敛效果。     

带结构刚度非线性的超音速弹翼颤振分析方法研究

提出了带结构刚度非线性的超音速弹翼的颤振分析方法。以非定常可压三维N-S方程为基础，利用NND-LU混合差分格式，Baldwin-Lomax紊流模型，C-H贴体运动网格计算非定常气动力，并和变系数颤振方程耦合，用时间推进求解的方法计算带非线性结构刚度的翼面颤振特性。通过一个带结构刚度非线性的导弹全动翼面的算例，计算研究了3种常见的结构刚度非线性(间隙型，转折型，三次型)对颤振特性的影响。计算结果表明，间隙型的刚度非线性使颤振速度下降。转折型和三次型刚度非线性能提高额振速度。在一定的速度范围内带结构刚度非线性的颤振系统可能出现极限环振荡。     

一种改进的四边界插值网格生成方法及其应用

通过改进一种四边界插值方法，并引入提出的网格正交控制及加权平均光顺措施，给出了一种非常有效地代数网格生成方法。应用该方法生成的网格贴体性，疏密性及均匀性都较好，且网格生成非常迅速，应用该方法对各种复杂气运外形进行了较高质量网格生成，表明了此方法的有效，实用性。     

基于变弹翼方式的非结构动网格建模方法研究

针对变弹翼导弹包含运动边界的复杂非定常流场的求解问题,基于弹簧变形原理,即在网格变形时,网格点与静止边界和运动边界的相对位置关系均保持不变,分别推导了旋转弹翼运动边界上的网格点和受运动边界扰动流场子域上网格点的变形算法,从而发展了一种新的非结构动网格方法,该方法在不降低计算效率的前提下,大大提高了网格的变形能力,并解决了现有研究通常忽略弹翼运动受弹身影响而产生的非刚性变形问题,是基于非结构网格框架下,求解具有复杂运动边界且网格变形较大的非定常流场的一种简单有效的方法。     

回转体局部空泡流的数值分析

通过在回转体和空泡表面分布源汇的方法，建立了求解回转体局部空泡流场的数值方法。依据建立的数学模型与数值方法，对两种具有典型头部线型的回转体局部空泡绕流场进行了数值计算。在数值计算过程中发现第一次迭代所得结果和最终的收敛结果相比误差仅为5％，并把计算结果与有关献的实测数据进行了比较，发现计算结果与实测数据十分接近。     

复杂流场数值模拟中的网格生成

采用求解二维及三维椭圆型偏微分方程方法生成复杂外形的贴体及与边界正交网格。并采用分区,局部修正等手段提高网格质量。所得的网格用于导弹、全机、全机带外挂等我种外形的流场计算,结果与实验吻合良好,证明所生成的网格具有很高的质量。     

悬停旋翼粘性绕流N-S方程数值模拟

通过求妥薄层Navier-Stokes方程数值模拟了悬停旋翼粘性绕流,为了方便实施旋转对称性边界条件和更为有效地捕捉尾涡系,采用了嵌套网格技术,空间离散方法采用中心有限体积法,时间推进格式为五步Runge-Kutta法,考虑到旋翼绕流的特点,对远场边界条件作了合理的修正,应用本方法对Caradonna模型旋翼的两个实验状态进行了数值模拟,计算结果与实验值吻合很好,说明本 方法 成功的     

一种参数化曲面网格生成方法

针对在计算机辅助设计(CAD)及流场数值计算中,经常遇到绘制曲面及曲面网格的问题,探讨了一种生成参数化曲面网格的新方法.介绍的参数化曲面网格生成方法是依据代数变换和椭圆变换得来的.通常所用的二维和三维域椭圆变换是基于在Laplace方程,而本文生成的参数化曲面网格,则基于Laplace beltrami方程.由代数变换和椭圆变换构成的复合变换服从Possion方程的椭圆系统,其控制函数由代数变换确定,且遵守Neumann边界条件,这就保证了网格边界上的正交化.该方法生成的网格最终由曲面本身的形状及边界上的网格点分布确定,减少了计算量,尤其对大型复杂型面可有效提高网格生成速度.     

圆柱体绕流的数值模拟

采用有限元方法离散求解了雷诺平均Navier-Stokes方程。为了提高流场的分辨率,湍流的模拟采用了亚格子模型的大涡模拟模型,对绕单个圆柱体的绕流问题进行了数值模拟。给出了随时间变化圆柱体周围的流动速度矢量图、柱体上的非定常力系数值计算结果与相关实验结果的比较。数值计算结果与相关实验结果具有合理的一致性。