Abstract

Abdominal aortic aneurysm (AAA) can lead to high mortality rates and further complications such as stroke or heart attack due to the risk of rupture and thrombosis. Wall mechanics play a crucial role in the development and progression of aneurysms. This study investigated the effects of wall mechanics on hemodynamic parameters in AAA to understand the risk of rupture and thrombosis. The impact of three aortic wall models (rigid, linear elastic, and hyperelastic) on structural and hemodynamic parameters was examined using CFD and FSI techniques. The blood was modeled using the Carreau non-Newtonian model, and the flow was simulated using the k-ω model. Physiological pulses were used for the velocity at the inlet and the pressure at the outlet. The results demonstrated close similarity between the predictions of the linear elastic and hyperelastic models, in contrast to the somewhat different results of the rigid model. The hyperelastic model predicted higher deformation and von Mises stress levels than the elastic model, although the difference in stress predictions was smaller than the difference in deformation predictions. The rigid model evaluated the time-averaged wall shear stress and oscillatory shear index higher than the other two models in the aneurysmal area but with a lower relative residence time. In general, the hyperelastic model predicted a higher risk of rupture than linear elastic models and a higher risk of thrombus formation than the other two models. The rigid model had the most optimistic prediction.

1. Introduction

Abdominal aortic aneurysm (AAA) is one of the prevalent cardiovascular disorders in the circulatory system and is commonly found in the infrarenal segment of the abdominal aorta, between the renal arteries and the iliac bifurcation. Aneurysm refers to the permanent and irreversible dilation of the arterial walls during which the arterial wall experiences significant degradation in its materials. In an aneurysm, at least three different stages are distinguishable as follows: genesis, progression, and rupture. The exact initiation process is yet to be understood; however, it is reported that the degradation of elastin and significant loss of the media layer in its smooth muscle cell are the underlying processes observed in the initiation of the aneurysm [1, 2]. In other words, the media layer in an aneurysmatic arterial wall approximately vanishes and the arterial wall becomes thinner. Therefore, any factor contributing to the weakening of the arterial wall can trigger the aneurysm. During the progression, an aortic aneurysm reaches at least 1.5 times the healthy aorta in diameter and can swell to four times its initial diameter [3]. The current guideline for clinical intervention is the maximum transverse diameter of 50 to 55 mm [4]. However, rupture is seen in cases below this threshold, and patients with quite larger geometrical characteristics than the threshold are reported to experience no rupture during their lifetime [4].

Computational modeling has become a valuable approach for studying abdominal aortic aneurysm (AAA) due to its numerous advantages over traditional experimental and clinical studies [5]. Computational methods offer a unique opportunity to simulate and analyze complex hemodynamic and structural processes occurring in the arterial wall and blood flow, providing insights into the mechanisms and pathophysiology of AAA [6]. Moreover, computational modeling allows for the evaluation of various wall mechanics and blood flow conditions, which are difficult to achieve in traditional studies [7]. Additionally, computational simulations can be conducted with reduced costs and time, without exposing patients to potential risks associated with invasive procedures or drug administration [8, 9].

In light of the critical role that advanced computational modeling plays in the analysis and treatment of cardiovascular disorders, particularly abdominal aortic aneurysm (AAA), recent research emphasizes the necessity of integrating state-of-the-art methodologies and findings. Innovative computational techniques, as detailed in recent studies, have shown promising advancements in simulating the hemodynamic effects of treatments and understanding the intricate dynamics within the aneurysm sac [10]. Moreover, these studies offer insights into optimizing treatment strategies through detailed analyses of stent design parameters, directly contributing to the precision and effectiveness of AAA interventions [11]. This integration of cutting-edge computational research underscores the potential for significant improvements in AAA management, highlighting the importance of incorporating these advancements into our study to advance the field further.

Blood is comprised of plasma, minerals, and red blood cells, and its viscosity decreases as the shear rate increases due to the shear-thinning behavior that is attributed to the deformation of red blood cells in the flow. Several studies have been conducted in recent years to establish an appropriate model for blood fluid behavior, with a focus on examining the WSS-dependent parameters for different rheological models [12, 13]. These investigations have revealed that the Carreau non-Newtonian model is better at predicting these parameters. Therefore, the Carreau model has been suggested as the most suitable model for blood in various studies [7, 14] and will also be used in this study.

