I'm a hydrogeophysicist interested in hydraulic well/packer testing, parameter estimation and uncertainty quantification, flow and transport in ultra-low permeability rocks (e.g., salt, shale, granite, and deep crystalline basement), and coupled hydrogeophysical processes (i.e., electrokinetics, poroelasticity, and thermohaline convection). I am interested in the numerical aspects of integral transform inversion (Laplace and Hankel transforms) and evaluation of analytical functions.

Kristopher L. Kuhlman Applied Systems Analysis & Research Dept. Sandia National Laboratories Albuquerque, New Mexico xxxxx@sandia.gov (replace "xxxxx" with "klkuhlm")

NB: website not endorsed by my employer

Journal Articles · Sandia Reports · Conference Proceedings · Conference Session Chairing · Software

Citation indexes and project/manuscript repositories:

Google,
arXiv,
ResearcherID,
ORCiD,
Scopus,
ResearchGate,
Publons,
PubMed,
Bitbucket, and
GitHub

- Thomas, M.A., K.L.
**Kuhlman**& A.L. Ward, 2017. Anthropogenic influences on groundwater in the vicinity of a long-lived radioactive waste repository,*Hydrological Processes*, in press. - Wang, Y., C.F. Jove Colón & K.L.
**Kuhlman**, 2016. Nonlinear dynamics and instability of aqueous dissolution of silicate glasses and minerals,*Scientific Reports, 6*, 30256. - Malama, B., K.L.
**Kuhlman**, R. Brauchler & P. Bayer, 2016. Modeling cross-hole slug tests in an unconfined aquifer,*Journal of Hydrology, 540*:784–796. (preprint at arXiv and software). **Kuhlman**, K.L., B. Malama & J.E. Heath, 2015. Multiporosity flow in fractured low-permeability rocks,*Water Resources Research, 51*(2):848–860. (preprint at arXiv and software).- Malama, B. & K.L.
**Kuhlman**, 2015. Unsaturated hydraulic conductivity models based on truncated lognormal pore-size distributions,*Groundwater, 53*(3):498–502. (preprint at arXiv, bitbucket) - Sturchio, N.C., K.L.
**Kuhlman**, R. Yokochi, P.C. Probst, W. Jiang, Z.-T. Lu, P. Mueller & G.-M. Yang, 2014. Krypton-81 in groundwater of the Culebra Dolomite near the Waste Isolation Pilot Plant, New Mexico,*Journal of Contaminant Hydrology, 160*:12–20. (draft) - Malama, B., K.L.
**Kuhlman**& S.C. James, 2013. Core-scale solute-transport model selection using Monte Carlo analysis,*Water Resources Research, 49*(6):3313–3147. (preprint at arXiv, bitbucket) **Kuhlman**, K.L., 2013. Review of inverse Laplace transform algorithms for Laplace-space numerical approaches,*Numerical Algorithms, 63*(2):339–355. (preprint at arXiv, bitbucket, citations)- Mishra, P.K. & K.L.
**Kuhlman**, 2013. Unconfined aquifer flow theory: from Dupuit to present*in*Mishra & Kuhlman (Eds.)*Advances in Hydrogeology*(185–202). New York: Springer. (preprint at arXiv, bitbucket) - Ji, S.-H., Y.-K. Koh, K.L.
**Kuhlman**, M.-Y. Lee & J.-W. Choi, 2013. Influence of pressure change during hydraulic tests on fracture aperture,*Groundwater, 51*(2):298–304. (draft) **Kuhlman**, K.L., 2013. Historic testing relevant to disposal of heat-generating waste in salt,*American Nuclear Society: Radwaste Solutions, 20*(4):22–28.- Mishra, P.K. & K.L.
**Kuhlman**(Eds.), 2013.*Advances in Hydrogeology*, New York: Springer. - Mishra, P.K., V.V. Vesselinov & K.L.
**Kuhlman**, 2012. Saturated-unsaturated flow in a compressible leaky-unconfined aquifer,*Advances in Water Resources, 42*:62–70. (preprint at arXiv) - Malama, B., K.L.
**Kuhlman**, W. Barrash, M. Cardiff & M. Thoma, 2011. Modeling slug tests in unconfined aquifers taking into account water table kinematics, wellbore skin and inertial effects,*Journal of Hydrology, 408*(1–2):113–126. (draft and software) - Bakker, M. & K.L.
**Kuhlman**, 2011. Computational issues and applications of line-elements to model subsurface flow governed by the modified Helmholtz equation,*Advances in Water Resources, 34*(9):1186–1194. (software) **Kuhlman**, K.L. & E. Pardo Igúzquiza, 2010. Universal cokriging of hydraulic heads accounting for boundary conditions,*Journal of Hydrology, 384*(1–2):14–25. (draft and software)- Malama, B., K.L.
**Kuhlman**& A. Revil, 2009. Theory of transient streaming potentials associated with axial-symmetric flow in unconfined aquifers,*Geophysical Journal International, 179*(2):990–1003. (software) **Kuhlman**, K.L. & S.P. Neuman, 2009. Laplace-transform analytic-element method for transient porous-media flow,*Journal of Engineering Mathematics, 64*(2):113–130. (draft and software)- Malama, B., A. Revil & K.L.
**Kuhlman**, 2009. A semi-analytical solution for transient streaming potentials associated with confined aquifer pumping tests,*Geophysical Journal International, 176*(3):1007–1016. (software) - Malama, B., K.L.
**Kuhlman**& W. Barrash, 2008. Semi-analytical solution for flow in a leaky unconfined aquifer toward a partially penetrating pumping well,*Journal of Hydrology, 356*(1–2):234–244. (software) **Kuhlman**, K.L., 2008.*Laplace transform analytic element method for transient groundwater flow simulation*, PhD Dissertation, University of Arizona. (software or official U of A pdf copy)**Kuhlman**, K.L. & A.W. Warrick, 2008. Quasilinear infiltration from an elliptical cavity,*Advances in Water Resources, 31*(8):1057–1065. (draft and software)**Kuhlman**, K.L., A.C. Hinnell, P.K. Mishra & T.-C.J. Yeh, 2008. Basin-scale transmissivity and storativity estimation using hydraulic tomography,*Groundwater, 46*(5):706–715. (draft)- Malama, B., K.L.
**Kuhlman**& W. Barrash, 2007. Semi-analytical solution for flow in leaky unconfined aquifer-aquitard systems,*Journal of Hydrology, 346*(1–2):59–68. (software)

