可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...     

基于映射函数的新型五阶WENO格式

研究高精度和高分辨率的差分格式对于复杂流场的数值模拟有重要意义。为了克服WENO-JS格式和WENO-Z在通量函数的一阶和二阶极值点处降阶的缺陷,基于重构权重系数的思想,设计一族映射函数并应用到五阶WENO格式中。近似色散关系表明,WENO-Pe的色散误差和数值耗散均小于WENO-JS、WENO-Z以及其他基于映射函数的WENO格式。新格式与其他格式数值模拟变形的高斯波问题,Sod激波管、Lax激波管、激波密度干扰问题等一维算例,Riemann问题、Rayleigh-Taylor不稳定性问题、双马赫反射问题等二维算例的结果表明：在精度阶相同的情况下,WENO-Pe格式拥有更良好的捕捉间断能力,分辨率更高,适合应用于复杂流场的数值模拟。     

一种求解激波问题的中心差分-WENO混合方法研究

为提高求解包含激波问题的计算精度和效率,发展了一种低耗散、高效率的中心差分-WENO(Weighted Essentially Non-Oscillatory)混合格式,Navier-Stokes方程的无黏项在光滑流场区域采用六阶中心差分格式离散,而间断附近采用五阶WENO格式求解;基于密度设计了一种新型格式开关实现两种格式在光滑区域和间断之间自动切换,确保数值解在间断附近基本无振荡。针对一维激波熵波作用问题验证了混合格式的低耗散特性;对二维Riemann问题的研究表明发展的基于密度的格式开关更为合理,混合格式的计算效率较WENO格式明显提高;将其应用到激波诱导燃烧问题中,混合格式能很好地捕捉定常流场中激波和化学反应锋面的位置以及非定常流场的振荡频率。     

激波绕流与激波碰撞的数值模拟

本文应用基于Poe平均构造的最小耗散的扩散与反扩散守恒差分格式，数值模拟单激波绕流三我形障碍物，及等强度激波碰撞和不等强度激波碰撞的二维非定常复杂流场。计算结果表明，该格式具有很高的单波精度及较好的激波与接触间断分辨度，可直接用于复杂流场的数值模拟。     

改进虚拟流体方法在气-水界面捕捉中的应用

提出了一种改进的Ghost方法,通过求解双介质的Riemann问题,将全部解参量外推得到虚拟流体状态参量,使单一的气-水间断问题转化为两个单介质的Riemann问题。流场控制方程求解采用二阶NND格式,LevelSet方程用五阶WENO格式离散,实现了对气-水界面稳定而准确的捕捉。激波与气-水界面作用的一维和二维算例验证了程序的有效性。对气液两相均运用可压模型,计算了单个水滴被马赫数为2的空气激波扫过后的二维变形过程。     

一类新型自适应反扩散近似Riemann求解器及其应用

对于包含激波、剪切层等复杂结构的流动问题，为了精确模拟剪切层等精细结构，且保证激波计算的稳定性，必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大，Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解，容易导致不稳定。针对这一问题，本文在Riemann求解器中通过合理设计反扩散矩阵，发展了一类具有自适应反扩散的新型Riemann求解器，并将其应用到高阶加权紧致格式，实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性，结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷，既能准确识别剪切层等精细结构，又能保证激波解的稳定性。     

适用于气动学问题的DRP类有限差分方法

研究了一类低耗散低色散的高阶精度有限差分方法,目的是直接计算非定常欧拉方程用于气动声学问题.采用的数值方法是空间四阶、时间三阶精度的色散关系保持(DRP:dispersion-relation-preserving)类有限差分格式,通过波动方程算例验证了格式模拟波动问题的能力;采用加变幅值的高波数人工耗散项来抑制高频数值振荡,使得该格式可推广用于含激波、初始间断或非线性波动的问题.通过含激波、初始间断条件的算例验证了所用变幅值人工耗散的有效性,最后作为初步应用计算了二维亚、超音速均匀流中简单点声源辐射问题,得到了很好的结果,表明此类方法在气动声学问题计算上深有潜力.     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。