Fully-Coupled FSI Computational Analyses in the Ascending Thoracic Aorta Using Patient-Specific Conditions and Anisotropic Material Properties

Effect of turbulence and viscosity models on wall shear stress derived biomarkers for aorta simulations

Ascending aorta simulations provide insight into patient-specific hemodynamic conditions. Numerous studies have assessed fluid biomarkers which show a potential to aid clinicians in the diagnosis process. Unfortunately, there exists a large disparity in the computational methodology used to model turbulence and viscosity. Recognizing this disparity, some authors focused on analysing the influence of either the turbulence or viscosity models on the biomarkers in order to quantify the importance of these model choices. However, no analysis has yet been done on their combined effect. In order to fully understand and quantify the effect of the computational methodology, an assessment of the combined effect of turbulence and viscosity model choice was performed. Our results show that (1) non-Newtonian viscosity has greater impact (2.9-5.0%) on wall shear stress than Large Eddy Simulation turbulence modelling (0.1-1.4%), (2) the contribution of non-Newtonian viscosity is amplified when combined with a subgrid-scale turbulence model, (3) wall shear stress is underestimated when considering Newtonian viscosity by 2.9-5.0% and (4) cycle-to-cycle variability can impact the results as much as the numerical model if insufficient cycles are performed. These results demonstrate that, when assessing the effect of computational methodologies, the resultant combined effect of the different modelling assumptions differs from the aggregated effect of the isolated modifications. Accurate aortic flow modelling requires non-Newtonian viscosity and Large Eddy Simulation turbulence modelling.

Numerical investigation on circular and elliptical bulge tests for inverse soft tissue characterization

The acquisition of insights concerning the mechanobiology of aneurysmatic aortic tissues is an important field of investigation. The complete characterization of aneurysm mechanical behaviour can be carried out by biaxial experimental tests on ex vivo specimens. In literature, several works proposed bulge inflation tests as a valid method to analyse aneurysmatic tissue. Bulge test data processing requires the adoption of digital image correlation and inverse analysis approaches to estimate strain and stress distributions, respectively. In this context, however, the accuracy of inverse analysis method has not been evaluated yet. This aspect appears particularly interesting given the anisotropic behaviour of the soft tissue and the possibility to adopt different die geometries. The goal of this study is to provide an accuracy characterization of the inverse analysis applied to the bulge test technique using a numerical approach. In particular, different cases of bulge inflation were simulated in a finite element environment as a reference. To investigate the effect of tissue anisotropic degree and bulge die geometries (circular and elliptical), different input parameters were considered to obtain multiple test cases. The specimen deformed shapes, resulting from the reference finite element simulations, were then analysed through an inverse analysis approach to produce an estimation of stress distributions. The estimated stresses were, at last, compared with the values from the reference finite element simulations. The results demonstrated that the circular die geometry produces a satisfactory estimation accuracy only under certain conditions of material quasi-isotropy. On the other hand, the choice of an elliptical bulge die was proven to be more suitable for the analysis of anisotropic tissues.

Impact of the Spatial Velocity Inlet Distribution on the Hemodynamics of the Thoracic Aorta

The impact of the distribution in space of the inlet velocity in the numerical simulations of the hemodynamics in the thoracic aorta is systematically investigated. A real healthy aorta geometry, for which in-vivo measurements are available, is considered. The distribution is modeled through a truncated cone shape, which is a suitable approximation of the real one downstream of a trileaflet aortic valve during the systolic part of the cardiac cycle. The ratio between the upper and the lower base of the truncated cone and the position of the center of the upper base are selected as uncertain parameters. A stochastic approach is chosen, based on the generalized Polynomial Chaos expansion, to obtain accurate response surfaces of the quantities of interest in the parameter space. The selected parameters influence the velocity distribution in the ascending aorta. Consequently, effects on the wall shear stress are observed, confirming the need to use patient-specific inlet conditions if interested in the hemodynamics of this region. The surface base ratio is globally the most important parameter. Conversely, the impact on the velocity and wall shear stress in the aortic arch and descending aorta is almost negligible.

Enabling supra-aortic vessels inclusion in statistical shape models of the aorta: a novel non-rigid registration method

