考虑分子转动自由度的离散统一气体动理学格式

离散统一气体动理学格式（DUGKS）是一种适用于全流域模拟的有限体积方法。之前的研究考虑了分子平动自由度,验证了DUGKS在多尺度问题中的准确性及稳定性。本文基于Rykov模型方程构造了离散统一气体动理学格式,并采用Landau-Teller-Jeans转动能量松弛模型,可用于双原子气体从连续流动到稀薄流动的多尺度问题计算。测试了激波结构、超声速平板绕流以及超声速圆柱绕流等非平衡流动问题,计算结果显示出双原子气体分子中存在平动自由度与转动自由度对应的能量交换过程,并与统一气体动理学格式（UGKS）、直接蒙特卡罗（DSMC）方法的解以及实验值吻合较好。     

激波结构内流动问题的气体运动论描述

从介观Boltzmann速度分布函数理论出发,发展计及分子粘性碰撞截面与扩散碰撞截面,可描述各流域一维气体流动问题的Boltzmann简化速度分布函数方程及其气体运动论数值计算方法。通过对不同Knudsen数下非定常激波管流动及不同马赫数定常正激波结构问题数值模拟,研究分析不同流区的激波突跃变化过程以及近连续流、稀薄过渡流特有的分子输运现象,揭示不同马赫数、不同分子模型的激波内流动与传热变化规律,证实基于Bo-ltzmann模型方程的气体运动论数值计算方法用于激波结构内流动研究的准确可靠性。     

分子气体稀薄效应的动理学建模

稀薄气体效应,是指气体在特征尺度与其分子平均自由程相当的系统中流动时出现的非平衡效应.相比于单原子气体,分子气体(每个气体分子包含两个及以上原子)流动因同时具有转动、振动等多种自由度的非平衡过程,其稀薄效应更为复杂.分子气体稀薄效应在航空航天、微机电系统和页岩气开采等民生、科技领域广泛存在,而描述该效应的动理学模型与数...     

振动温度的激波管测定

根据双原子分光带系及其强度理论，介绍了测定振动温度的原理及方法，同时，利用激波加热试验气体的方法，得到氧化铝分子绿带系（Ｂ＾２Σ＾＋－Ｘ＾２Σ＾＋）的全部光带，通过氧化铝分子绿带光带强度的测定，得到在激波马赫数为１０时，反射激波后气体的振动温度约为５２００Ｋ。     

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题.研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出.     

含转动非平衡效应Boltzmann模型方程统一算法与跨流域绕流问题模拟研究

通过对转动自由度松弛变化特性的研究,采用转动惯量描述气体分子自旋运动,利用分子总角动量守恒作为一个新的碰撞不变量,确立含转动非平衡效应各流域统一的Boltzmann模型方程。基于转动能量空间对分布函数进行守恒积分,得到计及转动非平衡效应的分子速度分布函数控制方程组,应用气体分子运动论离散速度坐标法对速度分布函数方程所依赖的速度空间离散降维,构造直接求解分子速度分布函数的气体动理论统一格式。基于物面质量流量通量守恒与能量平衡关系,发展计及转动非平衡影响的气体分子运动论边界条件数学模型及数值处理方法。由此提出模拟高稀薄自由分子流到连续流各流域转动非平衡效应的Boltzmann模型方程统一算法。通过对高、低不同马赫数1.5≤Ms≤25考虑转动非平衡效应的氮气激波结构与不同Knudsen数9×10-4≤Kn!≤10再入竖直平板、Ramp制动器、尖双锥外形跨流域高超声速绕流问题模拟研究,将计算结果与有关实验数据、广义退化Boltzmann与ES椭球统计模型等研究结果对比分析,验证含转动非平衡效应Boltzmann模型方程统一算法求解自由分子流到连续流跨流域绕流问题可靠性。     

十一组元高温平衡、非平衡空气光谱吸收系数计算

飞行器以很高的马赫数再入大气层时 ,头部激波层、尾迹的气体辐射产生的紫外、可见、红外特征信号 ,是地面监测和反导制导系统探测、识别的主要依据 ;同时 ,头部气体层热辐射也是再入体壁面加热热流的重要源项。激波层内高温气体吸收系数是目标紫外、可见光、红外辐射特性计算、热防护计算和气动流场辐射场耦合计算的基本参数。本文分析了再入过程中高超声速稀薄气体流的非平衡现象 ,采用三温度模型 (电子能温度 Te、振动能温度 TV、转动和平动能温度 TR)来表征激波层内非平衡态的气体组份各个能级占有数分布 ,并由原子分子辐射理论直接计算 N2 ,N+ 2 ,N,N+ ,O2 ,O+ 2 ,O,O+ ,NO,NO+和 e-等 1 1种主要空气组元各种能级跃迁对辐射的贡献 ,最终得到平稀和非平衡高温空气吸收系数计算模型。计算结果可见光区与试验结果符合较好 ,紫外、红外区结果稍差     

斜向环形激波聚焦诱导直接起爆的大涡模拟

基于大涡模拟(LES)数值方法,利用基元反应模型对斜向环形激波聚焦诱爆氢气-空气可燃混合气直接起爆进行了详细的数值模拟.研究结果表明:在进气马赫数为2.5的条件下,斜向环形激波聚焦产生的高能区可以在可燃气体中实现直接起爆形成爆震波,并且燃料能够高效完全反应;而垂直环形激波则没能诱导出爆震燃烧过程.可见,斜向环形入射比垂直环形入射更加有利于燃料混合气聚焦诱导起爆成爆震燃烧.      

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

管道系统激波运动的数值模拟

