901 research outputs found
Combination of WENO and Explicit Runge–Kutta Methods for Wind Transport in the Meso-NH Model
This paper investigates the use of the weighted essentially nonoscillatory (WENO) space discretization methods of third and fifth order for momentum transport in the Meso-NH meteorological model, and their association with explicit Runge–Kutta (ERK) methods, with the specific purpose of finding an optimal combination in terms of wall-clock time to solution. A linear stability analysis using von Neumann theory is first conducted that considers six different ERK time integration methods. A new graphical representation of linear stability is proposed, which allows a first discrimination between the ERK methods. The theoretical analysis is then completed by tests on numerical problems of increasing complexity (linear advection of high wind gradient, orographic waves, density current, large eddy simulation of fog, and windstorm simulation), using a fourth-order-centered scheme as a reference basis. The five-stage third-order and fourth-order ERK combinations appear as the time integration methods of choice for coupling with WENO schemes in terms of stability. An explicit time-splitting method added to the ERK temporal scheme for WENO improves the stability properties slightly more. When the spatial discretizations are compared, WENO schemes present the main advantage of maintaining stable, nonoscillatory transitions with sharp discontinuities, but WENO third order is excessively damping, while WENO fifth order provides better accuracy. Finally, WENO fifth order combined with the ERK method makes the whole physics of the model 3 times faster compared to the classical fourth-order centered scheme associated with the leapfrog temporal scheme
Spatially partitioned embedded Runge-Kutta Methods
We study spatially partitioned embedded Runge–Kutta (SPERK) schemes for partial differential equations (PDEs), in which each of the component schemes is applied over a different part of the spatial domain. Such methods may be convenient for problems in which the smoothness of the solution or the magnitudes of the PDE coefficients vary strongly in space. We focus on embedded partitioned methods as they offer greater efficiency and avoid the order reduction that may occur in non-embedded schemes. We demonstrate that the lack of conservation in partitioned schemes can lead to non-physical effects and propose conservative additive schemes based on partitioning the fluxes rather than the ordinary differential equations. A variety of SPERK schemes are presented, including an embedded pair suitable for the time evolution of fifth-order weighted non-oscillatory (WENO) spatial discretizations. Numerical experiments are provided to support the theory
Level Set Jet Schemes for Stiff Advection Equations: The SemiJet Method
Many interfacial phenomena in physical and biological systems are dominated
by high order geometric quantities such as curvature.
Here a semi-implicit method is combined with a level set jet scheme to handle
stiff nonlinear advection problems.
The new method offers an improvement over the semi-implicit gradient
augmented level set method previously introduced by requiring only one
smoothing step when updating the level set jet function while still preserving
the underlying methods higher accuracy. Sample results demonstrate that
accuracy is not sacrificed while strict time step restrictions can be avoided
A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws
We introduce a new methodology for adding localized, space-time smooth,
artificial viscosity to nonlinear systems of conservation laws which propagate
shock waves, rarefactions, and contact discontinuities, which we call the
-method. We shall focus our attention on the compressible Euler equations in
one space dimension. The novel feature of our approach involves the coupling of
a linear scalar reaction-diffusion equation to our system of conservation laws,
whose solution is the coefficient to an additional (and artificial)
term added to the flux, which determines the location, localization, and
strength of the artificial viscosity. Near shock discontinuities, is
large and localized, and transitions smoothly in space-time to zero away from
discontinuities. Our approach is a provably convergent, spacetime-regularized
variant of the original idea of Richtmeyer and Von Neumann, and is provided at
the level of the PDE, thus allowing a host of numerical discretization schemes
to be employed. We demonstrate the effectiveness of the -method with three
different numerical implementations and apply these to a collection of
classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the
Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a
classical continuous finite-element implementation using second-order
discretization in both space and time, FEM-C. Second, we use a simplified WENO
scheme within our -method framework, WENO-C. Third, we use WENO with the
Lax-Friedrichs flux together with the -equation, and call this WENO-LF-C.
All three schemes yield higher-order discretization strategies, which provide
sharp shock resolution with minimal overshoot and noise, and compare well with
higher-order WENO schemes that employ approximate Riemann solvers,
outperforming them for the difficult Leblanc shock tube experiment.Comment: 34 pages, 27 figure
Strong Stability Preserving Two-Step Runge-Kutta Methods
We investigate the strong stability preserving (SSP) property of two-step Runge– Kutta (TSRK) methods. We prove that all SSP TSRK methods belong to a particularly simple\ud
subclass of TSRK methods, in which stages from the previous step are not used. We derive simple order conditions for this subclass. Whereas explicit SSP Runge–Kutta methods have order at most four, we prove that explicit SSP TSRK methods have order at most eight. We present TSRK methods of up to eighth order that were found by numerical search. These methods have larger SSP coefficients than any known methods of the same order of accuracy, and may be implemented in a form with relatively modest storage requirements. The usefulness of the TSRK methods is demonstrated through numerical examples, including integration of very high order WENO discretizations
- …