On adaptive timestepping for weakly instationary solutions of hyperbolic conservation laws via adjoint error control

TL;DR

This paper introduces an adaptive timestep method for hyperbolic conservation laws using adjoint error control, optimizing computational efficiency by adjusting timestep sizes based on flow stationarity.

Contribution

It presents a novel approach that simplifies dual problem computation by focusing on the gradient of the dual solution, enabling efficient adaptive timestepping for weakly instationary flows.

Findings

Effective timestep adaptation in weakly instationary flows

Simplified dual problem computation using solution gradients

Demonstrated efficiency on scalar conservation law test problem

Abstract

We study a recent timestep adaptation technique for hyperbolic conservation laws. The key tool is a space-time splitting of adjoint error representations for target functionals due to S\"uli 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 {\phi} 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, {\phi}-{\phi}h . We can show that it is therefore sufficient to compute the spatial gradient of the dual solution, $w = {\nabla} {\phi}$. 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. We demonstrate the capabilities of the approach for a weakly instationary test problem for scalar conservation laws.