一种基于新型斜率限制器的理想磁流体方程的高分辨率熵相容格式

针对磁流体动力学方程, 通过分析数据重建所需的条件, 构造一种基于MUSCL(Monotone Upstream-Centred Scheme for Conservation Laws)型重建方法的斜率限制器, 获得了一种求解理想磁流体动力学方程的高分辨率熵相容格式。该格式在解的光滑区域具有高精度; 在解的间断区域可以合理地控制耗散, 可有效避免非物理现象的产生。采用熵稳定格式、熵相容格式和新的高分辨率熵相容格式对一维、二维理想磁流体动力学方程进行数值模拟。结果表明: 新格式能准确地捕捉解的结构, 且具有无振荡、高分辨、鲁棒等特性。     

高分辨率熵相容格式及其在激波流动中的应用

为解决熵守恒格式在激波附近出现数值振荡的问题,本文将熵相容格式与MUSCL格式相结合,提出一种既能适合于激波问题、又不依赖于传统人工黏性经验模型的高分辨率熵相容格式,通过对多个激波问题的数值计算,并对比二阶中心格式、熵守恒格式、熵相容格式和高分辨率熵相容格式的计算结果,发现:熵相容格式具有较好的激波捕捉能力,有效解决了熵守恒格式在激波附近的数值振荡问题;MUSCL重构格式进一步提高了熵相容格式的数值模拟能力,既能精确捕捉激波附近的流动细节,又在光滑区保持二阶精度;在对比的四种格式中,本文提出的高分辨率熵相容格式对激波问题的预测性能最佳。该项工作对发展激波湍流相互作用模型、提高跨/超音速叶轮机械流动预测精度具有理论价值和应用潜力。     

针对二维双曲守恒律方程求解方法的研究

对双曲守恒律方程进行数值求解是计算流体力学的重要研究内容。本文从物理概念出发,通过对计算流体力学和双曲守恒律方程研究现状及发展趋势进行引入,详细介绍了满足熵稳定条件的二维双曲守恒律方程的熵守恒、熵稳定、熵相容、高分辨率熵稳定格式,可将其格式应用于具体算例的数值求解中。     

求解双曲守恒律方程的高分辨率熵稳定格式

熵稳定格式从物理概念出发,保证总熵关于时间耗散,在计算过程中无需进行熵修正,有效避免如膨胀激波,负压力等非物理现象,显示出独特的优点.通过插入限制器和在单元交界面处进行高阶重构,得到一类高分辨率的熵稳定格式.算例结果表明,格式具有可靠性,高精度和基本无振荡性等特点.     

求解粘性Hamilton-Jacobi方程的高阶方法

提出了求解具有粘性项的Hamilton-Jacobi方程的二阶、四阶方法.该方法以加权基本无振荡(WENO)格式为基础,通过修正数值通量函数和构造右端粘性项的基于非线性限制器的二阶近似、基于Taylor展开的四阶近似,成功地求解了一维、二维的粘性Hamilton-Jacobi方程.给出的算例说明了本方法具有高分辨率、鲁棒性和无振荡特性.     

单个守恒型方程熵耗散格式中熵耗散函数的构造

对于一维单个守恒律方程,文[8]设计了一种非线性守恒型差分格式.此格式为二阶Godunov型的,用的是分片线性重构(reconstruction),重构函数的斜率是根据熵耗散得到的.格式满足熵条件.与传统的守恒格式不同的是此格式在计算过程中不仅用到了数值解还用到了数值熵.在此格式中一个所谓的熵耗散函数起到了很重要的作用,它在每一个网格的计算中耗散熵,以保证格式满足熵条件.文[8]中设计的熵耗散函数比较复杂,并且不是很完善.故数值地分析了在格式的构造中为何应给熵以一定的耗散,及应耗散多少.并且给出了一个新的以数值解的二阶差分作为基本模块的熵耗散函数.最后给出了相应的数值算例.     

采用Power限制器的PNND和PWNND格式

