一种利用正交流线坐标系求解跨音速管流的新方法

本文以正交流线坐标系为工具,采用位流模型,导出了正交流线坐标系下平面或轴对称流动的控制方程组,利用流函数与位函数的存在性,进而导出了两个首次积分。设计了一套求解匀熵管流的分凝迭代方案,模拟了管内无激波的平面或轴对称喷管的跨音速流动。数值实验表明该方法能够处理小喉道曲率半径平面或轴对称喷管跨音速流动,与已有的数值计算结果及实验结果符合较好。     

2000年前固体推进剂的发展趋势

固体推进剂的发展趋势与战术和战略导弹的发展有密切关系。战略导弹对推进剂的要求就战略导弹而言,研究推进剂的比冲和密度比冲是十分必要的,这项研究工作应与降低结构系数的其它技术(如采用高强度比的纤维缠绕壳体和先进的喷管和防热层材料)和核弹头的小型化研究相配合。装药比冲提高以后,可以针对不同任务,进行灵活     

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

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

含非均匀辐射气体空腔辐射特性计算研究

文章发展了求解积分辐射传输方程组的方法来计算充满非均匀辐射性气体封闭腔的辐射特性的方法,讨论了辐射空间分布及光谱分布的计算方法,并编制了相应的计算机程序.计算了某型发动机轴对称喷管的辐射分布,并与采用通常的辐射换热经验公式的计算结果进行了比较.     

进气旋流对推力喷管性影响的研究

介绍了模拟近代大涵道比发动机排气状况的旋流模型,和在喷管落压比1.2到6.5范围内,进气旋流对轴对称喷管,矩形收-扩喷管,楔形二元喷管,和单边膨胀斜板形喷管等性能和射流流态影响等的一些实验结果。实验数据指出:进气旋流可使喷管出口射流柱中的旋涡加强,射流与外界气流间的混合区扩大,中心核心流区缩小,但它对轴对称喷管出口流态的影响很小。进气旋流同时对喷管推力系数和流量系数均起不利影响,增大进气旋流角度会使喷管推力系数和流量都下降,特别是当旋流角超过某一“临界”值后二者降低很急剧,并且二元喷管降低的程度比轴对称喷管更严重,因之在具体应用情况中,保留多大发动机的旋流,要针对具体情况综合考虑各种因素来选定。     

自由分子流微电热推力器流动模拟与性能预示

为了精确的预示自由分子流推力器(FMMR)在不同工质、工况下的性能,采用直接模拟蒙特卡罗(DSMC)方法模拟FMMR内气体流动,计算了不同工质气体、工作压强和加热温度条件下的FMMR推力、比冲等性能,分析了工质、工作压强和加热温度对FMMR性能的影响。计算和分析结果表明,工质气体的选择不影响FMMR的推力性能,但随着工质气体相对分子质量的减小,推力器比冲迅速增加;在总压50~500Pa,温度300~600K条件下,推力器的比冲随总压和加热温度的增大而增加,推力随加热温度增大而减小,随总压增大而增大。     

火箭及导弹的底阻计算方法及其CRMBP程序

本文叙述了作者开发的基于Chapman-Korst理论的超音速底部压强计算方法及其计算程序CRMBP。该方法及程序计及了喷流、喷管和火箭后体多种参数的影响,以及羽流诱导分离情况。计算结果与实验结果符合得较好,特别是冷喷流情况。本文的方法及其程序可供设计部门作设计计算及参数分析使用。     

引射喷管的气动性能和传热计算

对某型飞机引射喷管装三种机尾罩时的推力性能和机尾罩温度进行了计算，主流和次流的特性分别用特征线法和一维方法计算，为了计算机尾罩的温度，建立了一种简化的一维模型，计算结果与试飞实测数据吻合良好，引射喷管装三种机尾罩的比较表明：切短的机尾罩可以在推力下降不太多的同时改善机尾罩温度。     

轴对称羽流流场数值模拟和实验验证

针对在轨航天器羽流效应对航天器整体寿命的影响,在"DFH"系列卫星姿控发动机喷管羽流实验研究的基础上,建立喷管的轴对称模型,选取矩形计算区域,设定一定的边界条件和计算条件,采用轴对称单组分直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法对相应实验状态的羽流流场进行数值模拟研究.通过对比数值模拟结果和实验测量结果,表明DSMC数值模拟羽流流场的计算结果符合实际流动规律,实验验证了所编制的轴对称DSMC数值模拟程序的正确性和有效性.      

谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.     

第二相质点在LC9超硬铝晶粒细化过程中的作用

本文研究了LC9超硬铝中三种质点,即析出相质点、杂质相质点和弥散相质点在晶粒细化过程中的作用。研究发现,析出相质点的存在是晶粒细化的关键;细小弥散质点对稳定细晶组织有决定性的作用;粗大杂质相对晶粒细化作用不大。定量金相分析表明,直径大于0.75μm的析出相质点的密度控制着再结晶晶粒的密度,0.75μm为析出相质点激发再结晶形核的临界尺寸。     

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

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

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

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

