相邻离心式喷嘴液膜撞击雾化过程仿真

基于Gerris软件建立的锥形液膜雾化破碎过程数值仿真方法,对相邻多个离心式喷嘴液膜撞击雾化过程进行了数值仿真,可视化展示了喷雾场三维形态和结构特征,并做了流场分析,讨论了喷嘴之间的相互干扰作用,获得了液滴空间分布,对相邻不同数目的喷嘴雾化效果进行了比较。结果表明:双喷嘴液膜撞击会形成一个类似互击式射流撞击形成的扇形区域,但又与其不同;锥形液膜的撞击效果主要依赖径向速度而不是切向速度,液滴平均粒径与撞击液膜间的旋向几乎无关。多喷嘴液膜撞击会使得液滴空间分布发生大的变化,液膜破碎长度会缩短。多喷嘴液膜撞击后的雾化特性是否改善与单喷嘴原先的雾化特性以及喷嘴间的距离密切相关。对于文中特定的喷嘴结构、喷嘴间距及排列方式,双喷嘴的液滴SMD均比单喷嘴的增大约4.8%~6.1%;四喷嘴的液滴SMD比单喷嘴的增大约7.3%,比双喷嘴的增大约2.4%。     

基于Gerris的离心式喷嘴锥形液膜破碎过程数值模拟

为了实现离心式喷嘴锥形液膜破碎过程的精细可视化和雾化特性的准确计算,基于自适应网格加密技术和VOF方法,建立了一种数值方法,可以捕捉到包含一次雾化和二次雾化的喷雾场结构特征,获得了全场液滴的空间分布和粒径分布,液滴捕捉效果逼真。针对给定结构尺寸的敞口型离心式喷嘴,计算了在不同流量和不同切向孔直径下的液膜锥角、液滴平均粒径SMD的变化,并与高速摄影拍摄的锥形液膜破碎过程图像和PDPA(Phase Doppler Particle Analyzer)测量的液滴SMD对比。结果表明,数值模拟与实验所获得的喷雾场吻合,获得的液膜锥角和液滴平均粒径最大相对误差分别为4.18%和11.82%,结果吻合较好,验证了数值方法的准确性,为离心式喷嘴的精细研究和设计应用提供一种新手段。研究表明,切向孔直径对液膜锥角和液滴平均粒径的影响较显著,在设计加工时,是一个比较重要的结构参数。     

内直外旋气液同轴式喷嘴流量及雾化特性

针对气体中心直流、液体涡流器离心式内缩进气液同轴喷嘴,采用高速摄影和PDPA测量了其流量和雾化特性。发现在常温常压下,喷嘴流量近似与喷注压降的0.477次方成正比。喷嘴喷雾场为空心锥形,雾化角约为110°,随液体压降增加,雾化锥角增大。液体压降不变,雾化锥角受中心气体卷吸作用随气体压降增加而减小。液膜破碎距离与雾化锥角有关,雾化锥角越大,液膜破碎距离越短,反之亦然。随径向距离增加,雾化液滴轴向速度、径向速度减小,SMD先增大然后略有减小,沿径向呈"马鞍"状双峰分布;随轴向距离增加,粒子轴向速度减小,SMD减小。     

气/液同轴离心式喷嘴流量及雾化特性实验

通过实验方法，研究气／液同轴离心式喷嘴的流量特性及气液之间的相互影响，研究压降及缩进长度对喷嘴雾化特征的影响。研究表明，液喷嘴流量受气喷嘴压降影响，流量增大或减小取决于缩进室内气喷嘴对液喷嘴产生的是引射还是堵塞作用；同一喷嘴，液膜破碎长度及液滴尺寸随压降增大而减小；液滴尺寸分布沿喷嘴横截面呈马鞍型分布，曲线的最高点和分布宽度随喷嘴几何特性及压降变化；提高气／液之间的相对运动速度，可提高雾化质量；随着缩进长度增大，雾化质量明显变好。     

离心式喷嘴在小尺度受限空间雾化特性的实验研究

为了了解小尺度受限空间中离心式喷嘴的雾化特性,设计了四旋流槽离心式喷嘴和小尺度模拟燃烧室,采用相位多普勒粒子动态分析仪(PDA),观测了喷射压力对模拟燃烧室内雾化性能的影响.引入离散系数cv对喷雾场雾化参数的周向分布规律进行了着重分析.结果表明:模拟燃烧室内,喷射压力越大,液滴索特尔平均直径D32的值越小;在z∈[41.5,94.2]的轴向距离范围内,喷射压力对Vz的影响不明显.由于壁面作用,液滴D32的周向分布均匀性随着径向距离的增加变差;距离喷嘴出口越远处,液滴D32的周向分布均匀性越好.在喷射压力相同条件下,同一喷嘴在模拟燃烧室内雾化参数的周向分布均匀性比在大气环境中波动更大.     

不同液压条件下锥形液膜破碎的实验研究

应用PIV激光粒子图像测速仪对由离心式压力雾化喷嘴喷出的锥形液膜的破碎进行了实验研究。澜得了不同压力下的扰动波波长、波速和振幅。研究结果表明，液膜半锥角随着液压的增加而增加，半波长的液膜段剥落后几乎立即破碎成液滴，看不到明显的液膜段收缩成液丝的过程。由实验拍摄的图像看，长、短扰动波同时存在，共同影响液膜的破碎。这些数据可用于发展锥形液膜初始雾化模型。     

切向孔单/双排布局对离心式喷嘴锥形液膜雾化特性影响

