Back to Search Start Over

Automated adjoints of coupled PDE-ODE systems

Authors :
Farrell, Patrick E.
Hake, Johan E.
Funke, Simon W.
Rognes, Marie E.
Publication Year :
2017

Abstract

Mathematical models that couple partial differential equations (PDEs) and spatially distributed ordinary differential equations (ODEs) arise in biology, medicine, chemistry and many other fields. In this paper we discuss an extension to the FEniCS finite element software for expressing and efficiently solving such coupled systems. Given an ODE described using an augmentation of the Unified Form Language (UFL) and a discretisation described by an arbitrary Butcher tableau, efficient code is automatically generated for the parallel solution of the ODE. The high-level description of the solution algorithm also facilitates the automatic derivation of the adjoint and tangent linearization of coupled PDE-ODE solvers. We demonstrate the capabilities of the approach on examples from cardiac electrophysiology and mitochondrial swelling.

Subjects

Subjects :
Mathematics - Numerical Analysis

Details

Database :
arXiv
Publication Type :
Report
Accession number :
edsarx.1708.07648
Document Type :
Working Paper