基于支板燃烧室的喷管化学非平衡效应

采用有限体积法全隐式格式和代用燃料C₁₂H₂₃的10组元13步化学反应Arrhenius有限速率模型研究煤油燃料超燃冲压发动机单边膨胀喷管(SERN, Single Expansion Ramp Nozzle)内的化学非平衡流动,通过建立支板燃烧室——喷管模型有效解决了单边膨胀喷管模型的"入口薄层"问题.计算结果表明,整个单边膨胀喷管内,流动呈现化学非平衡效应,喷管入口附近区域尤其显著;非平衡流动喷管性能明显高于冻结流动,随发动机当量混合比 ε增加,非平衡流动的喷管推力系数和升力系数相对冻结流动的百分比增量 δ不断升高,当 ε=0.8时,推力系数百分比增量 δ_CF达到9.41%,升力系数百分比增量 δ_CY达到16.39%,化学非平衡效应对煤油燃料超燃冲压发动机尾喷管性能的影响不可忽略.     

三维CT重建几何参数的非线性最小二乘估计

重建几何参数的精度是3D-CT(Three Dimensional Computed Tomography)重建中严格控制的指标,在实际的X射线扫描成像中,这些参数无法直接测量得到,从而难以保证其精度.介绍了利用非线性最小二乘估计Feldkamp三维重建几何参数的方法,根据空间质点质心的投影位置与其投影的质心位置重合的原理,通过计算空间质点(实验中采用近似质点的目标体)投影质心坐标,求解非线性方程组估计重建几何参数的最优解.计算机模拟结果表明,在参数向量初始值接近于真值的情况下,估计所得参数值有较好的精度和重复性.利用此方法估计所得实际扫描系统的几何参数值进行三维重建,也得到了很好的重建结果.     

使用GAO-YONG湍流方程组 对顶盖驱动方腔流的计算

运用SIMPLER方法、QUICK格式及交错网格技术求解GAO-YONG湍流方程组,对二维顶盖驱动方腔流进行数值模拟,得到方腔内的涡量分布和方腔中心水平剖线上的速度变化曲线,并与DNS(Directly Numerical Simulation)结果进行定性对比,有较好的一致性.结果表明,GAO-YONG湍流方程组中与湍流多尺度现象相对应的机械能方程对此湍流方程组能够合理地模拟出顶盖驱动方腔流中复杂的湍流粘性分布以及湍流能量逆转现象起到了关键作用.研究进一步证明了GAO-YONG湍流方程组能够得到复杂流动的真实粘性场.     

引射混合器内流场的数值分析

计算了以波瓣为主喷管的三维引射掺混流场；分析了包括复杂型面和引射二次流在内的边界条件处理；再现了流向涡的产生，计算结果与实验结果作了对比分析，吻合程度较好。     

轴对称推力矢量控制伺服机构空间运动解耦研究

针对由发动机摆动喷管矢量控制的两个伺服机构间的牵连效应导致的伺服机构运动耦合问题,提出了一种新型空间解耦计算方法。该方法利用欧拉角描述喷管姿态,通过空间齐次坐标变换矩阵原理对两路喷管伺服机构进行解耦计算,得到了两路伺服机构的伸长量解析解,确保了控制指令的精确,并通过ADAMS虚拟样机仿真验证了该计算方法的正确性。仿真实验结果表明,此计算方法与虚拟样机的仿真结果完全吻合,计算精度高,可以对喷管任意姿态进行运动精确控制。     

多单元直排塞式喷管发动机性能

为了了解优化设计的塞式喷管的性能及燃气流动中热力学参数变化对性能的影响,对比钟型喷管与塞式喷管的高度特性,从曲线坐标下的三维平均雷诺N-S方程和Euler方程出发,采用LU时间隐式格式、MUSCL空间离散方法,发展了模拟塞式喷管三维流场的数值程序.计算了从喉部圆转方内喷管的性能,比较了冻结和变化热力学参数对塞式喷管性能的影响及塞式喷管与相同面积比钟型喷管的高度特性曲线.计算结果可为塞式喷管的设计研制提供参考.      

单元二维直排塞式喷管冷流实验

为研究塞式喷管的高度特性和底部特性,采用高压空气为工作介质对单元直排塞式喷管进行实验.研究了底部盖板、底部二次流对性能的影响和塞锥壁面压强分布.实验结果表明:无底部盖板可提高塞式喷管低空性能2%~7%;底部二次流可提高底部压强,减少底部开闭状态转变过程带来的推力突降;底部二次流流量以1%~1.5%为宜,过大将引起喷管性能下降.本实验喷管设计点效率均超过99%,部分实验接近100%,高度补偿效果明显.