Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
MalcolmWoodruff`TerraPercolatio`
BoundaryFlux  |
|
| | ||||
▪
Integrates -(K.Grad[h]).n, where K = diag(kx,ky) and n is the outward normal, over the boundary segments selected by pred[x,y]. Positive values are outflow, negative are inflow.
▪
Use it to check the discharge across any boundary that is not already reported in "Q", such as a drain, a sub-length of a face, or a check that inflow equals outflow.
Examples
(1)
Basic Examples
(1)
Compute the discharge across the excavation floor:
In[1]:=
gUp=10.`;gExc=7.`;hUp=12.`;hDn=gExc;ytip=4.`;w=0.05`;L=20.`;kSand=0.00001`;excPoly={{-L,0.`},{L,0.`},{L,gExc},{w,gExc},{w,ytip},{-w,ytip},{-w,gUp},{-L,gUp}};excLayers={Association["Polygon"excPoly,"kx"kSand,"ky"kSand]};excHeadBCs=hUp,Function{x,y},y≥gUp-&&x≤-w,hDn,Function{x,y},y≥gExc-&&x≥w;pilePred=Function{x,y},Abs[x]≤w+&&y≥ytip-;outerPred=Function{x,y},y≤||Abs[x]≥L-;excSol=
[excLayers,excHeadBCs,{pilePred,outerPred},"MeshSize"0.4`];
excSol,Function{x,y},y≥gExc-&&x≥w
1
6
10.`
1
6
10.`
1
6
10.`
1
6
10.`
1
6
10.`
1
6
10.`
 | SeepageSolve |
 | BoundaryFlux |
1
3
10.`
Out[1]=
0.0000238332
 |
|
""