Compute steady drawdown in a confined aquifer - MATLAB Cody - MATLAB Central

Problem 57532. Compute steady drawdown in a confined aquifer

Difficulty:
(187)
Rate
A well extracting water from a confined aquifer will lower the piezometric head and create a cone of depression. In steady state, if the distance r to the point where the drawdown is wanted is smaller than the radius of influence R, then the drawdown s of a well pumping at rate Q0 in a confined aquifer of transmissivity T is
s = (Q0/2piT) ln(R/r)
If the distance r is greater than the radius of influence, then the drawdown is zero. If multiple wells are pumping, the drawdown at the requested point is the sum of the drawdowns from the individual wells.
Boundaries, such as no-flow and constant-head boundaries, can be modeled using the method of images, as described in Cody Problem 57497. Each real well will have a corresponding image, and the drawdown will be the sum of the drawdowns from all wells—real and image. Recall from the previous problem that for no-flow boundaries, the image wells pump in the same sense as the real wells, whereas for constant-head boundaries, the image wells pump in the opposite sense as the real wells.
Write a function to compute steady-state drawdown in a confined aquifer. Input to the function will be the x- and y-coordinates of the points where drawdown is requested, the x- and y-coordinates of the real wells, the pumping rates and radii of influence of the wells, and the transmissivity of the aquifer. If a boundary is present, it will be specified by two pairs of x- and y-coordinates as well as a character string (‘NF’ for no-flow, ‘CH’ for constant-head).
Solve

Solution Stats

28.57% Correct | 71.43% Incorrect
Last Solution submitted on Feb 23, 2025

Problem Comments

Solution Comments

Show comments
Primes and Rough Numbers, Basic ideas
What is a rough number? What can they be used...
John D'Errico in Tips & Tricks
2
4
View Post

Problem Recent Solvers2

Suggested Problems

More from this Author291

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Go to top of page