Search

Implicit-explicit Runge-Kutta (IMEX-RK) schemes are popular methods to treat multiscale equations that contain a stiff part and a non-stiff part, where the stiff part is characterized by a small parameter \varepsilon. In this work, we prove rigorously the uniform stability and uniform accuracy of a class of IMEX-RK schemes for a linear hyperbolic system with stiff relaxation. The result we obtain is optimal in the sense that it holds regardless of the value of \varepsilon and the order of accuracy is the same as the design order of the original scheme, i.e., there is no order reduction. [ABSTRACT FROM AUTHOR]

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

In this paper, we develop a family of third order asymptotic-preserving (AP) and asymptotically accurate (AA) diagonally implicit Runge-Kutta (DIRK) time discretization methods for the stiff hyperbolic relaxation systems and kinetic Bhatnagar-Gross-Krook (BGK) model in the semi-Lagrangian (SL) setting. The methods are constructed based on an accuracy analysis of the SL scheme for stiff hyperbolic relaxation systems and kinetic BGK model in the limiting fluid regime when the Knudsen number approaches 0. An extra order condition for the asymptotic third order accuracy in the limiting regime is derived. Linear von Neumann stability analysis of the proposed third order DIRK methods are performed to a simplified two-velocity linear kinetic model. Extensive numerical tests are presented to demonstrate the AA, AP, and stability properties of our proposed schemes. [ABSTRACT FROM AUTHOR]

Many hyperbolic and kinetic equations contain a non-stiff convection/transport part and a stiff relaxation/collision part (characterized by the relaxation or mean free time ε). To solve this type of problems, implicit-explicit (IMEX) multistep methods have been widely used and their performance is understood well in the non-stiff regime (ε = O(1)) and limiting regime (ε → 0). However, in the intermediate regime (say, ε =O(Δ t)), uniform accuracy has been reported numerically without a complete theoretical justification (except some asymptotic or stability analysis). In this work, we prove the uniform accuracy - an optimal a priori error bound - of a class of IMEX multistep methods, IMEX backward differentiation formulas (IMEX-BDF), for linear hyperbolic systems with stiff relaxation. The proof is based on the energy estimate with a new multiplier technique. For nonlinear hyperbolic and kinetic equations, we numerically verify the same property using a series of examples. [ABSTRACT FROM AUTHOR]

We study the stability of the semi-discrete central scheme for the linear damped wave equation with boundary. We exhibit a sufficient condition on the boundary to guarantee the uniform stability of the initial boundary value problem for the relaxation system independently of the stiffness of the source term and of the space step. The boundary is approximated using a summation-by-parts method and the stiff stability is proved using energy estimates and the Laplace transform. We also investigate if the condition is also necessary, following the continuous case studied by Xin and Xu (J. Differ. Equ.167 (2000) 388–437). [ABSTRACT FROM AUTHOR]

We consider the large time asymptotic behavior of the global solutions to the initial value problem for the nonlinear damped wave equation with slowly decaying initial data. When the initial data decay fast enough, it is known that the solution to this problem converges to the self‐similar solution to the Burgers equation called a nonlinear diffusion wave, and its optimal asymptotic rate is obtained. In this paper, we focus on the case that the initial data decay more slowly than previous works and derive the corresponding asymptotic profile. Moreover, we investigate how the change of the decay rate of the initial values affect its asymptotic rate. [ABSTRACT FROM AUTHOR]

This paper is concerned with modeling nonequilibrium phenomena in spatial domains with boundaries. The resultant models consist of hyperbolic systems of first-order partial differential equations with boundary conditions (BCs). Taking a linearized moment closure system as an example, we show that the structural stability condition and the uniform Kreiss condition do not automatically guarantee the compatibility of the models with the corresponding classical models. This motivated the generalized Kreiss condition (GKC)—a strengthened version of the uniform Kreiss condition. Under the GKC and the structural stability condition, we show how to derive the reduced BCs for the equilibrium systems as the classical models. For linearized problems, the validity of the reduced BCs can be rigorously verified. Furthermore, we use a simple example to show how thus far developed theory can be used to construct proper BCs for equations modeling nonequilibrium phenomena in spatial domains with boundaries.

We solve a nonlinear, one-dimensional initial boundary-value problem of thermoelasticity in generalized thermodynamics. A Cattaneo-type evolution equation for the heat flux is used, which differs from the one used extensively in the literature. The hyperbolic nature of the associated linear system is clarified through a study of the characteristic curves. Progressive wave solutions with two finite speeds are noted. A numerical treatment is presented for the nonlinear system using a three-step, quasi-linearization, iterative finite-difference scheme for which the linear system of equations is the initial step in the iteration. The obtained results are discussed in detail. They clearly show the hyperbolic nature of the system, and may be of interest in investigating thermoelastic materials, not only at low temperatures, but also during high temperature processes involving rapid changes in temperature as in laser treatment of surfaces. [ABSTRACT FROM AUTHOR]