**Kuhlman**, K.L., M.M. Mills & E.N. Matteo, 2017.*Consensus on Intermediate Scale Salt Field Test Design*, SFWD-SFWST-2017–000099, SAND2017–3179R, Albuquerque, NM: Sandia National Laboratories.- Sandia National Laboratories, 2016.
*Deep Borehole Field Test Laboratory and Borehole Testing Strategy*, FCRD-UFRD-2016–000072, SAND2016–9235R, Albuquerque, NM: Sandia National Laboratories. - Hansen, F.D., K.L.
**Kuhlman**& S. Sobolik, 2016.*Considerations of the Differences between Bedded and Domal Salt Pertaining to Disposal of Heat-Generating Nuclear Waste*, FCRD-UFRD-2016–000441, SAND2016–6522R, Albuquerque, NM: Sandia National Laboratories. **Kuhlman**, K.L., P.V. Brady, R.J. MacKinnon, J.E. Heath, C.G. Herrick, R.P. Jensen, M. Rigali, T. Hadgu, S.D. Sevougian, J. Birkholzer, B.M. Freifeld & T. Daley, 2016.*Conceptual Design and Requirements for Characterization and Field Test Boreholes: Deep Borehole Field Test*, FCRD-UFD-2015–000131 Rev 1, SAND2016–5692R, Albuquerque, NM: Sandia National Laboratories.- MacKinnon, R.J. & K.L.
**Kuhlman**, 2016.*A Control Variate Method for Probabilistic Performance Assessment: Improved Estimates for Mean Performance Quantities of Interest*. SAND2016–4679, Albuquerque, NM: Sandia National Laboratories. **Kuhlman**, K.L., 2015.*Deep Borehole Disposal Concept and Field Test*, SAND2015–6055O, American Rock Mechanics Association Newsletter, Fall:6–8(16).**Kuhlman**, K.L., P.V. Brady, R.J. MacKinnon, W.P. Gardner, J.E. Heath, C.G. Herrick, R.P. Jensen, T. Hadgu, S.D. Sevougian, J. Birkholzer, B.M. Freifeld & T. Daley, 2015.*Deep Borehole Field Test: Characterization Borehole Science Objectives*, FCRD-UFD-2015–000131, SAND2015–4424R, Albuquerque, NM: Sandia National Laboratories.**Kuhlman**, K.L., 2015.*Deep Borehole Field Test Site Characterization*, SAND2015–5507PE. Albuquerque, NM: Sandia National Laboratories. (Presentation at NWTRB briefing in Albuquerque)**Kuhlman**, K.L., 2015.*Deep Borehole: from Disposal Concept to Field Test*, SAND2015–3116PE. Albuquerque, NM: Sandia National Laboratories.**Kuhlman**, K.L., 2014.*Seeing into the Earth with Hydrogeophysics: Streaming Potential*, SAND2014–19145PE. Albuquerque, NM: Sandia National Laboratories.**Kuhlman**, K.L., 2014.*Summary Results for Brine Migration Modeling Performed by LANL, LBNL, and SNL for the Used Fuel Disposition Program*, FCRD-UFD-2014–000071, SAND2014–18217R. Albuquerque, NM: US Department of Energy, Used Nuclear Fuel Disposition Campaign. (bitbucket)**Kuhlman**, K.L. & B. Malama, 2014.*Assessment of Contaminated Brine Fate and Transport in MB139 at WIPP*, SAND2014–16153. Albuquerque, NM: Sandia National Laboratories.**Kuhlman**, K.L., 2014.*Technical Basis for Disposal of Spent Nuclear Fuel and High-Level Radioactive Waste in Salt*, SAND2014–1857C. Albuquerque, NM: Sandia National Laboratories. (Presentation at NWTRB March 19, 2014 public meeting; presentation transcript)**Kuhlman**, K.L., 2014.*Compliance Recertification Application 2014 for the Waste Isolation Pilot Plant: Appendix HYDRO-2014 Hydrological Investigations*, DOE/WIPP-14–3503. Carlsbad, NM: US Department of Energy, Carlsbad Field Office. (entire WIPP CRA 2014)**Kuhlman**, K.L., 2014.*Compliance Recertification Application 2014 for the Waste Isolation Pilot Plant: Appendix TFIELD-2014 Transmissivity Fields*, DOE/WIPP-14–3503. Carlsbad, NM: US Department of Energy, Carlsbad Field Office. (entire WIPP CRA 2014)

