Abstract: So-called corrected operator splitting methods are applied to a $1$-D scalar advection-diffusion equation of Buckley--Leverett type with \textit{general} initial data. Front tracking and a \textit{$2$nd} order Godunov method are used to advance the solution in time. Diffusion is modelled by piecewise linear finite elements at each new time level. To obtain correct structure of shock fronts independently of the size of the time step, a \textit{dynamically defined} residual flux term is grouped with diffusion. Different test problems are considered, and the methods are compared with respect to accuracy and runtime. Finally, we extend the corrected operator splitting to 2-D equations by means of dimensional splitting, and we apply it to a Buckley--Leverett type problem including gravitational effects.
mailto:conservation@math.ntnu.no Last modified: Mon Jan 19 09:47:08 1998