Physiologicallystructuredpopulationmodelsaretypicallyformulatedasapartial differential equation of transport type for the density, with a boundary condition describing the birth of new individuals. Here we develop numerical bifurcation methods by combin- ing pseudospectral approximate reduction to a finite dimensional system with the use of established tools for ODE. A key preparatory step is to view the density as the derivative of the cumulative distribution. To demonstrate the potential of the approach, we consider two classes of models: a size-structured model for waterfleas (Daphnia) and a maturity- structured model for cell proliferation. Using the package MatCont, we compute numer- ical bifurcation diagrams, like steady-state stability regions in a two-parameter plane and parametrized branches of equilibria and periodic solutions. Our rather positive conclusion is that a rather low dimension may yield a rather accurate diagram! In addition we show numerically that, for the two models considered here, equilibria of the approximating system converge to the true equilibrium as the dimension of the approximating system increases; this last result is also proved theoretically under some regularity conditions on the model ingredients.
Numerical Bifurcation Analysis of Physiologically Structured Population Models via Pseudospectral Approximation
Rossana Vermiglio;Dimitri BredaMembro del Collaboration Group
;
2021-01-01
Abstract
Physiologicallystructuredpopulationmodelsaretypicallyformulatedasapartial differential equation of transport type for the density, with a boundary condition describing the birth of new individuals. Here we develop numerical bifurcation methods by combin- ing pseudospectral approximate reduction to a finite dimensional system with the use of established tools for ODE. A key preparatory step is to view the density as the derivative of the cumulative distribution. To demonstrate the potential of the approach, we consider two classes of models: a size-structured model for waterfleas (Daphnia) and a maturity- structured model for cell proliferation. Using the package MatCont, we compute numer- ical bifurcation diagrams, like steady-state stability regions in a two-parameter plane and parametrized branches of equilibria and periodic solutions. Our rather positive conclusion is that a rather low dimension may yield a rather accurate diagram! In addition we show numerically that, for the two models considered here, equilibria of the approximating system converge to the true equilibrium as the dimension of the approximating system increases; this last result is also proved theoretically under some regularity conditions on the model ingredients.File | Dimensione | Formato | |
---|---|---|---|
VJOM-D-19-00208R1-2.pdf
non disponibili
Descrizione: versione accettata
Tipologia:
Documento in Pre-print
Licenza:
Non pubblico
Dimensione
488.36 kB
Formato
Adobe PDF
|
488.36 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.