高超声速马赫反射

通过实验和数值计算模拟了激波（马赫数７．６）绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项、粘性项和化学源项分别采用了迎风ＴＶＤ格式、中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

环形激波聚焦流场特性的间断有限元方法(英文)

针对环形激波聚焦时产生的高温、高压特性,采用间断有限元方法对环形激波在同轴圆柱形激波管内的聚焦流场进行了数值模拟。为了验证采用方法的有效性,首先对经典激波管问题进行了求解,计算结果与解析解吻合很好。其次,对入射激波马赫数为2.0到5.0的环形激波聚焦流场进行了计算,获得了环形激波在圆柱形激波管内聚焦过程的气动特性。模拟结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的绕射、反射和聚焦等主要流动特征,而且在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力。最后比较了网格尺度对聚焦点峰值压力的影响。     

由火焰聚心和激波聚焦诱导的爆震波研究

通过对几种不同结构形式的爆震发生器进行数值模拟,研究了激波聚焦和火焰聚心现象、气动特性及其机理.数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动Navier-Stokes方程来模拟流体动力学和化学动力学过程.数值计算表明用较低的点火能量对射流火焰燃烧器中可燃混合物点火,层流火焰在狭窄管壁作用下加速,射流火焰在轴线上汇聚过程有利于激波的加强,强激波加速火焰,在多重激波与火焰反复作用下,激波和火焰面之间出现热点,热点迅速放大并形成压力很高的过驱爆震波,而后衰减为稳定的爆震波,不同的激波聚焦腔爆震波的形成过程不同.通过对数值计算的结果进行分析,得到了起爆距离和稳定爆震距离,为进一步试验提供参考.     

三维管内激波/湍流附面层干扰流场的数值模拟

给出了可压缩性及曲率修正两方程湍流模型，用它对跨音速三维喷管管内激波／湍流附面层干扰流场进行了数值模拟；给出了不同纵向截面上的激波结构、等马赫线图和四周固壁上的摩擦力线谱。将计算所得的摩擦力线谱、激波结构与实验结果进行了比较。分析了横截面上的流动结构，计算与实验吻合较好，进一步证实可压缩性及曲率修正两方程湍流模型能够较好地模拟激波／湍流附面层干扰流场，满足工程设计和分析要求。     

激波绕楔形物体的流场分析

使用数值计算方法和数值干涉技术并结合实验干涉图片,分析了激波过拐角绕楔形物体运动的流场.并采用数值彩色干涉技术,将数值计算结果转换为彩色条纹编码的数值干涉图,并与实验干涉图比较,可以辨识过激波发生跳跃后的干涉条纹的位移数,从而确定实验干涉照片的密度场.数值计算采用显式二阶迎风TVD格式.采用数值干涉技术将计算密度场转化成数值干涉图与实验干涉图进行了比较,检验了数值结果的可靠性.数值计算结果显示了包括激波、接触间断和旋涡在内的详细的流场结构,与实验相比,更为细致地刻画了流场演变过程的细节.通过数值计算结果观察到激波传播到旋涡中心附近并没有发生明显耗散,旋涡上半部的激波得到减速,下半部分的激波得到加速.     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

超声速喷流激波噪声基础问题数值模拟研究进展

超声速喷流问题是一个包含激波、旋涡、湍流和声波的多尺度复杂流动问题,其数值模拟方法及激波噪声产生机制是相关研究的长期热点和难点。本文简要回顾了超声速喷流激波噪声研究进展,重点介绍了针对激波噪声计算方法的非物理噪声消除技术和光滑因子设计准则,针对超声速喷流激波噪声研究设计的模型问题（包括旋涡–旋涡相互作用、激波–旋涡相互作用和激波–剪切层相互作用等）,以及轴对称和三维超声速喷流的研究进展。本文还介绍了作者最近针对超声速喷流开展的三维直接数值模拟、实验验证工作和初步分析结果（包括轴对称模态定位、束缚波演化和摆动模态发展等）。     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

