This article has multiple issues. Please help improve it or discuss these issues on the talk page . (Learn how and when to remove these messages)
|
MODFLOW is the U.S. Geological Survey modular finite-difference flow model, which is a computer code that solves the groundwater flow equation. The program is used by hydrogeologists to simulate the flow of groundwater through aquifers. The source code is free public domain software, [1] written primarily in Fortran, and can compile and run on Microsoft Windows or Unix-like operating systems.
Since its original development in the early 1980s, [2] the USGS has made six major releases, and is now considered to be the de facto standard code for aquifer simulation. There are several actively developed commercial and non-commercial graphical user interfaces for MODFLOW.
MODFLOW was constructed in what was in 1980's called a modular design. This means it has many of the attributes of what came to be called object-oriented programming. For example, capabilities (called "packages") that simulate subsidence or lakes or streams, can easily be turned on and off and the execution time and storage requirements of those packages go away entirely. If a programmer wants to change something in MODFLOW, the clean organization makes it easy. Indeed, this kind of innovation is exactly what was anticipated when MODFLOW was designed.
Importantly, the modularity of MODFLOW makes it possible for different Packages to be written that are intended to address the same simulation goal in different ways. This allows differences of opinion about how system processes function to be tested. Such testing is an important part of multi-modeling, or alternative hypothesis testing. Models like MODFLOW make this kind of testing more definitive and controlled. This results because other aspects of the program remain the same. Tests become more definitive because they become less prone to being influenced unknowingly by other numerical and programming differences.
The governing partial differential equation for a confined aquifer used in MODFLOW is:
where
The finite difference form of the partial differential in a discretized aquifer domain (represented using rows, columns and layers) is:
where
This equation is formulated into a system of equations to be solved as:
where
or in matrix form as:
where
The USGS throughout the 1970s had developed several hundred models, written in different dialects of FORTRAN. At the time, it was common practice to rewrite a new model to fit the need of a new groundwater scenario. The concept for MODFLOW was originally designed in 1981 to provide a common modular groundwater model, which could be compiled on multiple platforms without major (or any) modification, and can read and write common formats. Different aspects of the groundwater system would be handled using the modules, similar to the idea of a "component stereo system". The original name of the code was "The USGS Modular Three-Dimensional Finite-Difference Ground-Water Flow Model", or informally as "The Modular Model". The name MODFLOW was coined several years after the initial code development, which started in 1981. [2]
The first version of MODFLOW [3] was published on December 28, 1983, and was coded entirely in FORTRAN 66. The source code for this version is listed in USGS Open File Report 83-875 referred to above.
This version of MODFLOW [4] was rewritten in FORTRAN 77, and was originally released on July 24, 1987. The current version of MODFLOW-88 is 2.6, released on September 20, 1996.
MODPATH, was initially developed in 1989 to post-process the steady-state MODFLOW-88 data to determine three-dimensional pathlines of particles. This innovation has been indispensable for the fields of contaminant hydrogeology. It is still used as a post-processor in recent versions of MODFLOW.
A separate program, MODFLOWP, was developed in 1992 to estimate various parameters used in MODFLOW. This program was eventually built into MODFLOW-2000.
MODFLOW-96 (version 3.0) was originally released on December 3, 1996, and is a cleaned-up and revised continuation of MODFLOW-88. [5] [6] There are three final releases of MODFLOW-96:
Several graphical interfaces were first developed using the MODFLOW-96 code.
MODFLOW-2000 (version 1.0; version numbering was reset) was released on July 20, 2000, which merged MODFLOWP and HYDMOD codes into the main program and has integrated observation, sensitivity analysis, parameter estimation, and uncertainty evaluation capabilities. [7] Many new packages and enhancements were also included, including new solvers, stream and saturated flow packages. The internal design concepts also changed from previous versions, such that packages, processes and modules are distinct. This version was coded in a mixture of FORTRAN 77, Fortran 90, and one solver was programmed in C. MODFLOW-2000 can also be compiled for parallel computing, which can allow multiple processors to be used to increase model complexity and/or reduce simulation time. The parallelization capability is designed to support the sensitivity analysis, parameter estimation, and uncertainty analysis capabilities of MODFLOW-2000.
The final version of MODFLOW-2000 (or MF2K) is version 1.19.01, released on March 25, 2010. There are four related or branched codes based on MODFLOW-2000:
MODFLOW-2005 [8] differs from MODFLOW-2000 in that the sensitivity analysis, parameter estimation, and uncertainty evaluation capabilities are removed. Thus, the support for these capabilities now falls to "clip on" codes that are supported externally to the MODFLOW support effort. In addition, the code was reorganized to support multiple models within one MODFLOW run, as needed for the LGR (Local Grid Refinement) capability. [9] MODFLOW-2005 is written primarily in Fortran 90 and C, with C being used for one solver.
The current version of MODFLOW-2005 is version 1.12.00, released on February 3, 2017. Related or branched codes include:
MODFLOW 6 (MF6), first released in 2017, is the sixth core version of MODFLOW to be released by the USGS. [16] This release is a rewrite of MODFLOW following an object-oriented programming paradigm in Fortran, and provides a platform that includes the capabilities from several previous MODFLOW-2005 versions, including MODFLOW-NWT, MODFLOW-USG, and MODFLOW-LGR. [17] MODFLOW 6 supports structured or unstructured grids, has full support for the Newton-Raphson formulation, and has a unique Water Mover Package that allows flows to be routed between the advanced packages, including the Streamflow Routing, Lake, Multi-Aquifer Well, and Unsaturated Zone Flow Packages. MODFLOW 6 also contains a Groundwater Transport (GWT) model that simulates transient three-dimensional solute transport on structured or unstructured grids and through the advanced flow and mover packages. An Application Programming Interface (API) is also available for MODFLOW 6, which allows the program to be coupled with other models or controlled with popular scripting languages, such as Python. While there are a few features lacking in the current release that are supported in MODFLOW-2005, most of the popular capabilities in previous MODFLOW versions are available in MODFLOW 6. The current version is 6.2.2, released July 30, 2021. [18]
The names in this table are the labels used to turn MODFLOW capabilities on and off via a key input file. Most capabilities have many alternatives or can be omitted, but the ones related to the BASIC Package are always required. Many of the capabilities introduced are supported in later versions, though the grid change enabled with MODFLOW-USG and MODFLOW 6 meant that such backward compatibility was rather selective.
Name | Long name | Version introduced |
---|---|---|
Basic Package and its Components | ||
BAS | Basic Package | original |
OC | Output Control | original |
DIS | Discretization | MODFLOW-2000 (1.0) |
DISU | Unstructured Discretization | MODFLOW-USG (1.0) |
DISV | Discretization by Vertices | MODFLOW 6 (1.00) |
IC | Initial Conditions | MODFLOW 6 (1.00) |
Groundwater flow packages | ||
BCF | Block-Centered Flow Package | original |
CLN | Connected Linear Network Process | MODFLOW-USG (1.0) |
GNC | Ghost Node Correction Package | MODFLOW-USG (1.0) |
HFB | Horizontal Flow Barrier Package | MODFLOW-88 |
HUF | Hydrogeologic Unit Flow Package | MODFLOW-2000 (1.1) |
LPF | Layer-Property Flow Package | MODFLOW-2000 (1.0) |
NPF | Node Property Flow | MODFLOW 6 (1.00) |
SWI2 | Seawater Intrusion Package | MODFLOW-2005 (1.11) |
UPW | Upstream Weighting Package | MODFLOW-NWT (1.0) |
UZF | Unsaturated-Zone Flow Package | MODFLOW-2005 (1.2) |
Conjunctive Use and Land Use Simulation | ||
FMP | Farm Process | MODFLOW-FMP |
SWO | Surface Water Operations | MODFLOW-OWHM (2.0) |
Specified Head boundary packages | ||
CHD | Constant-Head Boundary / Time-Variant Specified-Head | MODFLOW-88 |
FHB | Flow and Head Boundary Package | MODFLOW-96 (3.2) |
Specified flux boundary packages | ||
FHB | Flow and Head Boundary Package | MODFLOW-96 (3.2) |
RCH | Recharge Package | original |
WEL | Well Package | original |
Head-dependent flux boundary packages | ||
DAF | DAFLOW | MODFLOW-96 |
DRN | Drain Package | original |
DRT | Drain Return Package | MODFLOW-2000 (1.1) |
ETS | Evapotranspiration Segments Package | MODFLOW-2000 (1.1) |
EVT | Evapotranspiration Package | original |
GHB | General-Head Boundary Package | original |
LAK | Lake Package | MODFLOW-2000 (1.1) |
MAW | Multi-Aquifer Well | MODFLOW 6 (1.00) |
MNW | Multi-Node, Drawdown-Limited Well Package | MODFLOW-2000 (1.11) |
RES | Reservoir Package | MODFLOW-88 (2.6) |
RIP | Riparian Evapotranspiration Package | MODFLOW-OWHM (1.0) |
RIV | River Package | original |
SFR | Streamflow-Routing Package | MODFLOW-2000 (1.14.00) |
STR | Stream Package | MODFLOW-88 |
SWR | Surface-Water Routing Process | MODFLOW-NWT 1.08 |
UZF | Unsaturated-Zone Flow Package | MODFLOW-2005 (1.2) |
Solvers | ||
DE4 | Direct Solver Package | MODFLOW-88 (2.5) |
GMG | Geometric Multigrid Solver | MODFLOW-2000 (1.15.00) |
LMG | Link-AMG Package | MODFLOW-2000 (1.4) [Note 1] |
NWT | Newton-Raphson | MODFLOW-NWT (1.0) |
PCG | Preconditioned Conjugate-Gradient Package | MODFLOW-88 |
PCGN | Preconditioned Conjugate Gradient Solver with Improved Nonlinear Control | MODFLOW-2005 (1.9.0) |
SIP | Strongly Implicit Procedure Package | original |
SMS | Sparse Matrix Solver | MODFLOW-USG (1.0) |
SOR | Slice Successive Over-Relaxation Package | original |
Miscellaneous packages | ||
GAG | Gage | MODFLOW-2000 |
HYD | HYDMOD | MODFLOW-2000 (1.1) |
IBS | Interbed-Storage | MODFLOW-88 |
KDEP | Hydraulic-Conductivity Depth-Dependence Capability | MODFLOW-2000 (1.12) |
LMT | Link-MT3DMS | MODFLOW-2000 (1.5) |
LVDA | Model-Layer Variable-Direction Horizontal Anisotropy Capability | MODFLOW-2000 (1.12) |
MVR | Water Mover | MODFLOW 6 (1.00) |
STO | Storage | MODFLOW 6 (1.00) |
SUB | Subsidence and Aquifer-System Compaction | MODFLOW-2000 (1.12) |
SWT | Subsidence and Aquifer-System Compaction Package for Water-Table Aquifers | MODFLOW-2000 (1.18) |
CSUB | Skeletal Storage, Compaction, and Subsidence | MODFLOW 6 (6.1.0) |
UTL | Utility | original |
Observation process input files | ||
OBS | Input File For All Observations | MODFLOW-2000 |
HOB | Head-Observation | MODFLOW-2000 |
DROB | Drain Observation | MODFLOW-2000 |
DTOB | Drain Return Observation | MODFLOW-2000 |
RVOB | River Observation | MODFLOW-2000 |
GBOB | General-Head-Boundary Observation | MODFLOW-2000 |
CHOB | Constant-Head Flow Observation | MODFLOW-2000 |
ADV | Advective-Transport Observation | MODFLOW-2000 (1.0) |
STOB | Stream Observation | MODFLOW-2000 |
Obsolete packages | ||
GFD | General Finite-Difference | MODFLOW-88 to 96 |
TLK | Transient Leakage | MODFLOW-88 to 96 |
There are several graphical interfaces to MODFLOW, which often include the compiled MODFLOW code with modifications. These programs aid the input of data for creating MODFLOW models.
Non-commercial MODFLOW versions are free, however, their licensing usually limit the use to non-profit educational or research purposes.
Commercial MODFLOW programs are typically used by governments and consultants for practical applications of MODFLOW to real-world groundwater problems. Professional versions of MODFLOW are generally priced at a minimum of around $1000 and typically range upward to US$7000. This is a list of commercial programs for MODFLOW:
All current versions of these programs run only on Microsoft Windows, however previous versions of GMS (up to Version 3.1) were compiled for several Unix platforms.
Combustion, or burning, is a high-temperature exothermic redox chemical reaction between a fuel and an oxidant, usually atmospheric oxygen, that produces oxidized, often gaseous products, in a mixture termed as smoke. Combustion does not always result in fire, because a flame is only visible when substances undergoing combustion vaporize, but when it does, a flame is a characteristic indicator of the reaction. While activation energy must be supplied to initiate combustion, the heat from a flame may provide enough energy to make the reaction self-sustaining. The study of combustion is known as combustion science.
Hydrogeology is the area of geology that deals with the distribution and movement of groundwater in the soil and rocks of the Earth's crust. The terms groundwater hydrology, geohydrology, and hydrogeology are often used interchangeably, though hydrogeology is the most commonly used.
In power engineering, the power-flow study, or load-flow study, is a numerical analysis of the flow of electric power in an interconnected system. A power-flow study usually uses simplified notations such as a one-line diagram and per-unit system, and focuses on various aspects of AC power parameters, such as voltages, voltage angles, real power and reactive power. It analyzes the power systems in normal steady-state operation.
In fluid mechanics, materials science and Earth sciences, permeability is a measure of the ability of a porous material to allow fluids to pass through it.
Darcy's law is an equation that describes the flow of a fluid through a porous medium and through a Hele-Shaw cell. The law was formulated by Henry Darcy based on results of experiments on the flow of water through beds of sand, forming the basis of hydrogeology, a branch of earth sciences. It is analogous to Ohm's law in electrostatics, linearly relating the volume flow rate of the fluid to the hydraulic head difference via the hydraulic conductivity. In fact, the Darcy's law is a special case of the Stokes equation for the momentum flux, in turn deriving from the momentum Navier-Stokes equation.
Optical flow or optic flow is the pattern of apparent motion of objects, surfaces, and edges in a visual scene caused by the relative motion between an observer and a scene. Optical flow can also be defined as the distribution of apparent velocities of movement of brightness pattern in an image.
Large eddy simulation (LES) is a mathematical model for turbulence used in computational fluid dynamics. It was initially proposed in 1963 by Joseph Smagorinsky to simulate atmospheric air currents, and first explored by Deardorff (1970). LES is currently applied in a wide variety of engineering applications, including combustion, acoustics, and simulations of the atmospheric boundary layer.
Hydraulic head or piezometric head is a specific measurement of liquid pressure above a vertical datum.
Used in hydrogeology, the groundwater flow equation is the mathematical relationship which is used to describe the flow of groundwater through an aquifer. The transient flow of groundwater is described by a form of the diffusion equation, similar to that used in heat transfer to describe the flow of heat in a solid. The steady-state flow of groundwater is described by a form of the Laplace equation, which is a form of potential flow and has analogs in numerous fields.
Infiltration is the process by which water on the ground surface enters the soil. It is commonly used in both hydrology and soil sciences. The infiltration capacity is defined as the maximum rate of infiltration. It is most often measured in meters per day but can also be measured in other units of distance over time if necessary. The infiltration capacity decreases as the soil moisture content of soils surface layers increases. If the precipitation rate exceeds the infiltration rate, runoff will usually occur unless there is some physical barrier.
GMS is water modeling application for building and simulating groundwater models from Aquaveo. It features 2D and 3D geostatistics, stratigraphic modeling and a unique conceptual model approach. Currently supported models include MODFLOW, MODPATH, MT3DMS, RT3D, FEMWATER, SEEP2D, and UTEXAS.
Visual MODFLOW (VMOD) is a graphical interface (GUI) for the open source groundwater modeling engine MODFLOW. VMOD was developed by Waterloo Hydrogeologic and first released in 1994, the first commercially available GUI for MODFLOW. In May 2012 a .NET version of the software was rebranded as Visual MODFLOW Flex. The program includes proprietary extensions, such as MODFLOW-SURFACT, MT3DMS and a three-dimensional model explorer. Visual MODFLOW supports MODFLOW-2000, MODFLOW-2005, MODFLOW-NWT, MODFLOW-LGR, MODFLOW-SURFACT, and SEAWAT.
Groundwater models are computer models of groundwater flow systems, and are used by hydrologists and hydrogeologists. Groundwater models are used to simulate and predict aquifer conditions.
HydroGeoSphere (HGS) is a 3D control-volume finite element groundwater model, and is based on a rigorous conceptualization of the hydrologic system consisting of surface and subsurface flow regimes. The model is designed to take into account all key components of the hydrologic cycle. For each time step, the model solves surface and subsurface flow, solute and energy transport equations simultaneously, and provides a complete water and solute balance.
The Benjamin–Bona–Mahony equation is the partial differential equation
A hydrologic model is a simplification of a real-world system that aids in understanding, predicting, and managing water resources. Both the flow and quality of water are commonly studied using hydrologic models.
In finance, the T-model is a formula that states the returns earned by holders of a company's stock in terms of accounting variables obtainable from its financial statements. The T-model connects fundamentals with investment return, allowing an analyst to make projections of financial performance and turn those projections into a required return that can be used in investment selection.
In fluid dynamics, a cnoidal wave is a nonlinear and exact periodic wave solution of the Korteweg–de Vries equation. These solutions are in terms of the Jacobi elliptic function cn, which is why they are coined cnoidal waves. They are used to describe surface gravity waves of fairly long wavelength, as compared to the water depth.
The finite water-content vadose zone flux method represents a one-dimensional alternative to the numerical solution of Richards' equation for simulating the movement of water in unsaturated soils. The finite water-content method solves the advection-like term of the Soil Moisture Velocity Equation, which is an ordinary differential equation alternative to the Richards partial differential equation. The Richards equation is difficult to approximate in general because it does not have a closed-form analytical solution except in a few cases. The finite water-content method, is perhaps the first generic replacement for the numerical solution of the Richards' equation. The finite water-content solution has several advantages over the Richards equation solution. First, as an ordinary differential equation it is explicit, guaranteed to converge and computationally inexpensive to solve. Second, using a finite volume solution methodology it is guaranteed to conserve mass. The finite water content method readily simulates sharp wetting fronts, something that the Richards solution struggles with. The main limiting assumption required to use the finite water-content method is that the soil be homogeneous in layers.
In geology, numerical modeling is a widely applied technique to tackle complex geological problems by computational simulation of geological scenarios.
{{cite book}}
: CS1 maint: multiple names: authors list (link){{citation}}
: CS1 maint: multiple names: authors list (link)