水滴撞击特性的动态图像测量方法

用石膏脱模的方法制作圆柱实验件.实验采集圆柱试件表面不同位置染色过程的动态视频数据;依据圆柱驻点处色度值时序信号的增长曲线,将实验划分为欠饱和区、平衡区和过饱和区3个阶段;将平衡区RGB色度空间的色度值c与水滴撞击特性的局部水收集系数建立关联.结果表明:平衡区的比色分布符合局部水收集系数相对驻点处β/β_m分布;不同位置c值增长曲线具有相似性,缩比时间尺度1/t_b分布符合数值仿真欧拉法β/β_m分布的定义;实验结果与数值仿真结果对比误差率小于10%.     

发动机燃气喷流红外辐射场的数值模拟

在高温流动中热辐射与流场是耦合的,但在发动机燃气喷流流场的温度范围内可以在流场计算中忽略热辐射的影响,于是流场与辐射场的模型方程可解耦.为此提出了模拟飞行器绕流流场红外辐射的解耦模型以及相应的算法.该算法首先采用总变差减小TVD(Total Variation Diminishing)格式对流场进行模拟,再引入有限体积概念,应用所取得的流场参数,在同一组网格上计算流场中的红外辐射特征.以轴对称喷管内N₂,O₂,CO₂,H₂O,CO,HCl,H₂等7组分高度欠膨胀冻结流与均匀绕流干扰流场为算例进行了验证性的模拟计算,并与已有试验和计算数据进行了对比,表明该算法是可靠的,可在较大幅度地降低计算量的情况下给出满足工程需要的结果.     

改进SA模型对翼型分离流动的数值模拟

采用(CFD,Computational Fluid Dynamics)软件Fluent,选取(SA,Spalart-Allmaras)模型,对风力机翼型S825进行了二维数值模拟,并与实验数据进行比较.针对SA模型不能准确预测翼型尾部流动分离的情况,分析了分离区域内湍流的强非平衡输运特性.研究表明,SA模型中系数C_b1与湍流的非平衡输运特性密切相关,进而提出了修正模型系数C_b1改进SA模型的方法,结果表明该方法能够更准确地模拟风力机翼型分离流动.     

加热器喷管热-流耦合传热分析

针对加热器喷管中复杂的气、固、液多相流动传热问题,建立了三维热流耦合换热计算模型,分别对燃气、冷却剂和喷管室壁建立不同的控制方程,将辐射热量作为源项加入到方程中,进行流动和传热的耦合计算.采用此方法对美国AEDC(Arnold Engineering Development Center)喷管的流动传热过程进行了计算,数值计算结果与试验结果吻合较好.在此基础上对某超燃冲压发动机试验台水冷式加热器喷管的换热问题进行了三维数值模拟,并定量分析了辐射换热对加热器喷管壁面温度分布的影响.结果表明:冷却水流量取2.0 kg/s时,加热器喷管气壁最高温度为660K,膜温度为430 K,加热器能可靠冷却,其热效率满足试验要求,对于含有H2O和CO2这样的高温燃气,辐射热量对喷管壁面温度分布有较大影响,必须引入到温度场的求解之中.     

浮升力对竖直螺旋管中超临界CO2换热的影响及判别准则

基于航空发动机间冷器的冷却换热问题，进行竖直螺旋管内超临界压力CO₂换热数值研究。探究运行参数对沿流向和周向换热的影响机制。通过管截面温度场和流场分布阐述了浮升力和离心力引起的周向非均匀换热机制，评估二次流速度和强度。根据管道结构特性提出新的浮升力参数和浮升力影响判别准则，建立新的换热关联式。结果表明：管上游非均匀换热源于浮升力与离心力综合作用，管下游非均匀换热由离心力主导。当满足浮升力因子Bu≥1.6×10^-5时，浮升力在换热中起主导作用。新换热关联式可以较好地适用于螺旋管内换热预测。     

激发跃迁速率对热力学非平衡氮气紫外辐射的影响

