A Coupled Diffusion Approximation for Spatiotemporal Hemodynamic Response and Deoxygenated Blood Volume Fraction in Microcirculation

Maryam Samavaki,Santtu Söderholm,Arash Zarrin Nia,Sampsa Pursiainen
2024-07-24
Abstract:Background and Objective: This article concerns a diffusion-based mathematical model for analyzing blood flow and oxygen transport within the capillaries, emphasizing its significance in understanding the physiological and biochemical dynamics of the cerebrovascular system and brain tissue. The focus of this study is, in particular, on neurovascular coupling and the spatiotemporal aspects of blood flow and oxygen transport in microcirculation. Methods: By adopting a coupled modelling approach that integrates the hemodynamic response function (HRF) with Fick's law and the Navier-Stokes equations (NSEs), we provide a computational framework for the diffusion-driven transport of deoxygenated and total blood volume fractions (DBV and TBV), essential for understanding blood oxygenation level-dependent functional magnetic resonance imaging (fMRI) and near-infrared spectroscopy (NIRS) applications. Results: The applicability of the model is further demonstrated through numerical experiments utilizing a 7 Tesla magnetic resonance imaging (MRI) dataset for head segmentation, which facilitates the differentiation of arterial blood vessels and various brain tissue compartments. By simulating hemodynamical responses and analyzing their impact on volumetric DBV and TBV, this study offers valuable insights into spatiotemporal modelling of brain tissue and blood flow. Conclusions: By integrating spatiotemporal modelling within a realistic head model derived from high-resolution 7 Tesla-MRI, we analyze the complex interplay between blood flow, oxygen transport, and brain tissue dynamics. This inclusion of a realistic head model not only enriches the accuracy of our simulations but is also beneficial for understanding the physiological and hemodynamic responses within the human brain.
Numerical Analysis,Analysis of PDEs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a diffusion - based mathematical model for analyzing the phenomena of blood flow and oxygen transport in capillaries, with special attention to neurovascular coupling and the spatio - temporal characteristics of blood flow and oxygen transport in the microcirculation. This research is of great significance for understanding the physiological and biochemical dynamics of the cerebrovascular system and brain tissue. Specifically, by combining the hemodynamic response function (HRF), Fick's law, and Navier - Stokes equations, the paper provides a computational framework for simulating the diffusion - driven transport processes of deoxygenated blood and total blood volume fractions (DBV and TBV). These processes are crucial for the applications of functional magnetic resonance imaging (fMRI) and near - infrared spectroscopy (NIRS). ### Main problems: 1. **Develop a mathematical model**: Establish a mathematical model that can describe blood flow and oxygen transport in the microcirculation, especially focusing on the influence of neural activity on local blood flow regulation. 2. **Understand physiological and biochemical dynamics**: Through this model, gain an in - depth understanding of the physiological and biochemical dynamics of the cerebrovascular system and brain tissue, especially blood flow and oxygen transport under the influence of neural activity. 3. **Apply high - resolution data**: Use 7 - Tesla magnetic resonance imaging (MRI) data sets for numerical experiments to verify the validity and accuracy of the model. ### Key points of the model: - **Hemodynamic response function (HRF)**: Describes the changes in blood flow caused by neural activity. - **Fick's law**: Describes the diffusion process of oxygen in tissues. - **Navier - Stokes equations**: Describes the dynamic behavior of arterial blood flow. ### Research methods: - **Model construction**: By coupling HRF, Fick's law, and Navier - Stokes equations, construct a mathematical model that describes blood flow and oxygen transport. - **Numerical experiments**: Use 7 - Tesla MRI data sets for head segmentation, distinguish different brain tissue regions, simulate hemodynamic responses, and analyze their influence on DBV and TBV. ### Results and conclusions: - **Model applicability**: The validity of the model has been verified through numerical experiments, showing its potential in simulating the spatio - temporal characteristics of brain tissue and blood flow. - **Physiological significance**: The model not only improves the accuracy of simulation but also deepens the understanding of human brain physiology and hemodynamic responses. ### Formula examples: - **Diffusion equation**: \[ \frac{\partial \tilde{q}}{\partial t}=-\nabla\cdot(\mathbf{u}\otimes\tilde{q})-\nabla\cdot(\mathbf{J}_{\tilde{q}})+s\quad\text{in }\hat{\Omega}\times[0, T] \] where \(\mathbf{u}\) is the blood velocity field, \(\mathbf{J}_{\tilde{q}} = -\delta\nabla\tilde{q}\) is the flux density of the deoxygenated blood volume fraction, and \(s\) is the source or sink term. - **Boundary conditions**: \[ g(\delta\nabla\tilde{q},\mathbf{n})=-g(\mathbf{J}_{\tilde{q}},\mathbf{n})=(1 - h)f\quad\text{on }B \] - **Hagen - Poiseuille model**: \[ \delta=\lambda\frac{A_a\bar{p}}{8\pi\mu} \] where \(\lambda=\frac{\xi_a}{\bar{\xi}_a}\), \(\xi_a\) is the microvascular length density, \(\bar{\xi}_a\) is the average value, and \(A_a\) is the cross - section of a single arteriole.