Abstract
Motivated by the termination of undesirable arrhythmia, a time optimal control formulation for the monodomain equations is proposed. It is shown that, under certain conditions, the optimal solutions of this problem steer the system into an appropriate stable neighborhood of the resting state. Towards this goal, some new regularity results and asymptotic properties for the monodomain equations with the Rogers-McCulloch ionic model are obtained. For the numerical realization, a monolithic approach, which simultaneously optimizes for the optimal times and optimal controls, is presented and analyzed. Its practical realization is based on a semismooth Newton method. Numerical examples and comparisons are included.
Original language | English |
---|---|
Pages (from-to) | 381-414 |
Number of pages | 34 |
Journal | ESAIM: Mathematical Modelling and Numerical Analysis |
Volume | 50 |
Issue number | 2 |
DOIs | |
State | Published - Mar 1 2016 |
Externally published | Yes |
Funding
The third author gratefully acknowledges the Austrian Science Fund (FWF) for financial support under SFB F32, “Mathematical Optimization and Applications in Biomedical Sciences”. Kunisch Karl 1 2 Pieper Konstantin 3 Rund Armin 1 The first and second author gratefully acknowledge support from the International Research Training Group IGDK 1754, funded by the German Science Foundation (DFG) and the Austrian Science Fund (FWF).
Funders | Funder number |
---|---|
German Science Foundation | |
Deutsche Forschungsgemeinschaft | |
Austrian Science Fund | SFB F32 |
Keywords
- Asymptotic behavior
- Monodomain equations
- Reaction diffusion system
- Semismooth Newton method
- Time optimal control