Modified Fletcher-Reeves and Dai-Yuan Conjugate Gradient Methods for Solving Optimal Control Problem of Monodomain Model
- 1 Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, Johor Bahru, Malaysia
- 2 Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, Johor Bahru, Malaysia
Abstract
In this paper, we present the numerical solution for the optimal control problem of monodomain modelwith Rogers-modified FitzHugh-Nagumo ion kinetic. The monodomain model, which is a well-known mathematical model for simulation of cardiac electrical activity, appears as the constraint in our problem. Our control objective is to dampen the excitation wavefront of the transmembrane potential in the observation domain using optimal applied current. Various conjugate gradient methods have been applied by researchers for solving this type of optimal control problem. For the present paper, we adopt the modified Fletcher-Reeves method and modified Dai-Yuan methodfor computing the optimal applied current. Numerical results show that the excitation wavefront is successfully dampened out by the optimal applied current when the modified Fletcher-Reeves method is used. However, this is not the case when the modified Dai-Yuan method is employed. Numerical results indicate that the modified Dai-Yuan method failed to converge to the optimal solution when the Armijo line search is used.
- Z. J. Zheng, J. B. Croft, W. H. Giles and G. A. Mensah, “Sudden Cardiac Death in the United States, 1989 to 1998,” Circulation, Vol. 104, 2001, pp. 2158-2163. Hdoi:10.1161/hc4301.098254
- E. H. M. Ong, “Proposal for Establishment of a National Sudden Cardiac Arrest Registry,” Singapore Medical Journal, Vol. 52, No. 8, 2011, pp. 631-633.
- S. Zhang, “Sudden Cardiac Death in China,” Pacing and Clinical Electrophysiology, Vol. 32, No. 9, 2009, pp. 1159-1162. Hdoi:10.1111/j.1540-8159.2009.02458.x
- C. Nagaiah, K. Kunisch and G. Plank, “Numerical Solution for Optimal Control of the Reaction-Diffusion Equations in Cardiac Electrophysiology,” Computational Optimization and Applications, Vol. 49, No. 1, 2011, pp. 149-178. Hdoi:10.1007/s10589-009-9280-3
- S. M. Shuaiby, M. A. Hassan and M. El-Melegy, “Modeling and Simulation of the Action Potential in Human Cardiac Tissues Using Finite Element Method,” Journal of Communications and Computer Engineering, Vol. 2, No. 3, 2012, pp. 21-27.
- Y. Belhamadia, A. Fortin and Y. Bourgault, “Towards Accurate Numerical Method for Monodomain Models Using a Realistic Heart Geometry,” Mathematical Biosciences, Vol. 220, No. 2, 2009, pp. 89-101. Hdoi:10.1016/j.mbs.2009.05.003
- C. Nagaiah and K. Kunisch, “Higher Order Optimization and Adaptive Numerical Solution for Optimal Control of Monodomain Equations in Cardiac Electrophysiology,” Applied Numerical Mathematics, Vol. 61, 2011, pp. 53-65. Hdoi:10.1016/j.apnum.2010.08.004
- E. Polak and G. Ribière, “Note Sur la Convergence de Méthodes de Directions Conjuguées,” Revue Francaise d’Informatique et de Recherche Opérationnelle, Vol. 16, 1969, pp. 35-43.
- B. T. Polyak, “The Conjugate Gradient Method in Extreme Problems;” USSR Computational Mathematics and Mathematical Physics, Vol. 9, 1969, pp. 94-112. Hdoi:10.1016/0041-5553(69)90035-4
- W. W. Hager and H. Zhang, “A New Conjugate Gradient Method with Guaranteed Descent and an Efficient Line Search,” SIAM Journal on Optimization, Vol. 16, 2005, pp. 170-192. Hdoi:10.1137/030601880
- Y. H. Dai and Y. Yuan, “A Nonlinear Conjugate Gradient Method with a Strong Global Convergence Property,” SIAM Journal on Optimization, Vol. 10, No. 1, 1999, pp. 177-182. Hdoi:10.1137/S1052623497318992
- K. W. Ng and A. Rohanin, “Numerical Solution for PDE-Constrained Optimization Problem in Cardiac Electrophysiology,” International Journal of Computer Applications, Vol. 44, No. 12, 2012, pp. 11-15.