Fluid-Structure Interaction (FSI)

Physical Intuition of Fluid-Structure Interaction (FSI)

In Noise, Vibration, and Harshness (NVH) engineering, structures are rarely isolated. They are almost always in contact with a surrounding or enclosed fluid medium—be it the air inside a passenger cabin, the oil in an engine oil pan, the coolant in a water jacket, or the fuel in a tank. Fluid-Structure Interaction (FSI) is the physical phenomenon where the behavior of a structure affects the behavior of an adjacent fluid, and vice versa.

In classical vibration analysis, we often assume that the air surrounding a vibrating panel behaves as a negligible load. However, when the fluid is enclosed (like in a vehicle cabin or acoustic cavity) or when the fluid density is high (like water, oil, or fuel), this assumption breaks down. The interaction manifest in two primary ways:

Fluid Loading (Added Mass and Damping): The vibrating structure must physically displace the adjacent fluid. For heavy fluids like water or oil, this acts as a significant added mass (virtual mass effect) that lowers the structural natural frequencies. In addition, as the structure vibrates, it radiates acoustic energy into the fluid, which acts as radiation damping, drawing energy away from the structure.

Cavity Stiffening and Resonance Coupling: In an enclosed cavity (e.g., a car cabin or acoustic enclosure), the fluid is trapped. When the boundary walls vibrate, they compress the fluid, which acts as an acoustic spring, shifting the system's structural resonances upward. Furthermore, when a structural natural frequency is close to an acoustic cavity natural frequency, the two modes couple strongly, leading to split resonances and high-amplitude noise levels (such as cabin boom).

To design quiet vehicle cabins, quiet engine covers, or efficient acoustic enclosures, we must understand how structural bending waves couple with acoustic pressure waves. This coupling is governed by the boundary conditions at the interface, where the normal acceleration of the structure must equal the normal acceleration of the fluid particles.

Core Concepts

1. Lagrangian Structure vs. Eulerian Fluid Descriptions

To model FSI, we must reconcile two fundamentally different mathematical descriptions of mechanics:

•

Lagrangian Description (Structure): Used for solid structures. We track the motion of individual material points (or finite element nodes) as they deform through space. The coordinates move *with* the material. This is ideal for solids where displacements are relatively small but geometric boundaries are complex.

•

Eulerian Description (Fluid): Used for fluids. Instead of tracking individual fluid molecules, we focus on fixed points in space and observe changes in fluid properties (like pressure $p$, velocity $\vec{v}$, and density $\rho$) as the fluid flows past. This is ideal for acoustics and fluid dynamics where tracking individual particles is computationally prohibitive.

In elasto-acoustics (low-amplitude FSI), the fluid undergoes very small oscillations about its mean state. We can therefore linearize the Eulerian fluid equations, allowing us to couple the Lagrangian structure (in terms of displacement $u$) to the Eulerian fluid (in terms of acoustic pressure $p$) at their shared boundary.

2. The Acoustic Wave Equation and Interface Boundary Conditions

The acoustic fluid is governed by the linearized conservation of mass, conservation of momentum, and the thermodynamic equation of state. Under the assumptions of an inviscid, compressible, and barotropic fluid with small fluctuations, these combine to yield the classical acoustic wave equation in terms of acoustic pressure $p(\vec{x}, t)$:

$$\nabla^2 p - \frac{1}{c_0^2} \frac{\partial^2 p}{\partial t^2} = 0 \quad \text{in } \Omega_f$$

where:

•

$\nabla^2$ is the Laplacian operator (div of grad).

•

$c_0$ is the speed of sound in the fluid, related to the fluid bulk modulus $K_f$ and mean density $\rho_f$ by $c_0 = \sqrt{K_f / \rho_f}$.

•

$\Omega_f$ is the fluid domain.

