The Mathematical Institute, University of Oxford, Eprints Archive

Discrete adjoint approximations with shocks

Giles, M. B. (2002) Discrete adjoint approximations with shocks. Technical Report. Unspecified. (Submitted)



This paper is concerned with the formulation and discretisation of adjoint equations when there are shocks in the underlying solution to the original nonlinear hyperbolic p.d.e. For the model problem of a scalar unsteady one-dimensional p.d.e. with a convex flux function, it is shown that the analytic formulation of the adjoint equations requires the imposition of an interior boundary condition along any shock. A 'discrete adjoint' discretisation is defined by requiring the adjoint equations to give the same value for the linearised functional as a linearisation of the original nonlinear discretisation. It is demonstrated that convergence requires increasing numerical smoothing of any shocks. Without this, any consistent discretisation of the adjoint equations without the inclusion of the shock boundary condition may yield incorrect values for the adjoint solution.

Item Type:Technical Report (Technical Report)
Subjects:H - N > Numerical analysis
Research Groups:Numerical Analysis Group
ID Code:1218
Deposited By: Lotti Ekert
Deposited On:19 May 2011 06:38
Last Modified:29 May 2015 18:54

Repository Staff Only: item control page