Numerical schemes for solving and optimizing multiscale models with age of hepatitis C virus dynamics

Vladimir Reinharz, Harel Dahari, Danny Barash

Research output: Contribution to journalArticlepeer-review

9 Scopus citations

Abstract

Age-structured PDE models have been developed to study viral infection and treatment. However, they are notoriously difficult to solve. Here, we investigate the numerical solutions of an age-based multiscale model of hepatitis C virus (HCV) dynamics during antiviral therapy and compare them with an analytical approximation, namely its long-term approximation. First, starting from a simple yet flexible numerical solution that also considers an integral approximated over previous iterations, we show that the long-term approximation is an underestimate of the PDE model solution as expected since some infection events are being ignored. We then argue for the importance of having a numerical solution that takes into account previous iterations for the associated integral, making problematic the use of canned solvers. Second, we demonstrate that the governing differential equations are stiff and the stability of the numerical scheme should be considered. Third, we show that considerable gain in efficiency can be achieved by using adaptive stepsize methods over fixed stepsize methods for simulating realistic scenarios when solving multiscale models numerically. Finally, we compare between several numerical schemes for the solution of the equations and demonstrate the use of a numerical optimization scheme for the parameter estimation performed directly from the equations.

Original languageEnglish
Pages (from-to)1-13
Number of pages13
JournalMathematical Biosciences
Volume300
DOIs
StatePublished - 1 Jun 2018

Keywords

  • Adaptive stepsize integration
  • Age-structured model
  • Hepatitis C virus
  • Multiscale model
  • Numerical solution
  • Partial differential equations

ASJC Scopus subject areas

  • Statistics and Probability
  • Modeling and Simulation
  • Biochemistry, Genetics and Molecular Biology (all)
  • Immunology and Microbiology (all)
  • Agricultural and Biological Sciences (all)
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Numerical schemes for solving and optimizing multiscale models with age of hepatitis C virus dynamics'. Together they form a unique fingerprint.

Cite this