Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
MalcolmWoodruff`TerraPercolatio`
SeepageUnconfined  |
|
| | ||||
▪
spec is an : "Profile" -> counterclockwise polygon of the dam or slope cross-section; "UpstreamLevel" -> Hu and "DownstreamLevel" -> Hd (flow runs left to right); "Layers" -> {<|"Polygon"->..,"kx"->..,"ky"->..|>,...} for multiple materials, else give homogeneous "kx", "ky".
▪
The base (the minimum-y edge) is impermeable. For non-trivial profiles give "UpstreamFace" and "DownstreamFace" explicitly, as polylines ordered from the base upward; the free surface and the faces are assumed single-valued (x-monotone and y-monotone).
▪
Rainfall or surface infiltration is applied as a flux on the phreatic surface with "Rainfall" -> r (mm/day) or "Recharge" -> q (m/s, added to "Rainfall"); it raises the water table and pore pressures, or drains laterally to a boundary for drain-spacing design. The water table is capped at the ground surface, with the excess shed as runoff; "WaterTableDaylights" in the result reports whether that happened.
▪
The following options can be given:
"Name" | |||
"MeshSize" | Maximum finite-element area used at each phreatic-surface iteration; Automatic scales with the saturated-region area. | ||
"Points" | 25 | Number of stations used to represent the phreatic surface. | |
"Relaxation" | 0.65 | Under-relaxation factor applied to the phreatic-surface update at each iteration. | |
"MaxIterations" | 100 | Maximum number of phreatic-surface iterations before giving up. | |
"Tolerance" | Mean phreatic-surface movement below which the iteration is taken as converged; Automatic scales with the head difference. | ||
"Monitor" | True prints the iteration history (mean movement and exit-point elevation) as the free surface converges. | ||
"Rainfall" | 0. | Steady surface infiltration in millimetres per day, applied as a flux across the phreatic surface. | |
"Recharge" | 0. | The same infiltration flux as "Rainfall", given directly in m/s and added to it. | |
"UpstreamNoFlow" | True makes the upstream face a mid-spacing no-flow divide for parallel-drain design; "UpstreamLevel" is then only the initial guess for the water-table apex. |
▪
Returns the same keys as , plus "PhreaticSurface" (list of points), "ExitPoint", "SaturatedRegion", "GroundSurface", "UnsaturatedZone", "UnsaturatedThickness" ({min,max} gap to the ground), "WaterTableDaylights", "Rainfall_mmday", "Converged" and "Iterations".
Examples
(11)
Basic Examples
(1)
Solve for unconfined seepage through a homogeneous earth dam:
In[1]:=
damSpec=Association["Profile"{{0.`,0.`},{52.`,0.`},{28.`,12.`},{24.`,12.`}},"UpstreamLevel"10.`,"DownstreamLevel"0.`,"kx"1.`*^-6,"ky"1.`*^-6];damSol=
[damSpec,"Points"16];{damSol["Q"],damSol["ExitPoint"],damSol["Converged"]}
 | SeepageUnconfined |
Out[1]=
{1.57542×,{43.1792,4.41041},True}
-6
10
Draw the flow net, including the phreatic surface:
In[2]:=
 | FlowNet |
Out[2]=
Options
(10)
""