# RD-E: 0601 Fluid Structure Coupling

Sloshing inside a fuel tank by simulating the fluid structure coupling. The tank deformation is achieved by applying an imposed velocity on the left corners. Water and air inside the tank are modeled with the ALE formulation. The tank container is described using a Lagrangian formulation.

A numerical simulation of fluid-structure coupling is performed on sloshing inside a deformable fuel tank. This example uses the ALE (Arbitrary Lagrangian Eulerian) formulation and the hydrodynamic bi-material law (/MAT/LAW37) to model interaction between water, air and the tank container.

## Options and Keywords Used

- Fluid structure coupling simulation, and ALE formulation
- Shell and brick elements
- Hydrodynamic and bi-phase liquid gas material (/MAT/LAW37 (BIPHAS))
- ALE boundary conditions (/ALE/BCS)
- J. Donea Grid Formulation (/ALE/GRID/DONEA)
- Boundary conditions (/BCS)
- Gravity (/GRAV)
- Imposed velocity (/IMPVEL)
- ALE material formulation (/ALE/MAT)

Velocity
(ms^{-1}) |
0 | 5 | 0 | 0 |

Time (ms) | 0 | 12 | 12.01 | 50 |

- Material velocity
- Grid velocity

All nodes, except those on the border have grid (/ALE/BCS) and material (/BCS) velocities fixed in the Z-direction. The nodes on the border only have a material velocity (/BCS) fixed in the Z-direction.

The /ALE/GRID/DONEA option activates the J. Donea grid formulation to compute the grid velocity. See the Radioss Theory Manual for further explanations about this option.

## Input Files

Refer to Access the Model Files to download the required model file(s).

The model files used in this example are available in:

/radioss/example/06_Fuel_tank

## Model Description

A rectangular tank made of steel is partially filled with water, the remainder being supplemented by air. The initial distribution pressure is known and supposed homogeneous. The tank container dimensions are 460 mm x 300 mm x 10 mm, with thickness being at 2 mm.

**Material Properties**- Density
- 0.0078 $$\left[\frac{g}{m{m}^{3}}\right]$$
- Young's modulus
- 210000 $\left[\mathrm{MPa}\right]$
- Poisson's ratio
- 0.29
- Yield stress
- 180 $\left[\mathrm{MPa}\right]$
- Hardening parameter
- 450 $\left[\mathrm{MPa}\right]$
- Hardening exponent
- 0.5

The material air-water bi-phase is described in the hydrodynamic bi-material liquid-gas law (/MAT/LAW37). Material LAW37 is specifically designed to model bi-material liquid gas.

- Viscosity
- $${S}_{ij}=2\rho v{\dot{e}}_{ij}$$
- Liquid EOS
- $${P}_{l}={P}_{0}+{C}_{l}\mu $$
- Gas EOS
- $P{v}^{\gamma}={P}_{0}{v}_{0}^{\gamma}=constant$

The equilibrium is defined by: $$\text{\Delta}{P}_{l}=\text{\Delta}{P}_{g}$$

- $${S}_{ij}$$
- Deviatoric stress tensor
- $${e}_{ij}$$
- Deviatoric strain tensor

**Material Parameters - For Liquid**- Liquid reference density, $${\rho}_{l0}$$
- 0.001 $$\left[\frac{g}{m{m}^{3}}\right]$$
- Liquid bulk modulus, $${C}_{l}$$
- 2089 N/mm
^{2} - Initial mass fraction liquid proportion, $${a}_{l}$$
- 100%
- Shear kinematic viscosity ${\nu}_{l}$, $=\mu /{\rho}_{0}^{l}$
- 0.001 mm
^{2}/ms

**Material Parameters - For Gas**- Gas reference density, ${\rho}_{0}^{g}$
- 1.22x10
^{-6}$$\left[\frac{g}{m{m}^{3}}\right]$$ - Shear kinematic viscosity $${\nu}_{g}$$, $=\mu /{\rho}_{0}^{g}$
- 0.00143 mm
^{2}/ms - Constant perfect gas, $\gamma $
- 1.4
- Initial pressure reference gas, $${P}_{0}$$
- 0.1 N/mm
^{2}

**Properties**- Quadratic bulk viscosity/linear bulk viscosity
- 10
^{-20} - Hourglass bulk coefficient
- 10
^{-5}

### Model Method

Using the ALE formulation, the brick mesh is only deformed by tank deformation the water flowing through the mesh. The Lagrangian shell nodes still coincide with the material points and the elements deform with the material: this is known as a Lagrangian mesh. For the ALE mesh, nodes on the boundaries are fixed in order to remain on the border, while the interior nodes are moved.

## Results

### Curves and Animations

Density | |

Velocity |

Density | |

Velocity |

Density | |

Velocity |