Abstract
AbstractThe paper considers a time-adaptive finite element method for determination of drug efficacy in a parameter identification problem (PIP) for a system of ordinary differential equations (ODE) that describes dynamics of the primary human immunodeficiency virus (HIV) infection with drug therapy. Tikhonov’s regularization method, optimization approach and finite element method to solve this problem are presented. A posteriori error estimates in the Tikhonov’s functional and reconstructed parameter are derived. Based on these estimates a time adaptive algorithm is formulated and numerically tested for different scenarios of noisy observations of virus population function. Numerical results show a significant improvement of reconstruction of drug efficacy parameter using the local time-adaptive mesh refinement method compared to the gradient method applied on a uniform time mesh.
Funder
Vetenskapsrådet
Russian Foundation for Basic Research
University of Gothenburg
Publisher
Springer Science and Business Media LLC
Subject
Applied Mathematics,Analysis
Reference48 articles.
1. Albersmeyer, J.: Adjoint based algorithms and numerical methods for sensitivity generation and optimization of large scale dynamic systems, Ph.D. thesis, Heidelberg University (2010)
2. Almudevar, A., Arruda, E.F.: Optimal approximation schedules for a class of iterative algorithms, with an application to multigrid value iteration. IEEE Trans. Autom. Control 57(12) (2012)
3. Arruda, E.F., Dias, C.M., de Magalhaes, C.V., Pastore, D.H., Thomé, R.C.A., Yang, H.M.: An optimal control approach to HIV immunology. Appl. Math., 1115–1130 (2015). https://doi.org/10.4236/am.2015.66102
4. Bakushinskii, A.B., Kokurin, MYu.: Iterative Methods for Approximate Solution of Inverse Problems. Springer, New York (2004)
5. Bangerth, W., Joshi, A.: Adaptive finite element methods for the solution of inverse problems in optical tomography. Inverse Probl. 24, 034011 (2008)