一种改进混合迎风格式在翼型流场激波捕捉中的应用

波阻是飞行器超音速飞行的关键设计因素,精确捕捉激波在流场中的位置,是数值模拟含激波流场和精确计算波阻的一个重要研究内容.本文基于网格节点有限体积空间离散方法,采用AUSM格式与FVS格式的混合格式(MAUSM方法)计算对流通量,从而抑制在数值模拟流场出现的激波处振荡和过冲现象,确保AUSM准确捕获接触间断的特性和FVS格式捕捉激波的能力.本文使用MAUSM方法分别计算了在跨声速和超声速条件下的NACA0012翼型流场,并与中心差分格式的计算结果进行比较.结果表明,对于存在激波的翼型流场,MAUSM方法是有效的.     

可压缩流向涡与激波轴对称干扰的数值模拟1)

用NS方程数值模拟了可压缩流向涡和激波轴对称相互作用现象.数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较.结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征.提出了判断流向涡与运动激波相互作用中旋涡破碎的准则.     

抽吸激波管的性能分析与计算

为了减小由于反射激波和透射激波分叉引起的反射型激波风洞试验气体提前受到污染的现象,本文研究了一种新的具有抽吸的激波管流动,分析了抽吸缝的作用,给出了这种抽吸激波管性能参数的计算方法,同时还给出了反射激波与边界层相互作用引起的激波分叉的形状随抽吸量变化的计算公式。实验证实了边界层抽吸可以有效地减小激波与壁面边界层相互作用所产生的分离现象。计算与测量结果是一致的。     

两介质流界面-激波相互作用RKDG 方法应用分析

为精确模拟多介质流界面运动现象,采用RKDG方法结合虚拟流体方法对气-气、气-液和液-气等多种界面-激波相互作用问题展开研究。数值结果表明,RKDG方法的时空高精度特征使其能够精确、稳健地求解各种复杂界面运动问题。最后,对水下激波自由面折射问题用多种DG格式限制器进行了计算,对比了它们的间断捕捉能力。     

可压缩流向涡与激波轴对称干扰的数值模拟

用ＮＳ方程数值模拟了可压缩流向涡和激波轴对称相互作用现象．数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较．结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征．提出了判断流向涡与运动激波相互作用中旋涡破碎的准则     

可压缩流向涡与激波轴对称干扰的数值模拟

用ＮＳ方程数值模拟了可压缩流向涡和激波轴对称相互作用现象。数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较,结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征。提出了判断流向涡与运动激波相互作用中旋涡破碎的准则。     

激波和湍流相互作用的数值模拟

本文综述了关于激波和湍流相互作用数值模拟的近期研究进展, 主要包括激波和均匀各向同性湍流、激波和湍流边界层、激波和射流以及激波和尾迹的相互作用. 激波和湍流相互作用特性受到诸多因素的影响,如激波的强度、位置、形状和流动边界以及来流的湍流状态和可压缩性等. 激波和湍流的相互作用会引起流场结构、激波特性和湍流统计特性的显著变化. 最后简要讨论了激波和湍流相互作用数值研究需要关注的一些问题.      

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成; (2) U型空泡团头部溃灭; (3) U型空泡团腿部溃灭.在U 型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象, 抑制了下一周期的空化发展.      

真实虚拟流方法在多介质可压缩流动模拟中的应用

为了克服原始虚拟流方法(ghost fluid method,GFM)在处理激波与大密度比流体-流体（气-水）界面相互作用时遇到的困难,采用真实虚拟流法(real ghost fluid method,RGFM)处理流体界面附近的虚拟点,结合HLLC(Harten-Lax-Van Leer with contact discontinuities)格式求解Euler方程,采用五阶WENO(weighted essentially nonoscillatory)格式求解level set输运方程。通过一维和二维算例的物质界面捕捉研究,证明RGFM在处理小密度比界面问题时优于GFM,同时RGFM还可用于求解激波与大密度比物质界面相互作用问题。计算表明,将RGFM引入到本文算法中,可精确捕捉到激波与界面（气-气、气-水界面）相互作用的变化细节,包括大密度比界面的剧烈变形和破碎,并具有较高的计算分辨率。     

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成;(2) U型空泡团头部溃灭;(3) U型空泡团腿部溃灭.在U型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象,抑制了下一周期的空化发展.     

多管脉冲爆轰发动机管外复杂波系的研究

为了研究多管脉冲爆轰发动机外流场复杂波系变化特点,用CE/SE方法对7管脉冲爆 轰发动机内、外流场进行了计算. 7个爆轰管的排列方式是1个居中,另外6个绕其形成环 状排列,管内填充按化学当量比混合的氢气氧气混和物. 计算结果表明：CE/SE方法是一种 好的爆轰波模拟方法,能有效地捕捉激波等强间断;7管脉冲爆轰发动机外流场中存在着多 道相互作用相互影响的激波与膨胀波. 计算结果对多管脉冲爆轰发动机外流场复杂波系的特 性研究具有重要参考价值.     

可压缩流向涡与反向运动激波相互作用的实验

