在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

基于最优误差动力学的时间角度控制制导律

针对带有攻击时间和终端攻击角度约束的导弹制导问题,设计了一种基于最优误差动力学的时间角度控制制导律,并给出了明确的性能指标。推导广义最优角度控制制导律作用下的剩余飞行时间估算表达式,在广义最优角度控制制导律的基础上增加攻击时间误差反馈项,将攻击时间误差看做跟踪误差,设计的制导律使跟踪误差以最优模式在有限时间内收敛到零,最终实现攻击时间和终端攻击角度的共同控制。对不同参数下的情况进行仿真,验证了所提制导律的有效性。     

超燃冲压发动机燃烧室的准一维计算与分析

在考虑有限速率化学反应的准一维Euler方程基础上，通过增加截面面积变化、壁面摩擦和添质的源项，发展了适用于超燃燃烧室性能分析的准一维计算方法。依次以中国空气动力研究与发展中心（CARDC）和日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为验证算例，分别采用传统的一步反应模型和发展的有限速率反应模型，模拟了燃烧室流场，并基于NAL燃烧室，计算分析了不同当量比和进口压强对燃烧室流动特性的影响。结果表明：两种方法都能得到与实验数据吻合良好的结果；和一步反应模型相比，有限速率反应模型不仅可以更细致地捕捉流场细节，而且能够初步分析化学非平衡效应的影响；对于NAL燃烧室，当量比≥0.6时，压强随当量比的升高而增大，当达到1.0时，反压已推进隔离段，且推进速度随当量比增大而增加；进口压强不大于110.444kPa时，反压随进口压强增大而升高，且当反压不小于82.833kPa时，反压被隔离在等直段燃烧室入口处；过小的当量比和过大的进口压强均会导致燃烧室出口马赫数严重下降，甚至出现亚声速出流状态。     

流-固-推进耦合下的机体/推进一体化性能分析

吸气式空天飞行器的一体化性能随扰动变化的敏感性高,在高马赫数飞行条件下,有必要开展流-固-推进耦合性能分析。针对机体/推进一体化布局的吸气式飞行器,明确一体化部件之间的耦合关系和耦合问题,利用CFD、有限元和准一维流方法,结合本征正交分解(Proper Orthogonal Decomposition,POD)降阶手段,建立吸气式空天飞行器流-固-推进多物理场耦合快速分析方法,并开展多场耦合特性分析。结果表明:(1)进气道压缩面的流-固耦合导致出口静压的最大振荡振幅约为平均静压的21.6%,而出口马赫数的最大振幅约为平均马赫数的8.45%。(2)进气道出口性能的振荡会影响发动机的推力性能,导致推力振荡幅值可达平均值的31%,且随着时间的推移,会在进气道外压缩流场产生大量的气动涡,涡结构进入进气道后会导致进气道出口性能的持续下降,进一步削弱了发动机的平均推力性能。     

多机构比对融合的分布式InSAR编队星间基线确定

星间基线高精度确定是分布式干涉合成孔径雷达（InSAR）系统完成科学任务的重要保证,受星载全球定位系统（GPS）接收机连续跟踪弧段短、个别弧段共视GPS卫星个数少或模糊度固定成功率低、频繁轨道机动等因素影响,分布式InSAR高精度基线确定仍有不可靠的风险。通过多机构产品互比来识别基线精度较差的时间段,降低不可靠风险,并通过多机构产品融合进一步提高基线精度。选用重力反演与气候实验（GRACE）卫星数据进行实验,国防科技大学（NDT）和西安测绘研究所（CHS）采用不同的基线处理软件和简化动力学策略,保证了各自的基线产品具有一定的独立性。实验表明,多机构互比对可以有效识别基线精度较差的时间段,NDT和CHS的基线产品之间具有很好的一致性,互比对残差的均方根（RMS）在R、T、N方向分别为0.7、0.9、0.7 mm,二者之间没发现明显系统偏差,大约97.86%的基线三维互比对残差量级在2 mm以内。两个机构基线产品融合后发现可进一步降低基线产品中的随机波动误差,K/Ka波段测距（KBR）系统校核结果表明融合基线产品精度较NDT基线产品提高8.97%,较CHS基线产品提高29.21%。     

