非结构网格通量重构算法下三种紧致WENO限制器对比研究

在通量重构算法框架下对比研究了3种最新提出的用于间断伽辽金算法的紧致型WENO(weighted essentially non-oscillatory)限制器,即简单WENO限制器、p阶加权WENO限制器和多精度WENO限制器。这3种WENO限制器能够较高精度地模拟流场并捕捉激波,其紧致性在于问题网格单元中的重构只涉及问题网格单元及其相邻网格单元。同时,使用保精度的保正限制器用以避免可能出现的非物理的密度和压力负值。最后,采用非结构四边形网格在双马赫反射、激波与涡相互作用、激波反射、黏性弓形激波和激波与层流边界层相互作用等多个二维算例中,对这3种WENO限制器进行分析比较。结果显示:多精度WENO限制器与p阶加权WENO限制器能够高精度模拟流场并捕捉激波,同时一定程度抑制通量重构方法本身在激波附近的数值伪振荡;p阶加权WENO限制器与多精度WENO限制器相比,其稳态收敛性相对更好;简单WENO限制器则性能较差。根据以上研究,提出亟须发展能够使高阶WENO限制器在稳态问题中收敛的间断探测器。

基于简单WENO-间断Galerkin的Euler方程自适应计算

为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于"问题单元"的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有"悬挂节点"的非结构网格上的激波捕捉。

基于简单WENO-RKDG的溃坝问题数值模拟

对溃坝问题水流间断面的高精度、高分辨率数值模拟是水动力学的重要内容。简单加权本质无振荡(WENO)限制器由"问题单元"及其相邻单元的解重构"问题单元"的解,从而抑制数值解的非物理振荡,能够很好地模拟间断问题。本文详细介绍了简单WENO限制器的基本原理和过程。将简单WENO限制器-Runge-Kutta间断Galerkin方法应用于二维浅水控制方程的求解中,对二维矩形明渠中大坝瞬间全溃、局部溃塌所致的水流运动进行了数值模拟,并将数值计算结果与理论分析进行了比较。计算结果表明,方法能够清晰地捕捉到溃坝全过程中的间断,没有非物理的振荡现象发生,简单WENO限制器-RKDG方法能够很好地模拟溃坝波的演进过程。

一种求解黏弹塑性介质流动的全隐格式算法

在岩石圈动力学数值模拟中,现有的黏弹塑性数值模型通常在每个时间步先使用迎风间断Galerkin方法对偏应力张量进行旋转,然后使用Particle-In-Cell(PIC)方法或场方法求解对流方程,所构成的时间离散格式为显格式或半隐格式.我们将黏弹塑性介质的经典数值模型和非牛顿流体力学领域的黏弹性流体问题计算方法相结合,提出了一种基于有限单元法的求解黏弹塑性介质流动的全隐格式算法.本文通过数值实验将这种全隐格式算法与PIC方法和半隐格式算法进行了详细的对比,实验结果表明全隐格式算法的数值稳定性优于PIC方法,而当Deborah数较高时精度优于半隐格式算法.同时,我们在应力场引入三阶WENO(Weighted Essentially Non-Oscillatory)限制器,可以在保留数值解精度的同时有效消除应力集中引起的数值振荡.