超燃试验用加热器三维粘性流场数值模拟

超燃冲压发动机试验时,使用燃烧型加热器模拟高焓、高马赫数的来流条件,而加热器的喷注结构对改善加热器出口流场均匀性具有重要作用.使用三维数值模拟软件Fluent对2种喷注结构的燃烧型加热器进行了三维定常反应流场模拟,获得了设计工况下的流场.结果表明,采用多孔斜喷结构能使氢气在横向有较大的穿透深度,加上空气的引射作用在燃烧室头部产生大尺度横向涡,增强了氢气和空气的混合燃烧以及燃气的掺混,使加热器出口温度场不均匀性从53%降为21%,氧气摩尔分数不均匀性从63%降为17%,更有利于开展超燃冲压发动机试验.     

超音通流风扇叶栅流场数值模拟

用矢通量分裂法求解Ｎ－Ｓ方程和Ｅｕｌｅｒ方程，对超音通流风扇叶栅流场进行了数值模拟，由计算结果给出的物理图象表明，叶栅通道内没有强激波，亦未出现分离回流，全叶栅流场均为超音速，说明计算的两种叶型符合设计要求，其性能有待实验检验。     

可变桨距冲压空气涡轮混合型流场数值研究

冲压空气涡轮流场数值研究主要目的是模拟在实际高空加油工作状态下涡轮变桨距气动特性.数值风洞实验采用代数Baldwin-Lomax湍流模型,雷诺平均Navier-Stokes方程,模拟可变桨距冲压空气涡轮全三维混合型流场,并分析其流场主要气动特性,探讨桨叶载荷分布原因,然后将数值模拟结果与风洞实验数据进行了比较.改进桨叶近轮毂区域气动性能可进一步提高冲压涡轮载荷.      

进气道结构对固体冲压发动机补燃室燃烧及内壁流场的影响

为研究进气道结构对固体冲压发动机补燃室燃烧及内壁烧蚀的影响,采用标准k-ε湍流模型,单步涡耗散燃烧模型与KING硼粒子点火燃烧模型,开展了双下侧90°进气结构和双侧180°进气结构固体冲压发动机补燃室内燃气燃烧数值模拟,对比分析了补燃室燃气燃烧流场特征和内壁烧蚀环境特征。结果表明:双侧180°进气结构在补燃室中形成大漩涡,有利于燃气与空气的掺混燃烧,至补燃室出口位置,总燃烧效率超过90%,且该结构有效减少了粒子对内壁的冲刷侵蚀;在双下侧90°进气结构补燃室中,凝聚相粒子和燃气贴近补燃室一侧运动,导致氧气浓度和温度分布不均,不利于燃气的掺混燃烧,总燃烧效率为74%,在远离补燃室进气道一侧形成高温热烧蚀、高浓度粒子侵蚀、高速射流冲刷和热应力集中的综合破坏;双侧180°进气结构的固体冲压发动机补燃室总体性能优于双下侧90°进气结构的冲压发动机补燃室。      

超声速燃烧冲压发动机进气道起动性能研究

采用等激波强度设计方法,并考虑变比热、激波与附面层干扰等因素的影响,对唇口平直和唇口带有斜楔的超燃冲压发动机二维混压式前体/进气道进行了初步设计,比较分析了几种方案进气道的设计点和非设计点性能,研究表明,在低飞行马赫数(Ma)下,唇口带有斜楔的前体/进气道起动性能和总压恢复优于唇口平直的,在高飞行Ma下,唇口平直的前体/进气道冲压比高、外罩阻力小,而唇口带有斜楔的前体/进气道总压恢复系数高,外罩阻力相对较大。针对超声速燃烧冲压发动机燃烧室和进气道间非定常干扰的问题,计算研究了飞行Ma=4,6下,燃烧室压力升高对进气道/隔离段流场和起动性能的影响,结果表明,在低飞行Ma条件下,燃烧引起的压力扰动容易往上游传播,甚至引起进气道不起动;随着飞行Ma的增大,隔离段的抗扰动能力是增强的;当进气道进入不起动后,进气道的捕获流量和总压恢复系数急剧下降,高飞行Ma时的捕获流量的下降幅度比低飞行Ma时大。     

