Problem 59029. Compute the water table in a layered unconfined aquifer
Write a function to compute the flow per unit width and the elevation of the water table above the bottom of unconfined aquifer with two soil layers. The input to the function will be the positions where the water table elevation (WTE) is requested, the hydraulic conductivities of the two layers, the thickness of the bottom layer, the WTEs at and , and the length L of the aquifer.
Background
Cody Problem 52070 dealt with the effective conductivity of a heterogeneous confined aquifer. In that case, the effective conductivity for the case of flow through soil units in parallel could be computed by realizing that the head difference is same for the two soils and the total flow is the sum of the flows through the soils. The effective conductivity then is the weighted average of the individual conductivities, using the layer thicknesses as weights.
For an unconfined aquifer, the effective conductivity is computed in a similar way, but because the WTE h varies with x, so does :
Conservation of mass says that in steady one-dimensional flow, the flow Q past any section is constant. Then the WTE can be computed by inserting the effective conductivity, integrating, and using the boundary conditions and to determine the flow Q and the constant of integration.
Solution Stats
Problem Comments
Solution Comments
Show commentsProblem Recent Solvers2
Suggested Problems
-
Back to basics 20 - singleton dimensions
268 Solvers
-
286 Solvers
-
Find nearest prime number less than input number
705 Solvers
-
308 Solvers
-
convert matrix to single column
404 Solvers
More from this Author275
Problem Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!