幂律型流体雾化SPH方法数值分析

幂律型流体的雾化过程存在复杂的界面运动,用传统网格法很难精确追踪运动界面.为研究幂律型流体的雾化特性,运用SPH方法对典型的双股幂律型流体撞击雾化问题进行三维数值分析.根据文献实验,采用SPH方法模拟获得与实验条件相应的数值结果,对比后表明,二者雾化角相当吻合,数值结果还成功捕捉到液膜向液丝的破碎过程及网状的液丝分布状态,验证了方法的有效性.分析了射流速度和撞击角度对雾化角的影响,得到雾化角随着射流速度和撞击角度的增加而增大;对雾化后的速度场进行数值分析后表明,撞击点附近,惯性力的作用使速度场变化剧烈;在流体远离撞击点的过程中,粘性耗散作用使速度场趋于稳定,但速度大小小于初始撞击速度.     

直流互击式喷注单元雾化特性准直接数值模拟

为实现直流互击式喷注单元雾化过程的数值求解并探究结构及工作参数对雾化特性的影响规律,基于一种树形自适应加密算法与流体容积(Volume-of-Fluid,VOF)方法实现了雾化过程的准直接数值模拟。计算得到了两股射流由喷射到撞击形成液膜,液膜进一步破碎形成液丝、液滴的全过程,获得了液膜的破碎长度、液滴的Sauter平均直径、雾化频率等雾化特性参数。通过将典型算例的计算结果与试验数据进行对比验证了计算的有效性,给出了数值求解精度。探讨了撞击波的形成机理,分析了雾场液滴的尺寸分布规律,撞击夹角、孔径比、射流速度、动量比对雾化特性的影响规律。结果表明:所采用的算法可以实现多相、多尺度雾化过程的数值求解;撞击波的形成是由于两股射流撞击时惯性力的不完全对称导致的;雾场液滴的尺寸分布近似服从Rosin-Rammler分布;撞击夹角的增大与射流速度的提高导致液膜的破碎长度减小,液滴的平均粒径减小,撞击夹角增大雾化频率呈减小的趋势,射流速度提高雾化频率呈增大的趋势;动量比为1而孔径比不为1时会形成凹形液膜,雾场存在一定程度的偏斜;动量比主要影响雾场的偏斜程度。     

幂律型流体撞击射流破碎特性直接数值模拟

凝胶推进剂的双股射流撞击雾化广泛应用于液体火箭发动机的燃烧室中,其破碎特征及雾化效果直接影响燃烧效率。为探究雾化特性的发展规律,采用直接数值模拟DNS方法,对射流速度为100m/s的剪切稀化非牛顿液体正交撞击产生的雾化特征、液体表面积、表面波、涡特性以及非牛顿特性开展研究。结果表明,射流下形成的雾化流场迅速扩张形成液膜,液膜两侧边缘破碎成大量的液丝与液滴,核心部分产生撞击波后在气体力的作用下逐步发展为带有凸起和褶皱的不稳定表面波,其撞击波波长最大可达2.46倍射流直径。液体表面积不断增长,但无量纲表面积总体呈现先下降再上升的趋势。气体中的涡量分布则分为有序附着区和无序爆炸区两类,并且涡量主要集中分布于气相区域。此外,射流撞击时产生强剪切使该液体内部的粘性系数下降,最低仅为初始粘性系数的0.3倍。     

基于SIFT算法的含碳颗粒凝胶推进剂雾化场速度分析

雾化是凝胶推进技术的关键问题之一,雾化过程中的液膜、液丝及液滴等的速度分布、液膜厚度等参数对于判断雾化效果、揭示雾化机理具有重要作用。为了对含碳颗粒凝胶推进剂雾化场的速度进行定量分析,提出了一种基于SIFT关键点匹配的雾化场速度计算方法,并以雾化场速度分析为基础,提出了一种新的液膜厚度估算方法。结果表明:雾化场速度随射流速度的增大而增大,随撞击角度的增大而减小;雾化场平均速度与射流速度的比值vato/vjet在0.6~0.9,可用于表征雾化效果,vato/vjet越小,雾化效果越好;射流撞击形成液膜的厚度在0.04~0.13mm,液膜厚度随着射流速度及撞击角度的增大而减小;凝胶推进剂的雾化效果随着碳颗粒浓度的增大而降低,随着碳颗粒粒径的增大而略有改善。     

撞击角对撞击式喷嘴雾化特性影响研究