Appendix TFIELD-2014 Attachment A shows the 100 final calibrated realizations.*(70MB pdf)***Kuhlman**, K.L. & S.D. Sevougian, 2013.*Establishing the Technical Basis for Disposal of Heat-Generating Waste in Salt*, FCRD-UFD-2013–000233, SAND2013–6212P. Albuquerque, NM: US Department of Energy, Used Nuclear Fuel Disposition Campaign.**Kuhlman**, K.L. & B. Malama, 2013.*Brine Flow in Heated Geologic Salt*, SAND2013–1944. Albuquerque, NM: Sandia National Laboratories. (bitbucket)- Hansen, F.D., K.L.
**Kuhlman**, W. Steininger & E. Biurrun, 2013.*Proceedings of 3rd US/German Workshop on Salt Repository Research, Design and Operation*, FCRD-UFD-2013–000100, SAND2013–1231P. Albuquerque, NM: US Department of Energy, Used Nuclear Fuel Disposition Campaign. - Kuhlman, K.L., S. Wagner, D. Kicker, R. Kirkes, C. Herrick, D. Guerin, 2012.
*Review and Evaluation of Salt R&D Data for Disposal of Nuclear Waste in Salt*, FCRD-UFD-2012–000380, SAND2012–8808P. Carlsbad, NM: US Department of Energy, Used Nuclear Fuel Disposition Campaign. **Kuhlman**, K.L., 2010.*Analysis Report for the CRA-2009 PABC Culebra Flow and Transport Calculations*, ERMS-552951, Carlsbad, NM: Sandia National Laboratories.**Kuhlman**, K.L., 2010.*Culebra Water Level Monitoring Network Design*, ERMS-554054, Carlsbad, NM: Sandia National Laboratories.

