Ground flow formula (Water Overlay): Difference between revisions

From Tygron Support wiki
Jump to navigation Jump to search
mNo edit summary
mNo edit summary
Line 46: Line 46:
*[[Underground model (Water Overlay)#Horizontal flow and aquifers| Underground horizontal flow]]
*[[Underground model (Water Overlay)#Horizontal flow and aquifers| Underground horizontal flow]]


==Related models==
==Related==
* [[Underground model (Water Overlay)|Underground model]]
The following topics are related to this formula.
* [[Substance flow model (Water Overlay)|Substance flow model]]
; Features
: [[Aquifer (Water Overlay)|Aquifer]]
; Formulas
: [[Groundwater level formula (Water Overlay)|Groundwater level formula]]
: [[Underground infiltration formula (Water Overlay)|Underground infiltration formula]]
; Models
: [[Underground model (Water Overlay)|Underground model]]
: [[Infiltration model (Water Overlay)|Infiltration model]]
: [[Substance flow model (Water Overlay)|Substance flow model]]


==References==
==References==

Revision as of 08:21, 19 June 2019

Underground flow is different from surface flow, since it has to account for the slowdown and porousness of the medium. In general, horizontal underground flow is calculated using formulas described in Harbaugh 2005[1][2]. However, when an aquifer is present, the aquifer variant is applied.

Default horizontal flow

Two adjacent cells, where underground water level of cell 1 is larger than cell 2.

The flow between the two cells is calculated as:

Δw = w1 - w2
Bex = max ( B1 , B2 ) - Dground
Ac = Δx * ( wavg - Bex )
C = min(C1 , C2 )
q = Δw * C * Ac / Δx * Δt

where:

wn = The underground water level of cell n.
Bn = The surface height of cell n.
Cn = The ground conductance of cell, defined in GROUND_INFILTRATION_MD of the underground terrain.
Dground = The configured ground bottom distance, defined in GROUND_BOTTOM_DISTANCE_M of the Water Overlay.
Ac = Area of conductance.
Δw = Underground water level difference.
Δt = Computational timestep.
Δx = Size of grid cell.
wavg = Averaged underground water level, based on water levels in underground, WATER_STORAGE_PERCENTAGE and potentially the surface water level, when the underground is filled to the top.

Aquifer formula

When an aquifer is present, its hydraulic diffusivity is used to calculate the water flow.

First, the hydraulic diffusivity dictates the fraction of the water height difference which will flow.

F = 2 * sqrt( KD / WSPsource ) * sqrt( Δt ) * ( 1 / cell )

Based on this fraction, the actual amount of water flow is calculated.

Δw = ( (WLsource/WSPsource) - (WLtarget/WSPsource) ) * F

Where:

  • Δw = The underground flow which takes place.
  • Δt = Computational timestep.
  • cell = Cell size.
  • F = Fraction of water which flows between cells
  • KD = The AQUIFER_KD attribute of aquifer.
  • WSPsource = The WATER_STORAGE_PERCENTAGE attribute of the underground terrain type of the origin cell.
  • WSPtarget = The WATER_STORAGE_PERCENTAGE attribute of the underground terrain type of the target cell.
  • WLsource = The amount of water in the saturated zone of the source cell. The height of the water column if the equivalent amount of water was placed on the surface.
  • WLtarget = The amount of water in the saturated zone of the target cell. The height of the water column if the equivalent amount of water was placed on the surface.

See also

Related

The following topics are related to this formula.

Features
Aquifer
Formulas
Groundwater level formula
Underground infiltration formula
Models
Underground model
Infiltration model
Substance flow model

References

  1. Harbaugh, A.W., 2005, MODFLOW-2005, the U.S. Geological Survey modular ground-water model-the Ground-Water Flow Process: U.S. Geological Survey Techniques and Methods 6-A16, variously paginated.
  2. Langevin, C.D., Hughes, J.D., Banta, E.R., Niswonger, R.G., Panday, Sorab, and Provost, A.M. (2017) ∙ Documentation for the MODFLOW 6 Groundwater Flow Model: U.S. Geological Survey Techniques and Methods, book 6, chap. A55 ∙ p 31 ∙ found at: https://doi.org/10.3133/tm6A55 (last visited 2019-02-04)

Template:WaterOverlay nav