Parachute FSI
Theory and Physics
Parachute FSI Overview
What kind of physics is involved in parachute deployment simulation?
The parachute canopy is an extremely lightweight, highly deformable membrane structure. The mass ratio $m^* = \rho_s h / (\rho_f D)$ is very small ($m^* \ll 1$), leading to strong fluid-structure interaction. During the deployment process, the canopy inflates from a folded state and eventually generates stable drag.
Governing Equations
What is the structural model for the canopy like?
The canopy is represented by a combination of membrane elements and cable elements (suspension lines). The equation of motion for the membrane is:
Here, $\mathbf{T}$ is the membrane stress tensor, and $\Delta p$ is the pressure difference between inside and outside. The canopy fabric is modeled as a nonlinear orthotropic material.
The fluid side uses the incompressible Navier-Stokes equations. When considering canopy permeability, the flow through the membrane is expressed by Darcy's law.
$k$ is the permeability, and $C_2$ is the inertial resistance coefficient. Permeability significantly affects the drag coefficient and stability.
How are the dynamic loads during deployment handled?
A large instantaneous opening shock load occurs at the initial stage of deployment. The maximum load coefficient $C_x$ depends on the Mach number and Dynamic Pressure $q = \frac{1}{2}\rho V^2$. Predicting this transient load is the core of parachute design.
Parachute "Inflation" — The Theory of the Most Dangerous 0.5 Seconds
When a parachute deploys, the "inflation process" from the packed state to full deployment is the most intense FSI moment of the entire sequence. While the folded canopy captures air and rapidly inflates, the fabric surface experiences an "inflation load" where the dynamic pressure momentarily reaches 3 to 5 times the design load. In U.S. Air Force tests during the 1950s-60s, suspension line failure due to this inflation load was a primary cause of accidents. Theoretically, there is a relationship that "inflation load decreases inversely with the square of the opening time." Slow-openers (parachutes that deploy slowly) use this principle to mitigate shock. In FSI analysis, reproducing this inflation process is the theoretical core of parachute design.
Physical Meaning of Each Term
- Structure-Thermal Coupling Term: Thermal expansion due to temperature change induces structural deformation, and deformation affects the temperature field. $\sigma = D(\varepsilon - \alpha \Delta T)$. 【Everyday Example】Railroad tracks in summer where the gap narrows as the rail expands — temperature rise → Thermal Expansion → a classic example of stress generation. Warping of electronic circuit boards after soldering is also due to differences in thermal expansion coefficients between materials. Engine cylinder blocks develop thermal stress from temperature differences between hot and cold sections, potentially leading to cracks.
- Fluid-Structure Interaction (FSI) Term: Fluid pressure/shear forces deform the structure, and structural deformation changes the fluid domain — a bidirectional interaction. 【Everyday Example】Suspension bridge cables vibrating in strong wind (Vortex-Induced Vibration) — wind force shakes the structure, the shaken structure alters the airflow, further amplifying vibration. Blood flow in the heart and elastic deformation of vessel walls, 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: Joule heating $Q = J^2/\sigma$ causes temperature rise, and temperature change alters electrical resistance — a feedback loop. 【Everyday Example】Nichrome wire in an electric heater heats up (Joule heat) and glows red when current flows — temperature rise changes resistance, altering current distribution. Eddy current heating in IH cooking heaters, increased sag in power lines due to temperature rise are also examples of this coupling.
- Data Transfer Term: Interpolation resolves mesh mismatch between different physical fields. 【Everyday Example】When calculating "feels-like temperature" by combining "air temperature data" and "wind data" in weather forecasting, interpolation is needed if observation points differ — similarly in CAE coupling analysis, structural and CFD meshes generally don't match, so data transfer (Interpolation) accuracy at the interface directly impacts 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, strong temperature dependence in electromagnetic-thermal coupling.
- Time Scale Separation: When characteristic times of each physical field differ greatly, efficiency can be improved via subcycling.
- Interface Condition Consistency: Ensure energy/momentum conservation at the coupling interface is numerically satisfied.
- Non-applicable Cases: When three or more physical fields are strongly coupled simultaneously, monolithic methods may be required.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Thermal Expansion Coefficient $\alpha$ | 1/K | Steel: ~12×10⁻⁶, Aluminum: ~23×10⁻⁶ |
| Coupling Interface Force | N/m² (Pressure) or N (Concentrated Force) | Verify force balance between fluid and structure sides. |
| Data Transfer Error | Dimensionless (%) | Interpolation accuracy depends on mesh density ratio. Below 5% is a guideline. |
Numerical Methods and Implementation
Numerical Methods
What methods are used for parachute FSI analysis?
Since large deformation, folding, and contact of the canopy must be handled, the Immersed Boundary (IB) method and Space-Time FEM are mainstream.
| Method | Fluid | Structure | Features |
|---|---|---|---|
| SSTFSI (Space-Time FSI) | DSD/SST | Membrane/Cable | Developed by Tezduyar lab. For parachutes. |
| IB-FEM | Fixed-grid FVM | Membrane FEM | Strong for large deformations. |
| ALE + remeshing | FVM | FEM | High interface accuracy but deployment process is difficult. |
| Overset CFD + FEM | FVM | FEM (e.g., LS-DYNA) | Can handle complex shapes. |
How is the Space-Time method different from regular FEM?
It's a method that discretizes space and time simultaneously. Even if the mesh deforms due to structural movement, the formulation on the space-time slab avoids mesh compatibility issues. Tezduyar et al. have extensively reported application examples to parachutes in a series of papers.
Contact Handling
How is self-contact of the folded canopy handled?
Numerous contacts between canopy membranes occur during deployment. LS-DYNA's *CONTACT_AUTOMATIC_SINGLE_SURFACE is widely used. Setting the penalty stiffness is crucial; excessive penalty for soft canopy material causes numerical oscillation.
"Porosity" Changes the Calculation — Parachute FSI Methods Considering Air Permeability
Parachute canopy fabric is not a completely impermeable membrane; air passes through slightly (porosity). It has been reported that calculations ignoring this permeability and those considering it can yield a 10-30% difference in parachute drag coefficient. As an FSI method considering porosity, there is a method to set an "equivalent porous boundary condition" on the fabric surface. Specifically, the flow rate per unit area of fabric is expressed as a function of the pressure difference across the surface based on the Ergun equation or Darcy's law, and incorporated into the CFD grid. The difficulty is that parachute fabric porosity changes with deformation during deployment — stretched fabric becomes coarser and more permeable. Creating a coupling model that can accurately represent this "deformation-dependent porosity" is a cutting-edge challenge in parachute FSI methodology.
Monolithic Method
Solves all physical fields simultaneously as one system of equations. Stable for strong coupling but complex to implement and memory-intensive.
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 for 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 coupling iteration relaxation factor. An adaptive method that prevents divergence from over-relaxation and accelerates convergence.
Stability Condition
Beware of the added mass effect (in fluid-structure coupling when structural density ≈ fluid density). If unstable, apply Robin-type interface conditions or the IQN-ILS method.
Analogy for Aitken Relaxation
Aitken relaxation is like "balancing a seesaw." If one side pushes too hard, the other side flies up, and the recoil 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
Model Construction Procedure
What are the steps to start parachute deployment analysis?
1. Convert 2D canopy pattern (gore shape) to 3D initial shape
2. Model suspension lines (truss/beam elements)
3. Set fluid domain (area around canopy ≥5D)
4. Set initial folded state (FEM folding simulation or forced displacement)
5. Deployment simulation (FSI coupling)
6. Evaluate drag coefficient in steady descent state
How are canopy material parameters determined?
Example parameters for typical nylon fabric (MIL-C-7020 Type I) are shown.
| Parameter | Value |
|---|---|
| Areal Density | 40–60 g/m² |
| Young's Modulus (Warp) | 400–600 MPa |
| Young's Modulus (Weft) | 300–500 MPa |
| Poisson's Ratio | 0.1–0.3 |
| Permeability | $10^{-9}$–$10^{-10}$ m² |
How is opening shock load verified?
Compare with wind tunnel test data or drop test measurement data. In NASA's CPAS (Capsule Parachute Assembly System) program for Orion spacecraft parachute design, CFD-FSI results were validated against drop test data. The Drag Coefficient $C_D$ and the peak opening load value are the main validation metrics.
Mars Rover's "Seven Minutes of Terror" — Parachute FSI Supports Space Exploration
NASA Mars rover atmospheric entry is called "Seven Minutes of Terror." The most technically difficult part is the parachute deployment sequence. The Martian atmosphere is only about 1% of Earth's, so at the same speed, dynamic pressure is 1% of Earth's — normal parachutes provide no braking at all. For Curiosity and Perseverance (2021), a DGB (Disk-Gap-Band) type parachute is deployed at supersonic speeds (Mach 1.7). Since Martian atmosphere cannot be replicated in Earth wind tunnels, parachute design relies almost entirely on FSI calculations and mathematical modeling. In the 2014 LDSD (Low-Density Supersonic Decelerator) test, a parachute designed using FSI calculations failed during an actual supersonic test — an event that highlighted the difficulty of design due to the discrepancy between calculation and experiment.
Analogy for Analysis Flow
Have you ever inflated a balloon? At that moment, advanced 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 catchball at each calculation step is FSI analysis.
Common Pitfalls for Beginners
"One-way coupling should be enough, right?" — This misjudgment is the most dangerous in coupling analysis. If structural deformation is微小, one-way may indeed suffice. But in cases like heart valve opening/closing where deformation significantly alters the flow path, one-way coupling is completely inadequate. A rule of thumb is "whether deformation exceeds 1% of the characteristic length." If it does, bidirectional coupling is mandatory. If you settle for one-way, the result can be "plausible but actually completely wrong" — this is the scariest pattern.
Thinking About Boundary Conditions
Data exchange at the coupling interface is like "border control." Each country (physical field) has its own laws (governing equations), but if the exchange of people and goods (force, temperature, displacement) at the border (interface) is not managed accurately, the economies (energy balance) of both countries collapse. Interpolation when meshes don't match is like a "translator" — the smaller the mistranslation (interpolation error), the better the result.
Software Comparison
Tool Comparison
What tools are available for parachute FSI analysis?