- Sassani, D.C., P.V. Brady, K.L.
**Kuhlman**, C.F. Jove Colón & C.M. Lopez, 2017.*Fluid-Rock Processes Driving Isolation of Crustal Fluids in Crystalline Basement Systems*, Goldschmidt Conference, Paris, France, August 13–18, 2017. - Stein, E.R., J.M. Frederick, G.E. Hammond K.L.
**Kuhlman**, P.E. Mariner & S.D. Sevougian, 2017.*Modeling Coupled Reactive Flow Processes in Fractured Crystalline Rock*, International High-Level Radioactive Waste Management Meeting, Charlotte, NC, April 9–13, 2017. (presentation and paper) **Kuhlman**, K.L., E.L. Hardin, G.A. Freeze, D.C. Sassani & P.V. Brady, 2017.*Drilling and Testing in the Deep Borehole Field Test*, International High-Level Radioactive Waste Management Meeting, Charlotte, NC, April 9–13, 2017. (presentation and paper)**Kuhlman**, K.L., E.L. Hardin, G.A. Freeze, D.C. Sassani & P.V. Brady, 2016.*Site Characterization for a Deep Borehole Field Test*, SAND2016–12408C, American Geophysical Union, San Francisco, CA, Fall Meeting 2016, H13F-1456. (poster)- Broome, S., A. Cashion, J. Feldman, A.J. Sussman, E. Swanson, J. Wilson, J.E. Heath & K.L.
**Kuhlman**, 2016.*Laboratory Noble Gas Migration Experiments Through Rock*, American Geophysical Union, San Francisco, CA, Fall Meeting 2016, S31A-2693. (poster) **Kuhlman**, K.L., D.C. Sassani, G.A. Freeze, E.L. Hardin & P.V. Brady, 2016.*Site Characterization for the Deep Borehole Field Test*, SAND2016–9387C, Geological Society of America Annual Meeting, September 25-28, 2016. (presentation)**Kuhlman**, K.L., E.L. Hardin, G.A. Freeze, D.C. Sassani & P.V. Brady, 2015.*Site Characterization for a Deep Borehole Field Test*, SAND2015–10663C, American Geophysical Union, San Francisco, CA, Fall Meeting 2015, H13-M01. (presentation)- Hansen, F.D., C. Howard, K.L.
**Kuhlman**& J. Holland, 2015.*Capturing Early Evolution of Salt Openings*. American Rock Mechanics Association Annual Meeting, San Francisco, CA, June 28–July 1, 2015, ARMA15-300. (paper + poster) **Kuhlman**, K.L., B.W. Arnold, P.V. Brady, D.C. Sassani, G.A. Freeze & E.L. Hardin, 2015.*Site Characterization for a Deep Borehole Field Test*. International High-Level Radioactive Waste Management Conference, Charleston, SC, April 12–16. (paper + oral presentation)- Park, H.D., G. Hammond, T. Zeitler & K.L.
**Kuhlman**, 2015.*Development of Three Dimensional Models for WIPP Performance Assessment Using PFLOTRAN*. Waste Management, Phoenix, AZ, March 15–19, Annual Meeting 2015, 15073. (paper + oral) - Freeze, GA., B.W. Arnold, P.V. Brady, D.C. Sassani, K.L.
**Kuhlman**& R.J. MacKinnon, 2015.*Siting Considerations for a Deep Borehole Disposal Facility*. Waste Management, Phoenix, AZ, March 15–19, Annual Meeting 2015, 15403. (paper + oral) **Kuhlman**, K.L., J.E. Heath, W.P. Gardner & D.G. Robinson, 2014.*Multiporosity Flow of Gases in Fractured Shale Formations*. American Geophysical Union, San Francisco, CA, Fall Meeting 2014, H51B-0610. (poster)- Heath, J.E., W.P. Gardner, K.L.
**Kuhlman**, D.G. Robinson & S. Bauer, 2014.*Noble Gas Tracing of Fluid Transport in Shale Reservoirs*. American Geophysical Union, San Francisco, CA, Fall Meeting 2014, V51D-4803. (poster) **Kuhlman**, K.L., B. Malama, J.E. Heath, W.P. Gardner & D.G. Robinson, 2013.*Multi-Porosity Transport of Natural Tracers in a Fractured System*. American Geophysical Union, San Francisco, CA, Fall Meeting 2013, H53A-1405. (poster)- Heath, J.E., W.P. Gardner, K.L.
**Kuhlman**, B. Malama, 2013.*Characterization of Hydraulic Fractures and Reservoir Properties of Shale Using Natural Tracers*. American Geophysical Union, San Francisco, CA, Fall Meeting 2013, MR11A-2213. (poster) - Sturchio, N.C., K.L.
**Kuhlman**, R. Yokochi, P.C. Probst, W. Jiang, Z.-T. Lu, P. Mueller & G.-M. Yang, 2013.*Krypton-81 in Groundwater of the Culebra Dolomite Aquifer Near the Waste Isolation Pilot Plant, New Mexico*. American Geophysical Union, San Francisco, CA, Fall Meeting 2013, H13I-1490. (poster) **Kuhlman**, K.L., 2013.*Historic Testing Relevant to Disposal of Heat-Generating Waste in Salt*. American Nuclear Society, Albuquerque, NM, International High-Level Radioactive Waste Management Conference, 1:72–77. (paper 6915 and oral presentation)- Mishra, P.K. & K.L.
**Kuhlman**, 2012.*Radial flow towards well in leaky unconfined aquifer*. American Geophysical Union, San Francisco, CA, Fall Meeting 2012, H43D-1375. (poster) **Kuhlman**, K.L., B. Malama & S.C. James, 2011.*Parameter Uniqueness And Uncertainty Associated For Multirate Transport Models Applied To Core-Scale Test Data*. American Geophysical Union, San Francisco, CA, Fall Meeting 2011, abstract H31D-1179. (poster)**Kuhlman**, K.L. & K.S. Barnhart, 2011.*Hydrogeology Associated with the Waste Isolation Pilot Plant*. American Nuclear Society, Albuquerque, NM, International High-Level Radioactive Waste Management Conference, 2:1126–1132. (paper 3317 and oral presentation)**Kuhlman**, K.L., 2011.*Hydrogeology and WIPP Compliance*. Waste Management Symposia, Phoenix, AZ, Waste Management Conference 2011, paper 11658. (paper and oral presentation)- Clayton, D.J., R.C. Camphouse, S.C. Dunagan, J.W. Garner, A.E. Ismail, T.B. Kirchner, K.L.
**Kuhlman**, J.L. Long, M.B. Nemer, 2010.*Overview of Performance Assessment for the Waste Isolation Pilot Plant*. 10th International Conference on Probabilistic Safety Assessment & Management, Seattle WA, 4:3401–3411 (paper and oral presentation) **Kuhlman**, K.L. & B. Malama, 2010.*Computing the Electrokinetic Response with Simple Models via Eigenvalue Decomposition*. American Geophysical Union, San Francisco, CA, Fall Meeting 2010, abstract H13D-0992. (poster)- Ferré, T.P.A, M.S. Mirjat, D.A. Rose, M. Tuller & K.L.
**Kuhlman**, 2010.*Optimizing Flood Irrigation to Flush Accumulated Salts Through Zonal Flooding*. American Society of Agronomy/Crop Science Society of America/Soil Science Society of America, Long Beach, CA, joint International Annual 2010 meeting. (poster, and software) - Malama, B., K.L.
**Kuhlman**& A. Revil, 2009.*Aquifer Characterization using Transient Streaming Potentials Generated by Flow During Pumping Tests – New developments*. American Geophysical Union, San Francisco, CA, Fall Meeting 2009, abstract H43C-1046. (poster) - Hillesheim, M.B., K.L.
**Kuhlman**& R.L. Beauheim, 2009.*Significant Development in the Understanding of Culebra Hydrology at the WIPP Site*. Transactions of the American Nuclear Society, Atlanta, GA, 100:122–124. (oral presentation and extended abstract) **Kuhlman**, K.L. & S.P. Neuman, 2008.*Particle Tracking in Transient Groundwater Flow Fields with the Laplace Transform Analytic Element Method*. American Geophysical Union, San Francisco, CA, Fall Meeting 2008, abstract H31F-0949. (poster)**Kuhlman**, K.L. & S.P. Neuman, 2007.*Transient Analytic Element Solutions for Flexible Aquifer Test Analyses*. American Geophysical Union, San Francisco, CA, Fall Meeting 2007, abstract H23I-08. (oral presentation)- Malama, B., K.L.
**Kuhlman**& W. Barrash, 2007.*Extension of Leakage Theory to Unconfined Aquifer Flow*. American Geophysical Union, Acapulco, Mexico, Spring Meeting 2007, abstract H53E-04. (poster) **Kuhlman**, K.L. & S.P. Neuman, 2006.*Recent Advances in Laplace Transform Analytic Element Method (LT-AEM) Theory and Application to Transient Groundwater Flow*. American Geophysical Union, San Francisco, CA, Fall Meeting 2006, abstract H41B-0397. (poster)**Kuhlman**, K.L. & S.P. Neuman, 2006.*Recent Advances in LT-AEM Theory and Application to Transient Groundwater Flow*. Copenhagen, Denmark, Conference on Computational Methods in Water Resources XVI. (paper and oral presentation; download LT-AEM software)