In the complex domain of cardiovascular studies, particularly in examining abdominal aortic aneurysms (AAAs), the accurate simulation of flow phenomena becomes paramount. The determination of a turbulent flow regime within the AAA leverages the unique geometric characteristics and pulsatile nature of blood flow inherent to this condition. In particular, the dilation and irregular shapes common to AAAs foster conditions conducive to flow separation, recirculation, and vortex shedding, underscoring the necessity of employing turbulence models to simulate the intricate hemodynamics accurately [15]. These phenomena, emblematic of turbulent flows, alongside the pulsatile flow’s variable velocities and directions, substantiate the turbulent nature of blood flow within AAAs, thereby making turbulence modeling an indispensable tool in these simulations [16, 17]. Among the various turbulence models available, the k-omega (k-ω) model is distinguished by its adeptness at capturing boundary layer phenomena and its efficiency in integrating into the viscous sublayer without requiring additional damping functions. This capability renders the k-ω model particularly advantageous for cardiovascular applications, where a granular understanding of near-wall flow behavior exerts a significant influence on wall shear stress distributions and subsequent biological responses [13, 1820]. The choice of the k-ω model over others, such as the k-epsilon (k-ε), is informed by its demonstrated precision in depicting near-wall turbulence—a critical aspect in cardiovascular flows, including those altered by AAA. This model’s preference in computational studies on AAA and other cardiovascular conditions is further justified by its role in enhancing the interaction between blood components and endothelial cells, thereby facilitating platelet activation and the accumulation of blood lipids in the intima [21, 22], linking turbulent flow characteristics directly to the processes contributing to thrombosis formation [23].

Structural parameters such as deformation of the artery wall and von Mises stress are critical in evaluating the rupture risk of abdominal aortic aneurysm (AAA) in computational studies. High deformation and von Mises stress can lead to mechanical failure and rupture [24, 25]. Several studies have demonstrated a significant correlation between these structural parameters and an increased risk of rupture [2629]. However, the choice of model for the aortic wall, whether a linear elastic or hyperelastic model, can significantly affect the evaluation of rupture risk in computational studies of AAA [30]. Hyperelastic models can better capture the nonlinear behavior of the arterial wall and have been shown to lead to more accurate predictions of rupture risk compared with linear elastic models [25, 31]. Hyperelastic models generally predict higher deformation and von Mises stress values than linear elastic models, especially for larger aneurysms and higher pressures. Overall, hyperelastic models may provide a more accurate prediction of rupture risk in AAA [32, 33].

Thrombosis is a common complication associated with AAA, with a reported prevalence ranging from up to 70% depending on the aneurysm’s size and location. It can lead to blockage of the aneurysm sac and cause ischemia or embolization. To assess the risk of thrombosis occurrence in AAA, three hemodynamic parameters that depend on wall shear stress (WSS) are useful, including time-averaged WSS (TAWSS), the oscillating shear index (OSI), and the relative residence time (RRT). Low and high TAWSS values, as well as high OSI values, indicate an increased risk of thrombus formation, whereas high RRT values indicate stagnant flow that can contribute to thrombosis. The choice of wall mechanics model can affect the evaluation of thrombosis risk in AAA using hemodynamic parameters. Rigid wall models tend to overestimate WSS and provide inaccurate predictions of hemodynamic parameters, while elastic and hyperelastic models better capture complex wall motion, leading to more accurate predictions of hemodynamic parameters and thrombosis risk [34, 35]. The choice of wall mechanics model should be considered carefully when assessing thrombosis risk in AAA using hemodynamic parameters [25, 36].

Prior research in the domain of abdominal aortic aneurysm (AAA) has often relied on both rigid and elastic wall models, with many studies introducing simplifications in aspects such as geometry, flow dynamics, or material characteristics. These simplifications have historically constrained the depth of understanding attainable regarding the predictive accuracy of these models for rupture and thrombosis risks under conditions that closely mimic clinical realities. In addressing this gap, this study advances the field by employing a comprehensive approach that integrates a realistic AAA geometry, the k-omega turbulent flow model, the Carreau model for non-Newtonian blood behavior, and physiological pulse simulations as boundary conditions. This methodological rigor enabled a nuanced comparison between rigid and elastic wall models, elucidating their respective impacts on structural integrity and hemodynamic factors critical to assessing rupture and thrombosis risks. Significantly, the findings of the present study not only challenge simplifying modeling assumptions but also help achieve more accurate, personalized diagnostic and therapeutic strategies in cardiovascular care. By bridging these methodological divides, this study contributes a pivotal resource for future research and clinical practice, underscoring the critical role of sophisticated simulation techniques in enhancing our understanding and management of AAA.

