非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.     

二维非定常分离准则探讨

利用奇点理论和连续方程证明,当坐标系与非定常分离点固结时,分离点应该是鞍型奇点。由此自然地导出了非定常分离准则。进而还讨论了在非定常流场中如何确定分离点的几种方法     

评二维非定常分离准则

说明二维非定常与三维定常比拟准则存在着一些内在的矛盾,进而说明用Lagrange方法揭示的二维非定常分离应该用染色线而不应用轨线来解释。     

地效飞行器近波浪地面大迎角飞行分离流动数值研究

数值模拟了地效飞行器近波浪地面大迎角飞行的分离流场,研究飞行高度和地面波浪对气动性能的影响.使用Fluent软件求解非定常可压缩流动的质量加权平均NS方程和标准k-ε湍流模型,使用滑移网格处理地效飞行器和波浪地面间的相对运动.飞行高度降低,机翼上翼面的附着流动转化为分离流动.波浪地面飞行气动力平均值和水平地面飞行气动力随飞行高度变化的规律完全一致.     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

单边膨胀喷管内流动分离非定常特性

结合聚焦纹影和动态压力测量技术,对基于特征线法设计的单边膨胀喷管(SERN)不同落压比(NPR)条件下喷管内流场结构和壁面压力进行了试验测量,通过壁面压力时域和频域综合分析获得了喷管内流动分离非定常特性。结果表明:过膨胀状态下单边膨胀喷管内流场结构具有明显的非对称特征,喷管上壁面流动分离模态为受限激波分离(RSS),而下壁面流动分离模态为自由激波分离(FSS);相比于FSS模态,RSS模态下出口附近壁面压力振荡更剧烈。喷管上、下壁面压力标准差峰值均在分离点附近,且概率密度函数分布向一侧偏斜或出现双峰现象。RSS模态下,激波运动呈明显低频特性;FSS模态下,激波非定常特性不仅受回旋区压力扰动的影响,且受分离剪切层的影响。      

等离子体激励控制激波与边界层干扰流动分离数值研究

针对高超声速进气道激波与边界层干扰流动分离控制问题,提出了一种低功率重频非定常激励方式,并基于雷诺平均Navier-Stokes(N-S)方程,从唯象学的角度出发,将等离子激励简化为功率密度源项,对比研究了定常与低功率重频非定常等离子体气动激励的作用机理与控制效果。结果表明:定常激励的能量沉积作用对于激波控制非常有效,并可诱导出斜激波,但是对于流动分离控制而言,其能量沉积显然过于强大,反而会使流动分离更加严重,无法满足控制要求;当采用低功率重频非定常激励方式时,对于不同功率密度的情况均存在最佳激励时长与频率,当功率密度为5.0×109W/m3时,最大射流速度可以达到895m/s,并且可以在一定程度上减弱激波与边界层干扰流动分离。      

大负荷低压涡轮叶型分离转捩流动的大涡模拟

应用动力模式大涡模拟数值方法,对来流无扰动、不可压、雷诺数为5×104(基于进口速度和轴向弦长),定常来流条件下大负荷低压涡轮叶型(Pak B)叶型吸力面非定常分离转捩流动进行了三维数值模拟。在与相关实验数据的对比基础之上,对非定常流动物理信息进行了详细的分析讨论,揭示了计算来流状态下的Pak B叶型吸力面非定常分离转捩流动机理。结果表明,由无粘Kelvin-Helmholtz机制产生的、空间线性增长的初始二维不稳定性在分离剪切内诱导展向旋涡形成并脱落,脱落过程中的展向涡在非线性增长的三维不稳定性作用下发生变形并最终破碎成湍流。计算得到的Kelvin-Helmholtz不稳定性特征频率处于相关实验测量范围内。     

纳秒等离子体激励控制翼型流动分离机理研究