- Perry, F., K.L.
**Kuhlman**, D.C. Sassani & F.A. Caporuscio,*Geoscience and Engineering Research Applied to Mined Geologic Repositories and Deep Borehole Concepts for Disposal of Radioactive Waste*. Geological Society of America Annual Meeting, 25–27 September, 2016, Denver, CO; session advocate and co-chair (oral session T68, Science news article). **Kuhlman**, K.L.,*Percolation and Dihedral Angle Topics in Salt*, Sixth US/German Workshop on Salt Repository Research, Design, and Operations, September 7–9, 2016, Washington, DC; session moderator.- Jackson, R., D.C. Sassani, K.L.
**Kuhlman**& J. Birkholzer,*Deep Boreholes and Wellbores: Monitoring, Modeling, and Mitigation*, American Geophysical Union, San Francisco, CA, Fall Meeting 2015; one oral session (H13M) and one poster session (H11B). - Mao, D., K.L.
**Kuhlman**& P.A. Marker,*Advances in Hydrogeology and Hydrogeophysics: Innovations in Experimental Methods, Data Processing and Modeling*, American Geophysical Union, San Francisco, CA, Fall Meeting 2014; three oral sessions (H43O, H44A, H54B) and one poster session (H51B). - Mishra, P.K., K.L.
**Kuhlman**& B. Malama,*Modern Approaches in Hydrogeology: Conceptual and Numerical Model Advances in Cross-Disciplinary Approaches*, American Geophysical Union, San Francisco, CA, Fall Meeting 2012; one oral session (H41L) and one poster session (H43D). **Kuhlman**, K.L.,*Granular Salt Reconsolidation*, Third US/German Workshop on Salt Repository Research, Design, and Operations, October 8–10, 2012, Albuquerque, NM; session moderator.- Mishra, P.K., K.L.
**Kuhlman**& M. Hilpert,*Recent Advances in Hydrogeology: Innovative Solutions to Single and Multi-phase Flow and Transport Issues*, American Geophysical Union, San Francisco, CA, Fall Meeting 2011; two oral sessions (H33J, H34B) and one poster session (H31D). - Bridge J.W., K.L.
**Kuhlman**& C.B. Graham,*High-Resolution Hydrogeophysical Characterization of Soils and Aquifers From Microscale to Field Scale*, American Geophysical Union, San Francisco, CA, Fall Meeting 2010, oral session H21K.