2. Methodology

2.1. Governing Equations

The general form of the incompressible unsteady Reynolds-averaged Navier–Stokes equations, which are used to determine fluid motion, is as follows [20, 37, 38]:

For laminar flows, the term in (2) is excluded. However, for turbulent flows, this term is modeled as follows:where S denotes the modulus of the mean strain rate tensor, represents the turbulence viscosity of blood, and ρ refers to the density. The determination of in the k-ω turbulence model used for this study is expressed as follows:

The correction for the low Reynolds number is triggered by the coefficient , which functions to dampen turbulent viscosity [37]. The values of k and ω are computed by the following equations, respectively.

Reference [37] provides the constants and parameters used in (4)–(6).

The typical approach to describe the mechanical characteristics of the artery wall is to use the solid elastodynamic momentum equation (39):

In the equation, refers to the vector indicating the displacement of the solid, represents the Cauchy stress tensor, f is the vector representing external forces exerted on the solid, and is the volumetric mass of the solid.

2.2. Problem Statement

The present study involves a fluid-structure interaction (FSI) analysis of an aneurysmal aorta using a realistic geometry obtained from CT scan data. The primary objective of this research is to investigate the impact of different models for the artery wall on structural and hemodynamic parameters. The three models under consideration for the artery wall are solid, linear elastic, and hyperelastic.

The fluid and solid domains’ geometries are depicted in Figure 1.

Physiological pulses were employed at the inlet and the outlet to create realistic boundary conditions, and the inlet velocity pulse and the outlet pressure pulse are illustrated in Figures 2(a) and 2(b), respectively [40, 41].

The solid model was investigated using a computational fluid dynamics (CFD) study, whereas FSI studies were used for the elastic and hyperelastic models. This study assesses the performance of different wall models and compares their impact on the structural and hemodynamic parameters, which are vital for the proper diagnosis and treatment of aneurysmal aorta.

The analysis and visualization of flow were carried out using the ANSYS Fluent 2020 R1 commercial CFD software. The third-order MUSCL differencing scheme was employed for discretizing the equations. The pressure-velocity correction method was coupled to the simulations. A residual error convergence threshold of 1e − 6 was set for the momentum equations. To ensure reliable results and eliminate transient effects, unsteady simulations were conducted for three consecutive pulses as recommended in prior investigations [42], and the outcome was based on the third pulse. The maximum number of iterations per step was fixed at 40 for all simulations, while the time step was 0.005 seconds. Furthermore, in sensitivity testing, no significant alterations were observed while using a time step of 0.002 and maximum iterations of 80.

2.3. Material Properties

The density of the blood was ρ = 1060 kg/m3, and the Carreau non-Newtonian model was employed in this simulation. The Carreau model equation is used to predict the non-Newtonian behavior of blood, and it can be expressed as follows [43, 44]:where the values of A, n, , and are 10.976, −0.3216, 0.0035, and 0.056, respectively.

The AAA wall is assumed to be a nonlinear, isotropic, hyperelastic material with a density of 2200 kg/m3, Young’s modulus E = 2.7 MPa, and Poisson’s ratio υ = 0.45. For hyperelastic materials, the strain-energy function based on the Mooney–Rivlin material model is given as follows:where W is the strain-energy density, B is the left Cauchy–Green tensor, and is the first invariant of B; and α and β are the model parameters indicative of the mechanical properties of the aneurysmal wall. The material parameters α = 17.4 N/cm2 and β = 188.1 N/cm2 are based on the means of the best-fit material parameters of the data examined by Raghavan and Vorp [45]. The bulk modulus () for the Mooney–Rivlin model is determined as follows:

2.4. Mesh Study

A simulation was conducted using four different meshes containing elements of varying sizes, to verify that the outcomes were not affected by the quantity of computational cells. Table 1 provides details on all the meshes used in the study.