背负式S形进气道流场控制技术

通过在已有的S形进气道入口第一弯底部布置一系列类机翼的导流叶片,来改善出口流场分布,抑制旋流,减小出口畸变指数.分析了导流叶片的高度、弦长、攻角和安装位置(轴向位置和横向位置)这几个参数对S形进气道内流的影响,并利用计算流体力学(CFD)方法针对导流叶片的组数展开研究,得到最优结果.结果表明,采用优化后的导流叶片,可以有效抑制内流分离,降低S形进气道出口畸变指数,改善流场品质. 这些分析结果对于涡流发生器的设计有一定的指导意义.      

超燃冲压发动机燃烧室设计计算方法的研究

为寻找一种适当的计算超燃冲压发动机燃烧室性能的方法并评估现有模型的优劣,提出一维化学动力学模型,且通过几个算例验证了该模型的可靠性.为研究设计过程中的性能计算方法的适用性,针对一具体的燃烧室,采用目前通用的性能计算方法,即冲量分析法、一维化学动力学及二维化学动力学方法,计算得到燃烧室内各气动参数的分布曲线,并得到上述方法在超燃冲压发动机燃烧室设计过程中性能计算适用性的初步结论.计算结果表明:一维化学动力学方法与冲量分析法都能够考虑到摩擦、通道面积变化以及燃烧释热的影响因素,具有较好的适应性.对于冲量分析法,在考虑燃烧时,还需设定放热规律; 而一维化学动力学方法则可以利用化学反应模型,会自动计算释热规律,具有更大的独立性.和前两种模型相对比,二维化学动力学方法可以更细致地捕捉到流场中的一些细节,但此种模型需要较长的运算时间.对比这几个模型,各具有不同特点,在超燃冲压发动机的设计与性能计算过程中,需综合考虑上述区别,最充分地发挥各模型的优势.     

轴对称降落伞小迎角稳定下降时流场特性

根据降落伞的特点,通过伞衣零厚度假设、伞轴对称假设和流场定常假设,建立轴对称降落伞的流体力学计算模型。在0°~5°的小迎角范围内,求解RNG(Renormalization Group)k-ε湍流模型下的N-S方程组,获得与有关单位试验相吻合的数值模拟结果。分析发现0°迎角时轴对称降落伞截面流场中鞍点与鞍点直接连接的状态是不稳定的。随着降落伞迎角增大,奇点类型和奇点间的连接方式将发生改变,降落伞物面附近的涡环里将出现极限环,截面流场中将出现新的鞍点和结点。迎角继续增大,极限环的面积将会减小直至消失,结点将附着在物面上,成为半结点。结果表明小迎角范围内降落伞截面流场的拓扑结构符合拓扑规律。文章揭示了小迎角范围内轴对称降落伞稳定下降时流场特性的演变规律,为进一步研究降落伞流场的流动机理和流固耦合问题打下基础。     

降落伞流固耦合问题的数值模拟和流场分析

提出了基于压力耦合的半隐式算法,即SIMPLE (Semi-Implicit Method for Pressure-Linked Equations)算法和生成阶梯网格方法的对稳定下降阶段的降落伞进行数值模拟的新方法.此方法在降落伞流固耦合计算中具有较好的稳定性和较高的效率.对降落伞稳定下降阶段的流固耦合问题进行了数值模拟.降落伞模型分别采用了圆形伞和锥形伞稳定下降时的结构数据,并对锥形伞模型进行了不同来流攻角情形下的模拟.与平行有限元方法相比,采用SIMPLE算法的新方法得到了更准确和更合理的结果.此新方法可以作为传统数值模拟方法的一个可靠替代和重要补充. 　　      

肼电弧喷射推力器工作方式可替换性数值研究

在正确理解电弧喷射推力器工作过程物理机理的基础上建立了其简化的物理模型,并基于此物理模型利用耦合电磁源项的N-S方程组来表达其数学模型,对相同结构尺寸的低功率肼电弧喷射推力器加电弧和不加电弧时的工作性能进行了比较计算.采用二阶精度无波动、无自由参数的耗散差分格式(NND格式)和显式时间推进法求出数值解.计算结果给出了有电弧和无电弧时推力器的流场结构,预计了不同工作方式及不同工作条件下推力器的推力和比冲,证实了推力器工作方式的可替换性.      

