{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-input" ] }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%config InlineBackend.rc = {'figure.dpi': 300}\n", "%config InlineBackend.figure_format = 'svg'\n", "import sys\n", "sys.path.append(\"../py_figures\")\n", "import diagrams\n", "from myst_nb import glue\n", "import inter_fig_chap3" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "(content:chap3:recharge_boundary_and_river)=\n", "# Areal recharge between an impermeable boundary and a river\n", "Consider uniform recharge at a rate $N$ on an unconfined aquifer bounded on the left by an impermeable boundary and on the right by a river (see {numref}`Figure %s `). The problem is the same as in {numref}`content:chap1:recharge_wall_river`, but now the transmissivity is a function of the head $h$ in the aquifer rather than approximated as constant. The length of the aquifer is $L$, and the elevation of the aquifer bottom is $z_b$. The boundary conditions are\n", "\n", "$$\n", "\\left. \\frac{\\text{d}h}{\\text{d}x}\\right\\vert_{x=0} = 0\n", "$$(eq_boundary1_3.1)\n", "\n", "$$\n", "h\\vert_{x=L} = h_r\n", "$$(eq_boundary2_3.1)\n", "\n", "where $h_r$ is the head in the river at $x=L$. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n \n \n \n \n 2026-03-09T11:32:14.427805\n image/svg+xml\n \n \n Matplotlib v3.10.6, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "application/papermill.record/text/plain": "
" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_sol3_1" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:14.461667\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig_sol3_1 = diagrams.solution3_1()\n", "glue(\"fig_sol3_1\", fig_sol3_1, display=False)" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{glue:figure} fig_sol3_1\n", "---\n", "name: fig_recharge3.1\n", "---\n", "Areal recharge on an unconfined aquifer between an impermeable boundary and a river.\n", "```" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "     The boundary conditions in terms of heads (Eqs. {eq}`eq_boundary1_3.1` and {eq}`eq_boundary2_3.1`) are converted to boundary conditions in terms of the discharge potential for unconfined flow\n", "\n", "$$\n", "\\left. \\frac{\\text{d}\\Phi}{\\text{d}x}\\right\\vert_{x=0} = 0\n", "$$(eq_boco_discharge1_3.1)\n", "\n", "$$\n", "\\Phi\\vert_{x=L}=\\Phi_r=\\tfrac{1}{2}k(h_r-z_b)^2\n", "$$(eq_boco_discharge2_3.1)\n", "\n", "Application of the boundary conditions to the general solution ({eq}`eq_secdif3.0`) gives\n", "\n", "$$\n", "\\Phi=-\\frac{N}{2}(x^2 - L^2) + \\Phi_r\n", "$$(eq_appboundary3.1)\n", "\n", "The head is obtained from the potential with Eq. ({eq}`eq_head3.0`), which gives, using Eq. ({eq}`eq_boco_discharge2_3.1`)\n", "\n", "$$\n", "h = z_b + \\sqrt{2\\Phi/k} = z_b + \\sqrt{-N(x^2 - L^2)/k + (h_r-z_b)^2}\n", "$$(eq_headpotential3.1)\n", "\n", "The discharge vector $Q_x$ is obtained by differentiation of the potential as\n", "\n", "$$\n", "Q_x = -\\frac{\\text{d}\\Phi}{\\text{d}x}= Nx\n", "$$(eq_difpotential3.1)\n", "\n", "An example is presented below. The solution is compared to solution ({eq}`eq_sol3h1.3`), where the saturated thickness was approximated as $H=10$ m ({numref}`Figure %s `). The specific choice of the approximate $H$ in the example leads to an underestimate of the head in the aquifer at the impermeable boundary. The discharge vector, on the other hand, is the same for both problems (compare Eqs. ({eq}`eq_difpotential3.1`) and ({eq}`eq_sol3Qx1.3`)), because $Q_x$ is independent of the choice of the approximate thickness $H$." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [], "source": [ "# parameters\n", "L = 1000 # aquifer length, m\n", "H = 10 # approximate aquifer thickness used in solution of Section 1.3, m\n", "zb = -5 # aquifer bottom, m\n", "k = 10 # hydraulic conductivity, m/d\n", "n = 0.3 # porosity, -\n", "hr = 4 # specified head at the right boundary, m\n", "N = 0.001 # areal recharge, m/d" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [], "source": [ "# solution\n", "phir = 0.5 * k * (hr - zb) ** 2\n", "x = np.linspace(0, L, 100)\n", "phi = -N / 2 * (x ** 2 - L ** 2) + phir\n", "h = zb + np.sqrt(2 * phi / k)\n", "happrox = -N / (2 * k * H) * (x ** 2 - L ** 2) + hr\n", "Qx = N * x" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-output" ] }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:14.508148\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# basic plot\n", "plt.subplot(121)\n", "plt.plot(x, h, 'C0')\n", "plt.plot(x, happrox, 'C0--')\n", "plt.subplot(122)\n", "plt.plot(x, Qx, 'C1');" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n \n \n \n \n 2026-03-09T11:32:14.563391\n image/svg+xml\n \n \n Matplotlib v3.10.6, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "application/papermill.record/text/plain": "
" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_recharge3_1" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:14.613507\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot\n", "fig = plt.figure(figsize=diagrams.figdef2)\n", "plt.subplot(121)\n", "plt.title('(a)', loc='left')\n", "plt.plot(x, h, 'C0', label='variable $T$')\n", "H = 10\n", "T = k * H\n", "h13 = -N / (2 * T) * (x ** 2 - L ** 2) + hr\n", "plt.plot(x, h13, 'C0--', label='constant $T$\\n $H=10$ m')\n", "plt.plot([0, L], [zb, zb], 'k')\n", "plt.text(500, -4.8, 'aquifer bottom', ha='center')\n", "plt.yticks([-5, 0, 4, 8])\n", "plt.xlabel('$x$ (m)')\n", "plt.ylabel('head (m)')\n", "plt.legend(loc='center left')\n", "plt.grid(True)\n", "plt.subplot(122)\n", "plt.title('(b)', loc='left')\n", "plt.plot(x, Qx, 'C1')\n", "plt.xlabel('$x$ (m)')\n", "plt.ylabel('$Q_x$ (m$^2$/d)')\n", "plt.grid(True)\n", "plt.tight_layout()\n", "\n", "glue(\"fig_recharge3_1\", fig, display=False)" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{glue:figure} fig_recharge3_1\n", "---\n", "name: fig_recharge_un3.1 \n", "---\n", "Areal recharge on an unconfined aquifer between an impermeable boundary and a river; head (left) and $Q_x$ (right).\n", "```" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "flux at x = 0: -0.000 m^2/d\n", "flux next to river at x = L: 1.000 m^2/d\n" ] } ], "source": [ "print(f'flux at x = 0: {-Qx[0]:.3f} m^2/d')\n", "print(f'flux next to river at x = L: {Qx[-1]:.3f} m^2/d')" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "     A plot of the stream function may also be created, but now the top of the aquifer is formed by the water table ({numref}`Figure %s `)." ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [], "source": [ "# solution\n", "psi = np.zeros((2, len(x)))\n", "psi[1] = -Qx\n", "xg = np.zeros_like(psi)\n", "xg[:] = x\n", "zg = np.zeros_like(psi)\n", "zg[0] = zb\n", "zg[1] = h" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-output" ] }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:14.663538\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# basic streamline plot\n", "plt.subplot(111, aspect=25)\n", "plt.contour(xg, zg, psi, 10, colors='C1', linestyles='-')\n", "plt.plot(x, h, 'k');" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n \n \n \n \n 2026-03-09T11:32:14.709161\n image/svg+xml\n \n \n Matplotlib v3.10.6, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "application/papermill.record/text/plain": "
" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_streamlines3_1" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:14.740867\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot\n", "fig = plt.figure(figsize=diagrams.figdef2)\n", "aspect = 20\n", "ax = plt.subplot(111, aspect=aspect)\n", "cs = plt.contour(xg, zg, psi, np.arange(-0.9, 0, 0.1), colors='C1', linestyles='-')\n", "plt.clabel(cs, fmt='%1.1f', fontsize='smaller')\n", "#plt.contour(xg, zg, psi, np.arange(-0.65, 0.3, 0.1), colors='C1', linestyles='-')\n", "plt.plot(x, h, 'C0')\n", "plt.xlabel('$x$ (m)')\n", "plt.ylim(-5, 9)\n", "plt.yticks(np.arange(-5, 8, 2.5))\n", "for xa in np.arange(100, 1000, 100):\n", " diagrams.arrow((xa, 9), (xa, 10), arrow=\"->\")\n", "plt.ylabel(f'$z$ (m) - VE:{aspect}x')\n", "plt.plot([0, L], [9, 9], 'k:')\n", "ax.spines['top'].set_visible(False)\n", "plt.tight_layout()\n", "\n", "glue(\"fig_streamlines3_1\", fig, display=False)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-input" ] }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "
\n", " \n", "
\n", " \n", "
\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
\n", "
\n", " \n", " \n", " \n", " \n", " \n", " \n", "
\n", "
\n", "
\n", "\n", "\n", "\n", "
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "inter_fig_chap3.fig_3_4()" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{figure} ../py_figures/1x1_pixel.png\n", "---\n", "name: fig_streamlines3.1 \n", "---\n", "Streamlines in a vertical cross section for areal recharge on an unconfined aquifer between two rivers with $\\Delta\\psi=0.1$ m$^2$/d and travel time tick marks every 2 years.\n", "```" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "     The travel time is computed for water that reaches the water table at position $x$ and flows to the river at $x=L$. The horizontal component of the velocity vector is obtained as\n", "\n", "$$\n", "v_x = \\frac{Q_x}{n(h-z_b)} = \\frac{\\sqrt{Nk}}{n}\\frac{x}{\\sqrt{2\\Phi_r/N + L^2 -x^2}}\n", "$$(eq_hvelocity3.1)\n", "The travel time can be computed, as before, as\n", "\n", "$$\n", "t_\\text{tr}(x, L) = \\int_{x}^{L} \\frac{\\text{d}x}{v_x}\n", "$$(eq_traveltime13.1)\n", "The integration is conducted analytically and gives (omitting the constant of integration)\n", "\n", "$$\n", "F(x) = \\int \\frac{\\text{d}x}{v_x} = \\frac{n}{\\sqrt{Nk}} \\sqrt{a-x^2} -\n", "\\sqrt{a} \\: \\text{arctanh}\\left(\\frac{\\sqrt{a-x^2}}{\\sqrt{a}}\\right)\n", "$$(eq_integration3.1)\n", "where \n", "\n", "$$\n", "a = 2\\Phi_r/N + L^2\n", "$$(eq_a3.1)\n", "The travel time can now be computed as \n", "\n", "$$\n", "t_\\text{tr}(x, L) = F(L) - F(x)\n", "$$(eq_traveltime23.1)\n", "\n", "The travel time is plotted from point $x$ to the river in {numref}`Figure %s `. " ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [], "source": [ "# solution\n", "def integral(x):\n", " a = 2 * phir / N + L ** 2\n", " return np.sqrt(a - x ** 2) - np.sqrt(a) * \\\n", " np.arctanh(np.sqrt(a - x ** 2) / np.sqrt(a))\n", "\n", "def traveltime(x):\n", " return n / np.sqrt(N * k) * (integral(L) - integral(x))\n", "\n", "x = np.linspace(10, L, 100)\n", "trtime = traveltime(x)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-output" ] }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:33.260116\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# basic travel time plot\n", "plt.plot(x, trtime);" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n \n \n \n \n 2026-03-09T11:32:33.285152\n image/svg+xml\n \n \n Matplotlib v3.10.6, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "application/papermill.record/text/plain": "
" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_traveltime3_1" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2026-03-09T11:32:33.305789\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot\n", "fig = plt.figure(figsize=diagrams.figdef1new)\n", "plt.plot(x, trtime)\n", "plt.xlabel('starting location (m)')\n", "plt.ylabel('travel time to river (d)')\n", "plt.grid()\n", "plt.tight_layout()\n", "\n", "glue(\"fig_traveltime3_1\", fig, display=False)" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{glue:figure} fig_traveltime3_1\n", "---\n", "name: fig_traveltime3.1\n", "---\n", "Travel time from starting location to river at $x=L$.\n", "```" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{exercise}\n", ":label: ex_3_1_value_H\n", "\n", "Determine the value of the effective saturated thickness $H$ such that the head at the impermeable boundary is the same for the solutions with variable saturated thickness and constant saturated thickness. Plot the head vs. $x$ for both solutions.\n", "\n", "
\n", " Solution \n", "\n", "$H = 11.227$ m\n", "\n", "\n", "\n", "
\n", "```" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{exercise}\n", ":label: ex_3_1_traveltime\n", "\n", "Compare the travel time from $x=400$ to the river computed with the analytic solution derived above and numerical integration of the integral (Eq. {eq}`eq_traveltime13.1`) using the `quad` function of `scipy.integrate`.\n", "\n", "
\n", " Solution \n", "\n", "$t_\\text{tr}=3182.66$ d for both\n", "\n", "
\n", "\n", "```" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{exercise}\n", ":label: ex_3_1_traveltime2\n", "\n", "Compute the travel time for water that recharges the aquifer at position $x$ to reach the river at $x=L$ for the case that the saturated thickness is approximated as constant and equal to $H=10$ m, and add the computed travel time to {numref}`Figure %s `. \n", "\n", "
\n", " Solution \n", "\n", "\n", "\n", "
\n", "\n", "```" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "hideCode": false, "hidePrompt": false, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{exercise}\n", ":label: ex_3_1_compare\n", "\n", "Consider the case of unconfined flow between two parallel rivers as in {numref}`content:chap1:areal_recharge_two_rivers`, but now the saturated thickness of the unconfined aquifer is not approximated by an effective thickness $H$, but varies with the head instead. Compute the difference between the solution for unconfined flow with a varying saturated thickness and the solution for unconfined flow with an effective thickness $H$. Use the parameters specified in {numref}`content:chap1:areal_recharge_two_rivers`. Compare the location of the water divide and the maximum head between the rivers for the two solutions.\n", "\n", "
\n", " Solution \n", "\n", "Water divide is at same location, but the head at the water divide is slightly higher for the case with $H=10$ m.\n", "\n", "\n", "\n", "
\n", "\n", "```" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "celltoolbar": "Tags", "hide_code_all_hidden": false, "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.13.7" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autoclose": false, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 4 }