为研究纳秒介质阻挡放电(NSDBD)等离子体控制翼型流动分离的物理机理,采用已建立的NSDBD唯象学模型耦合非定常Navier-Stokes方程模拟纳秒等离子体对流场的作用。使用非定常雷诺平均NavierStokes方程(URANS)和大涡模拟(LES)两种求解方法,研究纳秒等离子体激励对NACA0015翼型流动分离控制。结果表明:NSDBD等离子体激励促使边界层提前转捩,转捩对控制流动分离起重要作用;NSDBD激励开始时在翼型前缘形成展向涡,展向涡促使分离剪切层失稳并最终进入尾迹,展向涡贴近壁面运动,将外区的高能气流带入近壁区,使上翼面流场结构发生变化,然后翼型前缘流动提前转捩促使流动经过一个小层流分离泡后发生湍流再附,最终在上翼面形成稳定的附着流动。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

TIME ACCURATE 2D NAVIER-STOKES CALCULATIONS WITH A DUAL-TIME STEPPING METHOD

ＴＩＭＥＡＣＣＵＲＡＴＥ２ＤＮＡＶＩＥＲＳＴＯＫＥＳＣＡＬＣＵＬＡＴＩＯＮＳＷＩＴＨＡＤＵＡＬＴＩＭＥＳＴＥＰＰＩＮＧＭＥＴＨＯＤＹａｎｇＹｏｎｇ，ＱｉａｏＺｈｉｄｅ（Ｆａｃｕｌｔｙ５０３，ＮｏｒｔｈｗｅｓｔｅｒｎＰｏｌｙｔｅｃｈｎｉｃａｌＵｎｉ...     

弹体穿越平面激波三维非定常流场的数值模拟

用数值方法求解三维非定常可压缩Navier-Stokes方程，计算了复杂弹体在相同马赫数、不同迎角下的流场和气动力特征，计算结果与实验结果吻合，证明本文采用的方法能够准确计算此类复杂流场；然后采用该方法模拟了弹体穿越正激波的非定常过程，得到了弹体周围流场及气动力特征随穿越过程的非定常变化，阻力特征穿越后急剧下降。     

球柱三维分离及涡流的数值模拟

对三维可压非定常Ｎ－Ｓ方程进行有限体积中心离散,并用Ｒｕｎｇｅ－Ｋｕｔａ显式格式进行数值求解。数值模拟了球柱分离及涡流的形态结构,给出了具有激波、边界层、大分离区、涡流及其间相互干扰的复杂流型     

细干扰线对后向台阶分离流的作用

应用氢泡流动显示技术研究了后向台阶层流边界层分离、再附整个过程中流动结构的发展变化。内容涉及基本流动和非定常控制。非定常控制通过在台阶上游不同位置布置的细干扰线产生的尾涡来实现。流动显示结果表明细干扰线在大多数情况下延缓了剪切层的再附,只有当xw/h=0,yw/h=0.5时可使再附距离缩短。     

吸气条件对圆柱非定常分离流影响的数值研究

本文根据二维可压非定常Ｎ－Ｓ方程，基于Ｂｅａｍ ａｎｄ Ｗａｒｍｉｇｎ及Ｓｔｅｇｅｒ的方法，在圆柱绕流出现周期分离后，施加吸气条件，研究气对于圆柱非定常绕流的影响，计算结果表明，在恰当地区域施加恰当的吸气条件，可以希望通过减小压差阻力来减小总阻力，吸气还有抑制尾迹中非对称流动的作用，从而使圆柱所受的垂直于来流方向的升力的变化幅度减小，变化频率减小，但是吸气会增加磨擦阻力，过份的吸气会导致摩阻增大太     

座舱盖抛放过程中的非定常气动特性研究

通过对N-S方程的数值求解模拟了座舱盖在静止和旋转两种状态下的气动力分布,分析了来流马赫数、旋转角速度以及迎角对气动系数的影响,研究了座舱盖自由抛放过程的非定常效应,并最终建立了预测座舱盖自由抛放过程气动系数的气动系数图。计算表明:(1)非定常效应是影响抛放过程座舱盖气动力分布的主要因素;(2)外形左右对称的座舱盖,均有非定常气动系数平均值与无量纲角速度成线性关系的规律,而且变化斜率几乎与舱盖形状尺寸等无关。