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.     

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.

     

以中心- 迎风数值格式研究具有外力项的p 系统

本文以半离散中心- 迎风数值格式研究具有外力项的p 系统. 中心型数值格式用来处理双曲型守恒律或系统的优势是快速且简单, 因为不需要使用近似Riemann 解, 也不需要做特征分解. 我们的数值模拟验证了理论研究结果: 具有外力项的p 系统的解的收敛及爆破行为, 同时也指出一些尚待理论研究的问题.     

Time discretizations for numerical optimisation of hyperbolic problems

We consider a class of numerical schemes for optimal control problems of hyperbolic conservation laws. We focus on finite-volume schemes using relaxation as a numerical approach to the optimality system. In particular, we study the arising numerical schemes for the adjoint equation and derive necessary conditions on the time integrator. We show that the resulting schemes are in particular asymptotic preserving for both, the adjoint and forward equation. We furthermore prove that higher-order time-integrator yields suitable Runge-Kutta schemes. The discussion includes the numerically interesting zero relaxation case.     

双曲型守恒律方程的熵稳定格式的一些讨论

本文讨论双曲型守恒律方程的熵稳定格式.对于给定的熵对,格式所满足的熵条件中的数值熵通量是不唯一的.Tadmor的充分条件可以唯一地确定标量方程的熵守恒通量,但不能唯一确定方程组的熵守恒通量,却可以给出方程组的空间一阶精度的熵守恒格式.也讨论了在熵守恒通量上添加数值粘性得到的显式熵稳定格式需要满足的条件及常见的时间离散对熵守恒和熵稳定的影响.     

Relaxation WENO schemes for multidimensional hyperbolic systems of conservation laws

We present a class of high‐order weighted essentially nonoscillatory (WENO) reconstructions based on relaxation approximation of hyperbolic systems of conservation laws. The main advantage of combining the WENO schemes with relaxation approximation is the fact that the presented schemes avoid solution of the Riemann problems due to the relaxation approach and high‐resolution is obtained by applying the WENO approach. The emphasis is on a fifth‐order scheme and its performance for solving a wide class of systems of conservation laws. To show the effectiveness of these methods, we present numerical results for different test problems on multidimensional hyperbolic systems of conservation laws. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Some techniques for improving the resolution of finite difference component-wise WENO schemes for polydisperse sedimentation models

Polydisperse sedimentation models can be described by a system of conservation laws for the concentration of each species of solids. Some of these models, as the Masliyah–Locket–Bassoon model, can be proven to be hyperbolic, but its full characteristic structure cannot be computed in closed form. Component-wise finite difference WENO schemes may be used in these cases, but these schemes suffer from an excessive diffusion and may present spurious oscillations near shocks. In this work we propose to use a flux-splitting that prescribes less numerical viscosity for component-wise finite difference WENO schemes. We compare this technique with others to alleviate the diffusion and oscillatory behavior of the solutions obtained with component-wise finite difference WENO methods.     

On discreteness of the Hopf equation

The principle aim of this essay is to illustrate how different phenomena is captured by different discretizations of the Hopf equation and general hyperbolic conservation laws. This includes dispersive schemes, shock capturing schemes as well as schemes for computing multi-valued solutions of the underlying equation. We introduce some model equations which describe the behavior of the discrete equation more accurate than the original equation. These model equations can either be conveniently discretized for producing novel numerical schemes or further analyzed to enrich the theory of nonlinear partial differential equations.     

Fully adaptive multiresolution finite volume schemes for conservation laws

The use of multiresolution decompositions in the context of finite volume schemes for conservation laws was first proposed by A. Harten for the purpose of accelerating the evaluation of numerical fluxes through an adaptive computation. In this approach the solution is still represented at each time step on the finest grid, resulting in an inherent limitation of the potential gain in memory space and computational time. The present paper is concerned with the development and the numerical analysis of fully adaptive multiresolution schemes, in which the solution is represented and computed in a dynamically evolved adaptive grid. A crucial problem is then the accurate computation of the flux without the full knowledge of fine grid cell averages. Several solutions to this problem are proposed, analyzed, and compared in terms of accuracy and complexity.

     

