High‐order asymptotic‐preserving schemes for linear systems: Application to the Goldstein–Taylor equations

In this work, we address the numerical approximation of linear systems with possibly stiff source terms which induce an asymptotic diffusion limit. More precisely, we are interested in the design of high‐order asymptotic‐preserving schemes. Our approach is based on a very simple modification of the numerical flux associated with the usual HLL scheme. This alteration can be understood as a numerical diffusion reduction technique and allows to capture the correct asymptotic behavior in the diffusion limit and to consider uniformly high‐order extensions. We more specifically consider the case of the Goldstein–Taylor model but the overall approach is shown to be easily adapted to more general systems.     

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     

Asymptotic preserving numerical schemes for a non‐classical radiation transport model for atmospheric clouds

We present numerical schemes for the P1‐moment and M1‐moment approximations of a non‐classical transport equation modeling radiative transfer in atmospheric clouds. In contrast to classical radiative transfer, the photon path‐length is introduced as an additional variable and serves as pseudo‐time in this model. Because clouds may have optically thick regions, we introduce a diffusive scaling and show that the diffusion limits of the moment models and the original equations agree. Furthermore, we show that the numerical schemes also preserve the diffusion asymptotics as well as the set of admissible and realizable states, both for the explicit and the implicit discretization of the pseudo‐time variable. A source iteration‐like method is proposed, and we observe that it converges slowly in the optical thick case, but a suitable initialization can help to overcome this problem. We validate our method in 1D and present simulation results in the 2D‐case for real cloud data. Copyright © 2013 John Wiley & Sons, Ltd.     

An asymptotic preserving scheme on staggered grids for the barotropic Euler system in low Mach regimes

We present a new scheme for the simulation of the barotropic Euler equation in low Mach regimes. The method uses two main ingredients. First, the system is treated with a suitable time splitting strategy, directly inspired from the previous study that separates low and fast waves. Second, we adapt a numerical scheme where the discrete densities and velocities are stored on staggered grids, in the spirit of MAC methods, and with numerical fluxes derived from the kinetic approach. We bring out the main properties of the scheme in terms of consistency, stability, and asymptotic behavior, and we present a series of numerical experiments to validate the method.     

Efficient local structure‐preserving schemes for the RLW‐type equation

The multisymplectic schemes have been used in numerical simulations for the RLW‐type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we propose three novel efficient local structure‐preserving schemes for the RLW‐type equation, which preserve the local energy exactly on any time‐space region and can produce richer information of the original problem. The schemes will be mass‐ and energy‐preserving as the equation is imposed on appropriate boundary conditions. Numerical experiments are presented to verify the efficiency and invariant‐preserving property of the schemes. Comparisons with the existing nonconservative schemes are made to show the behavior of the energy affects the behavior of the solution.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1678–1691, 2017     

Asymptotic preserving time‐discretization of optimal control problems for the Goldstein–Taylor model

We consider the development of implicit‐explicit time integration schemes for optimal control problems governed by the Goldstein–Taylor model. In the diffusive scaling, this model is a hyperbolic approximation to the heat equation. We investigate the relation of time integration schemes and the formal Chapman–Enskog‐type limiting procedure. For the class of stiffly accurate implicit–explicit Runge–Kutta methods, the discrete optimality system also provides a stable numerical method for optimal control problems governed by the heat equation. Numerical examples illustrate the expected behavior. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1770–1784, 2014     

Stability of contact discontinuity for the Boltzmann equation

The Boltzmann equation which describes the time evolution of a large number of particles through the binary collision in statistics physics has close relation to the systems of fluid dynamics, that is, Euler equations and Navier-Stokes equations. As for a basic wave pattern to Euler equations, we consider the nonlinear stability of contact discontinuities to the Boltzmann equation. Even though the stability of the other two nonlinear waves, i.e., shocks and rarefaction waves has been extensively studied, there are few stability results on the contact discontinuity because unlike shock waves and rarefaction waves, its derivative has no definite sign, and decays slower than a rarefaction wave. Moreover, it behaves like a linear wave in a nonlinear setting so that its coupling with other nonlinear waves reveals a complicated interaction mechanism. Based on the new definition of contact waves to the Boltzmann equation corresponding to the contact discontinuities for the Euler equations, we succeed in obtaining the time asymptotic stability of this wave pattern with a convergence rate. In our analysis, an intrinsic dissipative mechanism associated with this profile is found and used for closing the energy estimates.     