基于氮气的碰撞-辐射（CR）模型,计算了速度为6.2 km/s、初始压力为133 Pa的高超声速流动激波中N₂和N₂+分子电子能级的分布情况,分析了不同激发跃迁速率模型对电子能级分布及辐射光谱模拟的影响。针对流动中热力学非平衡区域和平衡区域,在300~440 nm的辐射光谱开展了逐线法的数值模拟,并与激波管实验测量光谱进行了对比。结果表明,目前的激发跃迁速率均存在偏差,综合Park模型的爱因斯坦系数和Johnston模型的碰撞激发速率可以得到与实验结果最为接近的辐射光谱。      

加力燃烧室NOx排放生成规律试验分析

实现航空发动机大推力通常采用加力燃烧室设计,燃烧室进气温度升高会同步导致排放尾气中NO_x浓度提高.采用典型加力燃烧室主稳定器的V型结构试验件,测试加力燃烧前后(进口600 ℃,出口600~1 200 ℃)燃烧室排放断面的NO_x浓度分布,通过相同流场4种试验工况的NO_x生成浓度组分和氧含量变化分析,验证了加力燃烧室NO_x生成以热力型NO为主的"高温、富氧、贫油"燃烧特性.结果表明加力燃烧室燃烧过程中消耗的氧含量大多贡献于CO₂增量,在800~1 000 ℃燃气排放温度范围内,不完全燃烧产物CO对NO的热力生成有明显抑制作用.     

飞行器进入火星大气的流场预测

针对火星着陆探测器进入-下降-着陆过程的高超声速进入阶段, 利用三维并 行程序求解流体动力学Navier-Stokes方程与化学反应动力学模型, 分析火星 科学实验室进入火星大气时探测器周围的流场结构、化学非平衡效应影响和气 动特性变化规律. 结果表明, 对于完全气体模型, 来流的热力学性质参数选 取影响激波位置和强度. 在化学非平衡效应影响下, 探测器头部激波脱体距离 大幅减小, 驻点压力变化不大, 波后温度显著降低. CO₂在激波后大量分解, 消耗相当能量. 流线结构显示, 探测器尾迹流动中存在复杂的旋涡运动等流动 分离现象.      

火星大气模型参数对MSL气动特性的影响

火星大气与地球大气截然不同,飞行器在进入火星时气动特性不同于地球再入. 大气模型的差异主要表现为气体组份、密度和温度等物理参数. 针对火星进入器MSL在进入-下降-着陆过程中的高超声速进入段,利用三维并行程序求解耦合真实气体模型的流体动力学Navier-Stokes方程,分析MSL进入火星大气时大气模型参数对进入器气动特性的影响. 结果表明,通过与海盗号飞行数据的对比,验证了所采用的火星气体模型和计算方法,且其与NASA的 LAURA代码气动特性计算结果也较为一致;大气模型气体性质,即CO₂环境对进入器阻力系数和俯仰力矩系数影响较大,利用空气得到的计算和实验数据必须考虑CO₂效应;密度增大促进了化学非平衡效应,但对进入器气动特性基本没有影响;温度升高大大增强了化学非平衡效应,而对进入器气动特性影响较小.      

高超声速飞行器后体/尾喷管优化设计

结合试验设计方法、替代模型技术和遗传算法构建一套改进的优化方法,并将其应用于高超声速飞行器后体/尾喷管的一体化设计.在构建优化方法时,针对替代模型采用渐近全局策略提高精度;针对遗传算法,采用实数编码、多目标定级排序和改进小生境技术.在后体/尾喷管的一体化设计应用中,结合高精度的计算流体力学(CFD,Computational Fluid Dynamics)求解,以两个设计点的推力和升力为目标,以力矩为约束,得到优化问题的Pareto最优前沿面,优化结果在综合性能方面有很大提升.该优化方法可进一步推广应用于更复杂的优化设计当中.      

二元喉道倾斜矢量喷管的数值模拟

