807 research outputs found
Preconditioned fully implicit PDE solvers for monument conservation
Mathematical models for the description, in a quantitative way, of the
damages induced on the monuments by the action of specific pollutants are often
systems of nonlinear, possibly degenerate, parabolic equations. Although some
the asymptotic properties of the solutions are known, for a short window of
time, one needs a numerical approximation scheme in order to have a
quantitative forecast at any time of interest. In this paper a fully implicit
numerical method is proposed, analyzed and numerically tested for parabolic
equations of porous media type and on a systems of two PDEs that models the
sulfation of marble in monuments. Due to the nonlinear nature of the underlying
mathematical model, the use of a fixed point scheme is required and every step
implies the solution of large, locally structured, linear systems. A special
effort is devoted to the spectral analysis of the relevant matrices and to the
design of appropriate iterative or multi-iterative solvers, with special
attention to preconditioned Krylov methods and to multigrid procedures.
Numerical experiments for the validation of the analysis complement this
contribution.Comment: 26 pages, 13 figure
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
Method of lines transpose: High order L-stable O(N) schemes for parabolic equations using successive convolution
We present a new solver for nonlinear parabolic problems that is L-stable and
achieves high order accuracy in space and time. The solver is built by first
constructing a single-dimensional heat equation solver that uses fast O(N)
convolution. This fundamental solver has arbitrary order of accuracy in space,
and is based on the use of the Green's function to invert a modified Helmholtz
equation. Higher orders of accuracy in time are then constructed through a
novel technique known as successive convolution (or resolvent expansions).
These resolvent expansions facilitate our proofs of stability and convergence,
and permit us to construct schemes that have provable stiff decay. The
multi-dimensional solver is built by repeated application of dimensionally
split independent fundamental solvers. Finally, we solve nonlinear parabolic
problems by using the integrating factor method, where we apply the basic
scheme to invert linear terms (that look like a heat equation), and make use of
Hermite-Birkhoff interpolants to integrate the remaining nonlinear terms. Our
solver is applied to several linear and nonlinear equations including heat,
Allen-Cahn, and the Fitzhugh-Nagumo system of equations in one and two
dimensions
Fourth-order time-stepping for stiff PDEs on the sphere
We present in this paper algorithms for solving stiff PDEs on the unit sphere
with spectral accuracy in space and fourth-order accuracy in time. These are
based on a variant of the double Fourier sphere method in coefficient space
with multiplication matrices that differ from the usual ones, and
implicit-explicit time-stepping schemes. Operating in coefficient space with
these new matrices allows one to use a sparse direct solver, avoids the
coordinate singularity and maintains smoothness at the poles, while
implicit-explicit schemes circumvent severe restrictions on the time-steps due
to stiffness. A comparison is made against exponential integrators and it is
found that implicit-explicit schemes perform best. Implementations in MATLAB
and Chebfun make it possible to compute the solution of many PDEs to high
accuracy in a very convenient fashion
Hyperbolic/parabolic development for the GIM-STAR code
Flow fields in supersonic inlet configurations were computed using the eliptic GIM code on the STAR computer. Spillage flow under the lower cowl was calculated to be 33% of the incoming stream. The shock/boundary layer interaction on the upper propulsive surface was computed including separation. All shocks produced by the flow system were captured. Linearized block implicit (LBI) schemes were examined to determine their application to the GIM code. Pure explicit methods have stability limitations and fully implicit schemes are inherently inefficient; however, LBI schemes show promise as an effective compromise. A quasiparabolic version of the GIM code was developed using elastical parabolized Navier-Stokes methods combined with quasitime relaxation. This scheme is referred to as quasiparabolic although it applies equally well to hyperbolic supersonic inviscid flows. Second order windward differences are used in the marching coordinate and either explicit or linear block implicit time relaxation can be incorporated
- …