Software is free for use and redistribution via the MIT license; direct any comments or problems to me (contact info at top of the page). I will try my best to help or accommodate you. Programs are command-line programs (i.e., no graphical interface), either compiled from fortran source using a simple makefile or as python scripts. If you cannot compile or run the program yourself, I may be able to help you compile it or set up the right environment to get it working. I can provide a Windows executable from fortran source.

For most of the programs below, the files are available both as a zip archive (ZIP) and as a directory for browsing (DIR). Most programs have links to a Mercurial (bitbucket) or Git (github) repository.

These repositories are typically obtained by cloning the public repository (hosted at github or bitbucket) using the git or hg clients on your workstation. See help documentation on github and bitbucket sites.

- (github) ·
**Laplace Transform Analytic Element Method**

A Fortran03 program implementing the LT-AEM for circular (including wells as a special case) and elliptical elements (including lines as a special case), used in Bakker & Kuhlman (2011), my dissertation, and Kuhlman & Neuman (2009). This program is a continual work in progress, but it is functional and useful. This program includes an implementation of the same modified Mathieu library also given as a stand-alone Python library below. - (github) ·
**Unconfined well-test simulator**

A Fortran03 program implementing the major unconfined well test solutions (Neuman, 1972; Neuman, 1974; Moench, 1995; Moench, 2001; Mishra & Neuman, 2010; Mishra & Neuman, 2011; Malama, 2011), and special degenerate cases (Theis, 1935; Hantush, 1961). The program executes in parallel on multiple processors using OpenMP. It can compute a time-series solution at a point (including the derivative) or a contour map of drawdown in space at a given time. The software output is easily used with common plotting and parameter estimation software. The git repository also includes data from three large-scale unconfined pumping tests. - (bitbucket) ·
**Multiporosity well test simulator**