为了研究切向孔单/双排布局对离心式喷嘴雾化性能的影响,采用高速摄影仪对两种结构喷嘴形成的锥形液膜进行了拍摄分析,测量了不同压降下液膜的破碎长度、喷雾锥角、出口轴向速度以及表面波破碎点波长等雾化性能参数,并使用这些性能参数作为初始条件,求解了锥形液膜色散方程。结果表明:双排切向孔布局使得液膜雾化性能变差,体现为相同压降下喷雾锥角减小1°~2°,破碎长度增大,增大幅度最大达8.2%;色散方程求得的破碎点波长理论值与试验值吻合较好,验证了所用液膜破碎模型的有效性;单排切向孔喷嘴形成液膜的主导扰动波增长率大于双排切向孔喷嘴的主导波增长率,从直观上解释了单排切向孔喷嘴破碎长度小于双排切向孔喷嘴的原因。     

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

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

基于移动粒子半隐式方法的旋流液膜破碎过程模拟

为实现旋流初始雾化中液膜破碎过程的直接数值模拟,基于GPU (Graphics processing unit)加速和移动粒子半隐式方法 (Moving particle semi-implicit method, MPS),开发了离心式喷嘴液膜破碎过程的并行加速数值模拟方法与程序,模拟了喷雾场三维形态特征和初始雾化破碎过程。模拟结果成功捕捉到了液膜形成、液膜破碎成液丝继而破碎成液滴的瞬态过程。模拟得到的雾化破碎过程与实验拍摄结果基本吻合。模拟了不同射流速度下的旋流液膜破碎过程,模拟得到的液膜初始破碎长度与经验公式计算结果趋势一致,二者吻合较好,最大误差为24.2%,模拟得到的液膜半锥角与实验值较为吻合,误差为10.6%。表明开发的模拟方法与程序的准确性,为后续离心式喷嘴的液膜雾化过程及雾化特性研究打下基础。     

自激振荡过程中喷嘴内部扰动的产生及传递机制

气液同轴离心式喷嘴是液体火箭发动机采用的极其重要的喷注单元类型。但其容易产生自激振荡，可能会诱发不稳定燃烧，严重影响发动机的研制，而目前对于自激振荡过程中喷嘴内部非定常流动过程的认识还十分有限。借助试验及数值仿真手段对液体中心式气液同轴离心式喷嘴自激振荡过程中液膜的运动形变轨迹、振荡的产生原因及压力扰动在喷嘴内部的传递过程进行了系统分析。研究发现：喷嘴内部流场的压力振荡是由于缩进室内部液膜的周期性堵塞作用产生的，且振荡最先发生于缩进室内部液膜外侧；此后，压力振荡会以一定的相位差通过气体环缝及内部离心式喷嘴向上游传递，从而引起环缝气体、中心气核压力发生相同频率的振荡；当自激振荡强度足够大时，便会引起集液腔压力振荡，且振荡频率与自激振荡频率一致。     

ν-gap度量及其在飞行控制律评估中的应用

 传统控制律评估方法主要用于单输入单输出（SISO）系统，且对模型参数摄动考虑不够全面,针对这些不足，研究了ν-gap度量方法。在介绍系统广义稳定裕度相关概念的基础上，给出了ν-gap度量的定义、特点和性质以及近似摄动模型的计算，提出ν-gap度量评估控制律的步骤。实例结果表明，该方法不仅克服了上述传统评估方法的缺陷，而且还有根据所求的各摄动影响情况忽略影响小的元素，以减少计算量及可以找到最坏情况下的参数摄动组合等优点。     

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

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

Integrated Entry Guidance for Reusable Launch Vehicle

A method for the implementation of integrated three-degree-of-freedom constrained entry guidance for reusable launch vehicle is presented. Given any feasible entry conditions, terminal area energy management interface conditions, and the reference trajectory generated onboard then, the method can generate a longitudinal guidance profile rapidly, featuring linear quadratic regular method and a proportional-integral-derivative tracking law with time-varying gains, which satisfies all the entry corridor constraints and meets the requirements with high precision. Afterwards, by utilizing special features of crossrange parameter, establishing bank-reversal corridor, and determining bank-reversals according to dynamically adjusted method, the algorithm enables the lateral entry guidance system to fly a wide range of missions and provides reliable and good performance in the presence of significant aerodynamic modeling uncertainty. Fast trajectory guidance profiles and simulations with a reusable launch vehicle model for various missions and aerodynamic uncertain-ties are presented to demonstrate the capacity and reliability of this method.     

二维翼段颤振的μ控制

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

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

Abnormal Shape Mould Winding

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

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

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

航空发动机管路测量数据分割方法

 为了提高航空发动机管路测量数据的反求建模效率,提出了一种区域增长分割算法。该方法主要是利用管路表面在曲面索引系数映射、主曲率映射、高斯映射上的特性,并基于均值漂移算法和遗传算法提出了区域分割算法中种子区域的选择策略。然后利用区域增长分割方法实现了对管路测量数据的分割。经过仿真和实测数据验证,所提管路分割算法具有较好的分割质量和效率。     

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

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

Interfacial Properties Modification of Carbon Fiber/ Polyarylacetylene Composites

This work was dedicated to performing surface oxidation and coating treatments on carbon fibers (CF) and investigating the changes of fiber surface properties after these treatments, including surface composition, relative volume of functional groups, and sur- face topography with X-ray photoelectron spectroscopy (XPS) and atom force microscopy (AFM) technology. The results show that, after oxidation treatments, interfacial properties between CF and non-polar polyarylacetylene (PAA) resin are remarkably modified by removing weak surface layers and increasing fiber surface roughness. Coating treatment by high char phenolic resin solution after oxida- tion makes interface of CF/PAA composites to be upgraded and the interfacial properties further bettered.