The Carreau model was used as the blood non-Newtonian model, and the k-ω model was used as the turbulent model for this steady-state simulation. The velocity and pressure at the time of the pressure peak were used, with an inlet flow rate of 6.4e 5 m3/s and an outlet pressure of 123 mmHg. The mean values of the wall shear stress (WSS) were calculated along the z direction and compared for different meshes in Figure 3. The figure shows that there was no change in the mean WSS even after increasing the number of cells from 263,766 to 349,681. Based on the mesh study, a mesh with a minimum cell length of 0.0008 m, a maximum cell length of 0.0016 m, and a total cell number of 263,766 were selected, as shown in Table 1. Polyhedral and hexagonal elements were used in this research, and the final mesh details are shown in Figure 4.

The artery wall underwent a similar mesh study, employing the von Mises stress parameter. Ultimately, a mesh consisting of 295823 total elements with an element size of 8e − 4 was chosen. Figure 5 depicts the details of this mesh.

2.5. Hemodynamic Parameters

In this study, three hemodynamic parameters that are commonly used in the study of thrombosis and atherosclerosis [46] were used to investigate the mechanics of the artery wall. These parameters are based on wall shear stress (WSS) and include time-averaged WSS (TAWSS), oscillating shear index (OSI), and relative residence time (RRT).

Earlier studies [47, 48] have used the following equation to calculate the TAWSS variable, which is used to determine the total shear stress on the artery wall throughout the entire cardiac cycle.

During a cardiac cycle, endothelial cells are constantly exposed to varying shear stresses. The oscillating shear index (OSI) is a measure used to quantify these cyclic fluctuations. It provides insight into the average duration of flow separation and reattachment [49]. Larger OSI values indicate regions where flow direction changes more frequently, whereas smaller values are associated with a lack of WSS cyclic variations [50]. The OSI value ranges from zero for pure unidirectional flows to 0.5 for entirely oscillatory flows [51]. For pure pulsatile flows, the OSI can be calculated using the following formula [52, 53]:

The distribution of frictional forces on the inner surface of the artery is evaluated by the relative residence time (RRT) parameter. An area with high RRT is subjected to both low and oscillating wall shear stresses simultaneously [54]. RRT is a measure of how closely the particles remain in proximity to the wall and can be calculated using the following formula [55]:

3. Results and Discussion

3.1. Deformation

Geometrical changes in the aneurysmal wall are significant and they cause considerable changes in the mechanics of the wall, which could lead to rupture. During one cardiac cycle, blood flows through the abdominal aorta, causing deformation of the aneurysmal aorta. The deformation of the aorta during a cardiac cycle is a complex process that is influenced by the pulsatile flow of blood, the elasticity of the arterial wall, and the geometry of the aorta. Understanding the deformation of the aneurysmal aorta during one cycle is important for assessing the risk of rupture and designing effective treatment strategies. In this research, the structural parameters, including deformation and von Mises stress, are investigated by presenting the contours of these parameters at the systolic peak of the pressure wave during the cardiac cycle.

The deformation contours for both linear elastic and hyperelastic models are shown in Figure 6. These contours correspond to the peak systolic pressure time. Due to the increase in pressure resulting from the increase in the cross-sectional area of the flow, the deformation level is higher in the area affected by the aneurysm. In addition, the difference in the prediction of deformation level by these two models is shown in Figure 7 for three different views. As can be seen in Figure 7, the hyperelastic model predicts slightly higher deformation levels than the linear elastic model. The highest deformation level is related to a relatively small area in the upper half of the lumen. The greatest difference between the two models is near the other side of the aortic wall and is about 0.6 mm. The difference between the two models in the area with the highest deformation is less than 0.2 mm.

To better understand the difference between the two models in predicting deformation, the average deformation in the z direction is shown in Figure 8(a). In addition, the difference between the mean predicted deformation by the two models is presented in Figure 8(b). As observed, the greatest difference between the two models in predicting the average deformation is about 0.33 mm, which seems negligible compared with the diameter of the aorta in this region.

