利用通量限制的高分辨Godunov型格式的熵条件

1．引言考虑非线性双曲型守恒律方程的Cauchy问题式中f（w）∈C2（R）f＂,（w）≥0,初值。u0∈BV（R）．此问题通常只存在弱解,且需附加熵条件以保证解的唯一性．方程（1.1）的数值方法研究发展很快,但一阶精度格式（如Godunov格式）分辨率很低,而二阶精度格式在间断附近存在振荡;TVD格式则是一种成功的高分辨率无振荡格式．此外,双曲型守恒律数值方法的收敛性取决于差分格式的总变差稳定和离散熵条件．文献[2]中给出了利用通量限制构造TVD格式的方法,[1]则讨论了SOR-TVD格式的熵条件．本文第2节回顾了问的方法,具体导出了…     

一类满足熵增条件的流体力学方程守恒型格式

Lax,Wandtoff曾经证明:对于与守恒律方程组相容的守恒型差分格式,如果其差分解几乎处处有界收敛,那么极限函数是原方程组的一个弱解,并且提出了二阶精度的L-W格式.但是,一些数值计算表明,用二阶守恒型格式(如L-W格式及Mac Corma-ck格式),可能得到非物理解的计算结果.通常称满足熵条件的弱解为物理解.对     

一维单个守恒型方程的二阶熵耗散格式

本文考虑一维单个守恒律方程，对其设计了一种非线性守恒型差分格式，此格式为二阶Godunov型的，用的是分片线性重构，重构函数的斜率是根据熵耗散得到的，格式满足熵条件，且数值实验表明格式具有非线性稳定性，在此格式中一个所谓的熵耗散函数起了很重要的作用，它在每个网格的计算中耗散熵，在文中我们给出了熵耗散函数应满足的条件，并给出了一种具体的构造形式，最后给出了一些数值算例，从中可看出熵耗散函数是如何抑制非物理振荡的，及格式对计算的有效性。     

求解二维Euler方程的旋转通量混合格式

为提高求解二维Euler方程数值结果的分辨率,提出了一种旋转通量混合格式.该算法采用旋转通量法的类一维处理思想,通量函数选用满足热力学第二定律的熵稳定数值通量和具有良好鲁棒性的HLL数值通量耦合的混合格式,时间方向采用三阶强稳定Runge-Kutta方法进行推进.该旋转通量混合格式具有结构简单、分辨率高的优点,数值结果表明了该算法的良好特性.     

求解二维浅水波方程的移动网格旋转通量法

为提高求解二维浅水波方程数值算法的分辨率,拟构造求解该方程的新算法:基于移动网格法,选用熵稳定数值通量函数,利用旋转不变性得到混合数值通量.该算法中,浅水波方程的数值求解和依据解的特性进行自适应疏密分布的网格计算过程交错进行.利用变分原理进行网格重构,新网格上的物理量采用二阶精度的守恒型插值公式计算,最终采用三阶强稳定Runge-Kutta法与满足热力学第二定律的熵稳定格式实现浅水波方程的数值求解.数值结果表明,新算法具有良好的间断捕捉能力,分辨率高.     

用WENO方法求解双曲型守恒律方程组的初(边)值问题

本文用WENO算法解决双曲型守恒律方程组初（边值）问题．给出一种满足熵条件、Sδ熵条件和边界熵条件的WENO算法．通过这个算法就能得到守恒律方程组的数值解，数值解和理论解是非常吻合的．     

求解二维浅水波方程的旋转混合格式北大核心CSCD

针对二维浅水波方程数值求解问题,构造了一种旋转通量混合格式.空间方向上,该算法利用浅水波方程通量函数的旋转不变性,在单元界面法线方向及单元界面切线方向上采用可消除红斑现象的HLL与满足热力学第二定律的熵稳定加权混合数值通量函数,时间方向上采用三阶强稳定Runge-Kutta法.数值结果表明,该混合格式对于二维浅水波方程数值求解具有分辨率高的良好特性.     

求解Euler方程的高精度熵相容格式

