shallowFoam

OpenFOAM-based solver for 2D shallow water equations.

This master branch has been tested under foam-extend-3.1, and should work for foam-extend-3.x in general. For OpenFOAM-4.1 see the branch OpenFOAM-4.1, which should work for OpenFOAM-4.x in general.

Authors:

• KM-Turbulenz GmbH (www.km-turbulenz.de), 2009
• Florian Mintgen, 2012

Description:

• Solves the depth-averaged 2D shallow water equations:

  dH/dt + d(HU_i)/d(x_i) = 0

  d(HU_i)/dt + d/dx_j (U_j * HU_i) = - g/2 * dH^2/dx_i - g * H * dz_b/dx_i - tau_bx_i / rho + d^2/dx_j^2 ( nu_t * HU )

  with:

  • H: flow depth
  • HU: specific discharge
  • U: depth averaged velocity
  • z_b: bottom elevation
  • tau_b: bottom stresses

• Bottom stresses are modeled via Strickler-equation

• Turbulence is captured by an eddy viscosity model

• Works in parallel

• Captures wet-dry fronts

• Mesh should have a height of 1 m in z-direction (see tutorials)

• Main advantages over shallowWaterFoam (the shallow water solver in the official OpenFOAM repository):

  • Explicit formulation of flow depth and bottom elevation
  • Bottom stresses / surface roughness taken into account
  • Custom boundary conditions well suited for river hydraulics