Computational Implementation of "Uncoupling electrokinetic flow solutions", published in Mathematical Geosciences

<p>This dataset includes Python and Mathematica scripts used to generate figures, and images used in the Mathematical Geosciences (MG) manuscript &quot;Uncoupling Electrokinetic Flow Solutions&quot; by Kuhlman and Malama (2020).</p> <p>Python scripts implementing eigenvalue uncoupling approach for differential equations governing 1D cylindrically symmetric electrokinetic flow problem (i.e., flow to a pumping well).</p> <ol> <li>mpmath python script (recombine-expint.py)&nbsp;implementing Theis &quot;type curve&quot; solution for an infinite domain&nbsp;(Figures 1-3 in MG manuscript)</li> <li>fipy python script&nbsp;(compare-via-fipy.py) and plotting&nbsp;script (plot_fipy_results.py) showing&nbsp;a finite-volume fully coupled solution for a similar finite domain for comparison against eigenvalue uncoupling approach (Figure 4 in MG manuscript). Also includes two shell scripts for driving python scripts for a variety of inputs.</li> </ol> <p>mathematica script (periodic-1D-steady-state-type-1.nb) for solving the algebra associated with the governing equations and plotting figures for&nbsp;analytical solution of&nbsp;periodically driven 1D solution (i.e., laboratory sinusoidal streaming potential and electroosmosis; Figures 4-9 in MG manuscript).</p> <p>&nbsp;</p>