熵相容格式相比于一般的熵稳定格式进一步控制了激波处的熵增量,一维情况下能有效消除膨胀激波及间断处的伪振荡等现象.对于Euler方程,可以通过对特征变量进行WENO重构以获得高阶熵相容格式的数值粘性项,然后与高阶熵守恒格式结合得到高精度熵相容格式,在WENO重构过程中的权重关于特征变量是非线性的,这导致了大量的向量内积运算.通过用压强和熵代替特征变量来计算权重,可以显著减少重构的计算量,并且数值算例表明这种权重的计算方式能很好地保持数值格式的高阶精度和基本无振荡的效果.     

线性传输方程带二次多项式重构的熵格式（英文）

最近几年来,茅德康等发展了一类有限体积格式计算偏微分方程[1,3-4,6-9].该类格式得到比较好的计算结果.在文[8]中王和茅提出一个满足两个守恒律和三个守恒律的熵格式计算线性发展方程,但是该格式是基于线性多项式重构.本文发展了一个基于二次多项式重构满足两个守恒律的熵格式.数值试验表明本文的格式在长时间计算方面优于文[8].     

一维单个守恒律的一类二阶、几乎熵耗散的格式

本文考虑一维单个守恒律方程，对其设计了一个基于熵耗散的非线性守恒型差分格式．本格式的数值流函数是Lax-Freidrichs格式和Lax-Wendroff格式数值流函数的凸组合，凸组合中的系数是由考虑耗散熵来决定的．这样在解的光滑区域内，格式几乎、甚至完全是Lax-Wendroff格式，而在解的间断处，格式几乎、甚至完全是Lax—Freidrichs格式．从而消除了间断附近的非物理振荡，实现了计算的非线性稳定性．理论分析表明本格式在解的非极值点处是二阶精度的，而在解的极值点处至少有一阶精度．数值试验表明格式是有效的．     

On entropy stable schemes for degenerate parabolic multispecies kinematic flow models

Entropy stable schemes for the numerical solution of initial value problems of nonlinear, possibly strongly degenerate systems of convection–diffusion equations were recently proposed in Jerez and Parés's study. These schemes extend the theoretical framework of Tadmor's study to convection–diffusion systems. They arise from entropy conservative schemes by adding a small amount of viscosity to avoid spurious oscillations. The main condition for feasibility of entropy conservative or stable schemes for a given model is that the corresponding first‐order system of conservation laws possesses a convex entropy function and corresponding entropy flux, and that the diffusion matrix multiplied by the inverse of the Hessian of the entropy is positive semidefinite. As a new contribution, it is demonstrated in the present work, first, that these schemes can naturally be extended to initial‐boundary value problems with zero‐flux boundary conditions in one space dimension, including an explicit bound on the growth of the total entropy. Second, it is shown that these assumptions are satisfied by certain diffusively corrected multiclass kinematic flow models of arbitrary size that describe traffic flow or the settling of dispersions and emulsions, where the latter application gives rise to zero‐flux boundary conditions. Numerical examples illustrate the behavior and accuracy of entropy stable schemes for these applications.     

Implicit kinetic schemes for scalar conservation laws

Based on kinetic formulation for scalar conservation laws, we present implicit kinetic schemes. For time stepping these schemes require resolution of linear systems of algebraic equations. The scheme is conservative at steady states. We prove that if time marching procedure converges to some steady state solution, then the implicit kinetic scheme converges to some entropy steady state solution. We give sufficient condition of the convergence of time marching procedure. For scalar conservation laws with a stiff source term we construct a stiff numerical scheme with discontinuous artificial viscosity coefficients that ensure the scheme to be equilibrium conserving. We couple the developed implicit approach with the stiff space discretization, thus providing improved stability and equilibrium conservation property in the resulting scheme. Numerical results demonstrate high computational capabilities (stability for large CFL numbers, fast convergence, accuracy) of the developed implicit approach. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 26–43, 2002     

Second order scheme for scalar conservation laws with discontinuous flux