At the fluid-structure interface boundary $\Gamma$, the structure and fluid must remain in contact. Therefore, the normal acceleration of the structural boundary must match the normal acceleration of the adjacent fluid particles. Applying the linearized momentum equation for the fluid (Euler's equation of motion: $\rho_f \frac{\partial \vec{v}}{\partial t} = -\nabla p$), we get the interface coupling boundary condition:

$$\frac{\partial p}{\partial n} = -\rho_f \ddot{u}_n \quad \text{on } \Gamma$$

where:

•

$n$ is the unit normal vector pointing *outward* from the fluid domain (into the structure).

•

$u_n = \vec{u} \cdot \vec{n}$ is the normal displacement of the structure at the boundary.

•

$\ddot{u}_n$ is the normal structural acceleration.

This boundary condition is the mathematical bridge: the structure's acceleration acts as a source term (or pressure gradient) for the fluid, while the fluid pressure $p$ acts as a distributed normal traction force on the structure.

3. Lagrangian-Eulerian Coupled FEM Formulation

When we discretize the coupled system using the Finite Element Method (FEM), we obtain a set of coupled algebraic equations. Let the structural displacements be represented by the vector $U_s$ and the acoustic pressures at the fluid nodes by the vector $P$. The coupled system of equations of motion is written as:

$$\begin{bmatrix} M_s & 0 \\ \rho_f L & M_f \end{bmatrix} \begin{bmatrix} \ddot{U}_s \\ \ddot{P} \end{bmatrix} + \begin{bmatrix} K_s & -L^T \\ 0 & K_f \end{bmatrix} \begin{bmatrix} U_s \\ P \end{bmatrix} = \begin{bmatrix} F_s \\ 0 \end{bmatrix}$$

where:

•

$M_s$ and $K_s$ are the structural mass and stiffness matrices.

•

$M_f$ and $K_f$ are the acoustic fluid mass and stiffness matrices (often called the pressure-form matrices).

•

$L$ is the fluid-structure coupling matrix (also denoted as $H$ or $Q$), which maps the interface pressures to structural forces and interface structural displacements to fluid volumetric acceleration sources.

•

$F_s$ is the external structural force vector.

Symmetry and Solver Challenges

Notice that the coupled mass and stiffness matrices are asymmetric. The term $-L^T$ appears in the structural equation (representing the force exerted by the fluid pressure on the structure), whereas $\rho_f L$ appears in the fluid equation (representing the acceleration boundary condition).

This asymmetry poses significant numerical challenges:

•

Standard symmetric eigensolvers (like Lanczos) cannot be used directly.

•

Calculating eigenvalues yields complex eigenvalues if damping is present, or requires asymmetric solvers, which are computationally expensive and less stable.

•

To circumvent this, engineers use symmetric reformulations (e.g., using fluid displacement potential $\phi$ where $p = -\rho_f \ddot{\phi}$, or fluid velocity potential $\chi$), or utilize reduced-order modal models (modal coupling).

4. Bending Waves in Plates and Wave Coincidence

A major aspect of FSI is how efficiently a vibrating structure can radiate noise into the surrounding fluid. The coupling strength depends heavily on the relationship between the wavelength of bending waves in the structure and the wavelength of acoustic waves in the fluid.

Bending Waves in Plates

The bending wave speed (phase velocity) $c_b$ of a thin, flat plate is frequency-dependent (dispersive):

$$c_b(f) = \sqrt{\omega} \left( \frac{D}{\rho_s h} \right)^{1/4} = \sqrt{2\pi f} \left( \frac{E h^2}{12\rho_s (1 - \nu^2)} \right)^{1/4}$$

where:

•

$D = \frac{E h^3}{12(1-\nu^2)}$ is the flexural rigidity of the plate.

•

$h$ is the plate thickness, $\rho_s$ is the structural material density.

•

$E$ is Young's modulus, and $\nu$ is Poisson's ratio.

Coincidence and Critical Frequency

•

Acoustic Wave Speed: The speed of acoustic waves in air $c_0$ is constant ($c_0 \approx 343 \text{ m/s}$ at room temperature) and does not vary with frequency (non-dispersive).

•

Coincidence Frequency ($f_c$): The frequency at which the structural bending wave speed equals the acoustic wave speed ($c_b = c_0$). At this frequency, the structural wavelength matches the acoustic wavelength projecting along the plate.

•

Critical Frequency ($f_{crit}$): The lowest coincidence frequency, occurring for grazing acoustic waves (waves traveling parallel to the plate surface). It is given by:

$$f_{crit} = \frac{c_0^2}{2\pi} \sqrt{\frac{\rho_s h}{D}} = \frac{c_0^2}{2\pi h} \sqrt{\frac{12\rho_s(1-\nu^2)}{E}}$$

Radiation Regions

The radiation efficiency $\sigma_{rad}$ is a dimensionless ratio comparing the actual radiated acoustic power $W$ to the power radiated by a piston of the same area $S$ vibrating with the same space-time mean-square velocity $\langle \overline{v}^2 \rangle$:

$$\sigma_{rad} = \frac{W}{\rho_f c_0 S \langle \overline{v}^2 \rangle}$$

The radiation efficiency varies dramatically across three frequency regions:

Frequency Region	Wave Speed Relation	Radiation Efficiency $\sigma_{rad}$	Physical Mechanism
:---	:---	:---	:---
Sub-Coincidence ($f < f_{crit}$)	$c_b < c_0$	$\sigma_{rad} \ll 1$	Acoustic short-circuiting: fluid slips back and forth between adjacent out-of-phase peaks without radiating far-field sound. Radiation occurs only from plate edges and corners.
Coincidence Peak ($f \approx f_{crit}$)	$c_b \approx c_0$	$\sigma_{rad} > 1$ (can peak to 2-5)	Strong phase matching. Bending waves match the acoustic projection angle, resulting in highly efficient radiation.
Super-Coincidence ($f > f_{crit}$)	$c_b > c_0$	$\sigma_{rad} \approx 1$	The plate behaves like an array of pistons. Acoustic waves propagate away at all angles, yielding full radiation efficiency.

High-Quality Schematic Diagrams

Diagram 1: Lagrangian-Eulerian Fluid-Structure Interface Mechanics

This diagram illustrates the interface boundary $\Gamma$ separating the Lagrangian structural domain (where displacement $\vec{u}$ and acceleration $\vec{\ddot{u}}$ are defined) and the Eulerian acoustic fluid domain (where pressure $p$ is solved). It highlights the continuity of normal acceleration and the transfer of pressure forces.

Diagram 2: Panel Radiation Efficiency & Wave Coincidence

This plot displays the structural radiation efficiency $\sigma_{rad}$ as a function of the frequency ratio $f / f_{crit}$ for a typical vehicle panel. It visually details the transition from acoustically short-circuited sub-critical radiation to the coincidence peak, stabilizing at unity in the super-critical zone.