A Fortran03 program implementing a radially symmetric solution to the multiporosity governing equations (Kuhlman et al., 2015). The multiporosity solution is a generalization of dual-porosity flow in low-permeability fractured rocks. The distribution of a potentially infinite number of porosities (i.e., an extension of double or triple porosity) degenerates to multiple published solutions (Warren & Root, 1963; Kazemi, 1969; Clossman, 1975). The multiporosity solution is the logical combination of the multirate solute transport solution of Haggerty & Gorelick (1995) with the double-porosity pressure-driven well test solution of Warren & Root. - (bitbucket) ·
**Moving water table unconfined slug-test solution**

A Fortran90 program implementing the unconfined (MWT solution of Neuman 1972) slug test solution, which accounts for linearized momentum balance in the slug test interval (i.e., oscillating slug responses) discussed in Malama et al. (2011), and includes the extensions for predicting drawdown during cross-hole slug tests, as discussed in Malama et al. (2016). This program computes the double numerical inverse Laplace and Hankel transforms (infinite Hankel transform in 2011 solution, and finite Hankel transform in 2016 solution) to compute the time-domain solution given in the paper. A makefile and source are currently included (email me if you cannot build it yourself).*NB: The extended or unified Hyder solution from Malama et al. (2011) (which extended the confined Hyder solution to unconfined conditions with a constant head boundary at the water table) is a MATLAB/Octave script that is also available. If interested, please contact me or the corresponding author.* - (ZIP)
(DIR) ·
**Leaky-unconfined aquifer flow solution**

A Fortran90 program for computing the leaky-unconfined solutions discussed in Malama, Kuhlman & Barrash (2007; 2008). This program does the double numerical inverse Laplace and Hankel transforms required to compute the time-domain solutions given in the papers. A Windows binary and documentation are included. - (ZIP)
(DIR)
(bitbucket) ·
**Electrokinetic pumping test response solution**

A Fortran90 program for computing the electrokinetic response from pumping a well in an unconfined system, as discussed in Malama et al. (2009a; 2009b). The numerical Laplace inverse transform uses the Stehfest algorithm, while a combination of tanh-sinh quadrature and accelerated Gauss-Lobatto quadrature (between the zeros of the J_{0}Bessel function) is used for the numerical Hankel inverse transform. These algorithms are much faster and more accurate than the algorithms used in the leaky-unconfined papers listed above. - (ZIP)
(DIR) ·
**Unsaturated flow from an elliptical pipe**