Burger, Karlsen, Torres and Towers in [9] proposed a flux TVD (FTVD) second order scheme with Engquist–Osher flux, by using a new nonlocal limiter algorithm for scalar conservation laws with discontinuous flux modeling clarifier thickener units. In this work we show that their idea can be used to construct FTVD second order scheme for general fluxes like Godunov, Engquist–Osher, Lax–Friedrich, … satisfying (A, B)-interface entropy condition for a scalar conservation law with discontinuous flux with proper modification at the interface. Also corresponding convergence analysis is shown. We show further from numerical experiments that solutions obtained from these schemes are comparable with the second order schemes obtained from the minimod limiter.     

On Wavewise Entropy Inequalities for High-Resolution Schemes. I: The Semidiscrete Case

We develop a new approach, the method of wavewise entropy inequalities for the numerical analysis of hyperbolic conservation laws. The method is based on a new extremum tracking theory and Volpert's theory of BV solutions. The method yields a sharp convergence criterion which is used to prove the convergence of generalized MUSCL schemes and a class of schemes using flux limiters previously discussed in 1984 by Sweby.

     

Entropy stable shock capturing space–time discontinuous Galerkin schemes for systems of conservation laws

We present a streamline diffusion shock capturing spacetime discontinuous Galerkin (DG) method to approximate nonlinear systems of conservation laws in several space dimensions. The degrees of freedom are in terms of the entropy variables and the numerical flux functions are the entropy stable finite volume fluxes. We show entropy stability of the (formally) arbitrarily high order accurate method for a general system of conservation laws. Furthermore, we prove that the approximate solutions converge to the entropy measure valued solutions for nonlinear systems of conservation laws. Convergence to entropy solutions for scalar conservation laws and for linear symmetrizable systems is also shown. Numerical experiments are presented to illustrate the robustness of the proposed schemes.     

非凸单个守恒律初边值问题的整体弱熵解的构造

本文研究具有两段常数的初始值和常数边界值的非凸单个守恒律的初边值问题.在流函数具有一个拐点的条件下,由相应的初始值问题弱熵解的结构和Bardos-Leroux-Nedelec提出的边界熵条件,给出初边值问题整体弱熵解的一个构造方法,澄清弱熵解在边界附近的结构.与严格凸的单个守恒律初边值问题相比,非凸单个守恒律初边值问题的弱熵解中包括下列新的相互作用类型:一个接触或非接触激波碰到边界,边界弹回一个非接触激波.     

ON THE CONVERGENCE OF IMPLICIT DIFFERENCE SCHEMES FOR HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are are interested in the fOllowing Cauchy problem for scalar conservation lawswhere the initial data uo E BV(R) and the flux function f 6 C'(n).It is well known that this problem may not always have a smooth global solution even ifthe i…     

High order accurate semi-implicit WENO schemes for hyperbolic balance laws

In this paper we propose a family of well-balanced semi-implicit numerical schemes for hyperbolic conservation and balance laws. The basic idea of the proposed schemes lies in the combination of the finite volume WENO discretization with Roe’s solver and the strong stability preserving (SSP) time integration methods, which ensure the stability properties of the considered schemes [S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Rev. 43 (2001) 89-112]. While standard WENO schemes typically use explicit time integration methods, in this paper we are combining WENO spatial discretization with optimal SSP singly diagonally implicit (SDIRK) methods developed in [L. Ferracina, M.N. Spijker, Strong stability of singly diagonally implicit Runge-Kutta methods, Appl. Numer. Math. 58 (2008) 1675-1686]. In this way the implicit WENO numerical schemes are obtained. In order to reduce the computational effort, the implicit part of the numerical scheme is linearized in time by taking into account the complete WENO reconstruction procedure. With the proposed linearization the new semi-implicit finite volume WENO schemes are designed.A detailed numerical investigation of the proposed numerical schemes is presented in the paper. More precisely, schemes are tested on one-dimensional linear scalar equation and on non-linear conservation law systems. Furthermore, well-balanced semi-implicit WENO schemes for balance laws with geometrical source terms are defined. Such schemes are then applied to the open channel flow equations. We prove that the defined numerical schemes maintain steady state solution of still water. The application of the new schemes to different open channel flow examples is shown.