为高精度捕捉激波等流场结构,引入一种Power限制器,对NND格式和WNND格式进行改进,分别得到二阶PNND(Power NND)格式和三阶PWNND(Power WNND)格式.该Power类型格式通过Power限制器对相邻待选模板上的一阶导数进行限制,改善了NND格式和WNND格式在间断附近的耗散效应.对各种格式的分析表明,在间断附近采用Power限制器的格式比原格式的表现要好,耗散小且捕捉间断精度高,其中PNND格式虽然只有二阶精度,但在所有算例中与三阶WNND格式的计算结果比较接近,在个别算例中甚至优于WNND格式.最后将PWNND格式应用到二维NACA0012翼型的强迫俯仰振动的数值模拟,计算结果与实验值、参考计算值吻合较好.     

反应流体的移动网格动理学格式

在移动网格上构造一种反应流的动理学格式.首先利用BGK模型推导含化学反应的流体力学方程组,并利用其积分形式构造移动网格上离散格式,再利用自适应移动网格方法得到网格速度,最后利用时间精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.一维与二维的数值实验表明这种格式同时具有高精度、高分辨率的特点.     

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工粘性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.     

低耗散非结构网格多维限制器

基于极值原理分析非结构网格多维限制器的构造准则，通过引入约束位置和约束模板的定义，将非结构网格多维限制器构造方法统一化，并以此为基础提出一种新型非结构网格多维限制器.通过典型算例计算与已有限制器进行比较.结果表明：对于包含激波、接触间断等复杂波系干扰流动，新型多维限制器在保证数值稳定的条件下具有更低的数值耗散，对于复杂波系结构具有更高的流动分辨率.     

Osher Flux with Entropy Fix for Two-Dimensional Euler Equations

We compare in this paper the properties of Osher flux with O-variant and P-variant (Osher-O flux and Osher-P flux) in finite volume methods for the two-dimensional Euler equations and propose an entropy fix technique to improve their robustness. We consider both first-order and second-order reconstructions. For inviscid hypersonic flow past a circular cylinder, we observe different problems for different schemes: a first-order Osher-O scheme on quadrangular grids yields a carbuncle shock, while a first-order Osher-P scheme results in a dislocation shock for high Mach number cases. In addition, a second-order Osher scheme can also yield a carbuncle shock or be unstable. To improve the robustness of these schemes we propose an entropy fix technique, and then present numerical results to show the effectiveness of the proposed method. In addition, the influence of grid aspects ratio, relative shock position to the grid and Mach number on shock stability are tested. Viscous heating problem and double Mach reflection problem are simulated to test the influence of the entropy fix on contact resolution and boundary layer resolution.     

Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations

One of the main challenges in computational simulations of gas detonation propagation is that negative density or negative pressure may emerge during the time evolution, which will cause blow-ups. Therefore, schemes with provable positivity-preserving of density and pressure are desired. First order and second order positivity-preserving schemes were well studied, e.g., [6], [10]. For high order discontinuous Galerkin (DG) method, even though the characteristicwise TVB limiter in [1], [2] can kill oscillations, it is not sufficient to maintain the positivity. A simple solution for arbitrarily high order positivity-preserving schemes solving Euler equations was proposed recently in [22]. In this paper, we first discuss an extension of the technique in [22], [23], [24] to design arbitrarily high order positivity-preserving DG schemes for reactive Euler equations. We then present a simpler and more robust implementation of the positivity-preserving limiter than the one in [22]. Numerical tests, including very demanding examples in gaseous detonations, indicate that the third order DG scheme with the new positivity-preserving limiter produces satisfying results even without the TVB limiter.     

Relativistic Hydrodynamics and Essentially Non-oscillatory Shock Capturing Schemes

Numerical solutions of relativistic hydrodynamic equations are obtained with essentially non-oscillatory (ENO) finite differencing schemes. The method is explicit, conservative, consistent with the entropy condition, and high order accurate in space and time. The present implementation is applicable to the most general, three-dimensional problems with an arbitrary equation of state. Numerical experiments, including computations of multi-dimensional flows, demonstrate that the method delivers sharp, non-oscillatory shock transitions without sacrificing high resolution of the smooth regions. This extends results already established for the Euler gas dynamics to the relativistic regime, suggesting the usefulness of ENO schemes for modelling relativistic nuclear collisions.     

