Waterway level change benchmark (Water Module)

Initially, the ground water table is the same as the water level in the waterway. Suddenly the level in the waterway rises with h0 meters. This causes water to infiltrate into the ground, raising the ground water table.

Formulas

When the water level in the water way changes, the water table level will also change over time. The relation between water level in the waterway and the ground water table at distance x from the water way is described by the following formulas : In case the water level is raised, the following formula should hold:

${\frac {h(x)}{h_{0}}}=erfc(x{\sqrt {\frac {S}{4kDt}}})$

In case the water level is lowered, the following formula should hold:

${\frac {h(x)}{h_{0}}}=erf(x{\sqrt {\frac {S}{4kDt}}})$

The formula erf is the complementary error function.

$erf(z)={\frac {2}{\sqrt {\pi }}}\int \limits _{0}^{z}e^{-t^{2}}dt$

$erfc(z)=1-erf(z)$

where:

$h(x)$: change in freatic groundwater table (m)

$h0$: change in water level in the waterway(m)

$x$: distance to the water way (m)

$t$: time since the water level change

$kD$: transmissivity of the aquifer (m2 / day)

$S$: freatic water storage fraction.

Setup

The test setup uses the same parameters used in the research results published in The grid chosen is 112 by 5. The ground water table comparisons will be conducted in the x-direction. The grid cell size is set to 25meters.

The terrain height is set up as followed:

0.0 for x <= 2

3.0 for x == 3

10.0 for x > 3

The initial water level is set to 2, for both surface and the groundwater table, by using a water level area set to 2 meters. The water will be raised to 4 meters, thereby simulating an $h0$ of 2. On x is ?, from y 1 to 3, a unique Inlet is placed, each with:

UPPER_THRESHOLD set to 4.

LOWER_THRESHOLD set to 4.

Inlet Q set to 0, such that is unlimited.

This ensures that the water level is instantly raised to level 4 and that it remains at 4 at all times.

The ground is setup with: Water storage fraction set to 0.25.

An aquifer is added for the region as well, initialized with the Aquifer KD set to 700 m2/day.

A ground water overlay is used for the simulation, with a simulation time T configured using the Weather's rain M parameter, set to [T, 0.0], resulting in a simulation of T seconds without rain. Note that T has to be converted to days to be used in the formula's described above.

References

Water Module

Tutorials

Overview

Theory

Benchmarks