A cell edge‐based linearity‐preserving scheme for diffusion problems on two‐dimensional unstructured grids

We propose a new finite volume scheme for 2D anisotropic diffusion problems on general unstructured meshes. The main feature lies in the introduction of two auxiliary unknowns on each cell edge, and then the scheme has both cell‐centered primary unknowns and cell edge‐based auxiliary unknowns. The auxiliary unknowns are interpolated by the multipoint flux approximation technique, which reduces the scheme to a completely cell‐centered one. The derivation of the scheme satisfies the linearity‐preserving criterion that requires that a discretization scheme should be exact on linear solutions. The resulting new scheme is then called as a cell edge‐based linearity‐preserving scheme. The optimal convergence rates are numerically obtained on unstructured grids in case that the diffusion tensor is taken to be anisotropic and/or discontinuous. Copyright © 2017 John Wiley & Sons, Ltd.     

Efficient structure‐preserving schemes for good Boussinesq equation

The good Boussinesq equation is endowed with symplectic conservation law and energy conservation law. In this paper, some new highly efficient structure‐preserving methods for the good Boussinesq equation are proposed by improving the standard finite difference method (FDM). The new methods only use and calculate values at the odd (or even) nodes to reduce the computational cost. We call this kind of methods odd‐even method (OEM). Numerical results show that the OEM and the standard FDM have nearly the same numerical errors under the same mesh partition. However, the OEM is much more efficient than the standard FDM, such as the consumed CPU time and occupied memory.     

The backward group preserving scheme for 1D backward in time advection‐dispersion equation

In this article, we propose a backward group preserving scheme (BGPS) on the advection‐dispersion equation (ADE) for tackling of the contamination problems. The BGPS has been successfully applied on the backward heat conduction problems as well as the backward in time Burgers equation, but it has never been applied to solve the ADE. The BGPS is able to recover the spatial distribution of groundwater contaminant concentration in this work. Several numerical examples are worked out, and we show, based on those numerical examples, that the BGPS is applicable to the ADE and the method can also handle the ADE with piecewise constant dispersion coefficients. When a steep gradient is appeared in the solution, several steps of the BGPS can be used to retrieve the desired initial data and its result is better than the marching‐jury backward beam equation (MJBBE) method as far as our examples are concerned. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

A vertex‐centered and positivity‐preserving scheme for anisotropic diffusion equations on general polyhedral meshes

We propose a new nonlinear positivity‐preserving finite volume scheme for anisotropic diffusion problems on general polyhedral meshes with possibly nonplanar faces. The scheme is a vertex‐centered one where the edge‐centered, face‐centered, and cell‐centered unknowns are treated as auxiliary ones that can be computed by simple second‐order and positivity‐preserving interpolation algorithms. Different from most existing positivity‐preserving schemes, the presented scheme is based on a special nonlinear two‐point flux approximation that has a fixed stencil and does not require the convex decomposition of the co‐normal. More interesting is that the flux discretization is actually performed on a fixed tetrahedral subcell of the primary cell, which makes the scheme very easy to be implemented on polyhedral meshes with star‐shaped cells. Moreover, it is suitable for polyhedral meshes with nonplanar faces, and it does not suffer the so‐called numerical heat‐barrier issue. The truncation error is analyzed rigorously, while the Picard method and its Anderson acceleration are used for the solution of the resulting nonlinear system. Numerical experiments are also provided to demonstrate the second‐order accuracy and well positivity of the numerical solution for heterogeneous and anisotropic diffusion problems on severely distorted grids.     

A Relaxation Scheme for Solving the Boltzmann Equation Based on the Chapman-Enskog Expansion