Understanding the deformation of the aneurysm during one cycle can help identify regions of the aneurysm that are particularly prone to rupture. For example, if certain areas of the aneurysm are subjected to higher levels of stress or strain during one cycle, they may be more likely to rupture in the future [49]. By identifying these regions, clinicians can develop more targeted treatment strategies to reduce the risk of rupture [50]. Figures 6 and 7 reveal that the wall deformation within the lumen of the aorta is greater than that in other regions, resulting in an evident increase in the risk of rupture. Additionally, a dilated section that has a higher intensity of aneurysm exhibits greater deformation, leading to a consequent increase in the risk of rupture. Moreover, the hyperelastic model predicts a higher deformation than the linear elastic model. This issue has also been observed in previous research studies [56]. This not only results in predicting a higher risk of rupture but also leads to differences in predicting other parameters based on wall shear stress, which are used in the analysis of thrombotic areas. The level of deformation is high on both sides of the wall, but due to the smaller area on the right side of the wall, there is likely a higher risk of rupture in this area. To obtain a more precise understanding, it is necessary to analyze the von Mises stress in this region.

3.2. Von Mises Stress

Von Mises stress is an important factor to consider in the evaluation and management of aortic aneurysms. By providing a more comprehensive understanding of the stress levels in the aortic wall, it can help guide treatment decisions and improve outcomes for patients with this condition.

Figure 9 displays the distribution of von Mises stress at various points on the artery wall during the peak pressure. The difference between the two models (i.e., hyperelastic and linear elastic) in terms of von Mises stress is shown in Figure 10. Stress in the aneurysmal region is higher than that in healthy areas, and points with high stress values are more prominent in the upper half of the aneurysm. As expected, the location that experienced the highest deformation also has the highest von Mises stress. Similar to the deformation parameter, hyperelastic models predict slightly higher von Mises stress than linear elastic models. Another notable point is that both models have almost the same prediction for the critical point locations.

The graphs in Figure 11 provide a clearer understanding of the difference between the two models in predicting the von Mises stress parameter. These graphs show the average value of von Mises stress on the aortic wall along the Z-axis. Additionally, the difference between the two models at different points is also shown in Figure 11(b). As can be seen, the difference between the two models in predicting stress values is lower than their difference in predicting deformation values.

Studies have shown that von Mises stress is a useful predictor of the risk of rupture in aortic aneurysms, with significantly higher levels found in ruptured aneurysms than in nonruptured ones [50]. As such, von Mises stress is a critical metric for identifying aneurysms that are at high risk of rupture, and its analysis can help clinicians make informed decisions about the best course of treatment. High von Mises stress levels may suggest that the aneurysm requires more aggressive treatment, such as surgery or endovascular repair, while lower levels may indicate that conservative management is more appropriate [57]. Therefore, von Mises stress analysis can guide treatment decisions for aortic aneurysms. As can be observed in Figures 9 and 10, the von Mises stress in the right wall of the aorta is significantly higher in a small area than at other points. In the deformation analysis section, it was observed that this area has a high deformation despite its small surface area. Therefore, this area can be considered the most susceptible point for aortic rupture. The predicted von Mises stress levels in both models are close to each other, although the hyperelastic model predicts slightly higher von Mises stress levels. Another notable point that has been emphasized in previous studies is the high average von Mises stress in the location where the aneurysm severity is higher.

3.3. Time-Averaged Wall Shear Stress (TAWSS)

Aneurysms occur when the arterial wall weakens and distends, causing a bulge or sac that protrudes from the normal arterial wall. This can alter the blood flow pattern and result in changes in wall shear stress (WSS) values. In particular, WSS in aneurysms tends to be lower than that in healthy arterial segments, indicating that the arterial wall is exposed to lower levels of shear stress [56].

In Figure 12, the TAWSS contours obtained from all three presented models are shown. Additionally, Figure 13 illustrates the differences between these models in predicting TAWSS values. As can be seen, the level of stress is very low in the region affected by the aneurysm. These low shear stress values are due to the sudden increase in cross-sectional area and the decrease in mean velocity in this region. Also, the presence of points with high shear stress at the iliac bifurcation is noteworthy. The reason for the high shear stress in this location may be the increase in pressure in the aneurysm sac, which leads to an increase in flow velocity in the region after it. By examining Figures 12 and 13, it is not possible to clearly determine which model predicts the shear stress in the aneurysm-affected region to be lower or higher. A notable point from this figure is that the linear elastic and hyperelastic models predict almost the same level of stress in the average time.

