Well near impermeable wall

import matplotlib.pyplot as plt
import numpy as np

import timflow.transient as tft

plt.rcParams["font.size"] = 8.0
plt.rcParams["figure.figsize"] = (8, 3)

ml = tft.ModelMaq(kaq=10, z=[10, 0], Saq=1e-4, tmin=0.01, tmax=10)
w = tft.Well(ml, 0, 20, rw=0.3, tsandQ=[(0, 100)], layers=0)
ml.solve()

self.neq  1
solution complete

t = np.logspace(-2, 1, 100)
h = w.headinside(t)

plt.plot(t, h[0])
plt.grid()

ml.plots.contour(win=[-30, 50, -10, 30], ngr=40, t=1, decimals=2);

ml2 = tft.ModelMaq(kaq=10, z=[10, 0], Saq=1e-4, tmin=0.01, tmax=10)
w2 = tft.Well(ml2, 0, 20, rw=0.3, tsandQ=[(0, 100)])
wall = tft.LeakyWallString(
    ml2, xy=[(-20, 0), (20, 0), (40, 20)], res="imp", order=5, layers=0
)
ml2.solve()

self.neq  13
solution complete

h2 = w2.headinside(t)

plt.plot(t, h[0], label="no wall")
plt.plot(t, h2[0], label="wall")
plt.legend()
plt.grid()

ml2.plots.contour(win=[-30, 50, -10, 30], ngr=40, t=1, decimals=2, parallel=True);