尾喷管内外超声速流场数值模拟

采用有限元法数值模拟带二次流可调收－扩喷管内外跨声速和超声速流场，在计算中采用四边形等参元确定位移函数，利用加权余量法中的Galerkin法建立有限元方程，数值研究三种型式喷管和四种飞行工况下尾喷管内外流场。计算结果与实验数据符合较好，说明可用此数值模拟技术帮助开展实验研究，进行尾喷管优化设计。     

前缘对进气道性能影响的数值模拟

应用质量平均的N-S方程和B-L代数湍流模型，在不同压缩面前缘半径条件下，以高超声速进气道的流场进行了数值模拟，对比研究了进气道压前缘半径对流场结构及性能的影响，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。     

尾喷管内部及其射流流场的数值模拟

通过求解三维Ｎ－Ｓ方程,对某型发动机尾喷管在亚临界及超临界状态下的性能进行了数值模拟,尤其对喷管射流流场及其对喷管内部流动的影响进行了研究。差分格式采用具有ＴＶＤ性质的三阶精度的Ｇｏｄｕｎｏｖ格式,射流场及内部流场简采用分区算法。数值试验表明,所采用的算法可以很好地处理来自于任意方向上的任意间断。     

燃料引射方式对超燃冲压燃烧室混合及燃烧效率的影响

应用含组分守恒方程的质量平均Navier Stokes方程和B L代数湍流模型,数值模拟了后台阶构型燃烧室在采用台阶上游支板引射和壁面垂直引射燃料时的内部流场。在计算过程中,对方程的对流项采用空间为二阶精度的TVD格式,扩散项则采用二阶中心差分离散。通过流场计算,对比研究了引射方式对燃料混合性能的影响。结果表明,台阶上游的支板在燃烧室的流场中产生了一对相对稳定的大尺度轴向旋涡,该旋涡不利于燃料的混合。采用壁面垂直引射时,在喷嘴下游的燃料流场中产生了小尺度轴向旋涡,该旋涡是提高燃料混合及燃烧效率的关键。     

支板构型超燃冲压燃烧室流场数值模拟

数值模拟了带凹腔的支板构型超燃冲压燃烧室内部的流场结构。在计算过程中，控制方程采用了含组分守恒方程，并与k-ω双方程湍流模型紧耦合的质量平均Navier—Stokes方程，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。采用在流体力学时间步内的当地积分法克服了非平衡源项的刚性问题，保持了LU-SSOR隐式求解算法的高效率。通过数值模拟，对比研究了凹腔构型的位置、长深比(S／H)对燃料混合性能的影响。结果表明，长深比大的凹腔增加了燃料在腔内的停留时间，相应地改善了燃料的混合。位于燃料喷嘴之后的凹腔比位于下游的凹腔更有助于提高燃料的混合效率。     

冲压燃烧室内的燃料扩散性能研究

应用含组分守恒方程的质量平均Navier-Stokes方程和B-L代数湍流模型,本文数值模拟了带支板引射器的模型燃烧室内部的流场结构,对比研究了支板引射器数量以及燃料注入角度对燃料混合性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的TVD和AUSMPW格式,扩散项则采用了二阶中心差分离散.采用流体力学时间步内的当地积分法克服了非平衡源项的刚性问题,保持了LU-SSOR隐式算法的高效率.研究结果表明,燃料有角度注入及采用多支板构型可提高燃料的混合性能.     

细长体大迎角湍流流场的数值模拟

通过数值方法对大迎角细长体超声速流场的模拟，探求正确模拟导弹大迎角绕流湍流流场的简单而有效的方法。数值模拟求解一般曲线坐标系下的三维可压缩Navier-Stokes方程，时间离散采用Euler向后差分，无粘项的空间离散采用二阶TVD格式，分别研究了Bald-win-Lomax代数模型及其修正形式(BLDS)对大迎角分离流动的模拟能力。数值试验表明修正的Baldwin-Lomax模型更精确地预测了流场的旋涡与分离情况，给出合理的表面压力分布，因而更准确地模拟了存在横向分离的导弹流场。     

带支板的冲压燃烧室的燃烧性能研究

应用含组分守恒方程的质量平均Navier-Stokes方程和B-L代数湍流模型,数值模拟了带支板引射器的模型燃烧室内部的流场结构,对比研究了燃料注入角度对燃烧及发动机性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的AUSMPW格式,对扩散项采用了二阶中心差分离散.研究结果表明,燃料有角度注入改善了燃料的混合性能,强化了燃料的燃烧效率,提高了发动机燃烧室的动量增量.此外,燃烧室的性能研究需考虑燃烧效率、总压损失及动量增量等因素.     

共用尾喷管多管脉冲爆震发动机数值模拟研究

采用基元反应模型和显式迎风TVD差分格式数值模拟方法,研究共用尾喷管多管脉冲爆震发动机(Pulse Detonation Engine,简称PDE)各爆震管之间相互影响.采用点隐算法解决化学反应引起的刚性问题,模拟了共用尾喷管内的流场,以及爆震波退化为激波遇到尾喷管收敛斜面反射及反射波向爆震管上游传播的过程.研究结果为多管PDE的共用尾喷管设计提供了理论基础.     

