分层Reynolds应力模式

本文根据切变湍流结构的特征提出一种分层Reynolds应力模式。在切变湍流中存在湍能生成项达到狭峰的极大值层次(本文称它为活跃层),根据这一特征,在活跃层中保留微分Reynolds应力模式,但它可以简化为常微分方程。在两层活跃层之间以及活跃层和壁面之间的Reynolds应力分布,用满足衔接相容条件的多项式近似。时均湍流场将由以上得到的Reynolds应力代入Reynolds方程求出。用本方法计算了两个比较复杂的湍流例子,结果表明,本方法较常用的k-ε模式有明显的优越性,它的计算时间并不比k-ε模式多,而计算结果与实测结果的吻合则较k-ε模式好得多。另一方面本方法较之完全的微分Reynolds应力模式可大大节省计算时间。     

轴对称喷管流场分析的有限体积法

本应用有限体积法对轴对称喷管的跨音速流场进行了分析，采用显式预校两步MacCormack格式进行计算。运用隐式残差平均技术及当地时间步长来加速计算的收敛，提高格式的稳定性。对三个不同几何参数的轴对称喷管的数值分析表明，计算结果与试验数据吻合良好。     

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。     

二元喷管跨音速流场的数值计算

本文采用有限体积法进行空间离散,Runge-Kutta法进行时间推进的办法,通过求解Euler方程来模拟二元收敛—扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡,在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解。计算结果与试验数据和其他数值方法吻合良好。      

来流参数对防热瓦横缝旋涡结构及热环境的影响

针对高超声速飞行器表面缝隙内部流动,通过求解可压缩Navier-Stokes方程,自主研发了一套能够较好模拟缝隙流动特性的计算流体力学(CFD)软件。利用该软件研究了来流参数对防热瓦横缝旋涡结构及热环境的影响。计算结果表明:随着来流雷诺数的增加,缝内旋涡结构呈现主涡个数增多形态趋于饱满的变化趋势,缝隙壁面绝对热流和无量纲热流增加;随着来流马赫数的增加,缝内主涡个数、形态基本不变,但主涡旋转速度增加,缝隙壁面绝对热流增加,无量纲热流基本不变;随着来流迎角的增加(迎角较小时),缝内旋涡结构和热流变化规律基本与增加来流雷诺数相同。由此分析可知,涡量向下传递并形成旋涡的距离,即形成所谓"死水区"的深度,主要由来流雷诺数和来流迎角决定。     

平衡空气模型的流动稳定性及转捩预测

研究高超声速平板边界层考虑真实气体效应的流动稳定性问题.采用7组元化学反应平衡模型,黏度和导热系数采用混合律,同时考虑组元浓度扩散引起的能量传递,在马赫数为10~20、壁面温度为500~3500K、飞行高度为20~30.5km等条件下,对平板边界层流动的稳定性进行了分析,给出了扰动演化相对增长的N值.计算结果表明:高马赫数飞行中不稳定扰动的第3模态将与第2模态合并,共同影响转捩;高温真实气体的流动稳定性特征,随着马赫数、壁面温度、飞行高度变化的基本趋势与完全气体的基本一致;与完全气体相比,真实气体的相对增长N值包络线较小,表明高温真实气体将抑制转捩发生.      

单级环境静叶等离子体流动控制机理仿真研究

在轴流压气机等离子体扩稳研究中,针对单转子压气机流动控制的研究较多,而针对单级环境下静叶流动控制的研究却很少.采用静叶轮毂轴向等离子体激励方式,通过数值模拟方法研究单级环境下静叶流场特性,揭示轴流压气机静叶等离子体流动控制扩稳机理.结果表明:等离子体激励器的轴向位置对单级轴流压气机的扩稳效果影响显著,越靠近叶片前缘,扩稳效果越好;布置在静叶通道后半部的等离子体激励器无法提高压气机的稳定性,而在静叶前缘施加轴向等离子体激励时,近轮毂区气流被诱导加速,主流的轴向速度提高,有效抑制了静叶近轮毂区吸力面的流动分离,静叶近轮毂区的堵塞减小,使得单级轴流压气机的稳定性提高.     

纳秒脉冲等离子体分离流控制频率优化及涡运动过程分析

将纳秒脉冲驱动的介质阻挡放电等离子体激励器应用到NASA SC（2）-0712翼型上,在迎角分别为15°和20°时,开展了在不同雷诺数下的分离流动控制研究。通过模型表面静压测量,得到了不同激励频率下的分离流动控制效果。对翼型表面压力进行分布积分,得到了在不同雷诺数和激励频率下的升力系数,表明分离流的控制效果有一个较宽的激励频率范围,只要激励频率落在相应的频带范围内,均能实现有效的分离抑制。流动显示结果表明,分离流的控制在瞬时表现为放电后可形成大尺度旋涡拟序结构。旋涡的周期性产生、运动和演化造成了分离剪切流动的动态变化过程,从而促进了高/低速气流的动态掺混。     

远场假设对喷流噪声预测中格林函数求解的影响

目前基于雷诺平均Navier-Stokes（RANS）的喷流噪声预测方法在格林函数求解时, 为简化求解过程,通常对喷流流动做平行流假设, 对观测点做远场假设。随着格林函数求解方法发展,近年来的研究表明平行流假设对下游观测点格林函数的计算会引起较大偏差,而目前远场假设对格林函数求解的影响仍不清楚。为研究远场假设对喷流格林函数求解的影响,以二维喷流为例, 采用计算气动声学方法（CAA）分别数值求解了观测点远场假设条件与实际条件下90°~150°方向喷流内伴随格林函数,进而分析远场假设对格林函数求解的影响。研究结果表明,对于不同方向的观测点,由观测远场假设导致的伴随格林函数求解偏差不尽相同,且对于越靠近喷流中心线方向的观测点,远场假设导致的偏差越大,其中150°方向观测点,采用远场假设后,格林函数计算结果最大偏差达到-15 dB以上。因此,对于靠近喷流中心线方向的噪声观测点而言,为避免预测偏差,应采用实际观测条件求解喷流格林函数。     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。