In this paper, we explore the Lax-Wendroff (LW) type time discretization as an alternative procedure to the high order Runge-Kutta time discretization adopted for the high order essentially non-oscillatory (ENO) Lagrangian schemes developed in [2, 4]. The LW time discretization is based on a Taylor expansion in time, coupled with a local CauchyKowalewski procedure to utilize the partial differe...

Abstract We consider numerical boundary conditions for high order finite difference schemes for solving convection-diffusion equations on arbitrary geometry. The two main difficulties for numerical boundary conditions in such situations are: (1) the wide stencil of the high order finite difference operator requires special treatment for a few ghost points near the boundary; (2) the physical bou...

We present a general technique to solve Partial Differential Equations, called robust stencils, which make them tolerant to soft faults, i.e. bit flips arising in memory or CPU calculations. We show how it can be applied to a two-dimensional Lax-Wendroff solver. The resulting 2D robust stencils are derived using an orthogonal application of their 1D counterparts. Combinations of 3 to 5 base ste...

A parallel Hopscotch algorithm for the simulation of semiconductor laser arrays which overlaps the computations and communications is presented. The method uses a Lax-Wendroff scheme for the hyperbolic part of the operators and allows to obtain near-to-perfect speed-ups on MIMD computers. Numerical experiments indicate that semiconductor laser arrays exhibit a rich dynamic behavior ranging from...

A hardly known and very important result of Kreiss is proven explicitly: Outflow boundary extrapolation, which complements stable dissipative schemes for linear hyperbolic initial value problems, maintains stability. In view of this result, the Lax-Wendroff and the Gottlieb-Turkel schemes are applied to a test problem. As expected from the rate-of-convergence theory by Gustafsson, global order ...

We study complex-valued symmetric matrices. A simple expression for the spectral norm of such matrices is obtained, by utilizing a unitarily congruent invariant form. Consequently, we provide a sharp criterion for identifying those symmetric matrices whose spectral norm does not exceed one: such strongly stable matrices are usually sought in connection with convergent difference approximations ...

The advection equation is the basis for mathematical models of continuum mechanics. In the approximate solution of nonstationary problems it is necessary to inherit main properties of the conservatism and monotonicity of the solution. In this paper, the advection equation is written in the symmetric form, where the advection operator is the half-sum of advection operators in conservative (diver...

Based on general partitions of unity and standard numerical flux functions, a class of mesh-free methods for conservation laws is derived. A Lax-Wendroff type consistency analysis is carried out for the general case of moving partition functions. The analysis leads to a set of conditions which are checked for the finite volume particle method FVPM. As a by-product, classical finite volume schem...