Transient Stream Drawdown and Depletion
バージョン 1.0.0 (9.29 KB) 作成者:
Bwalya Malama
A collection of 5 scripts for computing transient aquifer and stream drawdown induced by groundwater pumping from a well near a stream.
There are 5 scripts, four of which are functions.
- Prerun_sp.m: Frontend script in which most model parameters are specified as well as the optimization option (using lsqcurvefit) for parameter estimation. This script call the function streamdrawdown(parameters, times);
- streamdrawdown.m: The script is of the function streamdrawdown(params,times) where params are aquifer parameters, namely Kx=Hydraulic conductivity in x-direction, Ss= Aquifer specific storage, Cr = Stream storage coefficient, \beta = streambed conductance, and \kappa = anisotropy ratio, i.e., \kappa = Ky/Kx. In this script, the fixed variables are specified, namely the pumping rate Q, aquifer thickness b, stream channel width w, (x,y) positions of observation location, and R the perpendicular (shortest) distance from the pumping well to the stream bank. The script calls the function dehoog2sDr(t,x,y,tmax,N,params) to perform the Laplace transform inversion. Note: for stream ; pumping well side of stream is , and other side of stream is .
- dehoog2sDr.m Inverts Laplace transform to get the time solution. This can be parallelized.
- inverse_cosFourier.m: performs inversion of the Fourier cosine transform by simple integration relative to the transform variable \xi from 0 to inf to get the solution in y.
- fourierlaplace_sD.m: This is the solution in the Laplace-Fourier cosine space.
引用
Bwalya Malama (2024). Transient Stream Drawdown and Depletion (https://www.mathworks.com/matlabcentral/fileexchange/118155-transient-stream-drawdown-and-depletion), MATLAB Central File Exchange. 取得済み .
MATLAB リリースの互換性
作成:
R2022b
すべてのリリースと互換性あり
プラットフォームの互換性
Windows macOS Linuxタグ
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!バージョン | 公開済み | リリース ノート | |
---|---|---|---|
1.0.0 |