Patient-Specific Modeling Could Predict Occurrence of Pediatric Stroke

In silico analysis of embolism in cerebral arteries using fluid-structure interaction method

Ischemic stroke, particularly embolic stroke, stands as a significant global contributor to mortality and long-term disabilities. This paper presents a comprehensive simulation of emboli motion through the middle cerebral artery (MCA), a prevalent site for embolic stroke. Our patient-specific computational model integrates major branches of the middle cerebral artery reconstructed from magnetic resonance angiography images, pulsatile flow dynamics, and emboli of varying geometries, sizes, and material properties. The fluid-structure interactions method is employed to simulate deformable emboli motion through the middle cerebral artery, allowing observation of hemodynamic changes in artery branches upon embolus entry. We investigated the impact of embolus presence on shear stress magnitude on artery walls, analyzed the effects of embolus material properties and geometries on embolus trajectory and motion dynamics within the middle cerebral artery. Additionally, we evaluated the non-Newtonian behavior of blood, comparing it with Newtonian blood behavior. Our findings highlight that embolus geometry significantly influences trajectory, motion patterns, and hemodynamics within middle cerebral artery branches. Emboli with visco-hyperelastic material properties experienced higher stresses upon collision with artery walls compared to those with hyperelastic properties. Furthermore, considering blood as a non-Newtonian fluid had notable effects on emboli stresses and trajectories within the artery, particularly during collisions. Notably, the maximum von Mises stress experienced in our study was 21.83 kPa, suggesting a very low probability of emboli breaking during movement, impact, and after coming to a stop. However, in certain situations, the magnitude of shear stress on them exceeded 1 kPa, increasing the likelihood of cracking and disintegration. These results serve as an initial step in anticipating critical clinical conditions arising from arterial embolism in the middle cerebral artery. They provide insights into the biomechanical parameters influencing embolism, contributing to improved clinical decision-making for stroke management.

A parametric study of the effect of 3D plaque shape on local hemodynamics and implications for plaque instability

The vast majority of heart attacks occur when vulnerable plaques rupture, releasing their lipid content into the blood stream leading to thrombus formation and blockage of a coronary artery. Detection of these unstable plaques before they rupture remains a challenge. Hemodynamic features including wall shear stress (WSS) and wall shear stress gradient (WSSG) near the vulnerable plaque and local inflammation are known to affect plaque instability. In this work, a computational workflow has been developed to enable a comprehensive parametric study detailing the effects of 3D plaque shape on local hemodynamics and their implications for plaque instability. Parameterized geometric 3D plaque models are created within a patient-specific coronary artery tree using a NURBS (non-uniform rational B-splines)-based vascular modeling pipeline. Realistic blood flow features are simulated by using a Navier-Stokes solver within an isogeometric finite-element analysis framework. Near wall hemodynamic quantities such as WSS and WSSG are quantified, and vascular distribution of an inflammatory marker (VCAM-1) is estimated. Results show that proximally skewed eccentric plaques have the most vulnerable combination of high WSS and high positive spatial WSSG, and the presence of multiple lesions increases risk of rupture. The computational tool developed in this work, in conjunction with clinical data, -could help identify surrogate markers of plaque instability, potentially leading to a noninvasive clinical procedure for the detection of vulnerable plaques before rupture.

Image-guided subject-specific modeling of glymphatic transport and amyloid deposition

The glymphatic system is a brain-wide system of perivascular networks that facilitate exchange of cerebrospinal fluid (CSF) and interstitial fluid (ISF) to remove waste products from the brain. A greater understanding of the mechanisms for glymphatic transport may provide insight into how amyloid beta (A β ) and tau agglomerates, key biomarkers for Alzheimer's disease and other neurodegenerative diseases, accumulate and drive disease progression. In this study, we develop an image-guided computational model to describe glymphatic transport and A β deposition throughout the brain. A β transport and deposition are modeled using an advection-diffusion equation coupled with an irreversible amyloid accumulation (damage) model. We use immersed isogeometric analysis, stabilized using the streamline upwind Petrov-Galerkin (SUPG) method, where the transport model is constructed using parameters inferred from brain imaging data resulting in a subject-specific model that accounts for anatomical geometry and heterogeneous material properties. Both short-term (30-min) and long-term (12-month) 3D simulations of soluble amyloid transport within a mouse brain model were constructed from diffusion weighted magnetic resonance imaging (DW-MRI) data. In addition to matching short-term patterns of tracer deposition, we found that transport parameters such as CSF flow velocity play a large role in amyloid plaque deposition. The computational tools developed in this work will facilitate investigation of various hypotheses related to glymphatic transport and fundamentally advance our understanding of its role in neurodegeneration, which is crucial for the development of preventive and therapeutic interventions.

Rapid and Efficient Computation of Cell Paths During Ultrasonic Focusing

This biophysical analysis explores the first-principles physics of movement of white blood cell sized particles, suspended in an aqueous fluid and experiencing progressive or standing waves of acoustic pressure. In many current applications the cells are gradually nudged or herded toward the nodes of the standing wave, providing a degree of acoustic focusing and concentration of the cells in layers perpendicular to the direction of sound propagation. Here the underlying biomechanics of this phenomenon are analyzed specifically for the viscous regime of water and for small diameter microscopic spheroids such as living cells. The resulting mathematical model leads to a single algebraic expression for the creep or drift velocity as a function of sound frequency, amplitude, wavelength, fluid viscosity, boundary dimensions, and boundary reflectivity. This expression can be integrated numerically by a simple and fast computer algorithm to demonstrate net movement of particles as a function of time, providing a guide to optimization in a variety of emerging applications of ultrasonic cell focusing.