A pseudospectral method for investigating the stability of linear population models with two physiological structures

03/22/2022
by   Alessia andò, et al.
0

The asymptotic stability of the null equilibrium of a linear population model with two physiological structures formulated as a first-order hyperbolic PDE is determined by the spectrum of its infinitesimal generator. We propose an equivalent reformulation of the problem in the space of absolutely continuous functions in the sense of Carathéodory, so that the domain of the corresponding infinitesimal generator is defined by trivial boundary conditions. Via bivariate collocation, we discretize the reformulated operator as a finite-dimensional matrix, which can be used to approximate the spectrum of the original infinitesimal generator. Finally, we provide test examples illustrating the converging behavior of the approximated eigenvalues and eigenfunctions, and its dependence on the regularity of the model coefficients.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/21/2023

A numerical method for the stability analysis of linear age-structured models with nonlocal diffusion

We numerically investigate the stability of linear age-structured popula...
research
02/05/2021

Sharp stability for finite difference approximations of hyperbolic equations with boundary conditions

In this article, we consider a class of finite rank perturbations of Toe...
research
06/10/2022

Spectral analysis and fast methods for large matrices arising from PDE approximation

The main goal of this thesis is to show the crucial role that plays the ...
research
09/07/2021

Bivariate collocation for computing R_0 in epidemic models with two structures

Structured epidemic models can be formulated as first-order hyperbolic P...
research
08/18/2019

A linear system for pipe flow stability analysis allowing for boundary condition modifications

An accurate system to study the stability of pipe flow that ensures regu...
research
08/29/2023

Hyperbolicity of a semi-Lagrangian formulation of the hydrostatic free-surface Euler system

By a semi-Lagrangian change of coordinates, the hydrostatic Euler equati...
research
09/07/2020

Unfolding by Folding: a resampling approach to the problem of matrix inversion without actually inverting any matrix

Matrix inversion problems are often encountered in experimental physics,...

Please sign up or login with your details

Forgot password? Click here to reset