用k-ω紊流模型和Osher格式在混合网格上计算紊流流场

用有限体积法求解了雷诺平均 Navier-Stokes方程。空间离散采用了 Osher逆风格式 ,紊流模型采用了两方程模型 ,并对绕 NACA0 0 1 2翼型及双椭球外型的紊流流动进行了数值模拟。结果表明用上述方法求解紊流流动非常有效     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

分区对接网格下多段翼型N-S方程计算

对多段翼型的复杂粘性流场，各翼段物面及尾迹区生成C型网格．外场采用O型网格，中间用H型附加区过渡．自动生成高效的分区对接C-O-H型网格。采用时间精确有限体积高精度格式解可压缩N—S方程。层流、紊流的N—S方程求解结果，说明了本方法的有效性。     

壁面冷却的三股混合气流的数值计算

对双层壁扩压器与波瓣喷管组合的三股气流的引射混合进行了流场温度场的数值计算。计算域划分为各子域 ,通过求解拉普拉斯方程 ,分别生成三维贴体曲线坐标网格 ,边界网格加密且波瓣处正交 ,各子域组合成整个网格。为避免因波瓣造成的网格强烈的非正交而引起解的发散 ,本文采用了同位网格 SIMPLEC计算方法和Chen-Kim修正的 k-ε湍流模型。固体区域采用大粘性解决流固耦合 ,还结合了线性欠松弛和拟瞬态欠松弛。结果证明 :双层壁间有外界冷气流被引射进入 ,形成壁面的冷却气流 ,相对单层壁扩压器 ,双层壁扩压器的壁面温度明显降低。对于小扩压角 ( 0～ 5°) ,随着扩压角的增大 ,壁面平均温度降低 ;但当扩压角较大时 ( 5～ 1 0°) ,壁面平均温度则增加     

用预处理方法及可压Navier-Stokes方程计算低马赫数流动

用可压方法及预处理方法求解了低马赫数紊流流动,其中空间离散格式采用了Roe格式,紊流模型采用了DES模型.以绕NACA0012的紊流流动为例,分析了该方法对提高流动速度和精度的特性,并将其用到多段翼型及绕圆柱的非定常流动中.数值模拟结果表明,该方法能有效求解低速流动及混合流场问题,且具有较高的数值精度.     

Numerical Simulation of Aerodynamic Characteristics of Helicopter Rotor with Tip Slots

Effects of tip slots on the aerodynamic characteristics of helicopter rotor were investigated numerically by solving three-dimensional Navier-Stokes equations based on unstructured overset grids algorithm.Improved delayed detached eddy simulation (IDDES) based on the Spalart-Allmaras turbulence model and adaptive grid refinement technique were employed.Several slots in the rotor blade tip were designed on the base of Caradonna-Tung rotor to study the effect of tip slots.Numerical results show that tip slots are able to introduce the airflow from the leading edge and turn it in the spanwise direction to be ejected out of the face at the rotor blade tip,which can reduce the strength of the rotor blade tip vortex and accelerate the dissipation process.Although tip slots may lead to the decrease of airfoils' lift coefficient at the root of the rotor blade,it can increase the lift coefficient of airfoils at the rotor blade tip,so the lift of the rotor with tip slots is almost the same as that of the rotor without tip slots.In addition,tip slots can also reduce the intensity of the tip shock wave,which is beneficial to reduce the wave drag of the rotor.     

Y形进气道流场数值研究

Y形进气道在两侧进气布局的单发轻型战斗机和教练机上应用较为普遍。由于流动的复杂．不易对其内部的流动进行实测．为此本文基于Favre平均的N-S方程和B／L代数湍流模型．采用了Jameson格式和矩阵人工粘性耗散，对某型亚声速飞机的Y形进气道流场进行了数值模拟。提出了一种解决网格块之间插值问题的处理方法．研究了进气道两侧管道气流的掺混作用及二次涡．得出了进气道内外流场的马赫数、压力和速度分布，以及进气道的性能随反压和攻角的变化曲线。研究结果可为进气道的设计提供一定的参考。     

盒式机翼布局气动特性研究

采用数值模拟的方法研究了盒式机翼布局的气动特性。在数值模拟过程中,控制方程为雷诺平均N-S方程,采用基于非结构网格的有限体积方法进行求解。空间离散采用了Osher逆风格式,时间离散则采用五步Runge-Kutta格式。紊流模型采用了S-A方程模型。对盒式机翼布局外形在不同状态下的流场进行了数值模拟,给出了每个状态下的气动参数,并进行了分析,结果显示盒式机翼布局具有较好的气动特性。     

任意回转面叶栅C网格的数值生成

