Variational data assimilation for transient blood flow simulations: Cerebral aneurysms as an illustrative example.

Several cardiovascular diseases are caused from localised abnormal blood flow such as in the case of stenosis or aneurysms. Prevailing theories propose that the development is caused by abnormal wall shear stress in focused areas. Computational fluid mechanics have arisen as a promising tool for a more precise and quantitative analysis, in particular because the anatomy is often readily available even by standard imaging techniques such as magnetic resonance and computed tomography angiography. However, computational fluid mechanics rely on accurate initial and boundary conditions, which are difficult to obtain. In this paper, we address the problem of recovering high-resolution information from noisy and low-resolution physical measurements of blood flow (for example, from phase-contrast magnetic resonance imaging [PC-MRI]) using variational data assimilation based on a transient Navier-Stokes model. Numerical experiments are performed in both 3D (2D space and time) and 4D (3D space and time) and with pulsatile flow relevant for physiological flow in cerebral aneurysms. The results demonstrate that, with suitable regularisation, the model accurately reconstructs flow, even in the presence of significant noise.

Adjoint multi-start-based estimation of cardiac hyperelastic material parameters using shear data.

Cardiac muscle tissue during relaxation is commonly modeled as a hyperelastic material with strongly nonlinear and anisotropic stress response. Adapting the behavior of such a model to experimental or patient data gives rise to a parameter estimation problem which involves a significant number of parameters. Gradient-based optimization algorithms provide a way to solve such nonlinear parameter estimation problems with relatively few iterations, but require the gradient of the objective functional with respect to the model parameters. This gradient has traditionally been obtained using finite differences, the calculation of which scales linearly with the number of model parameters, and introduces a differencing error. By using an automatically derived adjoint equation, we are able to calculate this gradient more efficiently, and with minimal implementation effort. We test this adjoint framework on a least squares fitting problem involving data from simple shear tests on cardiac tissue samples. A second challenge which arises in gradient-based optimization is the dependency of the algorithm on a suitable initial guess. We show how a multi-start procedure can alleviate this dependency. Finally, we provide estimates for the material parameters of the Holzapfel and Ogden strain energy law using finite element models together with experimental shear data.

Poro-elastic modeling of Syringomyelia - a systematic study of the effects of pia mater, central canal, median fissure, white and gray matter on pressure wave propagation and fluid movement within the cervical spinal cord.

Syringomyelia, fluid-filled cavities within the spinal cord, occurs frequently in association with a Chiari I malformation and produces some of its most severe neurological symptoms. The exact mechanism causing syringomyelia remains unknown. Since syringomyelia occurs frequently in association with obstructed cerebrospinal fluid (CSF) flow, it has been hypothesized that syrinx formation is mechanically driven. In this study we model the spinal cord tissue either as a poro-elastic medium or as a solid linear elastic medium, and simulate the propagation of pressure waves through an anatomically plausible 3D geometry, with boundary conditions based on in vivo CSF pressure measurements. Then various anatomic and tissue properties are modified, resulting in a total of 11 variations of the model that are compared. The results show that an open segment of the central canal and a stiff pia (relative to the cord) both increase the radial pressure gradients and enhance interstitial fluid flow in the central canal. The anterior median fissure, anisotropic permeability of the white matter, and Poisson ratio play minor roles.

Computational fluid dynamics evaluation of flow reversal treatment of giant basilar tip aneurysm.

Therapeutic parent artery flow reversal is a treatment option for giant, partially thrombosed basilar tip aneurysms. The effectiveness of this treatment has been variable and not yet studied by applying computational fluid dynamics. Computed tomography images and blood flow velocities acquired with transcranial Doppler ultrasonography were obtained prior to and after bilateral endovascular vertebral artery occlusion for a giant basilar tip aneurysm. Patient-specific geometries and velocity waveforms were used in computational fluid dynamics simulations in order to determine the velocity and wall shear stress changes induced by treatment. Therapeutic parent artery flow reversal lead to a dramatic increase in aneurysm inflow and wall shear stress (30 to 170 Pa) resulting in an increase in intra-aneurysmal circulation. The enlargement of the circulated area within the aneurysm led to a re-normalization of the wall shear stress and the aneurysm remained stable for more than 8 years thereafter. Therapeutic parent artery flow reversal can lead to unintended, potentially harmful changes in aneurysm inflow which can be quantified and possibly predicted by applying computational fluid dynamics.

Computation of hemodynamics in the circle of Willis.

BACKGROUND AND PURPOSE: Wall shear stress (WSS) and pressure are important factors in the development of cerebral aneurysms. We aimed to develop a computational fluid dynamics simulator for flow in the complete circle of Willis to study the impact of variations in vessel radii and bifurcation angles on WSS and pressure on vessel walls. METHODS: Blood flow was modeled with Navier-Stokes equations as an incompressible newtonian fluid within rigid vessel walls. A model of the circle of Willis geometry was approximated as a network of tubes around cubic curves. Pulsatile inlet flow rates and constant outlet pressure were used as boundary conditions. RESULTS: The simulations confirmed that differences in vessel radii and asymmetric branch angles influence WSS magnitude and spatial distribution. High WSS occurred at locations where aneurysms are frequent and in anatomic variants known to be associated with an increased risk for aneurysm development. CONCLUSIONS: Computational fluid dynamics analysis can be applied to the complete circle of Willis and should be used to study the pathophysiology of this complex vascular structure, including risk factors for aneurysm development. Further development of the method should include simulations with flexible vessel walls.