利用数值模拟方法,对二元喉道倾斜矢量喷管进行了研究.研究了喉道单侧注气、扩张段辅助注气对喷管流场和性能的影响.研究结果表明:喉道单侧注气可以产生不对称的流动,产生矢量推力,但是推力矢量效率较低;扩张段辅助注气可以显著提高喷管的推力矢量性能;只有注气流量比较大时,才会出现典型的"喉道倾斜"现象;但是推力矢量控制效率最高的区域并不是出现在"喉道倾斜"之后,而是出现在弓形激波位置逐渐前移、扩张段注气口上游亚音速区域不断扩大的过程中.     

两种流体控制方案矢量喷管内流场计算及分析

在激波诱导和喉道倾斜2种流体控制方案下,对推力矢量喷管的二维内流场进行了数值模拟.用局部加密的四边形结构网格对流场进行网格划分之后,采用二阶迎风离散格式和二方程湍流模型求解强守恒形式的Navier-Stokes方程,通过对内流场的模拟获得流体控制推力矢量喷管方案的推力矢量特性.计算结果表明:对于激波诱导矢量控制方案,在小落压比大注气量且后缝注气的情况下所产生的矢量角最大;对于喉道倾斜矢量控制方案,在扩张片上以与主流成某一角度的方向注入气流产生的矢量角最大且没有带来较大的推力损失;喉道倾斜矢量控制方案因其推力损失较小且能够降低喷管的重量与造价而更具发展前景.      

吸气式高超声速飞行器冷流试验设计及验证

对于吸气式飞行器而言,地面冷流试验是检验其进气道性能及气动特性的一项重要手段.以二元混压式进气道、机体/推进系统耦合为基本特征,设计了采用超燃冲压发动机为推进系统的内外流一体化巡航飞行器,针对其高超声速特性开展了冷流风洞试验,来流速度范围Ma=5.0~7.0,攻角范围α=-4°~8°.测压试验结果表明,随着来流马赫数的增大,进气道的总压恢复系数下降;而流量系数先上升,在设计点达到最大值;在一定攻角范围内,进气道的总压恢复系数和流量系数提高,但当攻角增大至巡航攻角时,随着攻角的增大,进气道的总压恢复系数和流量系数逐渐下降.测力试验验证了数值算法的有效性,除轴向力系数以外,其余气动特性系数的发展规律及数值基本吻合,可通过修正试验值的方式外推出飞行器的气动特性数据.     

喷管两相流耦合数值模拟与粒子热增量预测

含有氧化铝粒子的两相流是固体火箭发动机喷管流场的重要特征.在有限体积方法框架下,采用基于热增量试验数据的粒子壁面反弹模型以及基于粒子轨道的单元内颗粒源(PSIC,Particle Source in Cell)两相流耦合算法,对喷管内两相流流场及粒子撞击产生的壁面热增量进行了计算和分析,研究了氧化铝粒子尺寸对粒子轨道分布和喷管壁面热增量分布的影响规律.研究结果表明:喷管扩张段内粒子稀疏区域范围随粒子直径增加而增大;粒子热增量只分布于喷管收缩段内,粒子直径越大,产生的壁面热增量越强.     

基于辅助进气门的进气道/发动机一体化控制

针对进气道与发动机的耦合问题,研究了低速大迎角状态下,基于辅助进气门的进气道/发动机一体化控制。首先,建立了飞行条件、迎角、辅助进气门开度与出口总压恢复系数和流量相关联的进气道实时模型,进而将进气道出口流量和发动机进口流量相匹配,建立了进气道/发动机一体化模型的控制仿真平台。其次,为了解决大机动过程中发动机进口流量不足和压力不均的问题,提出了一种带有辅助进气门调节的进气道/发动机一体化控制方法,即通过调节辅助进气门开度实现进气道出口总压恢复系数控制,在保证进气道出口性能稳定的情况下,基于H_∞鲁棒控制方法实现对发动机转速和压比的控制。研究结果表明,在整个大机动过程中,所提出的进气道/发动机一体化控制可以使得发动机各项性能保持稳定,在典型任务工况下,推力提高了16%,耗油率下降了6%。      

