气动谐振管内压力振荡数值计算

为了研究气动谐振管内的加热现象,采用二维轴对称Euler方程组和显式MacCormack格式对谐振管内的非定常振荡流场进行了数值模拟.针对不同压力的驱动气体进行了计算,得到了谐振管底端压力振荡时间历程以及谐振管内其它气动参数的变化规律,确认了谐振管内的压力振荡是由于激波和膨胀波交替进出谐振管引起的,而且与驱动气体的压力相关.计算所得的谐振管底端压力振荡曲线与实验结果相一致,表明了所采用的数值方法对分析谐振管内的振荡流动过程是适用的.     

气动谐振管加热机理数值研究

放置在高速喷流中的半开口气动谐振管在一定的条件下会出现加热现象.针对这一现象,利用数值计算方法计算了一维非定常非线性Euler方程组,采用MacCormark格式,对不同管口条件下谐振管内的振荡过程进行了数值研究.研究发现,当管口压力变化频率处在谐振管的线性谐振频率附近时,管内形成的激波以及激波的反射增强导致了谐振管端壁区域明显的加热效果.对不同的驱动气体和不同的谐振管型式(圆柱形和圆锥形)的数值研究也给出了有价值的结论.研究结果与实验事实或现象是相吻合的.      

环形喷嘴气动谐振加热性能

为了研究环形喷嘴的气动谐振加热性能,采用高精度高分辨率的NND格式差分求解二维轴对称雷诺平均Navier-Stokes方程,对环形喷嘴-谐振管系统气动谐振加热过程中谐振管内振荡流动过程进行了数值仿真,并对不同面积的环形喷嘴-谐振管系统进行了数值模拟.研究结果表明:环形喷嘴与圆喷嘴具有相似的谐振加热规律和流场特征,即利用环形喷嘴同样可在谐振管内产生强烈的高频激波振荡.随着环形喷嘴面积减小,每一个谐振周期中的谐振温升逐渐减小.而其能够产生强烈谐振的间距小于圆喷嘴-谐振管系统的间距,在小间距时可以在很宽的喷嘴入口压力范围内产生强烈的谐振,这为气动谐振点火器的结构小型化和工程实用性提供了理论依据.     

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

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

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

航天飞机气动加热计算

介绍了国内外预测航天飞机气动加热的工程计算方法.给出了轴对称、非轴对称和机翼前缘等各种外形驻点热流密度的计算,分析和比较了平板、锥体的层流和湍流热流密度计算的各种方法,概述了航天飞机迎风面中心线和离开中心线横向热流密度的计算.计算结果表明,本文方法用在航天飞机气动热环境初步设计中比较简单并且有足够精度.     

瞬态气动模拟加热控制中的快速高精度"E-T"转换

为了准确模拟飞行器在高速飞行时的瞬态气动加热状态,必须使用快速、高精度的计算机瞬态热能控制系统,对气动模拟试验的加热过程,实行快速、高精度的非线性动态控制.为此,传感器的快速、高精度"E-T"转换是一个必须解决的非常重要的问题.提出一种高速飞行器瞬态气动加热控制系统中传感器的快速、高精度"E-T"转换方法.该方法具有计算简单、转换速度快、校正精度高的优点,使用该方法实现了高速飞行器气动加热过程中温度场高速变化状态下的瞬态非线性动态控制.     

南京地区斜向加热电离层效应数值模拟

与垂直加热相比,斜向加热电离层具有灵活性高、影响范围广和便于实际操作等优势.在非偏区考虑欧姆吸收,基于电子能量方程和连续性方程构建电波斜向加热低电离层的物理模型,并利用IRI-2007和NRLMSISE-00经验模型提供的背景参数对南京地区斜向加热低电离层进行数值模拟,对比不同加热条件下电子温度和电子密度的扰动情况.研究结果表明：电波加热效果随入射仰角和有效辐射功率的增大而增大;电子温度和电子密度增幅随电波频率增大而减小;X波模比O波模造成的电子温度扰动幅度和电子密度扰动幅度更大,同时X波模比O波模能更快地使电子温度和电子密度达到稳定状态;一定范围内较小仰角、较低频率、较大有效辐射功率的电波能使电子密度更快达到稳定,后两者还能加快电子温度达到稳定的过程;电子温度达到稳定所需时间随入射电波仰角呈单峰变化,仰角为62°时达到最大.

     

低功率氮氢电弧加热发动机非平衡数值模拟

文章对低功率氮氢电弧加热发动机进行了双温度化学非平衡数值模拟研究,模型中包含总的能量方程和电子能量方程,等离子体组分包括分子、原子、离子和电子等7个组分,采用的化学动力学模型中包含了氮氢组分解离、电离等重要的动力学过程,气体的物性根据当地的组分和温度实时计算。通过计算获得了发动机内部气体温度及各组分数密度分布。结果表明,发动机轴线附近等离子体接近热力学平衡,而在发动机阳极壁面电弧贴附区域等离子体明显偏离热力学平衡;计算获得的组分分布表明电弧加热发动机内存在反混合过程,即发动机内各组分分布与入口浓度分布明显不同。氢组分由发动机中心到阳极壁面沿径向呈现先减小后增大的趋势;而氮组分的浓度分布趋势与氢组分相反;进一步的分析表明,发动机内各组分的扩散主要受到气体解离和电离过程引起的浓度梯度所驱动。     

Euler方程数值模拟绕悬停旋翼桨叶的流动

本研究应用格心格式有限体积龙格－库塔时间推进法求解三维欧拉方程,模拟绕旋翼桨叶的流动,民流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角,选择Ｏ－Ｈ型式生成绕桨叶贴体的弦向和展向网格,通过算为进一步寻求有效的Ｅｕｌｅｒ或ＮＳ方程数值解积累了经验。     

