细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

预估-校正LU-SGS的隐式算法

为了进一步提高Euler方程求解效率,在传统LU-SGS格式基础上提出了预估-校正LU-SGS。结合高分辨率迎风通量格式,发展了一套高效的Euler方程求解器。以NACA0012翼型和RAE2822翼型跨音速无粘流动作为算例,对比分析了预估-校正LU-SGS、传统LU-SGS和GMRES格式的计算效率。结果表明:预估-校正LU-SGS算法能显著提高求解效率,其效率接近GMRES算法,比传统LU-SGS方法提高了近三倍。     

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

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

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

大迎角振动三角翼跨声速粘性绕流数值分析

本文通过数值求解三维可压缩雷诺平均ＮＳ方程，研究发析了振动三角翼的三维复杂非定常流场，其中对流场项计算格式的空间离散在平行物物理面方向与物面法线方向的空间导数分别采用二阶ＮＮＤ格式和由ＮＮＤ格式高阶插值而得到的三阶迎风格式，时间推进为隐式ＬＵ分解方法，并使用了修正的代数ＢＬ湍流模型模拟高雷诺数的湍流流动。     

大涡模拟模型燃烧室燃烧性能计算

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 k 方程亚网格尺度模型模拟亚网格湍流黏性;亚网格EBU(eddy-break-up)燃烧模型预估化学反应速率;多维经验分析法计算燃烧性能;并在非交错网格体系下气相采用SIMPLE(semi-implicit method for pressure-linked equations)算法对控制方程进行求解,液相采用随机离散模型,两相之间的耦合采用PSIC(particle-source-in-cell)算法.通过大涡模拟瞬态及时均计算结果表明:与粒子图像测速仪(PIV)测量的瞬态速度场、出口温度分布试验数据吻合,表明在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法,数值模拟模型燃烧室两相喷雾燃烧流场,所采用的亚网格模型可以用于燃烧室气液两相喷雾燃烧流场的大涡模拟;燃烧性能计算结果与试验测量结果基本一致,说明所采用多维经验分析法可以用来数值模拟航空发动机燃烧室燃烧性能的计算,特别是污染物的预估,为设计低污染高性能航空发动机燃烧室提供有用的设计依据.      

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

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

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

三维机翼非定常数值模拟

本文研究了三维非定常Euler方程的显式有限体积解法,并编制了相应的计算程序。在程序中,网格生成部分采用的是代数方法生成贴体网格,再用求解椭圆型方程的方法对网格进行优化;在计算部分采用二步Lax-Wendroff格式的一种变体进行迭代,它具有二阶时间精度,并在预估步和修正步采用同样的处理过程,避免了对中间点的定义,使程序简洁有效率。全部计算结果均与有关实验和其它资料进行了比较。     

喷嘴损伤对环形回流燃烧室性能的影响

利用Fluent商用软件对模型环形回流燃烧室三维两相喷雾燃烧流场进行了数值模拟,研究了喷嘴损伤引起雾化效果变化对燃烧室性能的影响,采用可实现的k-ε模型模拟湍流黏性、离散相模型（DPM）通过添加UDF(user defined function)程序追踪燃油运动轨迹、正庚烷作替代燃料及层流小火焰模型.计算结果表明:采用的数值模拟方法可以预估实际燃烧室燃烧流场以及喷嘴损伤对其性能的影响,雾化性能变化导致燃烧室出口温度分布不均匀度升高,品质降低,并导致燃烧室燃烧效率降低;当燃油流量降低约19%时,燃烧室性能已不符合运行要求.      

H2引燃雾化煤油超燃混合的数值研究

为研究H2引燃的煤油超燃混合问题，采用流体模型描述两相流动，气相反应系统的全N－S方程风迎风TVD格式求解，液相扩散Euler方程用NND格式求解，化学反应源项采用点隐处理，对单喷嘴喷H2超声速燃烧、H2引燃的煤油超燃混合问题进行了数值研究。获得了压力、密度、温度和组元浓度场的分布。结果表明：与单喷嘴喷H2相,比下游喷嘴H2的射流穿透浓度和扩散范围增大，燃烧区也变大；虽然煤油穿透深度大，但其扩散较H2差，且下游出现无煤油区。     

H_2引燃雾化煤油超燃混合的数值研究

为研究H2 引燃的煤油超燃混合问题 ,采用流体模型描述两相流动 ,气相反应系统的全N S方程风迎风TVD格式求解 ,液相扩散的Euler方程用NND格式求解 ,化学反应源项采用点隐处理 ,对单喷嘴喷H2 超声速燃烧、H2 引燃的煤油超燃混合问题进行了数值研究。获得了压力、密度、温度和组元浓度场的分布。结果表明 :与单喷嘴喷H2 相比 ,下游喷嘴H2 的射流穿透深度和扩散范围增大 ,燃烧区也变大 ;虽然煤油穿透深度大 ,但其扩散较H2 差 ,且下游出现无煤油区。     

环形燃烧室两相燃烧数值研究

