A CFD-Based Investigation of Local Luminal Curvature as a Primary Determinant of Hemodynamic Environments in Cerebral Aneurysms

This study utilizes CFD simulations on 76 patient-specific cerebral aneurysms to demonstrate that local luminal curvature is a primary determinant of hemodynamic environments, where saddle-like patches correlate with high wall shear stress and spherical-like patches with low shear stress, offering a geometric framework for improved risk stratification and intervention planning.

The Gist

Imagine your brain's blood vessels as a complex network of garden hoses. Sometimes, a weak spot in the hose balloons out, forming a bubble called a cerebral aneurysm. If this bubble pops, it causes a devastating stroke. Doctors need to know: Which bubbles are about to pop, and which are safe?

Currently, doctors look at the size and shape of the bubble to guess the risk. But this study suggests there's a much simpler, more reliable way to predict danger: look at the local "curvature" of the bubble's surface.

Here is the breakdown of the research using simple analogies:

1. The Two Shapes of the Bubble

The researchers looked at 76 different aneurysms and realized the surface of these bubbles isn't smooth and uniform. Instead, it's made of two distinct types of terrain, like a landscape:

• The "Saddle" Shape (Hyperbolic): Think of a horse saddle or a Pringles chip. It curves up in one direction and down in another. These shapes usually appear near the neck of the aneurysm (where it connects to the main artery).
• The "Sphere" Shape (Elliptic): Think of a basketball or a dome. It curves the same way in all directions. These shapes usually appear at the top (dome) of the aneurysm.

2. The River of Blood and the "Wind" on the Wall

Blood rushing through these vessels creates friction against the walls, known as Wall Shear Stress (WSS). Think of this like the wind hitting a building. The study found that the shape of the wall dictates how the "wind" hits it, regardless of how fast the river is flowing overall.

• On the Saddle (Neck):

  • The Analogy: Imagine a fast-moving river hitting a curved rock that splits the water. The water rushes past quickly, creating intense, swirling eddies (vortices) right against the rock.
  • The Result: The wall feels a strong, steady, high-speed wind (High Shear Stress). It's like a constant, powerful gust that never changes direction.
  • The Danger: This constant, high-speed "wind" acts like sandpaper, slowly grinding down the wall, making it thin and weak. This is often where the aneurysm starts to grow or rupture.

• On the Sphere (Dome):

  • The Analogy: Imagine the river hitting a round boulder. The water crashes into the front, stops dead (stagnation), and then swirls around in a lazy, chaotic circle before moving on.
  • The Result: The wall feels a weak, flickering wind that changes direction constantly (Low Shear Stress, High Oscillation). It's like a breeze that blows one way, then the other, then stops.
  • The Danger: This chaotic, back-and-forth motion confuses the cells in the vessel wall, causing them to become thick, stiff, and clogged (like plaque), rather than thin and weak.

3. The Big Discovery: Shape is Destiny

The most exciting part of this paper is that the shape of the wall predicts the blood flow behavior better than the size of the aneurysm does.

• If you see a Saddle shape, you can almost guarantee it will have high-speed, grinding blood flow (dangerous for thinning).
• If you see a Sphere shape, you can almost guarantee it will have slow, chaotic, swirling flow (dangerous for thickening/clogging).

This holds true whether the aneurysm has already burst or not, and whether it's on the left or right side of the brain. The local geometry is the "boss" of the local blood flow.

4. Why This Matters for Patients

Currently, doctors have to run complex, expensive computer simulations to guess where an aneurysm might break. This study suggests a shortcut:

Just look at the shape.

If a doctor can map the "saddle" areas (the Pringles-chip shapes) on a standard MRI scan, they can instantly identify the "thin, red, dangerous" spots that are likely to rupture. If they see "sphere" areas, they know those are the "thick, yellow, clogged" spots.

The Takeaway

Think of an aneurysm not just as a balloon, but as a landscape with hills and valleys.

• Valleys (Saddles) get hammered by strong, steady winds that wear the ground down.
• Hills (Spheres) get battered by chaotic, swirling winds that pile up debris.

By simply mapping these shapes, doctors can predict which part of the aneurysm is most likely to fail, allowing for safer, more precise treatments without needing to simulate every single drop of blood. It turns a complex physics problem into a simple geometry lesson.

Technical Summary

1. Problem Statement

Cerebral aneurysms (CAs) are pathological dilations of brain arteries with high risks of rupture, leading to severe morbidity and mortality. Current clinical decision-making regarding intervention versus monitoring relies on a mix of patient demographics and qualitative geometric metrics (e.g., size, location, presence of "blebs"). However, these scoring systems often yield inconclusive results (up to 43% non-definitive cases).

While it is established that hemodynamic forces drive aneurysm initiation and rupture, the specific mechanisms linking local wall morphology to local hemodynamic environments remain poorly characterized. Existing studies often focus on global flow patterns or averaged metrics, failing to account for how specific local surface curvatures (e.g., saddle vs. spherical shapes) dictate wall shear stress (WSS) and oscillatory shear index (OSI). This study aims to bridge that gap by investigating whether local luminal curvature is a primary, independent determinant of hemodynamics, regardless of global flow topology or rupture status.

2. Methodology

Data and Geometry

• Cohort: 76 patient-specific intracranial aneurysm geometries were selected from the Aneurisk repository (anterior circulation).
• Exclusion Criteria: Cases with multiple aneurysms, insufficient parent artery length, or outlets too close to the sac were excluded.
• Preprocessing: An automated clipping procedure standardized the inlet position at the fourth internal carotid artery bend to eliminate operator bias.

