Closed form solution and equivalent equation approximation of linear advection by a non dissipative second order scheme for step initial conditions

Grégory Arbia (1) , Daniel Bouche (2)

1 REO - Numerical simulation of biological flows
2 CMLA - Centre de Mathématiques et de Leurs Applications

We analyse the solution of the linear advection equation on a uniform mesh by a non dissipative second order scheme for discontinuous initial condition. These schemes are known to generate parasitic oscillations in the vicinity of the discontinuity. An approximate way to predict these oscillations is provided by the equivalent equation method. More specifically, we focus on the case of advection of a step function by the leapfrog scheme. Numerical experiments show that the equivalent equation method fails to reproduce the oscillations generated by the scheme far from the discontinuity. Thus, we derive closed form exact and approximate solutions for the scheme that accurately predict these oscillations. We study the relationship between equivalent equation approximation and exact solution for the scheme, to determine its range of validity.

Mots clés: Advection equation; Leapfrog scheme; Equivalent equation; Parasitic oscillations; Fourier analysis

Domaines: Analyse numérique [cs.NA]

Grégory Arbia, Daniel Bouche. Closed form solution and equivalent equation approximation of linear advection by a non dissipative second order scheme for step initial conditions. Acta Applicandae Mathematicae, 2013, ⟨10.1007/s10440-013-9844-1⟩. ⟨hal-00918645⟩