转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

现代自然层流翼型的设计方法

本文介绍一种设计跨声速自然层流翼型的计算流体力学(CFD)方法。本方法采用“正反迭代、余量修正”设计原理,通过将跨声速翼型设计软件NPU-TD2D中的反设计程序进行改进,并与含有层、湍流混合边界层修正的跨声速层流翼型计算程序DLRBGKWALZ耦合,实现了在跨声速粘性流动条件下直接设计层流翼型。亚、超临界的设计实例和风洞验证表明,本方法可以在几个设计迭代内设计出压力分布、转捩位置及气动参数均准确收敛于设计目标的新翼型,是一种设计现代自然层流翼型的有效而实用的CFD方法。     

高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。     

离心压气机中进口导叶尾迹的形态、输运及与叶轮流场的相互影响研究

对跨声速离心压气机中进口导叶与叶轮的相互影响作用进行了数值模拟研究,研究模型包括三种几何间距模型及对同一间距几何模型使用了两种不同的转静子交界面位置设定。计算结果表明：转静子交界面位置的设定对于多排叶片混合平面法定常计算性能会产生较大影响,5%叶轮轴向长度的变化就可造成1%的效率、2%的压比差别;进口导叶与叶轮的叶排间距大小对离心压气机级气动性能影响很小,但间距较小时会造成流场内较大的压力波动。在近叶排间距时,导叶切割叶轮激波会在导叶压力面侧产生一个高损失区,并沿导叶表面向上游传播;该高损失区由激波压力波造成。此外,激波存在与否,影响到两排叶片流场的相互影响作用强度：激波使叶轮对导叶流场的影响几倍甚至十几倍大于导叶对叶轮流场的影响;激波的影响作用使导叶压力面、吸力面侧压力波动由不同的原因造成。最后,在近叶排间距时,初始进入叶轮通道内且靠近主叶片压力面的尾迹相对靠近吸力面的尾迹片会滞后;完全进入叶轮通道后,近压力面侧尾迹的滞后更显著。     

跨声速风扇叶片变形对气动性能的影响

为了研究跨声速风扇叶片变形对气动性能的影响,发展了流固耦合方法计算风扇叶片在气动力作用下变形的数值程序,耦合求解了三维可压缩雷诺平均的Navier-Stokes (N-S)方程和描述叶片振动的结构动力学方程,并模拟了跨声速风扇Rotor 67转子叶片在气动载荷下的变形,重点分析了由于叶片变形带来的气动性能的变化.      

固体火箭发动机燃烧室冷态突扩气固多相流动的数值模拟

基于欧拉-欧拉双流体方法,建立固体火箭发动机燃烧室轴对称突扩气固多相流动的统一二阶矩两相湍流模型,采用k-ε双方程模拟气相湍流流动和颗粒相湍流雷诺应力kpg方程描述大尺度颗粒湍流运动、小尺度颗粒各向异性弥散特性以及气固相间作用。数值模拟了突扩流动中的颗粒轴向、径向平均速度和脉动速度以及颗粒-气相脉动速度关联轴向分量的分布特性,计算结果与试验吻合较好。研究结果表明:颗粒轴向速度脉动大于径向速度脉动,轴向平均速度梯度影响脉动速度峰值的发生位置和大小;两相脉动速度关联轴向分量大约是径向分量的2倍。     

Revisiting the space-time gradient method: A time-clocking perspective,high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternativ...     

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.