Dam Break in 3D
Problem Description
As shown above, the right part of the tank is filled with water (density: 997 kg/m3, kinematic viscosity: 8.9*10 -7 m2/s), with a height H of the water column. The experimental investigation comprise two different heights H = 300 mm and H = 600 mm, from which the former is simulated here. The setup is exposed to a gravity field with a constant acceleration of g = 9.81 m/s2 in negative z-direction. After the removal of the gate in positive z-direction, the water column collapses.
Numerical Setup
Initially, the particles are distributed on an equidistant Cartesian grid with dx = 5e-3 m. As the air is assumed to have only a very limited influence on the collapse of the water column it is neglected in the conducted simulation. To prevent the water particles from flying away, the geometry of the tank is mirrored to the upper half to close the geometry. From the beginning of the simulation on, the gate (blue) is removed with the approximate velocity in the experimental investigations of 3.46 m/s as stated by Lobovský et al. 1 The gate has to be chosen twice as thick as in the experiment as a representation by at least three particles is required for a wall boundary condition (see Adami et al. 2).
Results
These differences could also be the reason for the discrepancies in the collapsing splash, giving a secondary splash in the simulation in contrast to the experiment. Further investigations including the effect of the air could possibly shed light on these effects.
The rapid increase of the water level at the arrival of the backsplash is reproduced with varying discrepancies depending on the location, but still with reasonable agreement. The most pronounced difference at position H3 is obviously due to differences in the collapsing splash as already seen in Table 1.