Dupuit Divide-River 1D#
Note
This page and its static assets are auto-generated by python -m tools.doc_gallery. The Sphinx build only reads committed PNG and JSON artifacts.
Steady synthetic groundwater-flow case used to validate the launcher workflow against the Dupuit analytical profile for a homogeneous unconfined aquifer with:
Case Setup#
Homogeneous strip, west divide, east fixed head, uniform recharge
Available solver variants: MODFLOW-NWT, MODFLOW 6, MODFLOW 6 irregular triangles, Boussinesq.
What It Shows#
validate one explicit no-flow side boundary through a simple 1D benchmark,
keep a clear analytical reference while extending beyond the two-head case.
One-sided no-flow behavior plus recharge and one imposed downstream head
Primary metrics: head-profile RMSE, max abs error, cross-row spread.
Solver-specific figures and metrics are shown side by side so the same benchmark can be read across backends.
Analytical Reference#
compared quantity: domain-averaged head profile along x
reference: steady Dupuit profile with west-side divide and east-side river head
Solver Coverage#
Default solver: MODFLOW-NWT
Available variants: MODFLOW-NWT, MODFLOW 6, MODFLOW 6 irregular triangles, Boussinesq
Metrics - Head-profile RMSE: 0.0380 m - Head-profile max abs error: 0.0465 m - Cross-row head spread: 2.34e-07 m
Config file:
validation_cases/analytical/steady/dupuit_divide_river_1d/config_modflownwt.tomlTolerances:
validation_cases/analytical/steady/dupuit_divide_river_1d/tolerances.tomlExpected shape: 5 x 40
python -m validation_cases.analytical.steady.dupuit_divide_river_1d.run_case --no-show --solver modflow_nwt
Metrics - Head-profile RMSE: 0.0339 m - Head-profile max abs error: 0.0418 m - Cross-row head spread: 2.41e-06 m
Config file:
validation_cases/analytical/steady/dupuit_divide_river_1d/config_modflow6.tomlTolerances:
validation_cases/analytical/steady/dupuit_divide_river_1d/tolerances_modflow6.tomlExpected shape: 5 x 40
python -m validation_cases.analytical.steady.dupuit_divide_river_1d.run_case --no-show --solver modflow6
Metrics - Head-profile RMSE: 0.0268 m - Head-profile max abs error: 0.0365 m - Cross-row head spread: 8.88e-16 m
Config file:
validation_cases/analytical/steady/dupuit_divide_river_1d/config_modflow6_irregular_tri.tomlTolerances:
validation_cases/analytical/steady/dupuit_divide_river_1d/tolerances_modflow6_irregular_tri.tomlExpected shape: 5 x 40
python -m validation_cases.analytical.steady.dupuit_divide_river_1d.run_case --no-show --solver modflow6_irregular_tri
Metrics - Head-profile RMSE: 0.0217 m - Head-profile max abs error: 0.0282 m - Cross-row head spread: 8.88e-16 m
Config file:
validation_cases/analytical/steady/dupuit_divide_river_1d/config_boussinesq.tomlTolerances:
validation_cases/analytical/steady/dupuit_divide_river_1d/tolerances.tomlExpected shape: 5 x 40
python -m validation_cases.analytical.steady.dupuit_divide_river_1d.run_case --no-show --solver boussinesq
Reproduce#
Run the underlying example or validation case with:
python -m validation_cases.analytical.steady.dupuit_divide_river_1d.run_case --no-show
Refresh the committed gallery artifacts with:
python -m tools.doc_gallery
Case Parameters#
Reference Parameters#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Minimum x coordinate of the analytical reference domain. |
0 |
|
|
Maximum x coordinate of the analytical reference domain. |
400 |
|
|
River head imposed by the reference solution. |
5 m |
|
|
Recharge rate used by the reference solution. |
1 mm/day |
|
|
Hydraulic conductivity used by the analytical or benchmark reference. |
0.0001 m/s |
|
|
Aquifer thickness used by the reference formulation. |
20 m |
|
|
Axis along which the validation profile is extracted. |
0 |
|
|
Side of the domain where the analytical no-flow divide is enforced. |
west |
|
Common Numerical Setup#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Case-specific configuration field workflow used by the validation benchmark. |
simulation |
|
|
Synthetic-domain length along the x axis. |
400.0 m |
|
|
Synthetic-domain length along the y axis. |
50.0 m |
|
|
Grid cell count along the x axis. |
40 |
|
|
Grid cell count along the y axis. |
5 |
|
|
Synthetic topography shape used by the benchmark. |
flat |
|
|
Base land-surface elevation of the synthetic topography. |
20 m |
|
|
Case-specific configuration field domain.zone_ids used by the validation benchmark. |
[] |
|
|
Depth model used to build the aquifer support. |
constant_thickness |
|
|
Aquifer or support thickness used by the benchmark. |
20.0 m |
|
|
External data families loaded by the benchmark. |
[recharge] |
|
|
Recharge data source mode used by the benchmark. |
synthetic |
|
|
Recharge values used by the benchmark forcing. |
[1 mm/day] |
|
|
Runoff ratio applied when converting recharge forcing to effective infiltration. |
0 |
|
|
Steady or transient flow regime used by the benchmark. |
steady |
|
|
Sink and source families activated in the benchmark. |
[recharge] |
|
|
Boundary-condition families activated in the benchmark. |
[east_side] |
|
|
Hydraulic parameter families explicitly configured by the benchmark. |
[K] |
|
|
Parameterization mode used for K. |
homogeneous |
|
|
Homogeneous K value used by the benchmark. |
1e-4 m/s |
|
|
Initial-condition policy used to start the benchmark. |
custom |
|
|
Scalar initial head or top offset used by the initial-condition policy. |
9.0 m |
|
|
Fixed head applied on the east side boundary. |
5.0 m |
|
Solver-Specific Overrides#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Planar support construction mode used by MODFLOW-NWT. |
resample_to_shape |
|
|
Planar support cell count along x used by MODFLOW-NWT. |
40 |
|
|
Planar support cell count along y used by MODFLOW-NWT. |
5 |
|
|
Planar support resampling mode used by MODFLOW-NWT. |
nearest |
|
|
Number of vertical layers used by MODFLOW-NWT. |
1 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Solver-specific override applied to MODFLOW 6. |
false |
|
|
Linear-solver complexity preset used by MODFLOW 6. |
COMPLEX |
|
|
Planar support construction mode used by MODFLOW 6. |
resample_to_shape |
|
|
Planar support cell count along x used by MODFLOW 6. |
40 |
|
|
Planar support cell count along y used by MODFLOW 6. |
5 |
|
|
Planar support resampling mode used by MODFLOW 6. |
nearest |
|
|
Number of vertical layers used by MODFLOW 6. |
1 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Solver-specific override applied to MODFLOW 6. |
false |
|
|
Linear-solver complexity preset used by MODFLOW 6. |
COMPLEX |
|
|
Number of vertical layers used by MODFLOW 6. |
1 |
|
|
Committed unstructured mesh file used by the irregular-mesh solver variant. |
../../../shared/mesh_bundles/sloping_substratum_irregular_tri_400x50/mesh_2d.msh |
|
|
Committed mesh-bundle directory used to recover support metadata for the irregular-mesh solver variant. |
../../../shared/mesh_bundles/sloping_substratum_irregular_tri_400x50 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Runtime backend selected for the in-house solver. |
petsc |
|
Acceptance Criteria#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Simulated observable compared against the reference solution. |
watertable_elevation |
|
|
Expected spatial output shape checked by the validation helper. |
[5, 40] |
|
Acceptance Criteria by Solver#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected shape: 5 x 40 |
|
|
Maximum accepted root-mean-square error for head profile. |
0.05 |
|
|
Maximum accepted absolute error for head profile. |
0.1 |
|
|
Maximum accepted cross-row spread for head profile. |
2e-05 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected shape: 5 x 40 |
|
|
Maximum accepted root-mean-square error for head profile. |
0.04 |
|
|
Maximum accepted absolute error for head profile. |
0.05 |
|
|
Maximum accepted cross-row spread for head profile. |
1e-05 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected shape: 5 x 40 |
|
|
Maximum accepted root-mean-square error for head profile. |
0.31 |
|
|
Maximum accepted absolute error for head profile. |
0.41 |
|
|
Maximum accepted cross-row spread for head profile. |
1e-09 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected shape: 5 x 40 |
|
|
Maximum accepted root-mean-square error for head profile. |
0.05 |
|
|
Maximum accepted absolute error for head profile. |
0.1 |
|
|
Maximum accepted cross-row spread for head profile. |
2e-05 |
|
Source Pointers#
validation_cases/README.mdvalidation_cases/analytical/steady/dupuit_divide_river_1d/README.mdvalidation_cases/analytical/steady/dupuit_divide_river_1d/reference.pyvalidation_cases/analytical/steady/dupuit_divide_river_1d/comparison.pyvalidation_cases/analytical/steady/dupuit_divide_river_1d/plotting.pyvalidation_cases/analytical/steady/dupuit_divide_river_1d/run_case.pyvalidation_cases/analytical/steady/dupuit_divide_river_1d/metadata.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/config_modflownwt.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/runtime_boussinesq.pyvalidation_cases/analytical/steady/dupuit_divide_river_1d/tolerances.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/tolerances_modflow6.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/tolerances_modflow6_irregular_tri.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/config_modflow6.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/config_modflow6_irregular_tri.tomlvalidation_cases/analytical/steady/dupuit_divide_river_1d/config_boussinesq.tomlvalidation_cases/shared/boussinesq_uniform_strip.py
Artifacts#
docs/source/_static/capability_gallery/validation/dupuit_divide_river_1d__modflownwt.pngdocs/source/_static/capability_gallery/validation/dupuit_divide_river_1d__modflow6.pngdocs/source/_static/capability_gallery/validation/dupuit_divide_river_1d__modflow6_irregular_tri.pngdocs/source/_static/capability_gallery/validation/dupuit_divide_river_1d__boussinesq.pngdocs/source/_static/capability_gallery/validation/dupuit_divide_river_1d_summary.jsonstores the displayed metrics plus source hashes used bypython -m tools.doc_gallery --check.