Solvers#

Scientific content behind the HydroModPy solver layer. Focuses on what each solver represents mathematically and which numerical choices are active, not on the runtime orchestration.

Operational pages live elsewhere

For backend selection and TOML snippets, see Solvers. For the package layout and runtime handoff, see Solver Architecture. The in-house Boussinesq theory is in Boussinesq Theory.

MODFLOW family#

MODFLOW flow family overview

Reading map for MODFLOW 6 and MODFLOW-NWT in HydroModPy: common contract, backend split, comparison discipline, and version-specific emphasis.

MODFLOW Flow Family
Governing equation and CVFD formulation

Cell-by-cell groundwater balance, control-volume finite difference statement, sign conventions, and the relationship to the solver-agnostic problem.

MODFLOW Governing Equation And CVFD Formulation
Package semantics and boundary conditions

How HydroModPy maps the [flow] payload to NPF, STO, RCHA, WEL, CHD, DRN, and EVT packages, with version differences flagged explicitly.

MODFLOW Package Semantics And Option Selection
MODFLOW 6 vs MODFLOW-NWT

Scientific contrast of the two MODFLOW backends: support, solver settings, package coverage, transport coupling, and reasons to keep both alive.

MODFLOW 6 Versus MODFLOW-NWT: Scientific Comparison
XT3D on irregular DISV meshes

Why XT3D matters when a DISV-style mesh leaves the strictly-orthogonal regime, with the method-choice evidence published in the gallery.

XT3D On Irregular DISV Meshes
MODFLOW 6 flow page

MODFLOW 6-specific notes for flow/modflow6: structured and runtime DISV support, XT3D policy, and the GWT downstream route.

MODFLOW 6 Flow
MODFLOW-NWT flow page

MODFLOW-NWT-specific notes for flow/modflow_nwt: legacy structured sgrid continuity, MODPATH and MT3DMS compatibility.

MODFLOW-NWT Flow
Comparison and method choice

Checklist to follow before attributing a difference between flow/modflow6 and flow/modflow_nwt to the numerical backend.

Comparison And Method Choice

Mesh and discretization#

Mesh and discretization strategies

The separation between physical problem, planar support, vertical layering, and solver-specific cell interpretation.

Mesh And Discretization Strategies
Field-to-cell parameter transfer

How heterogeneous fields (geology, hydraulic conductivity, forcing) are projected to solver cells across mesh families.

Field-To-Cell Parameter Transfer
Vertical representation and storage

How vertical layering and saturated thickness assumptions shape storage, transmissivity, and the unconfined response.

Vertical Representation And Storage Assumptions
Mesh quality and acceptance

Quality diagnostics that decide whether a discretization is acceptable before solver results are trusted.

Mesh Quality And Acceptance Criteria

Worked cases and reference matrix#

Worked MODFLOW cases

Executable MODFLOW-family examples that link assumptions to TOML inputs, package assembly, and inspected outputs.

Worked MODFLOW Cases
Dupuit fixed-head 1D

The simplest analytical-style benchmark: Dupuit unconfined flow against fixed-head boundaries in a 1D strip.

Worked MODFLOW Case: Dupuit Fixed Head 1D
Linearized unconfined periodic 1D

Linearized unconfined response to a periodic recharge signal, with explicit comparison against the analytical solution.

Worked MODFLOW Case: Linearized Unconfined Periodic Recharge 1D
Nancon transient NWT

The Nancon basin transient MODFLOW-NWT case with ETP routed through the EVT package, the canonical real-basin teaching reference.

Worked MODFLOW Case: Nancon Transient NWT With ETP To EVT
Transport coupling

Which downstream transport solver consumes each MODFLOW-family flow field: MODFLOW 6 GWT, MODPATH, or MT3DMS.

Transport Coupling
Solver capability matrix

Compact inventory of what each solver can represent today and where the active development edges sit.

Solver Capability Matrix