To further investigate the impact of different models of arterial walls on TAWSS values in different regions, 2D graphs in Figure 14 showing the average values of this parameter along the Z-axis were used. As observed, both linear elastic and hyperelastic models had nearly identical predictions for TAWSS. However, in the region with very low shear stress, the rigid wall model predicted higher shear stress values than the other two models, while in points with high shear stress, the shear stress values calculated by linear elastic and hyperelastic models were higher.

Low wall shear stress (WSS) can lead to thrombus formation, endothelial cell dysfunction and damage, and low oxygen tension. Platelets can adhere to the endothelial cells lining the vessel wall, leading to the formation of a platelet-rich thrombus [58]. Low WSS can cause endothelial cells to produce less nitric oxide and upregulate pro-inflammatory cytokines and adhesion molecules, promoting platelet aggregation and thrombus formation [59]. Additionally, low WSS can create areas of low oxygen tension, resulting in the release of more adenosine diphosphate (ADP) from red blood cells, which is a potent activator of platelets [60]. Therefore, as seen in Figures 12 and 13, the magnitude of shear stress in the region most affected by the aneurysm is much lower than that in other areas, indicating the susceptibility of this area to thrombus formation. The rigid wall model predicts higher shear stress in this area than linear elastic and hyperelastic models, which means a lower risk of thrombus formation in this area. Another notable finding in these results is the good agreement between the linear elastic and hyperelastic models in predicting the risk of thrombus formation, although the linear elastic model has a slightly more optimistic prediction than the hyperelastic models. However, it should be noted that low shear stress alone does not necessarily lead to the identification of a thrombosis area, and it depends on the amount of shear stress oscillations at that point [61].

However, high WSS can have adverse effects, such as causing damage to endothelial cells that line the arterial wall. This is due to the creation of regions with high shear stress along the arterial wall. As a result, inflammation and further weakening of the arterial wall can occur [62]. Additionally, high WSS values can change blood flow patterns, potentially leading to aneurysm growth and rupture. For instance, vortex flows may form within the aneurysm sac, contributing to thrombus formation and further weakening of the arterial wall [63]. Based on what is observable in Figures 12 and 13, high shear stresses are present in two identifiable regions. The first region is located before the lumen, where the difference in shear stress predicted by all three models is not significant. However, in the second region, located at the iliac bifurcation, the solid wall model predicts lower levels of shear stress than the other two models, thus having a more optimistic prediction. Furthermore, the hyperelastic model predicts a higher risk of high shear stress in this region than the linear elastic model.

3.4. Oscillatory Shear Index (OSI)

The oscillatory shear index (OSI) is a hemodynamic parameter that measures the magnitude of the oscillations in shear stress that occur during the cardiac cycle. High OSI values have been associated with the development and progression of aneurysms and can lead to the activation of mechanosensitive signaling pathways within the endothelial cells that line the arterial wall [64].

Figure 15 displays the predicted values of OSI by all three wall models at different locations along the artery. Figure 16 shows the differences among the models in evaluating the magnitude of shear stress fluctuations. Three locations with high OSI values are observed. The first location is before the onset of the aneurysm, which can be justified by the high curvature of the artery and changes in the direction of these fluctuations in shear stress. The second location is at the bottom of the aneurysm sac, where the flow is stagnant and rotational, and points with high OSI values can be predicted. The third location is after the bifurcation of the iliac, where the increase in velocity and turbulence leads to high OSI values. In contrast to TAWSS, the predicted 3D contours of OSI by the rigid wall model differ from the other two models. In the rigid wall model, a larger area is affected by oscillatory shear stress, while in the linear elastic and hyperelastic models, points with high shear stress fluctuations are concentrated at two points. Similar to the previous parameters, the linear elastic and hyperelastic models have almost similar predictions for the OSI parameter, despite a slight difference in the location of maximum shear stress fluctuations in the aneurysm sac.

Figure 17 shows the average value of the OSI parameter along the Z-axis. As can be seen, the average OSI predicted by the solid wall model is higher than the average predicted by the other two models at the location of the aneurysm sac. However, in the location of the iliac bifurcation, this model predicts fewer oscillations in shear stress compared with the other two models. The noticeable difference between the two linear elastic and hyperelastic models is not observed in calculating the average OSI.

