THE DISCRETE ADJOINT METHOD FOR LARGE-SCALE OPTIMIZATION PROBLEMS WITH LINEAR TIME-DEPENDENT PDE CONSTRAINTS

Kai Rothauge,E. Haber,U. Ascher
Abstract:We discuss the application of the discrete adjoint method to large-scale distributed parameterestimation and sensitivity analysis problems that are constrained by linear time-dependent PDEs. The PDEs are solved using an arbitrary-order linear multistep or Runge-Kutta method, making the resulting procedures for the computation of the action of the sensitivity matrix and of the Hessian well-suited to situations where the forward problem needs to be solved to high accuracy. These procedures will depend on the time-stepping method being used and we discuss their implementation in detail. Adjoint linear multistep methods are seen to be equivalent to the corresponding forward schemes apart from the handling of the source terms, and for adjoint Runge-Kutta methods we give a simple argument to show that their order of accuracy is the same as that of the corresponding forward methods. This property is illustrated with a numerical example. The approach used here can easily be adapted to nonlinear PDEs solved using other types of time-stepping methods.
What problem does this paper attempt to address?