Linking Aneurysmal Geometry and Hemodynamics Using Computational Fluid Dynamics

This study utilizes a large-scale, patient-specific computational fluid dynamics framework to demonstrate that specific abdominal aortic aneurysm geometric features reliably dictate hemodynamic patterns, suggesting these geometry-driven flow signatures can serve as valuable biomarkers for predicting aneurysm growth and rupture risk.

The Gist

The Big Picture: Why Do Aneurysms Rupture?

Imagine your body's main highway, the aorta, is a giant rubber hose carrying blood from your heart to the rest of you. Sometimes, a weak spot in this hose starts to bulge out, like a balloon. This is called an Abdominal Aortic Aneurysm (AAA).

Doctors have traditionally measured these "bulges" with a ruler. If the balloon gets bigger than a certain size (usually 5.5 cm), they say, "Cut it out before it pops!" But this study suggests that size isn't everything. Just like a balloon can be big but sturdy, or small but fragile, the shape of the bulge and how the blood flows inside it are actually the real clues to whether it will burst.

The Experiment: A Digital Wind Tunnel

The researchers didn't just look at X-rays; they built a digital wind tunnel (called Computational Fluid Dynamics, or CFD) for 74 different patients.

Think of it like this: Instead of testing one real car in a wind tunnel, they built 74 different virtual cars, each with a slightly different shape, and ran them through a supercomputer to see how the air (blood) moved around them. They used a "multiscale" approach, which is like connecting a map of the whole country (the whole body's circulation) to a detailed street view of the specific problem area.

The Key Players: The "Traffic Cops" of Blood Flow

The study looked at how blood pushes against the walls of the artery. They used four main "traffic cops" to measure the stress:

1. TAWSS (Time-Averaged Wall Shear Stress): Imagine the wind blowing against a wall.
   • High Stress: A strong, steady wind that keeps the wall clean and healthy.
   • Low Stress: A dead calm. When the wind stops, dust (blood clots) starts to pile up, and the wall gets weak. The study found that big bulges often have "dead calm" zones where the wall gets sick.
2. OSI (Oscillatory Shear Index): Imagine the wind blowing one way, then the other, then back again, like a pendulum.
   • If the wind is constantly changing direction, it confuses the cells on the artery wall, making them angry and prone to damage.
3. RRT (Relative Residence Time): This measures how long a drop of blood gets "stuck" in a corner.
   • If blood hangs around too long (like a car stuck in a traffic jam), it's more likely to crash (form a clot) and damage the road.
4. LNH (Helicity): This measures if the blood is spinning like a corkscrew or a tornado.
   • A smooth corkscrew spin is actually good! It keeps the blood moving efficiently. But if the shape of the artery is weird, the spin breaks apart, creating chaotic turbulence.

The Big Surprise: It's Not Just the Bulge!

The most shocking discovery in this paper is about where the trouble is happening.

• The Old Belief: Everyone thought the danger was only in the big bulge (the aneurysm sac) right below the kidneys.
• The New Discovery: The researchers found that the iliac arteries (the two branches that split off to go to the legs) are actually the "troublemakers."

The Analogy: Imagine a river that widens into a lake (the aneurysm). You'd think the lake is the dangerous part. But this study found that the narrow streams leading out of the lake (the iliac arteries) are actually where the water gets the most chaotic and stressed. The shape of the bulge sends shockwaves downstream, messing up the flow in the leg arteries even more than the bulge itself.

What They Found (The Results)

1. Bigger isn't always worse, but shape is key: Large bulges created huge "eddies" (swirls) where blood got stuck. This is a recipe for clots and wall weakening.
2. The "Traffic Jam" Effect: In the big bulges, the blood slowed down, stopped, and swirled around. This is bad news for the artery wall.
3. The Iliac Connection: The study showed a very strong link between the shape of the aneurysm and the stress in the leg arteries. If the aneurysm is weirdly shaped, the leg arteries suffer the most. This means doctors need to look at the whole system, not just the bulge.
4. The "Corkscrew" Breakdown: In healthy arteries, blood spins smoothly. In these aneurysms, that smooth spin broke apart into chaotic messes, especially when the heart was relaxing (diastole).

Why Does This Matter?

Currently, doctors decide when to operate based mostly on how wide the balloon is. This paper argues that we need to look at how the blood flows.

If we can measure the "traffic patterns" (the shear stress and swirls) inside a patient's specific artery, we might be able to predict a rupture before the balloon gets huge. It's like knowing a bridge is about to collapse not because it's too wide, but because the wind is hitting it at a weird angle that creates a dangerous vibration.

The Bottom Line

This study is a massive step forward because it used a huge number of patients (74) to prove that geometry dictates destiny. The shape of the aneurysm creates specific flow patterns that either heal the artery or destroy it.

The Takeaway: Don't just measure the size of the balloon. Look at the shape, look at the flow, and don't forget to check the "downstream" pipes (the leg arteries), because that's where the real stress might be hiding.

Technical Summary

1. Problem Statement

Abdominal Aortic Aneurysms (AAAs) are characterized by localized dilation of the abdominal aorta, leading to complex flow patterns and wall-shear-driven mechanobiological stimuli that drive disease progression and rupture.

• Current Limitations: Clinical decision-making for endovascular repair currently relies primarily on the maximum aneurysm diameter. However, literature suggests this metric is insufficient for predicting rupture risk, as aneurysms of similar diameters can exhibit vastly different hemodynamic behaviors.
• Knowledge Gap: While the relationship between geometry and hemodynamics is acknowledged, a quantitative, large-scale analysis linking specific geometric descriptors (curvature, torsion, diameter) to hemodynamic biomarkers (Wall Shear Stress, Oscillatory Shear Index, etc.) across distinct anatomical regions (infrarenal aorta vs. iliac arteries) remains poorly defined.
• Objective: To establish a comprehensive, statistically significant correlation between patient-specific aneurysmal geometry and hemodynamic parameters using a large cohort of 74 AAA models.

2. Methodology

The study employed a multiscale Computational Fluid Dynamics (CFD) framework to simulate blood flow in 74 patient-specific infrarenal AAA models (23 from the SAFE-AORTA project and 51 from the VASCUL-AID repository).

A. Modeling and Simulation Setup

• Governing Equations: Blood was modeled as an incompressible Newtonian fluid ($\rho = 1060 \, \text{kg/m}^3$, $\mu = 0.004 \, \text{Pa}\cdot\text{s}$) governed by the Navier-Stokes and continuity equations.
• Multiscale Boundary Conditions:
  • Inlet: A physiologically consistent flow waveform was generated using a coupled 0D–1D systemic circulation model. This model solved fluid-structure interaction (FSI) equations to derive realistic flow rates based on patient-specific diameters.
  • Outlet: A 3-element Windkessel (RCR) model was applied to the iliac outlets. Parameters were tuned iteratively to match mean arterial pressure (MAP) and distribute flow based on anatomical outlet areas.
• Numerical Solver: Simulations were performed using SimVascular with the svSolver, utilizing the Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov-Galerkin (PSPG) methods for stability.
• Mesh and Convergence: Tetrahedral meshes ($5 \times 10^5$ to $10^6$ elements) were generated using TetGen. Simulations ran for 8 cardiac cycles (8 seconds), with results extracted from the final cycle to ensure periodicity. A time step of $\Delta t = 0.00334$ s was used.
• Assumptions: The aortic wall was treated as rigid (a common approximation for elderly AAA patients where deformation is minimal). Intraluminal thrombus (ILT) was not explicitly modeled due to data limitations.

B. Hemodynamic and Geometric Parameters

• Hemodynamic Indices:
  • TAWSS: Time-Averaged Wall Shear Stress.
  • OSI: Oscillatory Shear Index (measuring directional changes).
  • RRT: Relative Residence Time (combining TAWSS and OSI to measure particle trapping).
  • LNH: Local Normalized Helicity (measuring helical flow intensity).
• Geometric Descriptors: Extracted via the Vascular Modeling Toolkit (VMTK) from the centerline, including Maximum Diameter, Mean/Max Curvature, and Mean/Max Torsion.
• Statistical Analysis: The cohort was split into two anatomical regions: the infrarenal aorta and the iliac arteries. Spearman's rank correlation was used to assess relationships between geometric and hemodynamic variables due to non-normal data distribution.

3. Key Results

A. Flow Field Characteristics

• Systole: High-velocity jet-like flows entered the aneurysm, following the inner curvature. In larger sacs, these jets fragmented, creating intense recirculation zones and vortices.
• Diastole: Flow decelerated significantly, leading to widespread recirculation and stagnation, particularly in the aneurysm sac and distal necks.
• Helicity: Local Normalized Helicity (LNH) showed fragmented helical structures, especially during late diastole, indicating a loss of efficient flow organization in pathological geometries.

B. Hemodynamic Indices Distribution

• Pathogenic Zones: Regions with Low TAWSS (< 4 dyne/cm²)**, **High OSI (> 0.3), and High RRT were consistently identified in the aneurysm sac and necks. These conditions are strongly associated with thrombus formation, endothelial dysfunction, and wall weakening.
• Iliac vs. Infrarenal: While the infrarenal sac showed disturbed flow, the iliac arteries exhibited surprisingly strong correlations between geometry and hemodynamics, often showing more pronounced disturbances than the aorta itself.

C. Statistical Correlations

• Infrarenal Region:
  • Maximum Diameter showed the strongest correlations.
  • Negative Correlation: Max Diameter vs. Mean TAWSS ($\rho = -0.45, p < 0.001$) and Max/Mean OSI.
  • Positive Correlation: Max Diameter vs. Mean RRT ($\rho = 0.25$).
  • Curvature and Torsion showed weak or non-significant correlations with hemodynamics in this region.
• Iliac Region:
  • Stronger Correlations: The relationship between geometry and hemodynamics was significantly stronger here than in the infrarenal region.
  • Maximum Diameter dominated all correlations.
  • RRT: A very strong positive correlation was found between Max Diameter and Mean RRT ($\rho \approx 0.9$).
  • OSI: Unlike the infrarenal region, the correlation between Diameter and OSI in the iliacs was positive.
  • Regression Analysis: The coefficient of determination ($R^2$) for Diameter vs. Mean RRT in the iliacs was 0.8080, compared to only 0.0244 in the infrarenal region, confirming that iliac geometry is a primary driver of downstream hemodynamic disturbances.

4. Key Contributions

1. Large-Scale Cohort: One of the largest systematically analyzed CFD datasets for AAAs to date (74 patient-specific models).
2. Multiscale Framework: Integration of 0D–1D systemic circulation models with 3D stabilized finite-element simulations to generate highly physiologically consistent boundary conditions.
3. Downstream Effect Discovery: The study challenges the traditional focus solely on the infrarenal aorta by demonstrating that iliac arteries are dominant contributors to disturbed hemodynamics, with stronger geometry-hemodynamics correlations than the aneurysm sac itself.
4. Quantitative Biomarkers: Established that specific geometric features (particularly diameter) reliably shape shear-stress patterns, suggesting these could serve as patient-specific biomarkers for risk assessment beyond simple diameter thresholds.

5. Significance and Limitations

• Clinical Impact: The findings suggest that rupture risk assessment models should incorporate downstream geometric descriptors (specifically iliac curvature and diameter) and not just the maximum aortic diameter. The strong correlation in the iliac region implies that AAA morphology induces flow disturbances that propagate significantly downstream.
• Future Directions: The authors propose incorporating Fluid-Structure Interaction (FSI) to account for wall deformation and integrating Intraluminal Thrombus (ILT) modeling, which was omitted in this study. Validation with 4D MRI data is also suggested.
• Limitations: The study assumes rigid walls and Newtonian blood properties. The dataset, while large for CFD, is moderate for generalizing statistical trends, and the lack of ILT data prevents analysis of its specific morphological impact.

Conclusion: The study provides robust evidence that aneurysm shape dictates blood flow behavior, with the iliac arteries playing a critical, underappreciated role in the hemodynamic environment of AAA. These geometry-driven flow signatures offer a pathway for improved, personalized rupture risk prediction.