二维跨音速流动计算的显式多步有限面积方法

本文按Jameson等的思路,采用新的人工粘性项系数,建立显式多步有限面积方法,给出绕NACA0012翼型和RAE2822翼型跨音速流动的计算结果。     

具有变截面驱动段及自由开口端的激波管流场计算

本文用随机选取法计算了一个具有变截面驱动段及自由开口端的激波管流场。在自由开口端,考虑了全部可能的出流及入流情况。计算结果表明,这种边界条件处理方法正确地预示了波系在自由开口端的相互作用及反射,与试验结果及BRL的一维流计算结果相比,三者吻合程度令人满意。     

平板上钝缘舵在超声速绕流中的三维分离特性研究

对于平板上的直立钝缘舵在M_∞=1.79、2.04和2.50的条件下,研究了该平板干扰区中的三维分离特性。着重讨论了钝缘直径、舵面迎角和来流马赫数对主分离线位置和形状的影响。分析了该类三维分离在性态上更强地依赖于无粘流动特性。此外,还对分离区内的二次分离发展过程和形态进行了初步讨论。     

2 AU以内的“渐进型”太阳高能粒子事件模拟

太阳高能粒子（Solar Energetic Particle,SEP）事件是影响地球空间以及深空辐射环境的主要因素之一。“渐进型”太阳高能粒子事件中的高能粒子主要来自于日冕物质抛射（Coronal Mass Ejection,CME）所驱动的激波扩散加速（Diffusive Shock Acceleration,DSA）过程。CME驱动的激波在行星际的传播过程中,其结构不断演化,进而影响到高能粒子的加速过程。本文利用二维太阳高能粒子加速和传播模型,对发生于2014年4月18日的太阳高能粒子事件实例进行了数值模拟。模型考察了黄道面上2 AU的距离以内包含地球所在位置的4个不同点,分别计算了每个点上高能粒子的通量。数值模拟的结果表明：黄道面内不同位置的观察点,与激波波前的磁力线连接不同,从而导致观察点处高能粒子的通量有着显著的差异。该模型的计算结果可以为深空探测计划开展辐射环境研究提供必要的输入。     

左行运动激波主导的管内流动特性研究

为了分析左行运动激波主导的管内流动特征，本文采用非定常数值仿真方法，对亚声速进口条件下等直管道内左行运动激波传播与演化特性、左行运动激波/边界层干扰特征开展研究。研究结果表明：在出口周期性强压力脉动干扰下管内存在连续的左行运动激波，该左行运动激波传播特征具有相似性，激波强度、传播速度按幂函数规律衰减。气流经过左行运动激波后总压、总温、静压阶跃式升高，随后受膨胀波影响气流总压、总温、静压下降；左行运动激波/边界层干扰诱发形成翼型回流区，该回流区随运动激波强度衰减逐渐减小。理论与数值分析表明存在左行运动激波后速度为零和运动激波两侧总压相等的两个临界状态。波前马赫数低于临界值或左行运动激波强度高于临界值时，左行运动激波后为倒流、波后总压高于波前。     

基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

超声速可调进气道内流双解现象及其节流特性

为研究超声速可调进气道喉道调节过程中的内流结构及节流特性，设计了工作马赫数范围为0～4的超声速可调进气道，在来流马赫数为2.9的风洞中借助高速纹影观测系统和动态压力测量系统开展了试验研究。结果表明当内收缩比（ICR）为1.79时，进气道通流流场存在设计和非设计流态的内流双解现象。其中，在设计流态下进气道内通道中的波系结构正常建立；而在非设计流态下进气道内收缩段中存在局部分离诱导的复杂波系结构，导致其整体压升高于设计流态，并存在宽频、低频的小幅振荡。此外设计和非设计流态下进气道的下游节流性能相当，其临界堵塞度分别为42.4%和41.7%，临界压比分别为15.8和16.0。在两类流态的下游节流过程中扰动均以结尾斜激波串的形式向上游传播，且临界状态下激波串头波刚好位于喉道附近，但两类流态的结尾激波串在空间分布特征和振荡特性上均存在明显区别。