Navier-Stokes Modelling of Non-Newtonian Blood Flow in Cerebral Arterial Circulation and its Dynamic Impact on Electrical Conductivity in a Realistic Multi-Compartment Head Model

Background and Objective: This study aims to evaluate the dynamic effect of non-Newtonian cerebral arterial circulation on electrical conductivity distribution (ECD) in a realistic multi-compartment head model. It addresses the importance and challenges associated with electrophysiological modalities, such as transcranial electrical stimulation, electro-magnetoencephalography, and electrical impedance tomography. Factors such as electrical conductivity's impact on forward modeling accuracy, complex vessel networks, data acquisition limitations (especially in MRI), and blood flow phenomena are considered. Methods: The Navier-Stokes equations (NSEs) govern the non-Newtonian flow model used in this study. The solver comprises two stages: the first solves the pressure field using a dynamical pressure-Poisson equation derived from NSEs, and the second updates the velocity field using Leray regularization and the pressure distribution from the first stage. The Carreau-Yasuda model establishes the connection between blood velocity and viscosity. Blood concentration in microvessels is approximated using Fick's law of diffusion, and conductivity mapping is obtained via Archie's law. The head model used corresponds to an open 7 Tesla MRI dataset, differentiating arterial vessels from other structures. Results: The results suggest the establishment of a dynamic model of cerebral blood flow for arterial and microcirculation. Blood pressure and conductivity distributions are obtained through numerically simulated pulse sequences, enabling approximation of blood concentration and conductivity within the brain. Conclusions: This model provides an approximation of dynamic blood flow and corresponding ECD in different brain regions. The advantage lies in its applicability with limited a priori information about blood flow and compatibility with arbitrary head models that distinguish arteries.Comment: 13 pages; 8 figures; 2 tabl

In silico study of the effects of cerebral circulation on source localization using a dynamical anatomical atlas of the human head

Objective. This study focuses on the effects of dynamical vascular modeling on source localization errors in electroencephalography (EEG). Our aim of this in silico study is to (a) find out the effects of cerebral circulation on the accuracy of EEG source localization estimates, and (b) evaluate its relevance with respect to measurement noise and interpatient variation. Approach. We employ a four-dimensional (3D + T) statistical atlas of the electrical properties of the human head with a cerebral circulation model to generate virtual patients with different cerebral circulatory conditions for EEG source localization analysis. As source reconstruction techniques, we use the linearly constraint minimum variance (LCMV) beamformer, standardized low-resolution brain electromagnetic tomography (sLORETA), and the dipole scan (DS). Main results. Results indicate that arterial blood flow affects source localization at different depths and with varying significance. The average flow rate plays an important role in source localization performance, while the pulsatility effects are very small. In cases where a personalized model of the head is available, blood circulation mismodeling causes localization errors, especially in the deep structures of the brain where the main cerebral arteries are located. When interpatient variations are considered, the results show differences up to 15 mm for sLORETA and LCMV beamformer and 10 mm for DS in the brainstem and entorhinal cortices regions. In regions far from the main arteries vessels, the discrepancies are smaller than 3 mm. When measurement noise is added and interpatient differences are considered in a deep dipolar source, the results indicate that the effects of conductivity mismatch are detectable even for moderate measurement noise. The signal-to-noise ratio limit for sLORETA and LCMV beamformer is 15 dB, while the limit is under 30 dB for DS. Significance. Localization of the brain activity via EEG constitutes an ill-posed inverse problem, where any modeling uncertainty, e.g. a slight amount of noise in the data or material parameter discrepancies, can lead to a significant deviation of the estimated activity, especially in the deep structures of the brain. Proper modeling of the conductivity distribution is necessary in order to obtain an appropriate source localization. In this study, we show that the conductivity of the deep brain structures is particularly impacted by blood flow-induced changes in conductivity because large arteries and veins access the brain through that region.Peer reviewe

L1-norm vs. L2-norm fitting in optimizing focal multi-channel tES stimulation : linear and semidefinite programming vs. weighted least squares

Background and Objective: This study focuses on Multi-Channel Transcranial Electrical Stimulation, a non-invasive brain method for stimulating neuronal activity under the influence of low-intensity currents. We introduce a mathematical formulation for finding a current pattern that optimizes an L1-norm fit between a given focal target distribution and volumetric current density inside the brain. L1-norm is well-known to favor well-localized or sparse distributions compared to L2-norm (least-squares) fitted estimates. Methods: We present a linear programming approach that performs L1-norm fitting and penalization of the current pattern (L1L1) to control the number of non-zero currents. The optimizer filters a large set of candidate solutions using a two-stage metaheuristic search from a pre-filtered set of candidates. Results: The numerical simulation results obtained with both 8- and 20-channel electrode montages suggest that our hypothesis on the benefits of L1-norm data fitting is valid. Compared to an L1-norm regularized L2-norm fitting (L1L2) via semidefinite programming and weighted Tikhonov least-squares method (TLS), the L1L1 results were overall preferable for maximizing the focused current density at the target position, and the ratio between focused and nuisance current magnitudes. Conclusions: We propose the metaheuristic L1L1 optimization approach as a potential technique to obtain a well-localized stimulus with a controllable magnitude at a given target position. L1L1 finds a current pattern with a steep contrast between the anodal and cathodal electrodes while suppressing the nuisance currents in the brain, hence, providing a potential alternative to modulate the effects of the stimulation, e.g., the sensation experienced by the subject.publishedVersionPeer reviewe

Multi-compartment head modeling in EEG: Unstructured boundary-fitted tetra meshing with subcortical structures.

This paper introduces an automated approach for generating a finite element (FE) discretization of a multi-compartment human head model for electroencephalographic (EEG) source localization. We aim to provide an adaptable FE mesh generation tool for EEG studies. Our technique relies on recursive solid angle labeling of a surface segmentation coupled with smoothing, refinement, inflation, and optimization procedures to enhance the mesh quality. In this study, we performed numerical meshing experiments with the three-layer Ary sphere and a magnetic resonance imaging (MRI)-based multi-compartment head segmentation which incorporates a comprehensive set of subcortical brain structures. These experiments are motivated, on one hand, by the sensitivity of non-invasive subcortical source localization to modeling errors and, on the other hand, by the present lack of open EEG software pipelines to discretize all these structures. Our approach was found to successfully produce an unstructured and boundary-fitted tetrahedral mesh with a sub-one-millimeter fitting error, providing the desired accuracy for the three-dimensional anatomical details, EEG lead field matrix, and source localization. The mesh generator applied in this study has been implemented in the open MATLAB-based Zeffiro Interface toolbox for forward and inverse processing in EEG and it allows for graphics processing unit acceleration

Mesh visualization: Ary sphere.

Quadrants of the downsampled surface grids (a-c) and the boundary-fitted tetrahedral mesh (e-g) for 3.0, 2.0 and 1.3 mm (millimeter) mesh sizes, respectively. Surface grids (d) and unfitted tetrahedral mesh (h) for 1.0 mm mesh size are included for comparison. The presented layers are (top-bot): scalp (brown), skull (white), and grey matter.</p

Mind map of the meshing process.

The unfitted mesh is obtained after the first labeling stage, while the boundary-fitting process includes additional processing stages for refinement, re-labeling, smoothing, inflation, and optimization via Delaunay turns. A graphics processing unit (GPU) can be applied to accelerate the solid angle labeling and re-labeling stages as well as the surface extraction stage which finds the compartment boundaries after labeling. The re-labeling process is run recursively as long as one or more compartment labels change their value. A pseudocode of this mind map is provided in data in S1 Appendix.</p