浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

基于滑移网格的螺旋桨滑流影响研究

基于滑移网格技术,通过采用准定常N—S方法,分析研究螺旋桨滑流对全机气动特性的干扰影响。以某带四台螺旋桨发动机的巡航构型为例,对比研究滑移交界面的影响,确立准确的滑移交界面建立方式,进而重点对某带后缘襟翼的螺旋桨飞机的流场特性进行计算分析,研究螺旋桨滑流对全机流场和气动性能的影响。数值计算表明：提出的用于研究螺旋桨滑流影响效应的滑移网格方法,是螺旋桨飞机设计评估的一种有效数值模拟手段,可简化动静网格的数据交换,具有一定的工程实用价值。     

静止可燃气热射流爆震起爆块结构自适应网格加密数值模拟

为了深入认识理解静止可燃气热射流爆震起爆的机理过程,采用块结构自适应网格加密方法的开源程序AMROC,进行二次开发,在LINUX系统中建立小型集群并行计算系统,应用于静止可燃气热射流爆震起爆过程高效精细数值模拟,并根据流场特性划分三个阶段进行分析。结果表明:对于一定条件下的热射流,静止可燃气中的热射流爆震起爆对壁面的依赖性很强,正是通过壁面的反射作用形成局部爆震燃烧;通过左壁面和上壁面的有效二次壁面反射,上壁面处的局部爆震燃烧进一步增强,强化的局部爆震燃烧逐渐向整个等直管道中传播;上壁面和下壁面的局部爆震燃烧同时形成之后,双三波点结构在管壁中来回碰撞碰,形成周期性的稳定爆震燃烧。     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

改进的弹簧比拟非结构化网格优化方法

提出了改进的弹簧比拟非结构化网格优化方法.分别从网格的稀密和单元质量两个方面对非结构化网格进行优化.首先采用非等刚度弹簧,实现了弹簧比拟方法对网格疏密的控制;其次,分别采用修正刚度系数的顶点弹簧模型和修正初始平衡长度的棱边弹簧模型,改善了原弹簧比拟光顺算法中存在的局部网格单元过渡不够平滑的问题.该方法保留了原弹簧比拟法易于实现、数据结构简单的优点,提高了原算法对网格优化能力,算例中最短边与最长边之比从0.24提高到了0.42(越接近1越好),且对于存在交错单元的初始网格也能优化.      

气助式直接喷射喷油器动态喷雾场数值模拟及试验

为了研究小型飞机发动机气助式直接喷射(air-assisted direct injection,AADI)喷油器喷雾形态及发展规律,采用数值模拟和高速摄影技术对其动态喷雾场进行研究.以锥阀运动的物理学模型为数值计算的移动边界条件,结合动网格技术对AADI喷油器动态喷雾场的速度分布、静动压分布和不同时刻雾化颗粒的空间分布进行动态数值模拟分析.以0.2ms为步长,对3.0ms内不同锥阀开启时刻的喷雾形态、雾化颗粒和贯穿距离进行了试验和仿真对比分析.结果表明:采用CFD仿真模拟技术可对AADI喷油器动态喷雾场的速度分布、压力分布、喷雾形态等参数进行定性分析,有利于揭示喷雾的发展规律,同时通过对喷雾粒径和贯穿距离的定量对比分析,表明了仿真计算的精确性,数值模拟计算与试验结果误差小于5%,满足工程应用的要求.      

模拟存在运动激波流场的无网格自适应算法

将无网格自适应算法发展用于求解存在运动激波的非定常流动问题。为了降低初始布点要求,同时提高对运动激波的分辨率,提出了一种结合点云加密和点云粗化的动态点云技术。采用该技术,在非定常计算的每一个时间步,首先通过基于压强梯度的流场探针识别出激波所在位置,然后对激波附近点云进行加密,从而提高对激波的分辨率,同时对远离激波的点云进行粗化,尽可能减少冗余节点,提高计算效率。采用所提出的无网格自适应算法求解了N-S方程,首先对一维运动激波问题进行了数值模拟,通过将数值模拟结果与精确解比较,验证了所提算法求解存在运动激波的非定常流动的可行性,在此基础上,将方法推广用于二维非定常流动问题的求解,分别对瞬时激波碰撞圆柱以及双弧翼型激波振荡流动进行了数值模拟,计算结果表明,方法可以显著提高对激波的分辨率,同时相对于直接采用较密布点的非自适应算法,自适应算法在计算效率上具有优势。     

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

水下燃气-蒸汽弹射气-液两相流场数值研究

基于计算流体动力学方法和均质多相流理论,采用汽化模型模拟燃气与冷却水的汽化过程,使用动态网格分层技术模拟导弹的运动,对燃气-蒸汽弹射气-液两相流场进行三维非定常数值研究.通过与固体火箭发动机尾焰喷水效应实验对比,验证了数值方法的有效性.研究了弹射过程中流场结构、二次流现象和变深度弹射载荷和内弹道变化规律.结果表明:导弹弹射过程中汽化现象发生在弯管壁面处,汽化后的燃气-蒸汽高温区随着发射筒体积的增加发生偏移;燃气-蒸汽在发射筒和弯管截面上存在隔离区和二次流动现象,并且形成不同结构形式的漩涡;弯管受到的最大冲击力出现在0.16s附近和30°~ 60°弯曲角范围内;弹射深度每增加10m,观测点压力增加0.25MPa,弹动时间延迟0.01s.