Shallow water motion

For incompressible fluids integration of the governing equations over the fluid depth and subsequent linearization leads to the following equation:

\nabla \cdot (- gH\boldsymbol{I} \cdot \nabla{\eta}) + \ddot{\eta} = 0, (347)

where g is the earth acceleration, H is the fluid depth measured from a reference level, \boldsymbol{I} is the unit tensor and \eta is the fluid height with respect to the reference level. Usually the fluid level at rest is taken as reference level. The derivation of the equation is described in [95]. The following assumptions are made:

Due to the integration process the above equation is two-dimensional, i.e. only the surface of the fluid has to be meshed. By comparison with the heat equation, the correspondence in Table  ([*]) arises. Therefore, shallow water motion can be simulated using the *HEAT TRANSFER procedure.


Table: Correspondence between the heat equation and the shallow water equation.
heat quantity shallow water quantity
T \eta
\boldsymbol{q} \frac{ \partial }{\partial t} (H \overline{ \boldsymbol{ v}})
q_n H \frac{ \partial }{\partial t} ( \overline{ { v}}_n)
\boldsymbol{\kappa} H g \boldsymbol{ I}
\rho h -
\rho c 1

The quantity \overline{ \boldsymbol{ v}} is the average velocity over the depth and \overline{ v}_n is its component orthogonal to the boundary of the domain. Due to the averaging the equations hold for small depths only (shallow water approximation). Notice that the equivalence of the heat conduction coefficient is proportional to the depth, which is a geometric quantity. For a different depth a different conduction coefficient must be defined.

There is no real two-dimensional element in CalculiX. Therefore, the two-dimensional Helmholtz equation has to be simulated by expanding the two-dimensional fluid surface to a three-dimensional layer with constant width and applying the boundary conditions in such a way that no variation occurs over the width.

Notice that, similar to the acoustic equations, the shallow water equations are of the Helmholtz type leading to a hyperbolic system. For instationary applications eigenmodes can be calculated and a modal analysis performed.