Statistical Shape Models (SSMs) are well-established tools for assessing the variability of 3D geometry and for broadening a limited set of shapes. They are widely used in medical imaging due to their ability to model complex geometries and their high efficiency as generative models. The principal step behind these techniques is a registration phase, which, in the case of complex geometries, can be a critical issue due to the correspondence problem, as it necessitates the development of correspondence mapping between shapes. The thoracic aorta, with its high level of morphological complexity, poses a multi-scale deformation problem due to the presence of several branch vessels with varying diameters. Moreover, branch vessels exhibit significant variability in shape, making the correspondence optimization even more challenging. Consequently, existing studies have focused on developing SSMs based only on the main body of the aorta, excluding the supra-aortic vessels from the analysis. In this work, we present a novel non-rigid registration algorithm based on optimizing a differentiable distance function through a modified gradient descent approach. This strategy enables the inclusion of custom, domain-specific constraints in the objective function, which act as landmarks during the registration phase. The algorithm's registration performance was tested and compared to an alternative Statistical Shape modeling framework, and subsequently used for the development of a comprehensive SSM of the thoracic aorta, including the supra-aortic vessels. The developed SSM was further evaluated against the alternative framework in terms of generalisation, specificity, and compactness to assess its effectiveness.

Extended law of laplace for measurement of the cloverleaf anatomy of the aortic root

The cross-sectional shape of the aortic root is cloverleaf, not circular, raising controversy regarding how best to measure its radiographic "diameter" for aortic event prediction. We mathematically extended the law of Laplace to estimate aortic wall stress within this cloverleaf region, simultaneously identifying a new metric of aortic root dimension that can be applied to clinical measurement of the aortic root and sinuses of Valsalva on clinical computerized tomographic scans. Enforcing equilibrium between blood pressure and wall stress, finite element computations were performed to evaluate the mathematical derivation. The resulting Laplace diameter was compared with existing methods of aortic root measurement across four patient groups: non-syndromic aneurysm, bicuspid aortic valve, Marfan syndrome, and non-dilated root patients (total 106 patients, 62 M, 44 F). (1) Wall stress: Mean wall stress at the depth of the sinuses followed this equation: Wall stress = BP × Circumscribing circle diameter/(2 × Aortic wall thickness). Therefore, the diameter of the circle enclosing the root cloverleaf, that is, twice the distance between the center, where the sinus-to-commissure lines coincide, and the depth of the sinuses, may replace diameter in the Laplace relation for a cloverleaf cross-section (or any shaped cross-section with two or more planes of symmetry). This mathematically derived result was verified by computational finite element analyses. (2) Diameters: CT scan measurements showed a significant difference between this new metric, the Laplace diameter, and the sinus-to-commissure, mid-sinus-to-mid-sinus, and coronal measurements in all four groups (p-value < 0.05). The average Laplace diameter measurements differed significantly from the other measurements in all patient groups. Among the various possible measurements within the aortic root, the diameter of the circumscribing circle, enclosing the cloverleaf, represents the diameter most closely related to wall stress. This diameter is larger than the other measurements, indicating an underestimation of wall stress by prior measurements, and otherwise provides an unbiased, convenient, consistent, physics-based measurement for clinical use. "Diameter" applies to circles. Our mathematical derivation of an extension of the law of Laplace, from circular to cloverleaf cross-sectional geometries of the aortic root, has implications for measurement of aortic root "diameter." The suggested method is as follows: (1) the "center" of the aortic root is identified by drawing three sinus-to-commissure lines. The intersection of these three lines identifies the "center" of the cloverleaf. (2) The largest radius from this center point to any of the sinuses is identified as the "radius" of the aortic root. (3) This radius is doubled to give the "diameter" of the aortic root. We find that this diameter best corresponds to maximal wall stress in the aortic root. Please note that this diameter defines the smallest circle that completely encloses the cloverleaf shape, touching the depths of all three sinuses.

In-silico investigations of haemodynamic parameters for a blunt thoracic aortic injury case

Accounting for 1.5% of thoracic trauma, blunt thoracic aortic injury (BTAI) is a rare disease with a high mortality rate that nowadays is treated mostly via thoracic endovascular aortic repair (TEVAR). Personalised computational models based on fluid-solid interaction (FSI) principals not only support clinical researchers in studying virtual therapy response, but also are capable of predicting eventual outcomes. The present work studies the variation of key haemodynamic parameters in a clinical case of BTAI after successful TEVAR, using a two-way FSI model. The three-dimensional (3D) patient-specific geometries of the patient were coupled with three-element Windkessel model for both prior and post intervention cases, forcing a correct prediction of blood flow over each section. Results showed significant improvement in velocity and pressure distribution after stenting. High oscillatory, low magnitude shear (HOLMES) regions require careful examination in future follow-ups, since thrombus formation was confirmed in some previously clinically reported cases of BTAI treated with TEVAR. The strength of swirling flows along aorta was also damped after stent deployment. Highlighting the importance of haemodynamic parameters in case-specific therapies. In future studies, compromising motion of aortic wall due to excessive cost of FSI simulations can be considered and should be based on the objectives of studies to achieve a more clinical-friendly patient-specific CFD model.

