921 research outputs found
Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows
We present a review of the semi-Lagrangian method for advection-diusion and incompressible Navier-Stokes equations discretized with high-order methods. In particular, we compare the strong form where the departure points are computed directly via backwards integration with the auxiliary form where an auxiliary advection equation is solved instead; the latter is also referred to as Operator Integration Factor Splitting (OIFS) scheme. For intermediate size of time steps the auxiliary form is preferrable but for large time steps only the strong form is stable
Flux form Semi-Lagrangian methods for parabolic problems
A semi-Lagrangian method for parabolic problems is proposed, that extends
previous work by the authors to achieve a fully conservative, flux-form
discretization of linear and nonlinear diffusion equations. A basic consistency
and convergence analysis are proposed. Numerical examples validate the proposed
method and display its potential for consistent semi-Lagrangian discretization
of advection--diffusion and nonlinear parabolic problems
The remapped particle-mesh advection scheme
We describe the remapped particle-mesh method, a new mass-conserving method
for solving the density equation which is suitable for combining with
semi-Lagrangian methods for compressible flow applied to numerical weather
prediction. In addition to the conservation property, the remapped
particle-mesh method is computationally efficient and at least as accurate as
current semi-Lagrangian methods based on cubic interpolation. We provide
results of tests of the method in the plane, results from incorporating the
advection method into a semi-Lagrangian method for the rotating shallow-water
equations in planar geometry, and results from extending the method to the
surface of a sphere
Conservative and non-conservative methods based on hermite weighted essentially-non-oscillatory reconstruction for Vlasov equations
We introduce a WENO reconstruction based on Hermite interpolation both for
semi-Lagrangian and finite difference methods. This WENO reconstruction
technique allows to control spurious oscillations. We develop third and fifth
order methods and apply them to non-conservative semi-Lagrangian schemes and
conservative finite difference methods. Our numerical results will be compared
to the usual semi-Lagrangian method with cubic spline reconstruction and the
classical fifth order WENO finite difference scheme. These reconstructions are
observed to be less dissipative than the usual weighted essentially non-
oscillatory procedure. We apply these methods to transport equations in the
context of plasma physics and the numerical simulation of turbulence phenomena
Stability analysis for Eulerian and semi-Lagrangian finite-element formulation of the advection-diffusion equation
The article of record as published may be located at http://dx.doi.org/10.1016/S0898-1221(99)00185-6This paper analyzes the stability of the finite-element approximation to the linearized two-dimensional advection-diffusion equation. Bilinear basis functions on rectangular elements are considered. This is one of the two best schemes as was shown by Neta and Williams [1]. Time is discretized with the theta algorithms that yield the explicit (theta = 0), semi-implicit (theta = 1/2), and implicit (theta = 1) methods. This paper extends the results of Neta and Williams [1] for the advection equation. Giraldo and Neta [2] have numerically compared the Eulerian and semi-Lagrangian finite-element approximation for the advection-diffusion equation. This paper analyzes the finite element schemes used there. The stability analysis shows that the semi-Lagrangian method is unconditionally stable for all values of a while the Eulerian method is only unconditionally stable for 1/2 < theta < 1. This analysis also shows that the best methods are the semi-implicit ones (theta = 1/2). In essence this paper analytically compares a semi-implicit Eulerian method with a semi-implicit semi-Lagrangian method. It is concluded that (for small or no diffusion) the semi-implicit semi-Lagrangian method exhibits better amplitude, dispersion and group velocity errors than the semi-implicit Eulerian method thereby achieving better results. In the case the diffusion coefficient is large, the semi-Lagrangian loses its competitiveness. Published by Elsevier Science Ltd
The semi-Lagrangian method on curvilinear grids
International audienceWe study the semi-Lagrangian method on curvilinear grids. The classical backward semi-Lagrangian method [1] preserves constant states but is not mass conservative. Natural reconstruction of the field permits nevertheless to have at least first order in time conservation of mass, even if the spatial error is large. Interpolation is performed with classical cubic splines and also cubic Hermite interpolation with arbitrary reconstruction order of the derivatives. High odd order reconstruction of the derivatives is shown to be a good ersatz of cubic splines which do not behave very well as time step tends to zero. A conservative semi-Lagrangian scheme along the lines of [2] is then described; here conservation of mass is automatically satisfied and constant states are shown to be preserved up to first order in time
A conservative semi-Lagrangian method for oscillation-free computation of advection processes
The semi-Lagrangian method using the hybrid-cubic-rational interpolation
function [M. Ida, Comput. Fluid Dyn. J. 10 (2001) 159] is modified to a
conservative method by incorporating the concept discussed in [R. Tanaka et
al., Comput. Phys. Commun. 126 (2000) 232]. In the method due to Tanaka et al.,
not only a physical quantity but also its integrated quantity within a
computational cell are used as dependent variables, and the mass conservation
is completely achieved by giving a constraint to a forth-order polynomial used
as an interpolation function. In the present method, a hybrid-cubic-rational
function whose optimal mixing ratio was determined theoretically is employed
for the interpolation, and its derivative is used for updating the physical
quantity. The numerical oscillation appearing in results by the method due to
Tanaka et al. is sufficiently eliminated by the use of the hybrid function.Comment: 17 pages, 8 figures, accepted for publication in Comput. Phys.
Commun., Some misprint correcte
- …