喷流干扰气动热数值模拟的若干影响因素

工程上对喷流气动热进行近似计算时,常常对真实的全化学反应流近似处理,这时需要了解喷流气体和喷管外形的近似方法对计算结果的影响规律。采用数值试验对比分析的方法,探讨了近似计算的2个影响因素:其一是使用量热完全的空气喷流近似燃气喷流进行数值模拟时,喷管几何形状对热流分布的影响;其二是使用空气喷流近似燃气喷流进行数值模拟的方法,分析不同的喷流热力学参数匹配方法对喷流形态及热流分布的影响。通过对不同方法进行对比计算,揭示了各个方法对喷流干扰气动热计算的影响规律,可以用于指导工程分析。      

集气腔总压对电弧喷射推力器工作过程的影响

为了研究集气腔总压对电弧喷射推力器工作过程的影响,在分析其工作机理的基础上采用化学非平衡流动和稳态电磁流体电磁场模型对不同集气腔总压下推力器工作过程进行了数值模拟.流体力学方程组和电磁场方程考虑了多种流动机理及电磁场与高温电离气体的相互作用,化学动力学模型考虑了各种碰撞反应.采用二阶精度NND格式求解流体力学方程组,采用有限速率化学反应模型计算组分生成率,采用交替方向隐式(ADI)超松弛迭代法求解电磁场离散方程.给出了不同集气腔总压下推力器内部参数分布及其宏观性能.研究表明,集气腔总压对推力器工作过程具有多方面的影响,在保证电弧稳定的情况下,适当提高集气腔总压可同时提高比冲和推进效率.     

高超声速热化学非平衡空间格式的扩展与改进

扩展、改进和分析了ROE,MAPS,LDFSS,HLLC,HCUSP和ECUSP6种高分辨率空间格式于热化学非平衡流的计算中.采用5组分17反应的Dunn-Kang化学反应模型和Park的双温模型计算了三维球头算例,得到如下结论:ROE和MAPS格式精度高于LDFSS,HLLC和CUSP类格式;基于ROE平均的界面音速可提高格式的鲁棒性及斜激波捕捉能力,同时对流场参量的计算无影响;ROE平均音速的引用改进了MAPS和ECUSP格式.为热化学非平衡流场计算格式的选择提供了参考.     

Quantifying the relationship between the measurement precision and specifications of a UV/visible sensor on a geostationary satellite

To investigate the feasibility of new satellite observations, including air quality (AQ) observations from geostationary (GEO) orbit, it is essential to link the measurement precision (ε) with sensor specifications in advance. The present study attempts to formulate the linkage between ε and specifications of a UV/visible sensor (signal-to-noise ratio (SNR), full width at half maximum (FWHM) of the slit function, and sampling ratio (SR)) on a GEO satellite. A sophisticated radiative transfer model (JACOSPAR) is used to calculate synthetic radiance spectra that would be measured by a UV/visible sensor observing the atmosphere over Tokyo (35.7°N, 139.7°E) from GEO orbit at 120°E longitude. The spectra, modified according to given sensor specifications, are analyzed by the differential optical absorption spectroscopy technique to estimate the ε for slant column densities of O₃ and NO₂. We find clear relationships: for example, the ε of the O₃ slant column density (molecules cm⁻²) and SNR at 330 nm are linked by the equation log(ε) = −1.06 · log(SNR) + 20.71 in the UV region, and the ε of the NO₂ slant column density and SNR at 450 nm are linked by log(ε) = −0.98 · log(SNR) + 18.00, at a FWHM = 0.6 nm (for the Gaussian slit function) and SR = 4. The relationships are mostly independent of other specifications (e.g., horizontal and temporal resolutions), as they affect ε primarily through SNR, providing constraints in determining the optimal SNR (and alternatively FWHM and SR) for similar UV/visible sensors dedicated for AQ studies.