High OSI values can cause regions of the arterial wall to experience bidirectional shear stress, leading to the activation of mechanosensitive signaling pathways within the endothelial cells. This activation can cause the production of inflammatory cytokines, growth factors, and enzymes that further weaken the arterial wall and contribute to aneurysm growth and rupture [65]. In particular, high OSI values can activate nuclear factor kappa B (NF-kB) in endothelial cells, leading to the production of pro-inflammatory cytokines and adhesion molecules, which promote inflammation and contribute to the growth and rupture of aneurysms [62, 66]. Furthermore, high OSI values can cause changes in the distribution and organization of the extracellular matrix (ECM) within the arterial wall. The shear stress associated with high OSI values can cause mechanical strain on the arterial wall, leading to the degradation of collagen and elastin fibers, which are responsible for maintaining the structural integrity of the arterial wall. The degradation of these fibers can lead to a loss of structural support and an increased susceptibility to aneurysm formation and rupture [67]. Based on Figures 15 and 16, it can be found that the area most affected by high oscillatory shear indices (OSIs) is at the end of the lumen and before the bifurcation. Considering the previously mentioned decrease in shear stress in this area, it can be concluded that the risk of thrombus formation is high. Other areas also show points with high OSI, which cannot be equated with a high risk of thrombus formation in these areas due to high shear stress. Linear elastic and hyperelastic models show similar results in terms of OSI, similar to shear stress. Interestingly, despite the analysis of shear stress, the rigid model predicts a higher risk of thrombus formation than the other two models based on the analysis of OSI values. Therefore, it cannot be conclusively determined which model is pessimistic or optimistic, and it is necessary to examine other hemodynamic parameters, such as RRT.

3.5. Relative Residence Time (RRT)

In Figure 18, the predicted values of RRT by all three models used for the aneurysmal wall are shown. Unlike OSI, the points with high RRT values are only observed in one region, which is located in the lower part of the sac. The TAWSS values in this region are very low, and the OSI values are very high; therefore, the high RRT values, which have a direct relationship with OSI and an inverse relationship with TAWSS, are justifiable. The rigid model predicts a wider region for high RRT values, while in the linear elastic and hyperelastic models, points with high RRT values are concentrated in two small regions. This difference in the prediction of the RRT parameter by the models used is more apparent in Figure 19. By looking at this figure, it can be concluded that the difference between the linear elastic and hyperelastic models is very small and only related to the maximum RRT region.

From Figures 18 and 19, it is not possible to determine which model predicts a higher or lower level of RRT. To better understand this issue, Figure 20 is useful. By examining Figure 20, which shows the average RRT values along the Z-axis, it can be concluded that in a wide region of the aneurysm location under the influence, the hyperelastic model evaluates the RRT value higher than the other two models, although the predicted values by the linear elastic model are close to those of the hyperelastic model. However, it should be noted that the rigid model has predicted a higher value for this parameter at the maximum RRT location compared to the other two models.

The duration for which particles remain attached to the wall is directly related to the combined effect of OSI and TAWSS, which also has a significant connection to the biological processes of blood clotting [68]. Increased levels of OSI often indicate low TAWSS, and both of these factors serve as indicators of vascular damage [69]. However, it should be noted that regions with high OSI values may not always correspond to regions with low TAWSS values [61]. In studies that examine blood flow patterns, RRT is a reliable measure used to identify areas where there is low wall shear stress and high oscillation [70]. The regions of the aorta that are affected by aneurysmal dilation have notably higher RRT than the nondilated regions. This positive correlation between RRT and aneurysm growth rate suggests that elevated RRT levels may indicate lower wall shear stress and higher oscillating shear stress, both of which can contribute to the development and progression of an aneurysm. [71]. As can be seen in Figures 18 to 20, in all three models, the RRT value is higher at two points in the lower half of the expanded region, and as a result, the probability of thrombus formation is much higher at these points than at other points. In contrast to what was observed in the analysis of TAWSS and high OSI values, a smaller region of the aneurysmal aorta is susceptible to thrombus formation. As expected, linear elastic and hyperelastic models predict similar values for the RRT parameter, but the situation is somewhat different for the rigid model. Generally, elastic models predict higher RRT values in a wide region of the aorta than the rigid model, but it is the rigid model that considers a higher risk of thrombus formation at critical points. Furthermore, the distance between critical points in the rigid model is greater than in the two elastic models. This issue is rooted in the absence of differentiation in this model, which heavily influences other hemodynamic parameters that are dependent on the shear stress, which is itself dependent on the differentiation of the aortic wall. Accurate knowledge of the amount and location of critical points is essential for clinical studies, and therefore, caution must be exercised when using the rigid model.