Impact of wall displacements on the large-scale flow coherence in ascending aorta

In the context of aortic hemodynamics, uncertainties affecting blood flow simulations hamper their translational potential as supportive technology in clinics. Computational fluid dynamics (CFD) simulations under rigid-walls assumption are largely adopted, even though the aorta contributes markedly to the systemic compliance and is characterized by a complex motion. To account for personalized wall displacements in aortic hemodynamics simulations, the moving-boundary method (MBM) has been recently proposed as a computationally convenient strategy, although its implementation requires dynamic imaging acquisitions not always available in clinics. In this study we aim to clarify the real need for introducing aortic wall displacements in CFD simulations to accurately capture the large-scale flow structures in the healthy human ascending aorta (AAo). To do that, the impact of wall displacements is analyzed using subject-specific models where two CFD simulations are performed imposing (1) rigid walls, and (2) personalized wall displacements adopting a MBM, integrating dynamic CT imaging and a mesh morphing technique based on radial basis functions. The impact of wall displacements on AAo hemodynamics is analyzed in terms of large-scale flow patterns of physiological significance, namely axial blood flow coherence (quantified applying the Complex Networks theory), secondary flows, helical flow and wall shear stress (WSS). From the comparison with rigid-wall simulations, it emerges that wall displacements have a minor impact on the AAo large-scale axial flow, but they can affect secondary flows and WSS directional changes. Overall, helical flow topology is moderately affected by aortic wall displacements, whereas helicity intensity remains almost unchanged. We conclude that CFD simulations with rigid-wall assumption can be a valid approach to study large-scale aortic flows of physiological significance.

Uncertainty Quantification in the In Vivo Image-Based Estimation of Local Elastic Properties of Vascular Walls

Introduction: Patient-specific computational models are a powerful tool for planning cardiovascular interventions. However, the in vivo patient-specific mechanical properties of vessels represent a major source of uncertainty. In this study, we investigated the effect of uncertainty in the elastic module (E) on a Fluid–Structure Interaction (FSI) model of a patient-specific aorta. Methods: The image-based χ-method was used to compute the initial E value of the vascular wall. The uncertainty quantification was carried out using the generalized Polynomial Chaos (gPC) expansion technique. The stochastic analysis was based on four deterministic simulations considering four quadrature points. A deviation of about ±20% on the estimation of the E value was assumed. Results: The influence of the uncertain E parameter was evaluated along the cardiac cycle on area and flow variations extracted from five cross-sections of the aortic FSI model. Results of stochastic analysis showed the impact of E in the ascending aorta while an insignificant effect was observed in the descending tract. Conclusions: This study demonstrated the importance of the image-based methodology for inferring E, highlighting the feasibility of retrieving useful additional data and enhancing the reliability of in silico models in clinical practice.

An image-based approach for the estimation of arterial local stiffness in vivo

The analysis of mechanobiology of arterial tissues remains an important topic of research for cardiovascular pathologies evaluation. In the current state of the art, the gold standard to characterize the tissue mechanical behavior is represented by experimental tests, requiring the harvesting of ex-vivo specimens. In recent years though, image-based techniques for the in vivo estimation of arterial tissue stiffness were presented. The aim of this study is to define a new approach to provide local distribution of arterial stiffness, estimated as the linearized Young's Modulus, based on the knowledge of in vivo patient-specific imaging data. In particular, the strain and stress are estimated with sectional contour length ratios and a Laplace hypothesis/inverse engineering approach, respectively, and then used to calculate the Young's Modulus. After describing the method, this was validated by using a set of Finite Element simulations as input. In particular, idealized cylinder and elbow shapes plus a single patient-specific geometry were simulated. Different stiffness distributions were tested for the simulated patient-specific case. After the validation from Finite Element data, the method was then applied to patient-specific ECG-gated Computed Tomography data by also introducing a mesh morphing approach to map the aortic surface along the cardiac phases. The validation process revealed satisfactory results. In the simulated patient-specific case, root mean square percentage errors below 10% for the homogeneous distribution and below 20% for proximal/distal distribution of stiffness. The method was then successfully used on the three ECG-gated patient-specific cases. The resulting distributions of stiffness exhibited significant heterogeneity, nevertheless the resulting Young's moduli were always contained within the 1-3 MPa range, which is in line with literature.

Experimental and Mouse-Specific Computational Models of the Fbln4SMKO Mouse to Identify Potential Biomarkers for Ascending Thoracic Aortic Aneurysm

PURPOSE

