The Role of One-Dimensional Model-Generated Inter-Subject Variations in Systemic Properties on Wall Shear Stress Indices of Intracranial Aneurysms.

Variations in systemic properties of the arterial tree, such as aging-induced vessel stiffness, can alter the shape of pressure and flow waveforms. As a consequence, the hemodynamics around a cerebral aneurysm change, and therefore, also the corresponding in- and outlet boundary conditions (BCs) used for three-dimensional (3D) calculations of hemodynamic indices. In this study, we investigate the effects of variations in systemic properties on wall shear stress (WSS) indices of a cerebral aneurysm. We created a virtual patient database by varying systemic properties within physiological ranges. BCs for 3D-CFD simulations were derived using a pulse wave propagation model for each realization of the virtual database. WSS indices were derived from the 3D simulations and their variabilities quantified. Variations in BCs, caused by changes in systemic properties, yielded variabilities in the WSS indices that were of the same order of magnitude as differences in these WSS indices between ruptured and unruptured aneurysms. Sensitivity analysis showed that the systemic properties impacted both in- and outlet BCs simultaneously and altered the WSS indices. We conclude that the influence of variations in patient-specific systemic properties on WSS indices should be evaluated when using WSS indices in multidisciplinary rupture prediction models.

Application of an Adaptive Polynomial Chaos Expansion on Computationally Expensive Three-Dimensional Cardiovascular Models for Uncertainty Quantification and Sensitivity Analysis.

When applying models to patient-specific situations, the impact of model input uncertainty on the model output uncertainty has to be assessed. Proper uncertainty quantification (UQ) and sensitivity analysis (SA) techniques are indispensable for this purpose. An efficient approach for UQ and SA is the generalized polynomial chaos expansion (gPCE) method, where model response is expanded into a finite series of polynomials that depend on the model input (i.e., a meta-model). However, because of the intrinsic high computational cost of three-dimensional (3D) cardiovascular models, performing the number of model evaluations required for the gPCE is often computationally prohibitively expensive. Recently, Blatman and Sudret (2010, "An Adaptive Algorithm to Build Up Sparse Polynomial Chaos Expansions for Stochastic Finite Element Analysis," Probab. Eng. Mech., 25(2), pp. 183-197) introduced the adaptive sparse gPCE (agPCE) in the field of structural engineering. This approach reduces the computational cost with respect to the gPCE, by only including polynomials that significantly increase the meta-model's quality. In this study, we demonstrate the agPCE by applying it to a 3D abdominal aortic aneurysm (AAA) wall mechanics model and a 3D model of flow through an arteriovenous fistula (AVF). The agPCE method was indeed able to perform UQ and SA at a significantly lower computational cost than the gPCE, while still retaining accurate results. Cost reductions ranged between 70-80% and 50-90% for the AAA and AVF model, respectively.

A constitutive model for developing blood clots with various compositions and their nonlinear viscoelastic behavior.

The mechanical properties determine to a large extent the functioning of a blood clot. These properties depend on the composition of the clot and have been related to many diseases. However, the various involved components and their complex interactions make it difficult at this stage to fully understand and predict properties as a function of the components. Therefore, in this study, a constitutive model is developed that describes the viscoelastic behavior of blood clots with various compositions. Hereto, clots are formed from whole blood, platelet-rich plasma and platelet-poor plasma to study the influence of red blood cells, platelets and fibrin, respectively. Rheological experiments are performed to probe the mechanical behavior of the clots during their formation. The nonlinear viscoelastic behavior of the mature clots is characterized using a large amplitude oscillatory shear deformation. The model is based on a generalized Maxwell model that accurately describes the results for the different rheological experiments by making the moduli and viscosities a function of time and the past and current deformation. Using the same model with different parameter values enables a description of clots with different compositions. A sensitivity analysis is applied to study the influence of parameter variations on the model output. The relative simplicity and flexibility make the model suitable for numerical simulations of blood clots and other materials showing similar behavior.

A benchmark study of numerical schemes for one-dimensional arterial blood flow modelling.

Haemodynamical simulations using one-dimensional (1D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. Recent interest in verifying 1D numerical schemes has led to the development of alternative experimental setups and the use of three-dimensional numerical models to acquire data not easily measured in vivo. In most studies to date, only one particular 1D scheme is tested. In this paper, we present a systematic comparison of six commonly used numerical schemes for 1D blood flow modelling: discontinuous Galerkin, locally conservative Galerkin, Galerkin least-squares finite element method, finite volume method, finite difference MacCormack method and a simplified trapezium rule method. Comparisons are made in a series of six benchmark test cases with an increasing degree of complexity. The accuracy of the numerical schemes is assessed by comparison with theoretical results, three-dimensional numerical data in compatible domains with distensible walls or experimental data in a network of silicone tubes. Results show a good agreement among all numerical schemes and their ability to capture the main features of pressure, flow and area waveforms in large arteries. All the information used in this study, including the input data for all benchmark cases, experimental data where available and numerical solutions for each scheme, is made publicly available online, providing a comprehensive reference data set to support the development of 1D models and numerical schemes.

A constitutive modeling interpretation of the relationship among carotid artery stiffness, blood pressure, and age in hypertensive subjects.

Aging has a profound influence on arterial wall structure and function. We have previously reported the relationship among pulse wave velocity, age, and blood pressure in hypertensive subjects. In the present study, we aimed for a quantitative interpretation of the observed changes in wall behavior with age using a constitutive modeling approach. We implemented a model of arterial wall biomechanics and fitted this to the group-averaged pressure-area (P-A) relationship of the "young" subgroup of our study population. Using this model as our take-off point, we assessed which parameters had to be changed to let the model describe the "old" subgroup's P-A relationship. We allowed elastin stiffness and collagen recruitment parameters to vary and adjusted residual stress parameters according to published age-related changes. We required wall stress to be homogeneously distributed over the arterial wall and assumed wall stress normalization with age by keeping average "old" wall stress at the "young" level. Additionally, we required axial force to remain constant over the cardiac cycle. Our simulations showed an age-related shift in pressure-load bearing from elastin to collagen, caused by a decrease in elastin stiffness and a considerable increase in collagen recruitment. Correspondingly, simulated diameter and wall thickness increased by about 20 and 17%, respectively. The latter compared well with a measured thickness increase of 21%. We conclude that the physiologically realistic changes in constitutive properties we found under physiological constraints with respect to wall stress could well explain the influence of aging in the stiffness-pressure-age pattern observed.