Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
MalcolmWoodruff`TerraPercolatio`
SeepagePlot  |
|
| | ||||
▪
field is "Head", "PorePressure", "Speed", or "Gradient".
▪
For an unconfined solution, "PorePressure" extends above the water table into the unsaturated zone, contouring negative pore pressure (suction) with a tension cut-off at "SuctionCap" metres of head (default 1 m ≈ 9.81 kPa); set "IncludeSuction" -> for the saturated zone only.
▪
The following options can be given:
"PlotPoints" | 60 | ||
"SuctionCap" | 1. | Positive metres of head at which negative pore pressure (suction) above the water table is clipped in the "PorePressure" map. | |
"IncludeSuction" | Whether the "PorePressure" map extends into the unsaturated zone above the water table; Automatic does so whenever sol has a phreatic surface. |
▪
Examples
(4)
Basic Examples
(1)
Plot the pore-pressure field for an unconfined dam solution:
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,"PorePressure","PlotPoints"25]
 | SeepageUnconfined |
 | SeepagePlot |
Out[1]=
Plot the total head field:
In[2]:=
 | SeepagePlot |
Out[2]=
Options
(3)
""