利用通量限制思想改进紧致格式

利用通量限制思想改进紧致格式计算有间断流场的性能,并设计出一种限制器,该限制器被运用在一系列3至8阶的紧致格式上.数值实验表明,通量限制型紧致格式不仅具有较高的精度和分辨率,而且还能有效地抑制非物理振荡,适用于各种高低Mach数的流动,捕捉到的流场间断所占网格点数少.     

Entropy considerations in numerical simulations of non-equilibrium rarefied flows

Non-equilibrium rarefied flows are encountered frequently in supersonic flight at high altitudes, vacuum technology and in microscale devices. Prediction of the onset of non-equilibrium is important for accurate numerical simulation of such flows. We formulate and apply the discrete version of Boltzmann’s H-theorem for analysis of non-equilibrium onset and accuracy of numerical modeling of rarefied gas flows. The numerical modeling approach is based on the deterministic solution of kinetic model equations. The numerical solution approach comprises the discrete velocity method in the velocity space and the finite volume method in the physical space with different numerical flux schemes: the first-order, the second-order minmod flux limiter and a third-order WENO schemes. The use of entropy considerations in rarefied flow simulations is illustrated for the normal shock, the Riemann and the two-dimensional shock tube problems. The entropy generation rate based on kinetic theory is shown to be a powerful indicator of the onset of non-equilibrium, accuracy of numerical solution as well as the compatibility of boundary conditions for both steady and unsteady problems.     

On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes

We construct uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for Euler equations of compressible gas dynamics. The same framework also applies to high order accurate finite volume (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO)) schemes. Motivated by Perthame and Shu (1996) [20] and Zhang and Shu (2010) [26], a general framework, for arbitrary order of accuracy, is established to construct a positivity preserving limiter for the finite volume and DG methods with first order Euler forward time discretization solving one-dimensional compressible Euler equations. The limiter can be proven to maintain high order accuracy and is easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the positivity property. Following the idea in Zhang and Shu (2010) [26], we extend this framework to higher dimensions on rectangular meshes in a straightforward way. Numerical tests for the third order DG method are reported to demonstrate the effectiveness of the methods.     

Modified kinetic flux vector splitting schemes for compressible flows

We investigate the traditional kinetic flux vector splitting (KFVS) and BGK schemes for the compressible Euler equations. First, based on a careful study of the behavior of the discrete physical variables across the contact discontinuity, we analyze quantitatively the mechanism of inducing spurious oscillations of the velocity and pressure in the vicinity of the contact discontinuity for the first-order KFVS and BGK schemes. Then, with the help of this analysis, we propose a first-order modified KFVS (MKFVS) scheme which is oscillation-free in the vicinity of the contact discontinuity, provided certain consistent conditions are satisfied. Moreover, by using piecewise linear reconstruction and van Leer’s limiter, the first-order MKFVS scheme is extended to a second-order one, consequently, a nonoscillatory second-order MKFVS scheme is constructed. Finally, by combing the MKFVS schemes with the γ-model, we successfully extend the MKFVS schemes to multi-flows, and propose therefore a first- and second-order MKFVS schemes for multi-fluid computations, which are nonoscillatory across fluid interfaces. A number of numerical examples presented in this paper validate the theoretic analysis and demonstrate the good performance of the MKFVS schemes in simulation of contact discontinuities for both single- and multi-fluids.     

On maximum-principle-satisfying high order schemes for scalar conservation laws

We construct uniformly high order accurate schemes satisfying a strict maximum principle for scalar conservation laws. A general framework (for arbitrary order of accuracy) is established to construct a limiter for finite volume schemes (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO) schemes) or discontinuous Galerkin (DG) method with first order Euler forward time discretization solving one-dimensional scalar conservation laws. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. It is straightforward to extend the method to two and higher dimensions on rectangular meshes. We also show that the same limiter can preserve the maximum principle for DG or finite volume schemes solving two-dimensional incompressible Euler equations in the vorticity stream-function formulation, or any passive convection equation with an incompressible velocity field. Numerical tests for both the WENO finite volume scheme and the DG method are reported.