Mathematical and numerical modeling of coupled oxygen dynamics and neuronal electrophysiology
Keywords
Advanced Numerical Methods for Scientific Computing
Living Systems and Precision Medicine
Code:
28/2026
Title:
Mathematical and numerical modeling of coupled oxygen dynamics and neuronal electrophysiology
Date:
Thursday 19th March 2026
Author(s):
Daniele, F.; Leimer Saglio, C. B.; Pagani, S.; Antonietti, P. F.
Download link:
Abstract:
Modeling and simulating how oxygen supply shapes neuronal excitability is crucial for advancing the understanding of brain function in pathological scenarios, such as ischemia. This condition is caused by a reduced blood supply, leading to the deprivation of oxygen and other metabolites; this energy deficit impairs ionic pumps and causes cellular swelling. In this work, this phenomenon is modeled through a volumetric variation law that links cell swelling to local oxygen concentration and the percentage of blood flow reduction. The swelling law links volume changes to local oxygen and the degree of blood-flow depression, providing a simple mechanistic pathway from hypoxia to tortuosity-driven transport impairment. The interplay between oxygen supply and excitability in brain tissue is described by coupling the monodomain model with specific neuronal ionic and metabolic models that characterize ion and metabolite concentration dynamics.
The numerical approximation of this coupled multiscale problem is particularly challenging, owing to the presence of sharp and fast-propagating wavefronts and complex geometrical domains.
To address these challenges, suitable space- and time-adaptive schemes are employed to capture the action potential dynamics accurately.
This multiscale model is discretized in space with the high-order p-adaptive discontinuous Galerkin method on polygonal and polyhedral grids (PolyDG) and integrated in time with a Crank-Nicolson scheme. We numerically investigate different pathological scenarios on a two-dimensional idealized square domain and on a realistic geometry, both discretized with a polygonal grid, analyzing how subclinical and severe ischemia can affect brain electrophysiology and metabolic concentrations.