Computational Fluid Dynamics (CFD)

• Solver: OpenFOAM® (v23.12).
• Physical Model: Blood was modeled as an incompressible, Newtonian fluid ($\rho = 1056$ kg/m³, $\mu = 3.5 \times 10^{-3}$ Pa·s) in a laminar regime with rigid walls.
• Boundary Conditions:
  • Inlet: Spatially varying pulsatile velocity profiles based on patient age ("young" <50 vs. "older" ≥50) and artery diameter.
  • Outlet: Resistance conditions for pressure and zero-gradient for velocity.
• Mesh & Discretization: Hexahedral meshes with prismatic boundary layers (5 layers). Cell density $\ge 3000$ cells/mm³. Second-order upwind schemes for advection and PISO algorithm for pressure-velocity coupling.
• Simulation: Three cardiac cycles simulated; the final cycle used for analysis.

Geometric and Hemodynamic Analysis

• Curvature Classification: The aneurysm surface was partitioned using Gaussian Curvature ($K$):
  • Saddle-like (Hyperbolic): $K < 0$ (typically at the neck).
  • Spherical-like (Elliptic): $K > 0$ (typically at the dome).
• Hemodynamic Metrics:
  • TAWSS: Time-Averaged Wall Shear Stress.
  • OSI: Oscillatory Shear Index (directional variation).
  • TAWSSG: Time-Averaged WSS Gradient.
  • Flow Topology: Identification of fixed points (saddles, nodes, foci) and flow separatrices using the cycle-averaged WSS vector field.
  • Vorticity: $\lambda_2$-criterion to identify near-wall vortical structures.
  • Divergence: $\text{divWSS}$ to identify flow impingement (positive) vs. separation (negative).

Statistical Analysis

• Surface-averaged metrics were compared across patch types, rupture status, and aneurysm types using t-tests and Wilcoxon signed-rank tests ($\alpha = 0.05$).

3. Key Contributions

1. Decoupling Curvature from Global Flow: The study demonstrates that local hemodynamics are primarily dictated by local surface curvature, independent of global flow patterns (e.g., lateral vs. bifurcation aneurysms) or rupture status.
2. Mechanistic Link to Pathology: It establishes a direct mechanistic link between geometric patch types and specific pathological wall phenotypes (thin/red vs. thick/white).
3. Flow Structure Characterization: It maps specific 3D flow structures (vortices, impingement lines, separation lines) to specific curvature types, explaining why certain hemodynamic environments exist.
4. Clinical Tool Proposal: Proposes curvature-based mapping as an objective, non-invasive tool for risk stratification and endovascular planning, potentially replacing or augmenting qualitative assessments.

4. Key Results

Geometric Distribution

• Saddle-like patches ($K < 0$) occupy 30–45% of the surface area and are predominantly located at the aneurysm neck.
• Spherical-like patches ($K > 0$) are larger (average 56.7% larger area) and dominate the aneurysm dome.
• These distributions showed no significant difference between ruptured and unruptured aneurysms.

Hemodynamic Correlations

• Saddle-like Patches:
  • High TAWSS: Average 55.7% higher than spherical patches.
  • Low OSI: Average 15.0% lower.
  • Mechanism: Driven by intense near-wall vortical activity (negative $\lambda_2$) and flow detachment. The flow sweeps across the neck, creating high shear.
• Spherical-like Patches:
  • Low TAWSS: Lower shear stress.
  • High OSI: Elevated oscillatory shear.
  • Mechanism: Dominated by flow impingement (positive $\text{divWSS}$) and stagnation. Flow hits the dome and spreads, causing directional fluctuations.

Correlation with Pathological Phenotypes

Using established thresholds for "high-flow" (TAWSS > 10 Pa, OSI < 0.001) and "low-flow" (TAWSS < 5 Pa, OSI > 0.01) environments:

• Saddle-like regions are strongly associated with Type-I (Red-Thin) phenotypes, which are thin, reddish, and highly susceptible to rupture.
• Spherical-like regions are associated with Type-II (White-Thick) phenotypes, which are thick, atherosclerotic, and less prone to immediate rupture.
• This correlation held true regardless of whether the aneurysm was ruptured or unruptured.

Robustness

The correlations between curvature and hemodynamics remained statistically significant across different aneurysm types (lateral vs. bifurcation) and rupture statuses, confirming that local geometry overrides global flow variations.

5. Significance and Implications

• Risk Stratification: The study suggests that local curvature is a more reliable predictor of local hemodynamic risk than global geometric indices. Clinicians could use curvature mapping to identify "vulnerable" saddle-like regions prone to thinning and rupture.
• Intervention Planning: Understanding that saddle-like necks experience high shear and spherical domes experience stagnation can guide the placement of flow diverters or coils to optimize hemodynamic modification.
• Pathological Insight: The findings support the "high-flow" pathway for aneurysm growth (driven by saddle-like curvature and vortices) and the "low-flow" pathway for wall thickening (driven by spherical curvature and impingement).
• Future Directions: The authors suggest that curvature-based metrics could be integrated into automated clinical scoring systems to reduce the ambiguity currently present in aneurysm management decisions.

Limitations Noted: The study assumed rigid walls (ignoring fluid-structure interaction) and used arbitrary thresholds for defining flow regimes, though sensitivity analyses confirmed the robustness of the curvature-hemodynamics relationship against these variables.