气体二次喷射推力向量控制数值仿真

采用基于Favre平均的三维N-S方程和k-ε湍流模型对固体火箭发动机气体二次喷射推力向 量喷管复杂干扰内流场进行数值模拟;采用三阶精度空间差分格式和隐式Jacobi时间点迭代 方法进行求解;借助数值模拟技术对气体二次喷射推力向量喷管方案进行了初步探索,研究 了不同喷射条件对气体二次喷射流场特征及侧向控制力的影响作用;研究表明,二次喷射气 体喷射孔位置、总温及总压等因素对侧向力影响具有耦合作用,合适的喷射孔位置和高温燃 气以及较大喷射总压都能有效增加侧向力及向量角;燃气二次喷射系统具有较高的效率;侧 向力随二次喷射流总压增加而线性增大.      

潜艇应急燃气吹除过程的数值仿真及实验验证

进行了潜艇应急燃气吹除系统的小比例模型原理实验,模拟了水下100 m深度时燃气吹除的排水性能与规律以及燃气吹除过程中的主要性能参数变化情况.建立了模拟小比例实验吹除过程的物理模型,应用结构化六面体网格对其两相流体流动区域进行了计算区域的离散化,并对100m海水深度下的燃气吹除过程进行了数值仿真计算,得到了对燃气泡成长规律的描述,深化了对排水过程机理的认识.仿真结果与实际情况的实验结果十分吻合,混合气体压力和排水总体积的相对误差均不超过5%,验证了数值仿真方法的可行性.     

双组元统一推进系统减压器稳定性仿真

针对某双组元统一推进气路系统采用的高压卸荷膜片式减压器,建立起有限体积瞬变仿真模型,并搭建起减压器特性研究系统的模块化仿真模型.针对试验过程中出现的振荡现象,借助数值手段阐明了各结构和控制参数对减压器工作稳定性的影响,得出了减弱振荡的各种措施:减小阀芯质量,增大高、低压腔体积,增大弹性元件材料的阻尼系数,减小反馈孔直径,增大膜片刚度.与试验和相关文献研究情况的对比证实了某些措施的有效性.     

俯仰阻尼导数分量的CFD数值模拟

提出一种直接求解直接阻尼导数的方法,该方法不仅适用于轴对称外形,也适用于非轴对称外形.数值模拟飞行器的非定常强迫沉浮运动和强迫角振动,并在Etkin的非定常气动力模型基础上,辨识得到飞行器的洗流时差导数和俯仰阻尼导数,研究飞行器俯仰通道各阻尼导数的数值计算方法.对弹道外形和基本带翼导弹标模外形及Hyflex升力体外形进行研究分析表明,将阻尼导数分量相加得到的俯仰阻尼导数与直接求解强迫角振动得到的俯仰阻尼导数与试验结果吻合很好,各阻尼导数分量随质心位置的变化趋势也与理论预测相符.对于带翼飞行器,超声速条件下,洗流时差导数在俯仰阻尼导数中占主导作用.      

上游双扰动对圆柱气动特性影响的数值模拟

通过用商用计算流体力学软件FLUENT 6.0中的SIMPLE方法求解二维非定常层流模型,对雷诺数Re=200时上游双扰动体对下游圆柱气动特性影响进行了数值模拟研究.计算结果表明当双扰动体相隔较近(H/d=1.1,1.33)时,其后形成单一的Karman涡街,它们对下游圆柱的干扰作用类似于单扰动体时的情况,但是流动结构发生变化的临界间距Lcri较之有效宽度相同的单扰动体的要大,因此可以推断出在该种布置下双扰动体和圆柱绕流会出现2种流动模式.在空穴流动模式下面,下游圆柱的阻力和升力脉动最多减小94%和79%,系统阻力(包括圆柱和扰动体)也能减小33%.而当双扰动体尾流为双Karman涡街时(双扰动体相隔较远H/d=3.33,5.0),圆柱迎风面部分区域仍受到来流的冲击作用,气动力减小有限.     

航空拖曳诱饵空气动力及动态特性

航空拖曳式诱饵与被保护载机的相对空间位置及投放使用方法是诱饵能否成功干扰导弹的关键因素.针对传统有限段法的不足,提出改进的拖缆有限段法,将拖曳式诱饵的动态运动特性问题转变为拖曳线与拖曳飞机航迹轴系的夹角动态计算问题.通过分析拖曳飞机的速度矢量与拖曳线的几何关系,建立了拖曳式诱饵对拖曳载机的动态跟随模型.仿真计算了拖曳飞机机动飞行时,拖曳式诱饵的动态运动特性,结果证明算法简单有效,较好地模拟了其动态跟随特性.      

贮箱内低温推进剂汽化过程的CFD数值仿真

为研究贮箱内低温推进剂相变对推进剂温度和贮箱压力的影响,对贮箱内的传热传质过程进行了仿真.仿真涉及的物理过程包括贮箱与外界环境的换热、推进剂的自然对流、推进剂与贮箱内壁面的换热以及低温推进剂的相变过程等.根据热力学平衡原理建立了低温推进剂相变模型,使用CFD(Computational Fluid Dynamic)方法对处于地面常压停放状态的液氢贮箱进行了450 s的仿真.研究表明随着贮箱壁面传热过程的稳定,推进剂的温度分布、流动状态以及相变情况会趋于稳定;通过仿真获得了推进剂单位时间的汽化量;影响相变的主要因素是贮箱壁面漏热以及推进剂自身的对流运动.