Weakly nonoscillatory schemes for scalar conservation laws

A new class of Godunov-type numerical methods (called here weakly nonoscillatory or WNO) for solving nonlinear scalar conservation laws in one space dimension is introduced. This new class generalizes the classical nonoscillatory schemes. In particular, it contains modified versions of Min-Mod and UNO. Under certain conditions, convergence and error estimates for WNO methods are proved.

     

Convergence of relaxation schemes for conservation laws

We study the stability and the convergence for a class of relaxing numerical schemes for conservation laws. Following the approach recently proposed by S. Jin and Z. Xin we use a semilinear local relaxation approximation, with a stiff lower order term, and we construct some numerical first and second order accurate algorithms, which are uniformly bounded in the L^∞ and BV norms with respect to the relaxation parameter. The relaxation limit is also investigated.     

Conformal structure-preserving schemes for damped-driven stochastic nonlinear Schrödinger equation with multiplicative noise

Since many physical phenomena are often influenced by dispersive medium, energy compensation and random perturbation, exploring the dynamic behaviors of the damped-driven stochastic system has becoming a hot topic in mathematical physics in recent years. In this paper, inspired by the stochastic conformal structure, we investigate the geometric numerical integrators for the damped-driven stochastic nonlinear Schrödinger equation with multiplicative noise. To preserve the conformal structures of the system, by using symplectic Euler method, implicit midpoint method and Fourier pseudospectral method, we propose three kinds of stochastic conformal schemes satisfying corresponding discrete stochastic multiconformal-symplectic conservation laws and discrete global/local charge conservation laws. Numerical experiments illustrate the structure-preserving properties of the proposed schemes, as well as favorable results over traditional nonconformal schemes, which are consistent with our theoretical analysis.     

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     

Asymptotic‐preserving Godunov‐type numerical schemes for hyperbolic systems with stiff and nonstiff relaxation terms

We devise a new class of asymptotic‐preserving Godunov‐type numerical schemes for hyperbolic systems with stiff and nonstiff relaxation source terms governed by a relaxation time ε. As an alternative to classical operator‐splitting techniques, the objectives of these schemes are twofold: first, to give accurate numerical solutions for large, small, and in‐between values of ε and second, to make optional the choice of the numerical scheme in the asymptotic regime ε tends to zero. The latter property may be of particular interest to make easier and more efficient the coupling at a fixed spatial interface of two models involving very different values of ε. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Explicit finite difference schemes with extended stability for advection equations

In this study an explicit central difference approximation of the generalized leap-frog type is applied to the one- and two-dimensional advection equations. The stability of the considered numerical schemes is investigated and the scheme with the largest stable time step is found. For the linear and nonlinear advection equations numerical experiments with different schemes from the considered class are performed in order to evaluate the practical stability of the designed schemes.     

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.     

一类时空二阶精度高分辨率MmB差分格式的构造及数值试验

1．引言考虑如下二维双曲型守恒律初值问题的数值解．H．M．Wu和S．L．Yang在文山中给出了MmB差分格式的定义如下：给定（．1）M差分格式定义．若则称格式（1．2）为MmB差分格式．这里BmB表示局部MaximumandminimumBounds．由定义可知，若差分格式（1．2）可写为形式且。＼P’三0，＞。：r’一1．则格式（1．4）为MmB差分格式．j＝l文山构造了二维双曲型守恒律的二类二阶精度的MmB差分格式，使构造二维高分辨格式有了新的突破，但他们是从标量线性双曲型守恒律出发，然后把结果推广到非线性情形．本文直接从二维非线性双曲型守恒律…     

Centred TVD schemes for hyperbolic conservation laws

New first- and high-order centred methods for conservation lawsare presented. Convenient TVD conditions for constructing centredTVD schemes are then formulated and some useful results areproved. Two families of centred TVD schemes are constructedand extended to nonlinear systems. Some numerical results arealso presented.     

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.     

Error and convergence of two numerical schemes for stochastic differential equations

We examine the convergence and error rate of two stochastic numerical schemes using the method of proof used by G. N. Mil'shtein 1 . © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006.