为了研究液体火箭发动机中撞击式喷嘴的雾化问题,采用CLSVOF(Coupled Level-Set and Volume-of-Fluid Method)方法对撞击式喷嘴的雾化过程进行了数值仿真模拟,重点考察了撞击角对撞击式喷嘴雾化特性的影响。结果表明,通过CLSVOF方法能对射流撞击雾化形态进行较好的捕捉;分析整个雾化区域的液相分布,将射流撞击雾化过程中的液相分布大致分为三个区域:撞击前两股独立射流区域、撞击后形成的液膜区域、液膜破碎后液丝和液滴的生成区域;射流撞击雾化除了气液界面上速度差引起的不稳定之外,还存在由射流湍流或撞击波引起的其它不稳定因素;液膜破碎长度随撞击角的增大而减小,液膜表面上的表面波幅值及液膜破碎程度则随撞击角的增大而增大;射流撞击在撞击点位置处所形成的速度差对液膜的破碎和液膜上表面波幅值的大小起到了关键作用。     

液体射流撞击液膜振荡行为的实验研究

液体射流撞击过程是液体火箭发动机常见的高效雾化方式。为了进一步探究液体射流撞击雾化过程,采用高速摄像仪和图像处理技术对两股液体射流撞击液膜的振荡行为进行了实验研究,考察了射流韦伯数(21≤We≤1356)和撞击角(2θ=60°,90°,120°)对液膜振荡的影响。结果表明,随着射流韦伯数增大(We250),撞击液膜会从竖直稳定模式转变为振荡模式,促进液膜破裂和雾化。增大撞击角能加剧液膜振荡和增加雾化角。2θ=60°和90°的液膜振荡区间为250We400,2θ=120°的液膜振荡区间为We250,无量纲振幅的时均极大值约为5。结果揭示了随着液体射流速度的增大,射流不稳定性的发展会引起撞击点的动量不平衡,进而形成液膜振荡。此外,液膜向下传播时与空气界面间的Kelvin-Helmholtz不稳定性也促进液膜振荡。     

倾斜射流撞壁实验研究及液膜几何参数建模

倾斜射流撞壁在液体火箭发动机液膜冷却、射流撞壁雾化等领域具有广泛的应用。为了研究倾斜射流撞击壁面后形成的液膜的关键特征，开展单股圆柱射流撞击壁面的实验研究。从实验中研究各射流参数对液膜外形的影响规律，继而开展理论建模，获取液膜外形的关键几何参数表达式。实验研究发现随着射流倾角增大，液膜长度减小而宽度增大，随着射流孔径和射流速度增大液膜长度和宽度均增大这一定性规律。理论分析得到了液膜最大宽度位置与液膜对称面的夹角近似等于射流倾角α，液膜的长宽比近似等于1+cot α。通过进一步分析得到液膜宽度、最大宽度距撞击点的距离、液膜长度这几个液膜关键参数的表达式。建立的液膜几何参数表达式预测结果与实验结果的误差均在20%以内。     

动量比对两股互击式喷嘴雾化特征的影响

基于CLSVOF（coupled level-set and volume of fluid）方法对两股互击式喷嘴雾化过程进行了数值模拟,仿真结果与实验结果进行对比,并详细考察了动量比对液膜特性的影响.结果表明:①CLSVOF方法能对射流撞击形成的液膜形态进行有效捕捉,计算获得的液膜形态与实验结果一致.②两股不同孔径的射流在撞击点撞击过程中,由于大孔径一侧的部分射流未参与撞击,而中心位置的液体在撞击之后速度迅速下降,周围的液体则保持撞击前的速度,使得射流撞击之后形成弧形状液膜.③在不同动量比条件下,射流撞击后轴线位置的合成角主要取决于参与撞击射流的动量比,数值计算结果与理论计算结果吻合较好,误差均小于10%.      

基于PIV的凝胶模拟液撞击雾化速度场实验研究

为了更好研究凝胶推进剂的雾化,采用时间分辨粒子图像测速(TR-PIV),研究了不同撞击角度(45°,60°,75°,90°和120°)和射流压差(0.4MPa~0.8MPa)对凝胶推进剂雾化速度的影响。实验结果表明:雾化液滴速度对于撞击轴线呈单峰对称分布,距离撞击点越远,雾化液滴速度越小且分布越均匀;增大撞击角和增大射流压差都可提高凝胶推进剂有效撞击速度,即增加撞击后液体动能转换液体破碎所需的能量,雾化质量提高;当有效撞击速度大于27.9m/s时,实验室配置的凝胶推进剂可充分雾化。     

射流自由长度对凝胶推进剂撞击雾化影响的实验

