Adaptive timestepping for conservation laws via adjoint error representation
- Adaptive Zeitschrittweitensteuerung für Erhaltungssätze mittels adjungierter Fehlerdarstellung
Steiner, Christina; Noelle, Sebastian (Thesis advisor)
Aachen : Publikationsserver der RWTH Aachen University (2008)
Dissertation / PhD Thesis
Aachen, Techn. Hochsch., Diss., 2008
We study a recent timestep adaptation technique for hyperbolic conservation laws. The core of the method is a space-time splitting of adjoint error representations for target functionals due to Süli  and Hartmann . It provides an efficient choice of timesteps for implicit computations of weakly instationary flows. The timestep will be very large in regions of stationary flow, and become small when a perturbation enters the flow field. Besides using adjoint techniques which are already well-established, we also add a new ingredient which simplifies the computation of the dual problem. Due to Galerkin orthogonality, the dual solution does not enter the error representation as such. Instead, the relevant term is the difference of the dual solution and its projection to the finite element space. We can show that it is therefore sufficient to compute the spatial gradient of the dual solution. This gradient satisfies a conservation law instead of a transport equation, and it can therefore be computed with the same algorithm as the forward problem, and in the same finite element space. For this new conservative approach we will derive boundary conditions. First we demonstrate the capabilities of the approach for a weakly instationary test problem for scalar, 1D conservation laws. Then we extend the computations to the 2D Euler equations, where we couple the adaptive time-stepping with spatial adaptation. For the spatial adaptation, we use a multiscale-based strategy developed by Müller , which we combine with the time adaptive method. The combined space-time adaptive method provides an efficient choice of timesteps for implicit computations of weakly instationary flows. The timestep will be very large in regions of stationary flow, and becomes small when a perturbation enters the flow field. The efficiency of the solver is investigated by means of an unsteady inviscid 2D flow over a bump.