453 research outputs found
Effective order strong stability preserving Runge–Kutta methods
We apply the concept of effective order to strong stability preserving (SSP) explicit Runge–Kutta methods. Relative to classical Runge–Kutta methods, effective order methods are designed to satisfy a relaxed set of order conditions, but yield higher order accuracy when composed with special starting and stopping methods. The relaxed order conditions allow for greater freedom in the design of effective order methods. We show that this allows the construction of four-stage SSP methods with effective order four (such methods cannot have classical order four). However, we also prove that effective order five methods—like classical order five methods—require the use of non-positive weights and so cannot be SSP. By numerical optimization, we construct explicit SSP Runge–Kutta methods up to effective order four and establish the optimality of many of them. Numerical experiments demonstrate the validity of these methods in practice
Development of a Three-Dimensional Air Blast Propagation Model Based Upon the Weighted Average Flux Method
Accurate numerical modeling of complex, multi-dimensional shock propagation is needed for many Department of Defense applications. A three-dimensional code, based upon E.F. Toro\u27s weighted average flux (WAF) method has been developed, tested, and validated. Code development begins with the introduction and application of all techniques in a single dimension. First-order accuracy is achieved via Godunov\u27s scheme using an exact Riemann solver. Adaptive techniques, which employ approximate solutions, are implemented to improve computational efficiency. The WAF method produces second-order accurate solutions, but introduces spurious oscillations near shocks and contact discontinuities. Total variation diminishing (TVD) flux and weight limiting schemes are added to reduce fluctuation severity. Finally, the fully developed one dimensional code is validated against experimental data, and extended into two and three dimensions via dimension-splitting technique
Strong stability preserving explicit Runge-Kutta methods of maximal effective order
We apply the concept of effective order to strong stability preserving (SSP)
explicit Runge-Kutta methods. Relative to classical Runge-Kutta methods,
methods with an effective order of accuracy are designed to satisfy a relaxed
set of order conditions, but yield higher order accuracy when composed with
special starting and stopping methods. We show that this allows the
construction of four-stage SSP methods with effective order four (such methods
cannot have classical order four). However, we also prove that effective order
five methods - like classical order five methods - require the use of
non-positive weights and so cannot be SSP. By numerical optimization, we
construct explicit SSP Runge-Kutta methods up to effective order four and
establish the optimality of many of them. Numerical experiments demonstrate the
validity of these methods in practice.Comment: 17 pages, 3 figures, 8 table
A realizability-preserving high-order kinetic scheme using WENO reconstruction for entropy-based moment closures of linear kinetic equations in slab geometry
We develop a high-order kinetic scheme for entropy-based moment models of a
one-dimensional linear kinetic equation in slab geometry. High-order spatial
reconstructions are achieved using the weighted essentially non-oscillatory
(WENO) method, and for time integration we use multi-step Runge-Kutta methods
which are strong stability preserving and whose stages and steps can be written
as convex combinations of forward Euler steps. We show that the moment vectors
stay in the realizable set using these time integrators along with a maximum
principle-based kinetic-level limiter, which simultaneously dampens spurious
oscillations in the numerical solutions. We present numerical results both on a
manufactured solution, where we perform convergence tests showing our scheme
converges of the expected order up to the numerical noise from the numerical
optimization, as well as on two standard benchmark problems, where we show some
of the advantages of high-order solutions and the role of the key parameter in
the limiter
High-order TVD and TVB linear multistep methods
We consider linear multistep methods that possess the TVD (total variation diminishing) or TVB (total variation bounded) properties, or related general monotonicity and boundedness properties. Strict monotonicity or TVD, in terms of arbitrary starting values for the multistep schemes, is only valid for a small class of methods, under very stringent step size restrictions. This makes them uncompetitive to the TVD Runge-Kutta methods. By relaxing these strict monotonicity requirements a larger class of methods can be considered, including many methods of practical interest. In this paper we construct linear multistep methods of high-order (up to six) that possess relaxed monotonicity or boundedness properties with optimal step size conditions. Numerical experiments show that the new schemes perform much better than the classical TVD multistep schemes. Moreover there is a substantial gain in efficiency compared to recently constructed TVD Runge-Kutta methods
- …