为研究射流自由长度对凝胶推进剂撞击雾化的影响,建立了撞击雾化实验台,制备了煤油凝胶和水基模拟液,测试了雾化装置的喷射系数及模拟液的黏性和稳定性。分析了3种射流速度不同射流自由长度条件下的凝胶撞击雾化特性,观测了射流和撞击喷雾图像。测量了液膜破碎长度、雾化液滴粒径分布和相应的SMD(索太尔平均直径)值。研究结果表明:低速时,随着射流自由长度增大,撞击液膜的喷雾形态会发生较大变化,而高速条件下,雾化形态则基本一致。3种射流速度下,破碎长度在45~9mm之间,并随射流自由长度逐渐减小。液滴分布服从Rosin Rammler规律,并具有较高的拟合精度。均匀度指数均在3~4之间,并随射流自由长度逐渐降低,粒径均匀度降低;较高射流速度下,SMD随自由长度逐渐增大。低射流速度时,SMD随射流长度先减小后增大,射流自由长度存在一个最优值,在设定研究条件下其值为25/3。因此,在设计撞击喷嘴时,根据射流速度选取适当的自由长度值可以获得更好的雾化效果。      

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

基于H∞性能指标的质量矩拦截弹鲁棒控制

 以所建立的质量矩拦截弹数学模型为基础，利用微分几何的反馈线性化理论，得到一个解耦线性系统，考虑到拦截弹的鲁棒性要求和3个滑块的协调控制问题，提出采用双回路的设计方法，内回路采用线性二次调节器（LQR），外回路采用考虑混合灵敏度问题的H_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

航天器返回地球的气动特性综述

航天器返回地球的飞行过程中,气动特性是实现将宇宙飞行速度减到落地前速度、保证再入飞行得到有效控制以及再入防热安全可靠的关键因素。针对简单旋成体气动外形、半弹道式再入控制、烧蚀防热类返回航天器,综述了返回地球过程中变化的空气流域特性、航天器周围的气体绕流环境、空气与航天器作用产生的动力学与热效应等。系统地给出了该类航天器的再入气动特性参数与飞行性能的共性规律,包括:气动阻力与再入减速、气动升力与再入轨迹控制、配平攻角与飞行稳定性、气动加热与防热,以及再入过程中不同气动特性航天器、气象条件变化等对再入飞行性能的影响规律。为航天器开展返回飞行过程的跨流域气动性能工程研制提供设计参考。     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

基于马赫数分布可控曲面外/内锥形基准流场的前体/进气道一体化设计

提出了一种高超声速飞行器乘波前体的外锥形基准流场设计方法,在锥面马赫数分布规律给定的条件下,通过有旋特征线法实现反设计,提高了基准流场设计的灵活性。该基准流场通过锥形"下凹"弯曲激波和波后等熵压缩波系压缩气流,可以在较短的长度内完成高效压缩。基于反正切马赫数分布外锥形基准流场设计的乘波前体具有较高的容积率,乘波特性良好且出口均匀,设计点时有黏升阻比为1.89。另外,基于该乘波前体和马赫数分布可控的内收缩进气道给出了一种双乘波的前体与进气道一体化设计方案,实现了内外流分别独立乘波,充分发挥了乘波前体和内收缩进气道的各自优势。     

Potentials of manufacture and repair of nickel base turbine components used in aero engines and power plants by laser metal deposition and laser drilling

IntroductionExpensive turbine parts like HPT(HighPressure Turine)blades or vanes are replaced bynew parts in case of damage.For example theburn through of the inner side of a blade or vane(Figure 1)is a frequently appearing damage,which cannot be repaired…     

Enhancing the efectiveness of film cooling

Advanced gas turbine stages are designed to operate at increasingly higher inlet temperatures to increase thermal efficiency and specific power output.To maintain durability and reasonable life,film cooling is needed in addition to internal cooling,especially for the first stage.Film cooling lowers material temperature by forced convection inside film-cooling holes and by forming a layer of coolant about component surfaces to insulate them from the hot gases.Unfortunately,each cooling jet forms a pair of counter-rotating vortices that entrains hot gas and causes the film-cooling jet to lift off from the surface that it is intended to protect.This paper gives an overview of efforts to enhance the effectiveness of film-cooling.This paper also describes two new design concepts.One design concept seeks to minimize the entrainment of hot gases underneath of film-cooling jets by using flow-aligned blockers.The other design concept shifts the interaction between the approaching hot gas and the cooling jet to occur further above the surface by using an upstream ramp.For both design concepts,computational fluid dynamics results are presented to examine their usefulness in enhancing film-cooling effectiveness.      

Guide for Authors

<正>About Journal Chinese Journal of Aeronautics(CJA)is a comprehensive academic journal dealing with the fields of aeronautics and astronautics.It reports researches concerning the two fields in China and abroad to promote the academic exchange.Founded in 1988 and sponsored by the Chinese Society of Aeronautics and Astronautics and Beihang University,CJA publishes papers bimonthly,with issues released in February,April,June,August,October and December.