间断探测器在间断Galerkin方法中的应用

在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明：如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。     

一种适用于热流计算的改进WENO格式

数值格式是影响高超声速热流数值计算准确度的主要因素之一.根据高超声速流场的特性合理的控制了数值黏性,对WENO(weighted essentially non-oscillatory)格式进行了改进,使其适用于高超声速热流计算.研究结果表明:改进后的格式在捕捉激波时具有良好的激波辨率,在壁面附近尽量减小了数值黏性,提高了格式的黏性分辨率.数值试验也表明了这种改进算法的有效性.      

双曲守恒律方程的一种熵相容格式

提出了一种求解双曲守恒律方程的熵相容数值通量。在熵守恒通量中添加一个二阶迎风项和一个三阶的差商项来保持熵稳定并且抵消解在跨过激波时所产生的激波强度立方倍的熵增,从而实现熵相容。新的数值通量能精确保持定常的接触间断、消除非物理的膨胀激波及负压力等现象。通过采用近年发展起来的WENO方法在单元交界面处进行高阶重构,得到高阶精度的熵相容格式。数值算例采用空间半离散格式,并结合显式三步三阶Runge-Kutta(RK3)方法进行时间推进。不同的算例结果表明,格式具有稳定性、高分辨率和无振荡性等特点。     

限制器对高超声速体表面热流数值模拟的影响

采用Roe的FDS（Flux Difference Splitting）差分格式结合min mod,van Leer和Osher-C限制器,对高超声速体外部流场进行数值模拟.研究了限制器对高超声速体气动力/热数值模拟的区别,对比分析了3种限制器对三维模型表面热流数值模拟结果的影响.研究表明,限制器对数值模拟精度有重要影响,且对气动热数值模拟的影响大于对气动力的影响;对简单或较复杂的气动外形,采用Roe的FDS格式结合min mod限制器获得的表面热特性比它结合van Leer和Osher-C限制器时要好,可为热防护设计提供可靠结果.      

适用于Mie-Grüneisen状态方程的一种通用激波限制器

为了求解含Mie-Grüneisen状态方程多介质可压缩流动问题,针对界面追踪方法,发展了一种通用的激波限制器。根据平均守恒误差及L1误差随参数的变化情况,确定了参数的范围。由于表达式不含有状态方程中的物性参数,所以它适用于任何一种形式的状态方程。通过具体算例分析了激波限制器的开关闭状态。数值结果表明文中构造的激波限制器可以用于计算含强激波的一般状态方程的多介质可压缩流动问题。     

r型网格自适应在间断Galerkin有限元激波捕捉中的应用

间断Galerkin(DG)有限元方法因计算精度高、适用于非结构网格等特点得到广泛研究和应用,其在数值模拟包含强间断流场时存在残差收敛性和计算鲁棒性差问题,均匀分布的网格加剧这一问题并影响激波分辨率。针对该问题,发展了r型网格自适应方法,实现间断Galerkin有限元数值模拟过程中网格自适应加密。基于网格点归一化的压力值作为r型网格自适应中网格点移动驱动力的重要权值,并将网格自适应后的网格点位移变化量与网格点之间的初始位移之比作为驱动力的另一重要权值,实现网格沿激波方向各向异性自适应加密,并且激波附近网格点的相邻网格点同步向激波方向移动。发展了适合间断Galerkin有限元方法的Venkatakrishnan限制器。并列NACA0012翼型超声速算例及三维并列圆柱相互干扰算例结果表明：基于r型网格自适应的间断Galerkin有限元方法能够清晰锐利捕捉激波,提高模拟精度,具有良好的收敛性和鲁棒性。     

一种基于间断Galerkin格式的新型限制器设计

使用高阶间断Galerkin格式求解守恒律方程组时,激波附近的Gibbs效应容易导致非物理解的产生。为抑制这一现象,必须构造合理的限制器对数值解进行处理。目前间断Galerkin格式中的限制器多源于有限体积法,在非结构网格上只对低阶导数项进行限制,对高阶导数项则很难给出普适判据。文章对间断Galerkin解进行广义Fourier展开,实现不同频域范围内的谱分解;在新的模板坐标系下描述各阶方向导数的变化规律;结合当前单元和相邻单元的信息,分层限制各阶方向导数,实现对非物理解的抑制。通过求解Euler方程,对二维Riemann问题、翼型周围的亚、跨声速流动问题、前台阶问题以及超燃冲压发动机内流场激波反射问题进行数值模拟,检验了新型限制器的可靠性以及向高阶格式推广的可行性。     

应用维数分裂方法推广MUSCL和WENO格式的若干问题

首先对有限差分法和有限体积法的差异进行了讨论,在已有文献观点的基础上补充了二者在边界条件处理及网格需求等方面存在差异的新论据。介绍了逐维推导的MUSCL和WENO格式计算控制体界面通量的过程,认为此类格式计算界面通量的方法直接应用于高斯积分型有限体积法不够严谨,从而得到了应用维数分裂方法构造的MUSCL格式和WENO格式不属于高斯积分型有限体积法的观点,“积分格式”这一定义更能准确反映这类格式的特点。此外,讨论了MUSCL格式和WENO格式在曲线坐标系下不能保证守恒的原因,并简单介绍了消除方法。     

一种求解双曲守恒律方程的新方法

近年来,在计算流体力学中,WENO(WeightedEssentiallyNon-oscillatory)格式得到了很大的发展,这也促进了基于WENO格式基础上的CWENO(CentralWENO)格式的发展。WENO格式的主要思想是通过凸组合的方式进行插值重构,以得到高阶逼近,并且在间断附近防止伪震荡的产生。本文构造了一种新的五阶紧凑CWENO格式,该格式所选模板节点少,并且保留了原有中心格式的优点。最后通过一些数值算例来验证该格式的高精度,本性无震荡性质。     

SMOOTHNESS INDICATOR OF WENO SCHEME FOR RESOLVING SHORT WAVE

Based on the traditional fifth-order weighted essentially non-oscillatory(WENO)scheme,a smoothness indicator is introduced to improve the capability of WENO schemes for resolving short waves.In the construction of the new smoothness indicator,the proportion of the first-order term in the original smoothness indicator is reduced by replacing the square of the first-order term with the product of the first-order and the third-order terms.To preserve the fifth-order of convergence rate,the smoothness indicator is combined with the method of Borges,et al.The numerical results show that the proposed schemes are more suitable for simulating turbulent flows or aeroacoustics problems than the previous fifth-order WENO schemes,thanks to its improved resolution on short waves.