热容燃烧室基于单点温度的热流测量方法研究 *

为了研究推力室内壁面与燃气间的传热特性,获得多喷嘴燃烧室内壁热载分布,从而为推力室设计及热结构分析提供参考,开展了热容燃烧室壁面热流测量方法研究。以位于距内壁一定距离位置的单点测量温度作为输入条件,文章提出了三种热流计算方法,并进行火箭发动机试验研究,得到了对应位置燃气对内壁的热流密度。可以看出,运用单点方法得到的热流密度与传统两点法结果符合较好,两者计算结果误差在10%以内。结果显示:不同点火时长工况下在同一位置相同时间点上重复性很好,在燃烧室圆柱段热流密度随着点火时长的增加而减小,而在喷管段变化规律则相反,热流密度随时间增加而增加。在提出的计算平均热流的三种方法中,热积累方法的计算结果最高,热平衡方法次之,瞬态法最低。在火箭发动机热容试验中,本文提出的方法可以用于测量燃烧室内壁面热流密度。     

再生冷却推力室热障涂层系统的热结构分析

为了设计和优化适用于液体火箭发动机推力室的热障涂层,应用ANSYS的热-结构分析功能,对再生冷却推力室-热障涂层系统进行了热结构有限元分析,得到在不同涂层覆盖下,推力室壁中的温度场和应变场,并通过对热障涂层中应变场的分析,研究不同涂层发生分层剥落的关键位置以及主要驱动力。结果表明,陶瓷层厚度较大的YSZ+Ni Cr Al Y涂层拥有更优异的性能,使推力室壁在热试阶段的最大应变量减少约36.1%;工作循环中,涂层与推力室壁的接触面上会产生较大的应变量,最终有可能导致涂层剥落失效;粘结层能缓解涂层与推力室壁间的热膨胀系数不匹配,使陶瓷层在热试阶段的最大应变量减少约80%。     

Numerical simulations of turbulent flows in aeroramp injector/gas-pilot flame scramjet

To uncover the internal flow characteristics in an ethylene-fueled aeroramp injector/gaspilot(ARI/G-P)flame scramjet,a Reynolds-averaged Navier-Stokes(RANS)solver is constructed under a hybrid polyhedral cell finite volume frame.The shear stress transport(SST)k-x model is used to predict the turbulence,while the Overmann’s compressibility corrected laminar flamelet model is adopted to simulate the turbulent combustion.Nonreactive computations for Case 1(G-P jet on),Case 2(ARI jets on),and Case 3(both ARI and G-P jets on)were conducted to analyze the mixing mechanism,while reactive Cases 4–7 at equivalent ratios of 0.380,0.278,0.199 and0.167 respectively were calculated to investigate the flame structure and combustion modes.The numerical results are compared well to those of the experiments.It is shown that the G-P jet plays significant role in both the fuel/air mixing and flame holding processes;the combustion for the four reactive cases takes place intensively in the regions downstream of the ARI/G-P unit;Cases 4 and 5are under subsonic combustion mode,whereas Cases 6 and 7 are mode transition critical and supersonic combustion cases,respectively;the mode transition equivalent ratio is approximately 0.20.     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

高温处理对中密度针刺C/C复合材料性能的影响

针对C/C复合材料脆性问题,对密度为1.60 g/cm3的碳布叠层针刺C/C复合材料进行了1800、2 000、2 200和2 500℃的高温处理,研究了不同热处理温度对C/C复合材料微晶结构、力学和抗热震等性能的影响.结果表明,高温处理使针刺C/C复合材料的层间剪切和面内拉伸强度出现不同程度的降低,但材料的断裂伸长率和抗热震性能得到大幅度提高.其中,1 800℃高温处理后的C/C复合材料具有优异的力学和抗热震性能.