首先尝试采用任意曲面椭圆方程，数值生成任意回转面上叶栅Ｃ网格，经过对边界条件的仔细处理，生成的网格可保证在周期外边界和重合内边界上网格线导数连续，在叶片表面网格线相互垂直，提出了结合初始网格线分布，控制环绕叶片网格线分布的方法。     

不对称喷管欠膨胀超声速射流的数值模拟

采用改进的MUSCL格式解三维可压缩平均雷诺纳维尔 斯托克斯方程组,湍流模型为Spalart Allmaras代数模型,数值模拟了四种调整片不对称喷管在喷口压比为2.0、2.8、3.4和4.0条件下的欠膨胀超声速射流场,获得了射流场流谱和参数分布,并给出了单调整片不对称喷管出口压强比与射流偏转角和扩张角的分布曲线,计算结果和实验数据符合良好。     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

Highly-Efficient Aerodynamic Optimal Design of Rotor Airfoil Using Viscous Adjoint Method

In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.     

大展弦比机翼模型设计对翼型流场气动特性的影响

采用SST两方程湍流模型,通过求解非定常Navier-Stokes方程,模拟了大展弦比机翼风洞模型振动条件下的翼型流场,总结了翼型不同振动状况下的流场和气动力特点,分析了模型设计中的不同振动情况对风洞试验结果的影响.研究结果表明:在大展弦比机翼风洞模型的设计中,将翼型的重心设计在机翼的弹性轴之后,对风洞试验的精度较为有利.此结论对大展弦比机翼的风洞实验模型设计有指导意义.     

飞翼无人机低速着陆状态抖振响应研究

基于RANS方程和SST湍流模型,采用飞翼无人机升力系数判据、俯仰力矩判据和表面极限流判据等对抖振始发迎角度进行了预测估计。基于较为详细的结构模型和气动模型,构造了结构与气动的耦合求解技术,同时采用弹簧近似光滑和局部重构组合方法对气动的动网格技术进行更新;然后分别在时频域内分析了飞翼无人机的抖振载荷响应。研究结果表明:基于CFD的RANS方程、SST湍流模型及各类判据预计刚性飞翼无人机在低速大迎角状态下的抖振始发迎角可以得到较为合理的结果;采用升力系数判据、俯仰力矩系数判据预测的抖振始发迎角与表面极限流判据预测的始发迎角相比要保守一些;而表面极限流判据方法能给出翼面气流流动的细节,采用此判据可以分析对比不同迎角下的翼面气流流动变化情况;与气动弹性及嗡鸣响应相比,抖振是一种强迫振动,其响应频率并不单一。     

用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。     

平面COUETTE流动的湍流斑研究

湍流斑是近壁剪切流动中的重要现象,它们的生成和发展与流动转捩和湍流的形成密切相关.采用Navier-Stokes方程直接数值模拟和壁面脉冲模型研究平面Cou-ette流动中湍流斑的生成和演化.高精度高分辨率的三维耦合紧致差分格式被用于数值计算.探讨了平面Couette流动中湍流斑的重要特征,包括高频脉动的形成、雷诺应力的产生、扰动的急剧增长和湍流斑形状的变化,特别是流向涡的复杂演化过程.结果表明层流流动中的湍流斑具有湍流流动中一些基本特征.     

基于格子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时通过加密网格得到了更好的结果.为了避免转捩问题采用了无粘边界条件进行模拟,但得到的压力系数与实验结果吻合良好,证明了该方法的可靠性并显示了其模拟可压缩的潜力.该方法在翼型绕流上的应用为湍流流动的数值模拟提供了基础.     

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

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

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格,该方法的特点是在生成单元的同时引入结点,与其他预布结点的方法相比,其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构,有效地提高了程序运行效果,最后给出了Ｍ６机翼的数值算例。     

桨尖下反对旋翼BVI噪声特性的抑制分析

基于计算流体力学(Computational fluid dynamics,CFD)方法,分析研究了新型下反桨尖对前飞状态旋翼桨/涡干扰(Blade-vortex interaction,BVI)噪声特性的影响。首先,通过求解Navier-Stokes方程获得准确的噪声声源信息,湍流模型采用一方程S-A模型。采用双时间法进行时间推进,为了提高流场的收敛速度,在伪时间方向上使用高效的隐式LU-SGS格式推进,并采用并行算法进行加速。然后,在获得的可靠声源信息基础上,采用基于FW-H方程的Farassat 1A公式求解噪声。对有试验结果的算例进行了对比计算分析,验证了本文噪声预测方法的有效性。在此基础上,针对具有典型BVI噪声特征的前飞斜下降状态,开展了不同下反角度桨尖新型旋翼噪声辐射特性的计算分析,通过噪声辐射球的对比结果表明,选择适当的下反桨尖,可以有效地降低前飞斜下降状态下旋翼BVI噪声,从而达到较好的旋翼噪声抑制效果。