Example large diameter well

Experimental and for radial flow only

import matplotlib.pyplot as plt
import numpy as np
from scipy.special import k0

import timflow.steady as tfs

rw = 5
ml = tfs.Model3D(kaq=10, z=np.arange(20, -1, -2), kzoverkh=50)
w = tfs.well.LargeDiameterWell(ml, Qw=100, layers=[0, 1, 2, 3, 4], rw=rw)
rf = tfs.Constant(ml, 100, 0, 20)
ml.solve()

Number of elements, Number of equations: 2 , 11
.
.

solution complete

xg = np.linspace(rw, 50, 100)
h = ml.headalongline(xg, 0)
for i in range(10):
    plt.plot(xg, h[i])
plt.xlim(0, 50)
plt.xticks(np.arange(0, 51, 5))
plt.grid()

qx = ml.disvec(rw, 0)[0]
print(qx * 2 * np.pi * rw)
print(np.sum(qx * 2 * np.pi * rw))

[-1.50634520e+01 -1.54980142e+01 -1.66235366e+01 -1.95088952e+01
 -3.33061020e+01  8.15289261e-13 -1.08995891e-13  1.39950723e-13
 -1.29138331e-12 -1.95538628e-13]
-100.00000000000067

ml.head(rw, 0)

array([19.75128373, 19.75128373, 19.75128373, 19.75128373, 19.75128373,
       19.76475129, 19.77020492, 19.7733766 , 19.77522787, 19.77609269])

k0(5 / ml.aq.lab[1:])

array([2.06829668e-03, 6.74900159e-06, 3.32024764e-08, 2.57897091e-10,
       3.46416156e-12, 8.82738725e-14, 4.63575933e-15, 5.37054606e-16,
       1.44345016e-16])

ml.aq.eigvec

array([[ 1.00000000e-01,  4.41707654e-01,  4.25325404e-01,
         3.98470231e-01,  3.61803399e-01, -3.16227766e-01,
         2.62865556e-01, -2.03030724e-01, -1.38196601e-01,
        -6.99596196e-02],
       [ 1.00000000e-01,  3.98470231e-01,  2.62865556e-01,
         6.99596196e-02, -1.38196601e-01,  3.16227766e-01,
        -4.25325404e-01,  4.41707654e-01,  3.61803399e-01,
         2.03030724e-01],
       [ 1.00000000e-01,  3.16227766e-01,  3.22981015e-16,
        -3.16227766e-01, -4.47213595e-01,  3.16227766e-01,
         1.47542912e-16, -3.16227766e-01, -4.47213595e-01,
        -3.16227766e-01],
       [ 1.00000000e-01,  2.03030724e-01, -2.62865556e-01,
        -4.41707654e-01, -1.38196601e-01, -3.16227766e-01,
         4.25325404e-01, -6.99596196e-02,  3.61803399e-01,
         3.98470231e-01],
       [ 1.00000000e-01,  6.99596196e-02, -4.25325404e-01,
        -2.03030724e-01,  3.61803399e-01, -3.16227766e-01,
        -2.62865556e-01,  3.98470231e-01, -1.38196601e-01,
        -4.41707654e-01],
       [ 1.00000000e-01, -6.99596196e-02, -4.25325404e-01,
         2.03030724e-01,  3.61803399e-01,  3.16227766e-01,
        -2.62865556e-01, -3.98470231e-01, -1.38196601e-01,
         4.41707654e-01],
       [ 1.00000000e-01, -2.03030724e-01, -2.62865556e-01,
         4.41707654e-01, -1.38196601e-01,  3.16227766e-01,
         4.25325404e-01,  6.99596196e-02,  3.61803399e-01,
        -3.98470231e-01],
       [ 1.00000000e-01, -3.16227766e-01, -2.22286669e-16,
         3.16227766e-01, -4.47213595e-01, -3.16227766e-01,
         3.01756509e-16,  3.16227766e-01, -4.47213595e-01,
         3.16227766e-01],
       [ 1.00000000e-01, -3.98470231e-01,  2.62865556e-01,
        -6.99596196e-02, -1.38196601e-01, -3.16227766e-01,
        -4.25325404e-01, -4.41707654e-01,  3.61803399e-01,
        -2.03030724e-01],
       [ 1.00000000e-01, -4.41707654e-01,  4.25325404e-01,
        -3.98470231e-01,  3.61803399e-01,  3.16227766e-01,
         2.62865556e-01,  2.03030724e-01, -1.38196601e-01,
         6.99596196e-02]])