899 research outputs found
A new approximate matrix factorization for implicit time integration in air pollution modeling
Implicit time stepping typically requires solution of one or several linear systems with a matrix I−τJ per time step where J is the Jacobian matrix. If solution of these systems is expensive, replacing I−τJ with its approximate matrix factorization (AMF) (I−τR)(I−τV), R+V=J, often leads to a good compromise between stability and accuracy of the time integration on the one hand and its efficiency on the other hand. For example, in air pollution modeling, AMF has been successfully used in the framework of Rosenbrock schemes. The standard AMF gives an approximation to I−τJ with the error τ2RV, which can be significant in norm. In this paper we propose a new AMF. In assumption that −V is an M-matrix, the error of the new AMF can be shown to have an upper bound τ||R||, while still being asymptotically . This new AMF, called AMF+, is equal in costs to standard AMF and, as both analysis and numerical experiments reveal, provides a better accuracy. We also report on our experience with another, cheaper AMF and with AMF-preconditioned GMRES
Krylov implicit integration factor discontinuous Galerkin methods on sparse grids for high dimensional reaction-diffusion equations
Computational costs of numerically solving multidimensional partial
differential equations (PDEs) increase significantly when the spatial
dimensions of the PDEs are high, due to large number of spatial grid points.
For multidimensional reaction-diffusion equations, stiffness of the system
provides additional challenges for achieving efficient numerical simulations.
In this paper, we propose a class of Krylov implicit integration factor (IIF)
discontinuous Galerkin (DG) methods on sparse grids to solve reaction-diffusion
equations on high spatial dimensions. The key ingredient of spatial DG
discretization is the multiwavelet bases on nested sparse grids, which can
significantly reduce the numbers of degrees of freedom. To deal with the
stiffness of the DG spatial operator in discretizing reaction-diffusion
equations, we apply the efficient IIF time discretization methods, which are a
class of exponential integrators. Krylov subspace approximations are used to
evaluate the large size matrix exponentials resulting from IIF schemes for
solving PDEs on high spatial dimensions. Stability and error analysis for the
semi-discrete scheme are performed. Numerical examples of both scalar equations
and systems in two and three spatial dimensions are provided to demonstrate the
accuracy and efficiency of the methods. The stiffness of the reaction-diffusion
equations is resolved well and large time step size computations are obtained
Recommended from our members
A HYBRID METHOD FOR STIFF REACTION-DIFFUSION EQUATIONS.
The second-order implicit integration factor method (IIF2) is effective at solving stiff reaction-diffusion equations owing to its nice stability condition. IIF has previously been applied primarily to systems in which the reaction contained no explicitly time-dependent terms and the boundary conditions were homogeneous. If applied to a system with explicitly time-dependent reaction terms, we find that IIF2 requires prohibitively small time-steps, that are relative to the square of spatial grid sizes, to attain its theoretical second-order temporal accuracy. Although the second-order implicit exponential time differencing (iETD2) method can accurately handle explicitly time-dependent reactions, it is more computationally expensive than IIF2. In this paper, we develop a hybrid approach that combines the advantages of both methods, applying IIF2 to reaction terms that are not explicitly time-dependent and applying iETD2 to those which are. The second-order hybrid IIF-ETD method (hIFE2) inherits the lower complexity of IIF2 and the ability to remain second-order accurate in time for large time-steps from iETD2. Also, it inherits the unconditional stability from IIF2 and iETD2 methods for dealing with the stiffness in reaction-diffusion systems. Through a transformation, hIFE2 can handle nonhomogeneous boundary conditions accurately and efficiently. In addition, this approach can be naturally combined with the compact and array representations of IIF and ETD for systems in higher spatial dimensions. Various numerical simulations containing linear and nonlinear reactions are presented to demonstrate the superior stability, accuracy, and efficiency of the new hIFE method
Runge-Kutta-Gegenbauer explicit methods for advection-diffusion problems
In this paper, Runge-Kutta-Gegenbauer (RKG) stability polynomials of
arbitrarily high order of accuracy are introduced in closed form. The stability
domain of RKG polynomials extends in the the real direction with the square of
polynomial degree, and in the imaginary direction as an increasing function of
Gegenbauer parameter. Consequently, the polynomials are naturally suited to the
construction of high order stabilized Runge-Kutta (SRK) explicit methods for
systems of PDEs of mixed hyperbolic-parabolic type.
We present SRK methods composed of ordered forward Euler stages, with
complex-valued stepsizes derived from the roots of RKG stability polynomials of
degree . Internal stability is maintained at large stage number through an
ordering algorithm which limits internal amplification factors to .
Test results for mildly stiff nonlinear advection-diffusion-reaction problems
with moderate () mesh P\'eclet numbers are provided at second,
fourth, and sixth orders, with nonlinear reaction terms treated by complex
splitting techniques above second order.Comment: 20 pages, 7 figures, 3 table
An iterative semi-implicit scheme with robust damping
An efficient, iterative semi-implicit (SI) numerical method for the time
integration of stiff wave systems is presented. Physics-based assumptions are
used to derive a convergent iterative formulation of the SI scheme which
enables the monitoring and control of the error introduced by the SI operator.
This iteration essentially turns a semi-implicit method into a fully implicit
method. Accuracy, rather than stability, determines the timestep. The scheme is
second-order accurate and shown to be equivalent to a simple preconditioning
method. We show how the diffusion operators can be handled so as to yield the
property of robust damping, i.e., dissipating the solution at all values of the
parameter \mathcal D\dt, where is a diffusion operator and \dt
the timestep. The overall scheme remains second-order accurate even if the
advection and diffusion operators do not commute. In the limit of no physical
dissipation, and for a linear test wave problem, the method is shown to be
symplectic. The method is tested on the problem of Kinetic Alfv\'en wave
mediated magnetic reconnection. A Fourier (pseudo-spectral) representation is
used. A 2-field gyrofluid model is used and an efficacious k-space SI operator
for this problem is demonstrated. CPU speed-up factors over a CFL-limited
explicit algorithm ranging from to several hundreds are obtained,
while accurately capturing the results of an explicit integration. Possible
extension of these results to a real-space (grid) discretization is discussed.Comment: Submitted to the Journal of Computational Physics. Clarifications and
caveats in response to referees, numerical demonstration of convergence rate,
generalized symplectic proo
A class of high-order Runge-Kutta-Chebyshev stability polynomials
The analytic form of a new class of factorized Runge-Kutta-Chebyshev (FRKC)
stability polynomials of arbitrary order is presented. Roots of FRKC
stability polynomials of degree are used to construct explicit schemes
comprising forward Euler stages with internal stability ensured through a
sequencing algorithm which limits the internal amplification factors to . The associated stability domain scales as along the real axis.
Marginally stable real-valued points on the interior of the stability domain
are removed via a prescribed damping procedure.
By construction, FRKC schemes meet all linear order conditions; for nonlinear
problems at orders above 2, complex splitting or Butcher series composition
methods are required. Linear order conditions of the FRKC stability polynomials
are verified at orders 2, 4, and 6 in numerical experiments. Comparative
studies with existing methods show the second-order unsplit FRKC2 scheme and
higher order (4 and 6) split FRKCs schemes are efficient for large moderately
stiff problems.Comment: 24 pages, 5 figures. Accepted for publication in Journal of
Computational Physics, 22 Jul 2015. Revise
- …