根据环形燃烧室设计需要,运用计算机模拟技术研究模型环形燃烧室三维有旋液雾燃烧流场,研究两相流动与混合,液雾蒸发对燃烧过程的影响。两相流动模型采用颗粒群轨道模型,气相紊流输运方程求解采用双方程k-ε模型和EBU-Arrhenius紊流燃烧模型,热辐射采用热通量法模型,本文提供计算方法可以用来研究不同性能参数对环形燃烧室两相燃烧特性影响。     

Review of atomization mechanism and spray characteristics of a liquid jet in supersonic crossflow

The injection and atomization process of a liquid fuel jet is critical for an ignition start of a scramjet engine. Airwall-mounted crossflow injection strategy is widely used in scramjet combustors, avoiding high total pressure loss and allowing the liquid fuel to rapidly undergo atomization, mixing, and evaporation. In this review, research progress on a liquid jet in supersonic crossflow was evaluated from aspects of atomization mechanism and spray characteristics. When a liquid jet is injected into a supersonic crossflow, primary and secondary breakups occur successively. The surface instability of liquid can significantly affect the breakup process. This review discusses the current understanding of the breakup process and spray characteristics of a liquid jet in supersonic crossflow including the mechanism of atomization and the characteristics of distribution and atomization. The development of windward Rayleigh-Taylor (R-T) unstable waves is the main factor in column breakup. The development of Kelvin-Helmholtz (K-H) unstable waves along the circumferential direction of the jet or droplets is the main factor of surface and droplet breakups. The liquid–gas momentum ratio is the most important factor affecting the penetration depth. The span width of the liquid jet is affected by the windward area. Breakup and coalescence lead to a transformation of the size distribution of droplets from S- or C-shaped to I-shaped, and the velocity distribution of the droplets on the central symmetry plane has a mirrored S-shape. The droplet distribution on the spanwise cross-section retains a structure similar to an “Ω” shape. At last, some promising recommendations have been proposed, namely a theoretical predictive model which can describe the breakup mechanism of a liquid jet, the distribution characteristics and droplets size distribution of a liquid jet under a cavity combustion chamber, especially for enthalpy flows with complex wave structures.     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

基于火焰面模型的超声速燃烧混合LES/RANS模拟

为了明晰超燃冲压发动机燃烧室内部燃烧过程的细节,建立了超声速湍流燃烧稳态火焰面亚格子模型,并采用混合LES/RANS方法对氢燃料超燃冲压发动机进行算例验证.控制方程对流项用五阶精度WENO格式离散,时间方向采用二阶Runge-Kutta方法.研究表明:(1)冷流流场中燃料分布与大尺度结构分布相似,说明混合过程受大涡控制;(2)燃烧流场中涡的尺寸明显变大,且仅存在于火焰面上,另外温度分布和主要生成物分布与涡量云图基本相同,说明燃烧过程也由大涡控制;(3)时均计算结果与实验阴影基本符合,速度剖面和温度剖面与实验测量值定性一致,说明本文的数值模拟方法和燃烧模型可以较好地描述和预测超声速流动燃烧过程.     

后掠角度对超声速混合的影响

通过采用隐式方法求解三维、粘性、湍流 Navier-Stokes方程及组分方程 ,对以氢气为燃料的超燃冲压发动机燃烧室模型进行了数值模拟 ,分析了不同后掠角度对于超声速混合的影响。计算结果表明后掠斜坡产生的流向涡的强度以及混合增强效果优于无后掠斜坡 ,并随着后掠角度的增加而增强 ,但总压恢复系数随后掠角度的增加逐渐降低      

横流中喷雾掺混流场结构研究

为了获得有限空间内喷雾在横流影响下的掺混发展及机理,应用PIV系统测量了单个旋流雾化喷嘴产生的喷雾在横流中的掺混流场,实验在矩形通道内实施。获得了3种喷嘴入射角度和3种液气动量比下的掺混截面流场结构。结果表明,由于射流撞击、剪切和壁面约束作用,流场中形成前缘涡和反旋涡对,反旋涡对对掺混起了主要作用。比较不同喷嘴入射角度和液气动量比下的流场结构发现,随喷雾入射角度减小,两相作用提前发生,且反旋涡对尺度小,更利于掺混均匀;随液气动量比减小,气相作用增强,旋涡强度和尺度变小,更小的涡尺度更利于液滴均匀分布。     

二维横向喷流的N-S方程数值解

让一股射流横向注入超音速流产生所需要的干扰效应是一种飞行器飞行控制的常用手段,如航天飞机和机动导弹上的反作用控制系统(RCS).所以,对射流与主流相互干扰的某些流动特性应该有个基本的了解.本文通过应用MacCormack显格式和Baldwin-Lomax修正的代数湍流模型求解二维RANS方程对带横向射流的绕后台阶的M_∞为2.19和6.00的超音速外干扰流场进行了数值计算.与无喷射流的情形相比,回流涡旋区由2个增加到4个,且底部压力有大的回升,喷口宽度越大,回升效应越显著.另外,数值计算表明,在压力梯度变化比较平缓的区域附加一四阶人工粘性能有效地抑制数值振荡,加快收敛.