强磁场结构在行星际空间膨胀的数值模拟

本文研究了强磁场结构在18—240R_s的行星际空间的膨胀效应。模拟结果表明, 强磁场结构的膨胀能够形成典型的磁云剖面结构;膨胀速度明显依赖于磁场强度的大小;强磁场结构运动的速度直接影响到它的动力学特性。     

卫星表面电场尖端效应的数值计算

本文使用模拟电荷法和扁椭球体的电势系数法较好地解决了尖端效应的数值计算问题,并用这一方法,计算了卫星表面尖端的效应.     

热和振动复合环境中人体热调节的数值模拟

本文针对热和振动复合环境对人体的影响特点，建立了该复合环境中二维人体热调节系统的数学模型．采用数值方法对人体热调节系统数学模型进行求解，计算了不同热和振动复合环境条件下人体温度分布及动态响应，计算结果与实验结果相符较好，相对误差小于５％．     

航空发动机一维数值仿真

分析了航空发动机一维气动热力数值仿真模型,此模型能够用于航空发动机整机数值仿真.进行了整机性能数值仿真,求解过程考虑了粘性效应以及气体流动损失对于整机性能的影响;采用了具有TVD性质的高阶精度Godunov格式,此格式具有良好的数值稳定性和激波捕捉能力;应用了时间推进法,交替使用显/隐格式.计算了某型带加力的混排涡扇发动机的速度特性以及节流特性并且与设计数据进行了比较分析,数值仿真结果与设计数据比较吻合.分析了典型工作点的计算过程,模型收敛速度较快.     

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

无振荡、无自由参数格式在三维太阳风流动数值模拟中的应用

从快速、有效地进行空间天气数值预报的需要出发,针对1998年5月份行星际太阳风暴事件,采用全三维流体力学(HD)模型进行数值试验,考察了无振荡、无自由参数(NND)格式在三维太阳风流动数值模拟中的应用.计算中内边界的密度分布由观测K日冕数据来确定,由此得到的源表面密度分布具有和源表面电流片相似的结构.数值试验表明: 虽然该格式不需要加入人工粘性,但是具有很好的计算稳定性.在扰动计算时激波间断耗散小,在三维任何方向上间断所占的网格数比较少,没有数值色散现象.      

A CLASS OF TWO-STEP TVD MACCORMACK TYPE NUMERICAL SCHEME FOR MHD EQUATIONS

In this paper, a new numerical scheme of Total Variation Diminishing （TVD） Mac-Cormack type for MagnetoHydroDynamic （MHD） equations is proposed by taking into account of the characteristics such as convergence, stability, resolution. This new scheme is established by solving the MHD equations with a TVD modified MacCormack scheme for the purpose of developing a scheme of quick convergence as well as of TVD property. To show the validation, simplicity and practicability of the scheme for modelling MHD problems, a self-similar Cauchy problem with the discontinuous initial data consisting of constant states, and the collision of two fast MHD shocks, and two-dimensional Orszag and Tang’s MHD vortex problem are discussed with the numerical results conforming to the existing results obtained by the Roe type TVD, the high-order Godunov scheme, and Weighted Essentially Non-Oscillatory （WENO） scheme. The numerical tests show that this two-step TVD MacCormack numerical scheme for MHD system is of robust operation     

无人机雷达散射截面的数值模拟及突防分析

用数值模拟的方法获得了无人机雷达散射截面的空间分布，建立了无人机突防的数学模型，并对其突防能力进行了分析。     

行星际强磁场结构中的亚阿尔文波速流

1982年11月23日开始的极高速太阳风流追赶强磁场结构的事件,形成了高密度前鞘-强磁场本体-高密度后鞘-高速流的组合结构。本文分析了此组合结构的磁流体动力学特征,发现在强磁场本体内出现行星际亚阿尔文波速流,并推论了其可能的成因。