大動脈弁FSI解析
Theory and Physics
The aortic valve opens and closes with each heartbeat, undergoing approximately 100,000 load cycles per day. The opening and closing behavior of the valve leaflets is determined by the strong interaction between blood flow and valve tissue, so analysis of fluid alone or structure alone is insufficient. This technology is essential for design optimization of prosthetic valves (mechanical valves, bioprosthetic valves, TAVI valves), predicting the progression of valvular disease, and surgical planning.
What is the mechanics like when the valve opens and closes?
When left ventricular pressure exceeds aortic pressure during systole, the valve opens and blood is ejected. During diastole, the valve closes due to the reverse pressure gradient. The leaflets are thin tissue about 0.5mm thick, undergoing repeated large deformations. Vortex flow forms in the Valsalva sinus downstream of the valve, and Leonardo da Vinci already observed that this vortex assists in valve closure.
Governing Equations
What equations are solved in valve FSI?
The fluid side uses the incompressible Navier-Stokes equations. The Reynolds number reaches about $Re \approx 5000$ to $8000$ at peak systole, and turbulent transition should sometimes be considered.
The structural side models the leaflets as a hyperelastic shell. Constitutive laws like the Fung-type or Lee-Sacks type, which consider fiber reinforcement, are used.
Here, $E_{11}$ is the Green-Lagrange strain in the fiber direction (circumferential), and $E_{22}$ is the strain in the orthogonal direction.
How is leaflet contact handled?
The three leaflets contact each other (coaptation) during valve closure. This is the biggest technical challenge in aortic valve FSI. The structural side needs to handle contact mechanics, and the fluid side needs to handle the limit where the gap becomes zero. The advantage of the Immersed Boundary method is that it can handle contact without changing the fluid mesh topology.
The Aortic Valve Endures 100,000 Openings and Closings Per Day
The aortic valve leaflets are tissue only about 0.5mm thick, opening and closing about 100,000 times a day while withstanding a pressure difference of about 120mmHg (16kPa) during systole. The design life is about 30 years. The anisotropic structure, where collagen fibers are arranged in a cross-hatched pattern, supports this incredible durability. To reproduce this fiber structure in FSI theory, a simple isotropic elastic model is insufficient, and hyperelastic models like the Fung-type are necessary. The complexity of the theory is a reflection of the amazing nature of biological tissue.
Physical Meaning of Each Term
- Structure-Thermal Coupling Term: Thermal expansion due to temperature changes induces structural deformation, and deformation affects the temperature field. $\sigma = D(\varepsilon - \alpha \Delta T)$. 【Everyday Example】Railroad tracks expand in summer and the gaps narrow—a typical example of temperature rise→Thermal Expansion→stress generation. Warping of electronic circuit boards after soldering is also due to differences in thermal expansion coefficients between materials. Thermal stress occurs in engine cylinder blocks due to temperature differences between hot and cold parts, potentially leading to cracks.
- Fluid-Structure Interaction (FSI) Term: Bidirectional interaction where fluid pressure/shear forces deform the structure, and structural deformation changes the fluid domain. 【Everyday Example】Suspension bridge cables vibrating in strong wind (Vortex-Induced Vibration)—wind force shakes the structure, the shaken structure changes the airflow, further amplifying vibration. Blood flow in the heart and elastic deformation of blood vessel walls, and aircraft wing flutter (aeroelastic instability) are also typical FSI problems. One-way coupling may suffice in some cases, but bidirectional coupling is essential for large deformations.
- Electromagnetic-Thermal Coupling Term: A feedback loop where Joule heating $Q = J^2/\sigma$ causes temperature rise, and temperature change alters electrical resistance. 【Everyday Example】The nichrome wire in an electric stove heats up (Joule heat) and glows red when current flows—temperature rise changes resistance, altering current distribution. Eddy current heating in IH cooking heaters and increased sag of power lines due to temperature rise are also examples of this coupling.
- Data Transfer Term: Interpolation resolves mesh mismatches between different physical fields. 【Everyday Example】When calculating "feels-like temperature" by combining "temperature data" and "wind data" in weather forecasting, interpolation is needed if observation points differ—in CAE coupled analysis, structural and CFD meshes generally don't match, so data transfer (Interpolation) accuracy at the interface directly affects result reliability.
Assumptions and Applicability Limits
- Weak Coupling Assumption (One-Way Coupling): Valid when one physical field affects the other but the reverse is negligible.
- Cases Requiring Strong Coupling: Large deformations in FSI, cases with strong temperature dependence in electromagnetic-thermal coupling.
- Time Scale Separation: When characteristic times of each physical field differ greatly, efficiency can be improved with subcycling.
- Interface Condition Consistency: Ensure energy/momentum conservation at the coupling interface is satisfied numerically.
- Non-Applicable Cases: Cases where three or more physical fields are strongly coupled simultaneously may require monolithic methods.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Thermal Expansion Coefficient $\alpha$ | 1/K | Steel: ~12×10⁻⁶, Aluminum: ~23×10⁻⁶ |
| Coupled Interface Force | N/m² (Pressure) or N (Concentrated Force) | Check force balance between fluid and structural sides. |
| Data Transfer Error | Dimensionless (%) | Interpolation accuracy depends on mesh density ratio. Below 5% is a guideline. |
Numerical Methods and Implementation
There are three main approaches.
| Method | Advantages | Disadvantages |
|---|---|---|
| ALE-FEM | High interface accuracy | Mesh breakdown during valve closure |
| Immersed Boundary (IB) | Strong for Large Deformation and contact | Interface smearing |
| Immersogeometric | IGA accuracy + IB flexibility | Complex implementation |
Immersogeometric? That's unfamiliar.
A method proposed by Kamensky, Hsu, Bazilevs (2015), which embeds a NURBS-based leaflet model into a fixed fluid mesh. It overcomes the weakness of the IB method (delta function smearing) with interface conditions using the Nitsche method. The Bazilevs lab at UT Austin leads this field.
IBAMR/IBFE Implementation
What specific software is used?
IBAMR (Immersed Boundary Adaptive Mesh Refinement) is the leading open-source example. Developed by Professor Griffith (UNC), it embeds the FEM structure of leaflets into the fluid on a SAMRAI adaptive mesh using the IBFE (Immersed Boundary Finite Element) method.
The fluid solves the Navier-Stokes equations with a penalty-based IB method.
The leaflet elastic force $\mathbf{F}$ is calculated by FEM and spread to the fluid grid via the delta function.
Time Step and CFL Condition
How small a time step is needed?
Valve opening/closing occurs in about 30ms (from systole onset to full opening). With the CFL Condition $\Delta t \leq h/|\mathbf{u}_{max}|$, peak velocity 1.5 m/s, and minimum grid spacing 0.1mm gives $\Delta t \leq 67\mu$s. In practice, calculations use $\Delta t = 10$ to $50\mu$s. Simulating one heartbeat (0.8s) requires 16,000 to 80,000 steps.
The computational cost is enormous.
Yes. Using AMR (Adaptive Mesh Refinement) to refine only near the valve is essential. In IBAMR, the region within a few mm of the leaflet surface uses the finest grid, while areas farther away use coarser grids, reducing computational cost to less than 1/10.
Choosing Between IBM and ALE—Decision Points for Aortic Valve Analysis
The choice between "Immersed Boundary Method (IBM) or ALE method" is always debated in numerical methods for aortic valve FSI. IBM is advantageous for handling large deformation and contact, but boundary accuracy tends to suffer. ALE has high accuracy but requires measures against mesh collapse during complete closure. In practice, some research groups adopt a hybrid strategy of "switching to ALE during opening and IBM during closure." Which one you choose greatly changes code complexity.
Monolithic Method
Solves all physical fields simultaneously as one system of equations. Stable for strong coupling, but implementation is complex and memory consumption is high.
Partitioned Method (Partitioned Iterative Method)
Solves each physical field independently and exchanges data at the interface. Easy to implement and can utilize existing solvers. Suitable for weak coupling.
Interface Data Transfer
Nearest neighbor (simplest but low accuracy), projection (conservative), RBF interpolation (robust to mesh mismatch). Balance between conservation and accuracy is important.
Sub-iteration
Performs sufficient iterations within each coupling step to ensure interface condition consistency. Residual criteria are scaled based on typical values of each physical field.
Aitken Relaxation
Automatically adjusts the relaxation factor for coupling iterations. An adaptive method that prevents divergence from over-relaxation and accelerates convergence.
Stability Conditions
Beware of the added mass effect (in fluid-structure coupling when structural density ≈ fluid density). Apply Robin-type interface conditions or IQN-ILS method if unstable.
Analogy for Aitken Relaxation
Aitken relaxation is like "balancing a seesaw." If one side pushes too hard, the other side flies up, and the reaction causes it to push too hard again—Aitken relaxation automatically adjusts the pushing force to suppress this oscillation. It's an adaptive method that automatically adjusts the next correction amount based on the previous correction when coupling iterations oscillate and fail to converge.
Practical Guide
A typical workflow is as follows.
1. Valve Geometry Model Creation: Create CAD from parametric shape (leaflet height, coaptation height, annulus diameter) or reconstruct patient-specific shape from echocardiography/CT.
2. Aortic Sinus (Valsalva Sinus) Model: Construct fluid domain from annulus to ascending aorta.
3. Mesh Generation: For leaflets (shell or solid) and fluid domain. For IB method, leaflet mesh is independent of fluid mesh.
4. Material Parameter Setting: Hyperelastic constants for leaflets. Identified from biaxial tensile test data.
5. Boundary Conditions: Left ventricular pressure waveform at inlet, aortic pressure waveform or three-element Windkessel at outlet.
6. Calculation Execution: Three or more heartbeats (to remove initial transients).
7. Postprocessing: Effective Orifice Area (EOA), pressure drop, regurgitant volume, leaflet stress.
What is EOA (Effective Orifice Area)?
Abbreviation for Effective Orifice Area, the most important performance indicator for prosthetic valves. Defined by the Gorlin formula.
$Q_{rms}$ is RMS flow rate (mL/s), $\Delta p_{mean}$ is mean pressure drop (mmHg). For TAVI valves, EOA > 1.0 cm² is considered good.
Material Parameter Identification
How are the leaflet material constants determined?
Fit Fung-type or Lee-Sacks models to biaxial tensile test data. For native valves, $c_0 = 2$ to $10$ kPa is typical, and the fiber direction stiffness parameter $c_1$ is typically 5 to 20 times that of the orthogonal direction $c_2$.
| Valve Type | $c_0$ (kPa) | $c_1$ | $c_2$ | Source |
|---|---|---|---|---|
| Native Aortic Valve | 2–10 | 10–50 | 1–5 | Billiar & Sacks (2000) |
| Bovine Pericardium (TAVI Valve) | 5–20 | 30–80 | 5–15 | Varies by product |
| Porcine Valve (Bioprosthetic) | 3–15 | 15–60 | 2–10 | Stella & Sacks (2007) |
How are calcified valves handled?
In aortic stenosis, calcification progresses on the leaflets. A common approach is to extract calcified regions from CT and set the elastic modulus of those parts 10 to 100 times higher. The degree of calcification directly affects paravalvular leakage (PVL) after TAVI valve deployment, so accurate modeling is crucial.
The On-Site Struggle of Reconstructing Valve Geometry from CT Images
In the practical work of aortic valve FSI analysis, "patient-specific geometry reconstruction" takes the most time. Segmenting leaflet shape from CT or MRI images is a difficult task, taking even experienced researchers several days per case. Recently, deep learning-based automatic segmentation has emerged, reducing time to a few hours. Still, it's a world where "80% of analysis quality is determined by geometry quality."
Analogy for the Analysis Flow
Have you ever blown up a balloon? At that moment, a sophisticated fluid-structure interaction is actually occurring. Internal air pressure (fluid) pushes and expands the rubber wall (structure)→the expanded wall changes the internal pressure distribution→the changed pressure further deforms the wall... Repeating this back-and-forth at each calculation step is FSI analysis.
Related Topics
なった
詳しく
報告