4. Conclusions

The employment of linear elastic and rigid wall models is prevalent in the numerical simulation of abdominal aortic aneurysms, yet their accuracy remains underexplored. This investigation delineates the distinctions between these conventional models and the more complex, realistic hyperelastic models. Findings indicate that while linear elastic models align closely with hyperelastic models in terms of evaluated results, this concordance diminishes when compared to the rigid model.

Notably, hyperelastic models predict slightly higher strain differentials and assess von Mises stress within the wall more extensively than linear elastic models, suggesting a more conservative estimation regarding rupture risk by the latter, albeit by a marginal difference. In regions subjected to high shear stress, rigid models estimate lower time-averaged wall shear stress (TAWSS) than their counterparts, implying a diminished thrombosis risk assessment under conditions of extreme shear. Conversely, hyperelastic models forecast higher shear stress levels in these areas than linear elastic models, although the variance is minimal. Under the aneurysm, where shear stress is notably low, rigid models predict higher TAWSS values than both linear elastic and hyperelastic models, with linear elastic models showing substantial agreement with hyperelastic outcomes. Despite the rigid model’s higher shear stress oscillations in this region, analysis of the relative residence time (RRT) parameter reveals a lower thrombosis risk prediction across a broad segment of the aortic wall, in comparison with the alternative models. The linear elastic model also suggests a reduced thrombosis risk here, relative to hyperelastic models, but the difference is not markedly significant.

Considering hyperelastic models as a benchmark, the accuracy of linear elastic models in forecasting aortic rupture and thrombosis formation risk appears reasonable. However, caution is advised with the rigid model’s application. The findings substantiate the linear elastic model’s applicability within certain scenarios, simultaneously emphasizing the imperative for hyperelastic models to achieve a deeper, more nuanced comprehension of AAA dynamics. Through methodical model comparison, this investigation contributes substantially to theoretical knowledge and bears profound clinical practice implications, particularly in refining the precision of rupture and thrombosis risk evaluations. Addressing a critical research gap, the study advocates for the progressive embrace of advanced modeling techniques in AAA analysis, spotlighting their importance in both future research and enhanced clinical outcomes.

5. Limitations

This study critically assessed the accuracy of simplified wall models (rigid and linear elastic) versus the hyperelastic model in predicting AAA rupture and thrombosis risk. While simplified models are favored for their computational efficiency, they inherently lack the fidelity to fully represent the aortic wall’s biomechanics. The hyperelastic model, despite its closer approximation to the true mechanical behavior of arterial walls, also presents limitations primarily due to the reliance on literature-based material parameters. These parameters, while informed by existing studies, may not fully encapsulate the individual variability present in human tissues, possibly affecting the model’s precision in predicting rupture and thrombosis risk.

Additionally, the study’s boundary conditions were derived from literature data rather than patient-specific clinical measurements. This approach, while necessary in the absence of accessible patient data, introduces a degree of generalization that may not accurately reflect the unique physiological conditions of individual patients. Such generalizations could limit the model’s applicability for predicting patient-specific outcomes.

Despite these limitations, this study provides valuable insights into the comparative accuracy of simplified and hyperelastic models in predicting key hemodynamic parameters relevant to the assessment of rupture and thrombosis risk in AAAs. By highlighting the trade-offs between model complexity and computational efficiency, it contributes to the ongoing discourse on optimizing AAA modeling approaches for both research and clinical applications.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares that there are no conflicts of interest related to this research.

Acknowledgments

The author thanks Embodi3D website (https://www.embodi3d.com) for providing the CT scan data, which served as the foundation for the geometric preparations. The author also wishes to express thanks to ChatGPT, whose assistance was invaluable in polishing the text to enhance the manuscript’s clarity and readability.