An explicit stable numerical scheme for the 1D transport equation

Yohan Penel (1, 2)

1 LETR - Laboratoire d'Etudes Thermiques des Réacteurs
2 LAGA - Laboratoire Analyse, Géométrie et Applications

We derive in this paper a numerical scheme in order to calculate solutions of 1D transport equations. This 2nd-order scheme is based on the method of characteristics and consists of two steps: The first step is about the approximation of the foot of the characteristic curve whereas the second one deals with the computation of the solution at this point. The main idea in our scheme is to combine two 2nd-order interpolation schemes so as to preserve the maximum principle. The resulting method is unconditionaly stable and designed for classical solutions but turns out to remain valid when shocks occur.

Mots clés: Method of characteristics; Linear advection equation; Numerical scheme; MOC2

Domaines: Analyse numérique [math.NA]; Modélisation et simulation

Yohan Penel. An explicit stable numerical scheme for the 1D transport equation. Discrete and Continuous Dynamical Systems - Series S, 2012, 5 (3), pp.641-656. ⟨10.3934/dcdss.2012.5.641⟩. ⟨hal-00523197⟩