To use computational methods to explore geometric, mechanical, and fluidic biomarkers that could correlate with mouse lifespan in the Fbln4SMKO mouse. Mouse lifespan was used as a surrogate for risk of a severe cardiovascular event in cases of ascending thoracic aortic aneurysm.

METHODS

Image-based, mouse-specific fluid-structure-interaction models were developed for Fbln4SMKO mice (n = 10) at ages two and six months. The results of the simulations were used to quantify potential biofluidic biomarkers, complementing the geometrical biomarkers obtained directly from the images.

RESULTS

Comparing the different geometrical and biofluidic biomarkers to the mouse lifespan, it was found that mean oscillatory shear index (OSImin) and minimum time-averaged wall shear stress (TAWSSmin) at six months showed the largest correlation with lifespan (r2 = 0.70, 0.56), with both correlations being positive (i.e., mice with high OSImean and high TAWSSmin tended to live longer). When change between two and six months was considered, the change in TAWSSmin showed a much stronger correlation than OSImean (r2 = 0.75 vs. 0.24), and the correlation was negative (i.e., mice with increasing TAWSSmin over this period tended to live less long).

CONCLUSION

The results highlight potential biomarkers of ATAA outcomes that can be obtained through noninvasive imaging and computational simulations, and they illustrate the potential synergy between small-animal and computational models.

Introduction of a Novel Image-Based and Non-Invasive Method for the Estimation of Local Elastic Properties of Great Vessels

Background: In the context of a growing demand for the use of in silico models to meet clinical requests, image-based methods play a crucial role. In this study, we present a parametric equation able to estimate the elasticity of vessel walls, non-invasively and indirectly, from information uniquely retrievable from imaging. Methods: A custom equation was iteratively refined and tuned from the simulations of a wide range of different vessel models, leading to the definition of an indirect method able to estimate the elastic modulus E of a vessel wall. To test the effectiveness of the predictive capability to infer the E value, two models with increasing complexity were used: a U-shaped vessel and a patient-specific aorta. Results: The original formulation was demonstrated to deviate from the ground truth, with a difference of 89.6%. However, the adoption of our proposed equation was found to significantly increase the reliability of the estimated E value for a vessel wall, with a mean percentage error of 9.3% with respect to the reference values. Conclusion: This study provides a strong basis for the definition of a method able to estimate local mechanical information of vessels from data easily retrievable from imaging, thus potentially increasing the reliability of in silico cardiovascular models.

The Hemodynamic Effect of Modified Blalock–Taussig Shunt Morphologies: A Computational Analysis Based on Reduced Order Modeling

The Modified Blalock Taussig Shunt (MBTS) is one of the most common palliative operations in case of cyanotic heart diseases. Thus far, the decision on the position, size, and geometry of the implant relies on clinicians’ experience. In this paper, a Medical Digital Twin pipeline based on reduced order modeling is presented for fast and interactive evaluation of the hemodynamic parameters of MBTS. An infant case affected by complete pulmonary atresia was selected for this study. A three-dimensional digital model of the infant’s MBTS morphology was generated. A wide spectrum of MBTS geometries was explored by introducing twelve Radial Basis Function mesh modifiers. The combination of these modifiers allowed for analysis of various MBTS shapes. The final results proved the potential of the proposed approach for the investigation of significant hemodynamic features such as velocity, pressure, and wall shear stress as a function of the shunt’s morphology in real-time. In particular, it was demonstrated that the modifications of the MBTS morphology had a profound effect on the hemodynamic indices. The adoption of reduced models turned out to be a promising path to follow for MBTS numerical evaluation, with the potential to support patient-specific preoperative planning.

On the Role and Effects of Uncertainties in Cardiovascular in silico Analyses

The assessment of cardiovascular hemodynamics with computational techniques is establishing its fundamental contribution within the world of modern clinics. Great research interest was focused on the aortic vessel. The study of aortic flow, pressure, and stresses is at the basis of the understanding of complex pathologies such as aneurysms. Nevertheless, the computational approaches are still affected by sources of errors and uncertainties. These phenomena occur at different levels of the computational analysis, and they also strongly depend on the type of approach adopted. With the current study, the effect of error sources was characterized for an aortic case. In particular, the geometry of a patient-specific aorta structure was segmented at different phases of a cardiac cycle to be adopted in a computational analysis. Different levels of surface smoothing were imposed to define their influence on the numerical results. After this, three different simulation methods were imposed on the same geometry: a rigid wall computational fluid dynamics (CFD), a moving-wall CFD based on radial basis functions (RBF) CFD, and a fluid-structure interaction (FSI) simulation. The differences of the implemented methods were defined in terms of wall shear stress (WSS) analysis. In particular, for all the cases reported, the systolic WSS and the time-averaged WSS (TAWSS) were defined.