This paper is concerned with the large time behavior of solutions to the initial value problem for the damped wave equations with nonlinear convection in one-dimensional whole space. In 2007, Ueda and Kawashima showed that the solution tends to a self similar solution of the Burgers equation. However, they did not mention that their decay estimate is optimal or not. Under this situation, the aim of this paper was to find out the sharp decay estimate by studying the second asymptotic profile of solutions. The explicit representation formula and the decay estimates of the solution for the linearized equation including the lower order term play crucial roles in our analysis. [ABSTRACT FROM AUTHOR]

We propose a conservation-dissipation formalism (CDF) for coarse-grained descriptions of irreversible processes. This formalism is based on a stability criterion for non-equilibrium thermodynamics. The criterion ensures that non-equilibrium states tend to equilibrium in long time. As a systematic methodology, CDF provides a feasible procedure in choosing non-equilibrium state variables and determining their evolution equations. The equations derived in CDF have a unified elegant form. They are globally hyperbolic, allow a convenient definition of weak solutions, and are amenable to existing numerics. More importantly, CDF is a genuinely nonlinear formalism and works for systems far away from equilibrium. With this formalism, we formulate novel thermodynamics theories for heat conduction in rigid bodies and non-isothermal compressible Maxwell fluid flows as two typical examples. In these examples, the non-equilibrium variables are exactly the conjugate variables of the heat fluxes or stress tensors. The new theory generalizes Cattaneo's law or Maxwell's law in a regularized and nonlinear fashion. [ABSTRACT FROM AUTHOR]

Abstract: We investigate stability of multidimensional planar shock profiles of a general hyperbolic relaxation system whose equilibrium model is a system, under the necessary assumption of spectral stability and a standard set of structural conditions that are known to hold for many physical systems. Our main result, generalizing the work of Kwon and Zumbrun in the scalar relaxation case, is to establish the bounds on the Greenʼs function for the linearized equation and obtain nonlinear asymptotic behavior/sharp decay rate of perturbed weak shock profiles. To establish Greenʼs function bounds, we use the semigroup approach in the low-frequency regime, and use the energy method for the high-frequency bounds, separately. For the system equilibrium case, the analysis of the linearized equation is complicated due to glancing phenomena. We treat this difficulty similarly as in the inviscid and viscous systems, under the constant multiplicity condition. [Copyright &y& Elsevier]

We are concerned with the initial value problem for a damped wave equation with a nonlinear convection term which is derived from a semilinear hyperbolic system with relaxation. We show the global existence and asymptotic decay of solutions in W1,p (1 ≤ p ≤ ∞) under smallness condition on the initial data. Moreover, we show that the solution approaches in W1,p (1 ≤ p ≤ ∞) the nonlinear diffusion wave expressed in terms of the self-similar solution of the Burgers equation as time tends to infinity. Our results are based on the detailed pointwise estimates for the fundamental solutions to the linearlized equation. [ABSTRACT FROM AUTHOR]

Abstract: We study the large time behavior of solutions of a one-dimensional hyperbolic relaxation system that may be written as a nonlinear damped wave equation. First, we prove the global existence of a unique solution and their decay properties for sufficiently small initial data. We also show that for some large initial data, solutions blow-up in finite time. For quadratic nonlinearities, we prove that the large time behavior of solutions is given by the fundamental solution of the viscous Burgers equation. In some other cases, the convection term is too weak and the large time behavior is given by the linear heat kernel. [Copyright &y& Elsevier]

This paper is concerned with modeling nonequilibrium phenomena in spatial domains with boundaries. The resultant models consist of hyperbolic systems of first-order partial differential equations with boundary conditions (BCs). Taking a linearized moment closure system as an example, we show that the structural stability condition and the uniform Kreiss condition do not automatically guarantee the compatibility of the models with the corresponding classical models. This motivated the generalized Kreiss condition (GKC)—a strengthened version of the uniform Kreiss condition. Under the GKC and the structural stability condition, we show how to derive the reduced BCs for the equilibrium systems as the classical models. For linearized problems, the validity of the reduced BCs can be rigorously verified. Furthermore, we use a simple example to show how thus far developed theory can be used to construct proper BCs for equations modeling nonequilibrium phenomena in spatial domains with boundaries. [ABSTRACT FROM AUTHOR]