Mathematical and numerical modeling of coupled oxygen dynamics and neuronal electrophysiology
This work addresses brain ischemia modeling for neuroscience and medical research, but it is incremental as it builds on existing models with specific adaptations.
The study developed a multiscale model to simulate how oxygen dynamics affect neuronal excitability under ischemic conditions, linking cell swelling to local oxygen concentration and blood flow reduction, and demonstrated its application on idealized and realistic brain geometries.
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.