NumericalEnvironmental

Walt McNab's repositories

Pumping_Test_Interpretation_with_Python

This is a Python (version 2 or 3) script to aid in the interpretation of aquifer pumping tests using the Theis (1935) and Hantush and Jacob (1955) solutions for confined and leaky confined aquifers, respectively, as well as a simple numerical method-of-lines solution for an unconfined aquifer with wellbore storage.

Language:Python24 20

Coupled_gas_water_flow_and_transport_porous_media

This Julia language (version 0.5) script solves the coupled partial differential equations associated with (1) compressible gas flow in partially water-saturated porous media, (2) partially-saturated flow of water (essentially amounting to a solution of the Richards equation using the Van Genuchten equations for capillarity and relative permeability), and (3) advective-dispersive transport of both aqueous-phase solutes and gas components in the respective media.

Language:Julia10 20

Coupled_gas_water_transport_with_geochemistry

A two-phase (water + gas) reactive transport model, written in Julia, that also simulates geochemical reactions between water-gas-mineral phases via a link with PHREEQC.

Language:Julia10 10

Unsaturated_Flow_2D_Steady-state

This script, written in Julia (v. 1.0), solves the steady-state form of the Richards equation for variably-saturated flow in porous media in a two-dimensional vertical cross section.

Language:Julia9 10

RBF-based_correlated_random_field_generator

This is a streamlined python 3 script for generating spatially-correlated random fields in 2-D or 3-D using a radial basis function interpolator.

Language:Python7 10

random_field

Python-based script for generating realistic-looking correlated random fields in 2-D or 3-D; includes a tkinter-based user interface

Language:Python6 10

Theis_Well_Field_Drawdown

This python script uses scipy functionality (special function + spatial classes) to calculate the drawdown associated with a pumping well field according to the Theis equation.

Language:Python6 10

MFP

3-D multiphase flow model, written in D-language, plus sample input files

Language:D400

Partially_saturated_groundwater_flow_model_in_julia

This julia script is a complete groundwater flow model based on the integral finite difference method.

Language:Julia4 2 1

Simplified_surface_water-groundwater_screening_model

This python script models the transient flow of water through a network of connected stream reaches of differing slopes and widths in accordance with the Manning equation, also accounting for discharge through the stream bed to underlying groundwater.

Language:D4 10

vadose_multiphase

A python script for modeling 1-D multiphase flow, based on a method-of-lines approach to solving a modified form of the Richards equation

Language:D4 20

Facilitator_scripts_for_reactive_flow_through_fractured_rock_with_PHREEQC

This repository contains python scripts for pre- and post-processing PHREEQC simulations for flow through a dual-porosity column, emulating reactive transport through fractured porous media.

Language:Python3 1 1

Numerical_Solver_in_Python_for_Soil_Infiltration_Tests

This is a short Python 3 script that uses the SciPy package’s integrate.odeint method to solve a set of coupled ordinary differential equations approximating the infiltration of a wetting front into unsaturated soil.

Language:Python3 10

Python_Script_for_Statistical_Resampling_and_Modeling_of_Multivariate_Water_Quality_Data

This special-purpose Python script reads a multivariate water quality data set, generates a larger synthetic data set with matching parameter probability distributions and pairwise correlations, and runs PHREEQC speciation simulations on the synthetic data.

Language:Python3 10

Diamond-Square_Terrain_and_3D_Lattice_Generation_with_Python_and_Pandas

This is a relatively short python script (works with versions 2 or 3) for generating fractal surfaces or volumes, starting with an arbitrary initial lattice.

Language:Python2 10

Heuristic_model_for_aquifer-aquitard_structures_in_fluvial_environments

This julia script models the formation of aquifer-aquitard structures in a fluvial setting.

Language:Julia2 10

A_Prototype_Lake_Package_for_Solute_Transport_Modeling_with_PHAST

This Python script provides a limited groundwater-lake water interaction modeling capability to the U.S. Geological Survey’s PHAST reactive transport model.

Language:Python1 20

Heuristic_simulation_of_erosion_and_hydrologic_feature_formation

This julia language script reads an initial N x N grid of elevations across an x-y plane and simulates the formation of erosional and hydrologic features in response to user-supplied input parameters.

Language:Julia1 10

Leaky_Aquifer_Analytical_Model_for_Gas_Flow

This Python script models the steady-state pressure distribution and gas mass fluxes associated with operation of a soil vapor extraction (SVE) system in a leaky aquifer setting.

Language:Python1 10

Numerical_Toolkit_for_Approximating_Groundwater_Flow_Field_from_Scattered_Data

This Julia script imports and interpolates scattered hydraulic head and hydraulic conductivity data (for a single-layer aquifer) using ordinary kriging.

Language:Julia1 10

VOC_Plume_Meta-analysis_with_Python

This is a python script I have written to read and parse a large groundwater quality database (tens of thousands of wells, with hundreds of thousands of sample events) into individual groundwater plumes, subject to subsequent spatial analyses to develop a simple set of plume metrics.