对可压缩流向涡与反向运动激波相互作用的现象进行了实验研究．实验在９４ｍｍ×９４ｍｍ的方截面激波管中进行．在实验段上游安装了一个有限翼展平直机翼．当入射激波通过机翼后，波后２区气流在模型翼尖诱导出一条流向涡．入射激波在激波管端壁反射后，形成的反射激波在观察窗处和流向涡发生作用．实验中拍摄了激波与流向涡作用全过程的纹影照片，观察到了一些和定常激波与旋涡相互作用不同的现象，并与数值计算结果进行了初步比较     

再冲击载荷作用下流动混合的数值模拟

基于多介质的流体体积分数VOF（volume of fluid）方法和PPM（piecewise parabolic method）方法,给出和发展了可用于多介质粘性流体动力学的数值计算方法和计算代码MVPPM（multi-viscosity-fluid piecewise parabolic method）。为了检验和验证此计算代码,对某激波管实验再冲击载荷作用下的流体动力学不稳定性及其导致的混合过程进行了数值模拟,计算结果与实验结果一致。同时还研究了激波反射冲击作用下流体混合区的演化情况,在反射激波和混合区相互作用的瞬间,混合区的宽度明显减小,之后又迅速增大;另外,混合增长率与初始扰动的频谱有很大关联。通过对有粘性（分子动力学粘性）和无粘性结果的对比,发现粘性对混合区的影响很小。     

脉冲激光等离子体与正激波相互作用的PIV实验研究

脉冲激光等离子体与超声速流场相互作用在飞行器减阻隔热、点火助燃等方面具有重要的应用价值.纹影实验方法只能定性或半定量地反映流动状态.为定量研究速度分布和旋涡结构,针对激光等离子体及其与正激波相互作用过程开展粒子图像测速PIV实验研究.在激波管实验平台上建立了纳秒脉冲激光能量沉积系统和PIV测量系统,通过定量测量,探明了激光等离子体引致的激光空气泡以及热核的流动特性,揭示了激光等离子体在正激波冲击下的流动特性与演化规律,并给出了激光能量大小和位置对相互作用过程的影响.结果表明:激光空气泡内的速度分布在激光入射方向上并不关于击穿点对称,而是在靠近激光入射方向一侧的流速略大于远离激光入射方向一侧;斜压导致热核在演化初期产生涡环,后期则由剪切主导;正激波与激光空气泡界面、热核界面相互作用时,产生斜压涡量,当激光能量为87.8 mJ、正激波马赫数1.4时,热核在正激波作用下产生的涡量比在静止空气中演化时大1个数量级;激光与正激波相互作用的关键过程是热核在正激波冲击下演化成涡环,在激波波前注入激光能量能够获得更加显著的涡环.     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.      

基于level set的Eulerian-Lagrangian耦合方法及其应用

提出了一种基于水平集的Eulerian-Lagrangian耦合方法,其中Lagrangian方法采用相容显式有限元拉氏方法,Eulerian方法采用基于近似Riemann解的有限体积Eulerian方法,多介质界面处理采用新的水平集和Ghost方法计算. 给出了若干数值算例,包括激波管问题以及金属和气体的运动界面及其大变形问题,并分别与精确解和相容显式有限元拉氏方法的计算结果进行了对比. 数值结果表明,该方法计算结果正确,精度较高,能够准确捕捉物质界面,适用于处理大变形问题.      

激波/边界层相互作用诱导的激波风洞气体污染问题

应用多组分NS方程和频散可控耗散格式(DCD)计算了爆轰驱动激波风洞中反射激波/ 边界层/接触面的相互作用过程,分析了驱动气体与试验气体在壁面射流作用下的掺混机制及 其对风洞试验时间的影响. 为了延长风洞的试验时间,提出在风洞贮室内增设环形隔板, 以隔离壁面射流,延长风洞试验时间的方法. 计算结果表明：环形隔板确实可以限制驱动气 体与试验气体的过早掺混,显著增加激波风洞的有效试验时间.     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.     

磁流体动力学斜激波控制数值模拟分析

高超声速飞行器MHD(磁流体动力学)斜激波控制应用的关键在于理解等离子体斜激波流场与磁场的相互作用规律,这里发展了全MHD数值模拟方法对其进行研究,数值方法基于八波方程附加源项形式,进行有限体积离散,采用了Roe求解器、OC-TVD空间格式和LU-SGS方法,且采用投影方法降低磁场伪散度误差.考察外加均匀磁场的马赫10无粘导电拐角流动,压缩角为10°.结果中散度误差较低,并且通过激波参数验证了结果的准确性.流场显示,磁场使得激波角增大,部分情况下出现了快、慢激波结构,其中快激波变化更明显;壁面压强根据磁场的不同出现了不同程度的降低.最后采用群速度图方法进行了快慢激波形式分析,解释了磁场影响下流场形式变化机理.     

平面激波绕45°山峰后的马赫效应

本文通过激波管二维光测模型实验和采用三角形网格的FLIC方法进行数值计算两种手段,针对一平面激波S绕45°山峰后与地面相互作用而产生的、从规则反射RR到马赫反射MR转换性质等问题进行了分析和研究,并给出了RR→MR转换临界点L_k随入射激波马赫数Ma_i以及激波入射角ω_i、反射角ω_r,随距离L~*变化的规律。文中考虑了不同Ma_i的影响,并将二维光测模型实验和数值计算的结果同现有的理论结果做了比较。