对激波在管道，活塞腔系统内运动的两种数值模型－－一维非定常气体运动模型的二维轴对称气体运动模型，进行了数值模拟，其计算曲线与相应的实验曲线符合较好，表明基于Ｇｏｄｕｎｏｖ格式的激波运动计算方法能成功地处理含有激波且边界条件复杂的管道系统内气体运动问题，因而能较精确地预测管道系统内详细的和变化过程，特别是预测激波在管道系统内的传播过程。     

平面运动激波在定常超音速劈表面马赫反射的研究

本文介绍了在电控双驱动激波风洞中进行的平面运动激波在有定常超音速绕流的尖劈表面马赫反射的实验。这种反射现象属于准定常的,在实验中观察到了四种反射(RR,SMR,CMR和DMR)。本文还计算了三波点Τ的迹线及弯折点Κ的迹线分别和劈面的夹角χ和χ',以及各种反射相互转变的边界。发现运动激波波前定常超音速气流马赫数Μ_0会影响χ和χ'角,使(θ_ω+χ,Μ_s)平面上各种反射的转变边界发生变化,但不影响(θ_w,Μ_s)平面上的各条转变边界线。     

STUDY OF SWEPT SHOCK WAVE AND BOUNDARY LAYER INTERACTIONS

Ｔｈｅｉｎｔｅｎｓｉｖｅｒｅｓｅａｒｃｈｏｎｔｈｅｓｕｂｊｅｃｔｏｆｓｗｅｐｔｓｈｏｃｋｗａｖｅ／ｔｕｒｂｕｌｅｎｔｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎｓｈａｓｂｅｅｎｃｏｎｄｕｃｔｅｄａｎｄｓｕｐｐｏｒｔｅｄｆｏｒｍａｎｙｙｅａｒｓｂ...     

跨音叶栅三维激波损失的改进模型

在低展弦比、低轮毂比的轴流压气机中激波结构常常是三维的。叶尖区域的三维激波损失在整个流场的激波损失中占有十分重要的比例。Ｗｅｎｎｅｒｓｔｒｏｍ和Ｐｕｔｅｒｂａｕｇｈ于１９８４年推出的三维激波损失模型忽略了激波在叶尖区域的流动结构。通过分析叶尖区域的激波结构及其与机匣附面层流动的相互作用,提出了考虑叶尖区域激波与附面层相互作用所引起的激波结构及强度变化的跨音叶栅三维激波损失的改进模型,应用改进模型计算得出的结果和实验结果符合的很好。     

球面激波对平面的马赫反射的简化处理

本文利用Whitham绕射理论,给出了球面激波对平面绕射时马赫杆的近似解,同时求得了马赫反射三波点迹线及壁面超压。另外还导出了新的面积函数表达式。在壁面超压的计算上,这些理论结果与用Whitham方程求得的数值解相符合。     

用CSCM格式模拟二维收缩管道内的流动

 采用Lombard等人提出的CSCM(Conservative Supra-Characteristics Method)差分方法数值模拟二维收缩管道内波系干扰及激波边界层干扰等问题。所得结果同已往计算结果进行了比较。通过压力场和速度场的计算结果可看到Ma数对流场的影响。在Ma数较小时激波角较大,可看到激波经两次相交一次反射逐渐变弱的情况。Re数改变对流场也有明显影响,且给出由于激波干扰引起边界层分离的计算结果。最后还给出收缩管道内激波由正激波逐步形成斜激波的非定常过程。     

有限空间中激波运动的数值模拟

采用一种有限体积的TVD格式模拟喷流产生的启动激波在有限空间内的运动过程。计算获得的激波运动过程与激波在大气中传播时的实验结果一致。在传播过程中出现涡环和多次激波。计算表明,在真空箱为正方体的情况下,激波管出口附近的压力脉冲特性主要取决于有限空间的轴向长度。所得结果对毫秒级激波管喷流模拟装置的真空箱设计有参考价值。     

锥—柱体跨声速绕流的物面侧体激波系

本文在实验结果基础上讨论了锥-柱体跨声速飞行时物面上产生的侧体激波系。分析了波系分布随来流M_∞数的变化;讨论了波系形成与流场内部流动机制的关系,进而解释了波系形成原因;建立了波系的一种模拟方法。最后,阐述了侧体波系研究在实际应用中的重要性。     

在激波边界层干扰中后掠柱形区的独立性原理研究

在后掠压缩角引起的三维激波边界层干扰中,柱形的上游影响区表现出准二维的流动特性。实验结果表明,在柱形区中通常的边界层独立性原理不再适用。但是,如果考虑由后掠引起的横流效应修正,柱形区中的上游影响距离的尺度和相应的二维干扰上游影响区尺度之间可以得到很好的相关关系。     

凹面腔内二维激波会聚特性研究

为了研究凹面腔内二维激波会聚的特性,设计、建立了二维激波会聚冷态实验系统,进行了二维激波会聚冷态实验,结果表明:当射流入口宽度增大时,相同射流入射压力下,喷入凹面腔内的射流强度增大,导致凹面腔内的激波会聚强度增强,从而压力脉动幅值增大,频率也随之增大.数值模拟了二维冷态激波会聚过程,从模拟结果中的压力和频率上看,模拟结果与实验结果基本上是吻合的.     

振荡管内气柱谐振的研究

首次提出振荡管气柱谐振概念。对气柱谐振机理、条件、有关因素对谐振激励频率的影响及谐振状态下管内振荡流特性进行了探讨。结果表明：气柱谐振与管内激波运动有关;谐振状态下管内出现最强的压力波并产生最强的冷效应和热效应。