Hemodynamic Analysis of Blood Flow in Vessels

Hemodynamic analysis within the cardiovascular system serves as a cornerstone for understanding the complex mechanical interactions between flowing blood and vessel walls. Driven by advancements in medical imaging and computational fluid dynamics (CFD), patient-specific hemodynamic simulations have evolved into critical diagnostic and prognostic tools. They provide deep insights into the pathogenesis of vascular pathologies—such as atherosclerosis development and aneurysm rupture—and assist clinicians in optimizing endovascular interventions.
Accurate numerical simulation relies on a rigorous understanding of the underlying physical and mechanical properties of the cardiovascular environment.

• Rheological Behavior of Blood

  • Under physiological conditions, blood exhibits non-Newtonian characteristics due to the aggregation and deformation of red blood cells. However, in major arteries characterized by high shear rates ($>100\text{ s}^{-1}$), the fluid can be reasonably approximated as a Newtonian fluid with a constant dynamic viscosity of approximately $3.5\text{ mPa}\cdot\text{s}$.
  • The density of human blood is conventionally set at $1060\text{ kg/m}^3$.

• Vessel Wall Mechanics

  • Rigid-wall assumptions are computationally efficient and suitable for short-term steady-state analyses or in anatomically stiff regions like the healthy adult aortic arch.
  • Fluid-Structure Interaction (FSI) models account for arterial compliance, capturing the dynamic coupling between pulsatile pressure waves and vessel expansion. These models are indispensable for high-deformation domains such as cerebral aneurysms and heart valves.

• Flow Regimes and Pulsatility

  • Blood flow alternates between laminar and transitional states, governed by the dimensionless Reynolds number ($\text{Re} = \frac{\rho U D}{\mu}$). In large arteries, $\text{Re}$ typically ranges from 500 to 2000, necessitating advanced turbulence or transition models where appropriate.
  • The cardiac cycle introduces inherent unsteadiness, which is quantitatively characterized by the Womersley number ($\alpha = D\sqrt{\frac{\omega\rho}{\mu}}$).

Governing Equations

The mathematical formulation of blood flow is governed by the principles of continuum mechanics.

1. Conservation of Mass (Continuity Equation)
   $$\nabla \cdot \mathbf{u} = 0$$

2. Conservation of Momentum (Navier-Stokes Equations)
   $$\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u}\cdot\nabla\mathbf{u}\right) = -\nabla p + \mu \nabla^{2}\mathbf{u} + \mathbf{f}$$

3. Constitutive Laws for Vessel Walls

   • Linear elasticity ($\sigma = E \varepsilon$) is often applied for small displacements.
   • Hyperelastic material models (e.g., Ogden or Mooney-Rivlin formulations) are preferred for large deformations typical of soft biological tissues.

Computational Modeling Strategies

Selecting an appropriate modeling approach involves balancing computational cost with physical fidelity.

Standard CFD Workflow (e.g., OpenFOAM)

1. Image Segmentation & Reconstruction: Convert DICOM datasets from CTA or MRA scans into watertight surface meshes (.STL format) using segmentation tools like ITK-SNAP or 3D Slicer.
2. Computational Mesh Generation:
   • Generate volumetric grids using robust algorithms (e.g., snappyHexMesh).
   • Implement boundary layer prism cells to resolve near-wall velocity gradients, ensuring $y^{+} < 1$ for accurate wall shear stress calculations.
3. Boundary Condition Assignment:
   • Inlet: Time-varying physiological velocity profiles or volumetric flow rates.
   • Outlet: Multi-scale Windkessel models (RCR circuits) or pressure-outlet conditions to mimic downstream vascular beds.
   • Walls: No-slip conditions for rigid models, or coupled kinematic interfaces for FSI frameworks.
4. Solver Selection: Utilize transient pressure-based solvers (such as pimpleFoam) paired with suitable closure models.
5. Post-Processing: Extract primary hemodynamic metrics—including Wall Shear Stress (WSS), Oscillatory Shear Index (OSI), and pressure drops—using visualization platforms like ParaView.

Key Hemodynamic Metrics in Clinical Research

Numerical simulations yield quantitative markers that correlate strongly with the initiation and progression of cardiovascular diseases.

• Wall Shear Stress (WSS)

  • Regions exposed to persistently low WSS ($<0.4\text{ Pa}$) are strongly correlated with endothelial dysfunction and the focal accumulation of atherosclerotic plaques.
  • Conversely, abnormally elevated WSS peaks often coincide with structural weakening at aneurysm necks, signaling heightened rupture risk.

• Oscillatory Shear Index (OSI)
  $$\text{OSI} = \frac{1}{2}\left(1-\frac{|\int_{0}^{T}\mathbf{\tau}w dt|}{\int{0}^{T}|\mathbf{\tau}_w| dt}\right)$$

  • OSI quantifies the directional changes of WSS throughout a cardiac cycle. High OSI values denote disturbed, multidirectional flow patterns that promote atherogenesis.

• Pressure Gradients and Fractional Flow Reserve (FFR)

  • Pressure drops across stenotic lesions provide functional assessments of hemodynamic significance, guiding clinical interventions such as percutaneous coronary intervention (PCI).

Conclusion

Hemodynamic analysis bridges clinical medicine and fluid mechanics, offering unprecedented windows into cardiovascular pathophysiology. By integrating patient-specific imaging, sound numerical principles, and robust validation, researchers and clinicians can generate reliable physiological insights. Future developments will likely witness the integration of machine learning for real-time surrogate modeling and multi-scale frameworks that couple localized 3D models with systemic circulatory networks.