Ground flow formula (Water Overlay): Difference between revisions
Jump to navigation
Jump to search
No edit summary |
No edit summary |
||
| Line 1: | Line 1: | ||
Underground flow is calculated differently from surface flow, to account for the slowdown and porousness of the medium. | Underground flow is calculated differently from surface flow, to account for the slowdown and porousness of the medium. | ||
In general, Darcy's law is used. When an [[Hydraulic features (Water Overlay)# | In general, Darcy's law is used. When an [[Hydraulic features (Water Overlay)#Aquifer|aquifer]] is present a variant is applied. | ||
=====Darcy's law===== | =====Darcy's law===== | ||
| Line 51: | Line 51: | ||
* F = Fraction of water which flows between cells | * F = Fraction of water which flows between cells | ||
* KD = The | * KD = The [[Aquifer kd (Water Overlay)|AQUIFER_KD]] attribute of aquifer. | ||
* WSP<sub>source</sub> = The | * WSP<sub>source</sub> = The WATER_STORAGE_PERCENTAGE attribute of the underground terrain type of the origin cell. | ||
* WSP<sub>target</sub> = The | * WSP<sub>target</sub> = The WATER_STORAGE_PERCENTAGE attribute of the underground terrain type of the target cell. | ||
* WL<sub>source</sub> = 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. | * WL<sub>source</sub> = 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. | ||
* WL<sub>target</sub> = 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. | * WL<sub>target</sub> = 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. | ||
Revision as of 12:24, 4 April 2019
Underground flow is calculated differently from surface flow, to account for the slowdown and porousness of the medium.
In general, Darcy's law is used. When an aquifer is present a variant is applied.
Darcy's law
Underground flow between cells is calculated using Darcy's law[1].
- Cd = Δt * Iund * width * A * ( (WLsource - WLtarget) / distance )
Since both width and distance are directly related to the cell size, and the result should be in water height rather than volume, the formula can be rewritten as follows:
- Cd = Δt * Iund * A * (Wsource - Wtarget) / cell
Because the underground may have a different porousness, the maximum amount of water that can flow to another cell has to take into account the relative water storage capacities of the underground.
- Cwsp = (WLsource - WLtarget) * (WSPsource / (WSPsource + WSPtarget) )
The amount of water which flows from the source to the target cell is calculated as follows:
- Δw = max( 0 , min( Cd , Cwsp ) )
Where:
- Δw = The underground flow which takes place.
- Δt = Computational timestep.
- cell = Cell size.
- Cd = The capacity for water flow possible based on the relative water heights.
- Cwsp = The capacity for water flow possible based on the relative water storage percentages.
- 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.
- A = Contact area of the underground cells
- Iund = The GROUND_INFILTRATION_MD of the underground terrain type of the origin cell.
- WSPsource = The WATER_STORAGE_PERCENTAGE of the underground terrain type of the origin cell.
- WSPtarget = The WATER_STORAGE_PERCENTAGE of the underground terrain type of the target cell.
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.
References
- ↑ 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)