旋转爆震波发动机二维数值模拟

采用一种改进的化学非平衡流解耦方法对Euler反应流方程解耦处理,对流项采用五阶WENO格式离散,化学反应源项的刚性采用简化的隐式方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/Air预混气旋转爆震发动机内流场进行了二维数值模拟。模拟结果给出了不同发动机尺寸条件下的详细起爆过程,结果表明当发动机尺寸小于临界直径时无法成功起爆;详细分析了流场结构和爆震波形状,旋转爆震波的传播速度与理论预测值吻合;性能分析结果表明在喷注总压低于燃烧室平均压力时仍可实现推进剂喷注,没装尾喷管的情况下发动机比冲达到176.5s。     

湍流方程扩散项离散对某径向叶轮流场计算的影响

在相对旋转坐标系下采用Harten,Lax and van Leer contact （HLLC） 格式离散对流项,自行开发了基于多块结构化网格的有限体积程序,实现了对叶轮机械内部流场的数值求解.分别对半开式径向叶轮、闭式后弯叶轮展开数值模拟,程序和商业软件计算得到的不同叶高处表面压力数据,其相对差异不超过1%,验证了算法的正确性.针对湍流方程的扩散项,分别使用完全离散和略去交叉导数项离散,通过对湍流黏度等值线、气动轴向力和力矩的比较表明:在网格正交性较好的情况下,略去交叉导数项的离散对计算结果的影响小于1%,显著地减小湍流方程离散的计算量.      

叶片掠对跨声转子进口流动的影响

选择NASA Rotor 37跨声转子作为原型叶片,在叶片90%展高以上通过移动叶型基元的弦向位置构造叶尖掠特征,共有包括原型转子、2个前掠转子和2个后掠转子等5种不同叶型几何的算例.利用数值模拟软件计算得到5种转子的总体特性并获得三维流场.结果表明:叶尖前掠减小了叶片尖部进口的平均攻角,而后掠则使攻角增大.引入周向平均无黏动量方程,计算出由周向平均所产生的周向不均匀源项及径向平衡方程的各项数值,并对比其在5个转子进口的分布规律及其与攻角变化之间的关联.通过改变进口的流动平衡,叶尖前掠使轴向速度增大、周向速度减小,进而使得攻角减小;后掠则在叶尖起到了相反的效应.采用这种准三维方法研究了叶片掠对叶片进口流场的影响.      

气膜冷却涡轮导向叶片流场数值模拟

本文采用三维定常数值方法计算了某气膜冷却涡轮导向叶片流场。叶片的气膜冷却流动运用了源项模拟方法。将每排气膜孔用槽缝代替。计算的导叶能量损失系数和出口气流角与试验数据吻合良好。     

多管脉冲爆震发动机流场数值研究

针对二维带有收敛扩张喷管的双管脉冲爆震发动机(PDE),采用并行的时空守恒元和解元(CE/SE)方法对单管点火、双管同时点火和不同喉道面积情况下的爆震流场进行了数值模拟。化学反应模型采用H2和O2的8组分34基元反应。计算发现,单管点火情况下传入旁通爆震管的流场可以引爆其中的可燃气体;多管PDE内流场之间的相互作用使管内的压强处于比较高的状态;增大喉道面积削弱了多管PDE流场之间的相互影响,同时减小了封闭端面的压力,从而导致发动机性能下降。     

贴壁方柱湍流场TR—PIV实验研究

采用时变粒子图像速度场测试技术（TR—PIV），对低速循环水槽中贴壁二维方柱绕流湍流场进行了细致的测量。通过对实验得到的30000个连续瞬态速度场进行分析，得到了时均速度场和流线图谱，以及流向、法向速度分量的脉动强度场和涡量场。此外，对脉动速度时序信号进行谱分析，得到了流场中低频大尺度相干结构脱落频率。通过对瞬态速度场的分析，揭示了贴壁方柱绕流的旋涡脱落及其发展过程。     

轴对称弯道流动的N—S方程数值解

利用涡量—流函数方法求解离心压气机进气弯道中的轴对称粘性流场。本文推导出了q—τ湍流模型方程在正交曲线坐标系下的表达式,并用以预测湍流流场。文中在涡量方程及q—τ方程中引进非定常项,采用时间推进法求解。通过对二维直通道及复杂边界的轴对称弯道内的层流、湍流流动的计算,获得了令人满意的结果。      

近水面飞行气动特性数值计算方法研究

根据对水面波的数学描述，给出了其对流场影响的精确解析解。针对地(水)面效应飞行器的特点，分析了近水面波飞行时的流场特征，利用线性叠加原理，提出了飞行器与水面波耦合影响的非定常速度势修正项概念。分析了非定常气动特性的数值计算方法。对某实用型地效飞行器进行了计算，并对计算结果与试验结果进行了分析总结。