Abstract In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jinand Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weaklyparabolic, has a linearly hyperbolic convection part, and is endowed with a generalized eotropy inequality. Itagrees with the solution of the Boltzmann equation up to the Burnett order via the Chapman-Enskog expansion. We develop a one-dimensional non-oscillatory numerical scheme based on the relaxed Burnett system forthe Boltzmann equation. We compare numerical results for stationary shocks based on this relaxation scheme,and those obtained by the DSMC (Direct Simulation Monte Carlo), by the Navier-Stokes equations and bythe extended thermodynamics with thirteen moments (the Grad equations). Our numerical experiments showthat the relaxed Burnett gives more accurate approximations to the shock profiles of the Boltzmann equationobtained by the DSMC, for a range of Mach numbers for hypersonic flows, th     

A positivity‐preserving nonstandard finite difference scheme for the damped wave equation

A positivity‐preserving nonstandard finite difference scheme is constructed to solve an initial‐boundary value problem involving heat transfer described by the Maxwell‐Cattaneo thermal conduction law, i.e., a modified form of the classical Fourier flux relation. The resulting heat transport equation is the damped wave equation, a PDE of hyperbolic type. In addition, exact analytical solutions are given, special cases are mentioned, and it is noted that the positivity condition is equivalent to the usual linear stability criteria. Finally, solution profiles are plotted and possible extensions to a delayed diffusion equation and nonlinear reaction‐diffusion systems are discussed. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

Acoustic limit for the Boltzmann equation in the whole space

This paper is devoted to the following rescaled Boltzmann equation in the acoustic time scaling in the whole space

(0.1)     

一类非线性偏微分方程的四阶格子Boltzmann模型

通过Chapman-Enskog展开技术和多尺度分析,建立了一种新的D1Q4带修正项的四阶格子Boltzmann模型,一类非线性偏微分方程从连续的Boltzmann方程得到正确恢复．统一了KdV和Burgers等已知方程类型的格子BGK模型,还首次给出了组合KdV-Burgers,广义Burgers—Huxley等方程...     

Conservative numerical method for solving the averaged Boltzmann equation

A method is proposed for averaging the Boltzmann kinetic equation with respect to transverse velocities. A system of two integro-differential equations for two desired functions depending only on the longitudinal velocity is derived. The gas particles are assumed to interact as absolutely hard spheres. The integrals in the equations are double. The reduction in the number of variables in the desired functions and the low multiplicity of the integrals ensure the computational efficiency of the averaged equations. A numerical method of discrete ordinates is developed that effectively solves the gas relaxation problem based on the averaged equations. The method is conservative, and the number of particles, momentum, and energy are conserved automatically at every time step.     

Nonstandard methods for the convective‐dispersive transport equation with nonlinear reactions

A new nonstandard Eulerian‐Lagrangian method is constructed for the one‐dimensional, transient convective‐dispersive transport equation with nonlinear reaction terms. An “exact” difference scheme is applied to the convection‐reaction part of the equation to produce a semi‐discrete approximation with zero local truncation errors with respect to time. The spatial derivatives involved in the remaining dispersion term are then approximated using standard numerical methods. This approach leads to significant, qualitative improvements in the behavior of the numerical solution. It suppresses the numerical instabilities that arise from the incorrect modeling of derivatives and nonlinear reaction terms. Numerical experiments demonstrate the scheme's ability to model convection‐dominated, reactive transport problems. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 617–624, 1999     

Macroscopic regularity for the Boltzmann equation

The regularity of solutions to the Boltzmann equation is a fundamental problem in the kinetic theory. In this paper, the case with angular cut-off is investigated. It is shown that the macroscopic parts of solutions to the Boltzmann equation, i.e., the density, momentum and total energy are continuous functions of(x, t) in the region R3×(0, +∞). More precisely, these macroscopic quantities immediately become continuous in any positive time even though they are initially discontinuous and the discontinuities of solutions propagate only in the microscopic level. It should be noted that such kind of phenomenon can not happen for the compressible Navier-Stokes equations in which the initial discontinuities of the density never vanish in any finite time, see [22]. This hints that the Boltzmann equation has better regularity effect in the macroscopic level than compressible Navier-Stokes equations.