非定常欠膨胀射流的间断有限元模拟

采用间断有限元方法对非定常欠膨胀射流进行了数值模拟,将二维守恒方程的间断有限元方法发展到轴对称Euler方程,对欠膨胀射流问题进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉到包括激波、滑移线、涡环和多级马赫盘等复杂流场结构,与实验照片反映的流动特征吻合较好。此外,两个马赫盘的位置与其他数值结果和实验测量结果吻合很好,表明该方法具有很好的鲁棒性,对非定常欠膨胀燃气射流的数值模拟是有效的。     

瞬时爆轰TNT装药爆炸计算方法

对高精度TNT装药爆炸进行了理论计算,在使用二阶显式差分格式的基础上,详细地描述了使用样条函数进行激波装配、波形处理的全过程,收到了预期的效果     

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

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

高超声速侧压式模型进气道不起动特性分析

对某典型高超声速侧压式进气道三维流场进行了数值分析，就壁面压力分布、波系结构和近壁流谱图与实验结果进行了比较，计算反映了流动的基本特征。分析了在不同来流马赫数下的流动特征，随着来流马赫数的减小，激波角增大，压缩波在通道内的反射次数增加，而强度逐渐减弱，总压恢复系数逐渐增大。当马赫数减小到一定数值时，在等直隔离段入口出现喉道截面，进一步减小来流马赫数，流量阻塞，引起压力升高，波系向进口方向移动，波后出现亚声速流场，进气道不起动。同时还发现当不起动现象发生时，由于波后分离包的存在，在进气道的进口前形成向后倾斜的激波而不是正激波。本文还提出了一种确定不起动马赫数的方法。     

氮气辐射强度的激波管测量与验证

超高速飞行器的热防护设计必须考虑激波层内高温气体发射与吸收的辐射能量,需要有效的辐射加热评估手段。相应飞行条件下的光谱辐射强度地面实验测量是验证数值模型和方法、理解高焓流动的重要手段。基于燃烧驱动激波管,发展辐射强度标定技术,针对富氮气环境,开展高温气体光谱辐射强度的高分辨定量化测试,掌握辐射特征,为数值验证提供基础数据。实验获得了激波速度5.70和6.20km/s条件下的气体光谱辐射强度精细结构,数据表明激波波后的非平衡过程对辐射强度存在很大影响。通过求解耦合化学反应动力学模型的Navier-Stokes方程和辐射特性模型,得到对应实验条件下的流场参数和辐射强度,计算结果和实验数据符合很好,验证了数值模拟方法。     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构与在ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

光学干涉技术定量显示双马赫反射的密度场

本文介绍了对激波管流动采用两次曝光全息干涉、M-Z 干涉和差分干涉三种技术所获得的双马赫反射二维密度场的定量测量结果,并与数值计算结果作了比较,两者十分吻合。     

CO2中速度5~7 km/s自由飞圆球流场参数的实验测量和数值模拟

为研究火星进入条件下的非平衡流动特性，在中国空气动力研究与发展中心超高速空气动力研究所弹道靶上测量了CO₂中针对火星探测器进入速度范围5~7 km/s条件下的自由飞圆球的激波脱体距离。实验数据基于阴影法测量，并将其与数值计算结果进行对比，进一步计算了实验流场温度和组分分布等流场参数。一般认为激波脱体距离随来流速度升高而呈单调减小趋势，但研究结果表明:实验状态下，圆球飞行速度约5.5~7.0 km/s的范围内，圆球激波脱体距离随飞行速度升高而增大；采用Park的双温度非平衡模型和5组分6反应的CO₂化学反应动力模型可基本再现本文自由飞圆球激波脱体距离的实验测量数据；根据计算结果推测，本实验状态下自由飞圆球波后靠近激波一侧区域的流场主要处于热化学非平衡状态；当来流速度在约5.5~7.0 km/s的范围内时，流场组分CO开始发生显著离解，是引起圆球激波脱体距离在该速度范围内随速度升高反而增大的可能原因。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。