A set of MATLAB/Octave scripts for computing Quasilinear infiltration solutions from Kuhlman & Warrick (2008). These scripts compute the required ratios of Mathieu functions, compute the solution, and plot the figures shown in the paper. - (ZIP)
(DIR)
(bitbucket)
·
**Modified Mathieu functions of complex parameter**

A Python module for computing both angular and radial modified Mathieu functions (and their derivatives with respect to argument) of integer order and general complex Mathieu parameter (Re(q)<0), see a pdf summary of Mathieu function library output, and a numerical comparison to a published C++ Mathieu function library for real q (Alhargan, 2000). My Python module requires numpy and scipy. This same functionality is also implemented in Fortran in my LT-AEM program available above. - (ZIP)
(DIR) ·
**Cokriging hydraulic heads with boundary conditions**

A set of MATLAB/Octave scripts for cokriging heads taking boundary conditions into account, without requiring the user to explicitly come up with the cross-variograms, since they are given explicitly by the physics of the problem (Kuhlman & Pardo-Igúzquiza, 2010). These scripts set up and solve the cokriging equations using either the true derivative or finite-difference approximations, and can make all the figures and examples given in the paper. - (ZIP)
(PDF) ·
**Scrabble rhyming word lists**

A rhyming dictionary sorts words alphabetically starting at the last letter of the word, working back towards the front of the word. This type of list is useful when you are trying to come up with legal Scrabble words, for example, ending in “–work” (webwork, cribwork, beadwork, yardwork, needlework, etc.) or “–horse” (seahorse, redhorse, racehorse, packhorse, warhorse, etc.). This is not a true rhyming dictionary, as homophones are not put together. This list would really only be useful for Scrabble players. The zip archive includes the list-generating Python script, the OSPD4 word lists, and the LaTeX-generating Python script (e.g., see 11-letter word list). - (DIR) ·
**2D steady-state finite-difference solution to Richards' equation**

A MATLAB/Octave script modified and improved from an Excel-spreadsheet developed by Ty Ferré for his HWR-505 vadose zone hydrology class. It is not highly optimized, but it is faster than the spreadsheet version. It has no loops except for iteration. See fipy simulations below for transient Richards solution. - (bitbucket) ·
**Transient 2D Thermohaline & 1D Richards flow simulators**

Mercurial repository includes implementation of solutions to the transient 1D Richards' equation (richards1D.py) and the 2D fully implicit thermohaline convection problem (density_flow_2D.py – see animation of Elder problem) using the fipy finite-volume PDE solver framework. These are not finished products, but I find them useful. Other scripts are also in this repository that solve other problems I have had (e.g., non-linear heat conduction in salt, electrokinetic response of aquifer to pumping, thermoporoelastic response due to heating a borehole, and double-porosity flow and solute transport). - (github) ·
**Open-source version of nSIGHTS**

The MS-Windows graphical interface n-dimensional Statistical Inverse Graphical Hydraulic Test Simulator (nSIGHTS) is a hydraulic testing simulator developed originally at INTERA Canada (now Geofirma). Software development was originally funded by the Department of Energy through my employer (SNL) for the WIPP project. A 2012 New Mexico Small Business Assistance grant funded the conversion of it to open source (BSD license). The nSIGHTS GitHub repository includes the C++ source, documentation, and compiled installers. Building it requires the Microsoft Visual Studio 2010 compiler. The C++ source is free to download, build, use, and modify.*I was not involved in the original development of this software*, but I have recently worked on maintaining and bug-fixing the open-source version of it. - (github) ·
**Numerical Inverse Laplace Transform in mpmath or numpy/scipy**

I contributed numerical inverse Laplace transform functionality to the python-based multiprecition library mpmath, specifically the de Hoog, Knight & Stokes, Stehfest, and fixed Talbot algorithms. The documentation for this feature gives more details regarding the capabilities, limitations, and some simple test cases. These same three algorithms are also available in a bitbucket repository of a numpy/scipy version of the mpmath library that uses fixed precision.