Modeling fluid structure interaction
Fluid structure interaction (FSI) computations account for motion or deformation of the solid structure boundaries that are interacting with the adjacent fluid flow. The flow solver handles both transient and static fluid structure interactions.
Transient fluid structure interaction
The flow solver uses the arbitrary-Lagrange-Eulerian (ALE) form to discretize the governing equations for FSI modeling. Integration of the governing Navier-Stokes equations in the ALE form over a moving control volume V(t) [57] is given as:
where:
- U is a conserved quantity (mass, momentum and energy).
- F is the convective flux term.
- G is the viscous flux term.
- ∂V(t) is the control volume boundary.
- S is the control volume boundary surface.
- is the velocity of the control volume boundary.
- is the normal vector to the control volume boundary.
- t is the time.
The discrete form of the governing equations in the ALE form is given as:
where the discrete convective flux in the ALE form is:
and the discrete viscous flux is:
and x(t) are the control volume boundary positions.
The following source term is added to the right hand side of the discrete form of the governing equation, to satisfy the geometric conservation law, by assuming that the state of U is an exact constant solution for the governing equation.
Static fluid structure interaction
The flow solver uses the following discrete form of ordinary differential equation (ODE) to compute FSI modeling with static mesh.
where:
- U is a conserved quantity (mass, momentum and energy).
- is the grid face velocity.
- V is the control volume.
The grid face velocity is defined as:
where is the surface normal.
In static FSI modeling, the velocity of the control volume boundary, , is zero and the control volume depends on the control volume grid coordinates, x(t) as:
Knowing the grid coordinate positions and velocities as a function of time, for a first order time accurate scheme, the flow solver uses the following discretization of the ODE for the static FSI modeling:
where Δt = tn + 1 - tn is the current time step.