13,995 research outputs found
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 Improved Split-Step Wavelet Transform Method for Anomalous Radio Wave Propagation Modelling
Anomalous tropospheric propagation caused by ducting phenomenon is a major problem in wireless communication. Thus, it is important to study the behavior of radio wave propagation in tropospheric ducts. The Parabolic Wave Equation (PWE) method is considered most reliable to model anomalous radio wave propagation. In this work, an improved Split Step Wavelet transform Method (SSWM) is presented to solve PWE for the modeling of tropospheric propagation over finite and infinite conductive surfaces. A large number of numerical experiments are carried out to validate the performance of the proposed algorithm. Developed algorithm is compared with previously published techniques; Wavelet Galerkin Method (WGM) and Split-Step Fourier transform Method (SSFM). A very good agreement is found between SSWM and published techniques. It is also observed that the proposed algorithm is about 18 times faster than WGM and provide more details of propagation effects as compared to SSFM
An almost symmetric Strang splitting scheme for nonlinear evolution equations
In this paper we consider splitting methods for the time integration of
parabolic and certain classes of hyperbolic partial differential equations,
where one partial flow can not be computed exactly. Instead, we use a numerical
approximation based on the linearization of the vector field. This is of
interest in applications as it allows us to apply splitting methods to a wider
class of problems from the sciences.
However, in the situation described the classic Strang splitting scheme,
while still a method of second order, is not longer symmetric. This, in turn,
implies that the construction of higher order methods by composition is limited
to order three only. To remedy this situation, based on previous work in the
context of ordinary differential equations, we construct a class of Strang
splitting schemes that are symmetric up to a desired order.
We show rigorously that, under suitable assumptions on the nonlinearity,
these methods are of second order and can then be used to construct higher
order methods by composition. In addition, we illustrate the theoretical
results by conducting numerical experiments for the Brusselator system and the
KdV equation
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
- …