Linearized Unconfined Boundary Piecewise 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.
Transient synthetic groundwater-flow case used to validate the launcher workflow against the linearized unconfined analytical response for:
Case Setup#
Homogeneous strip, fixed east head, piecewise west-head series from CSV, no recharge
Available solver variants: MODFLOW-NWT, MODFLOW 6, MODFLOW 6 irregular triangles, Boussinesq.
What It Shows#
validate multi-step lateral boundary forcing through the launcher-managed CSV forcing path,
exercise analytical superposition beyond the single-step benchmark,
keep one compact transient 1D case where the reference remains closed-form.
Multi-step boundary forcing and superposition behavior beyond the single-step case
Primary metrics: space-time RMSE, space-time max abs error, final-profile RMSE, 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: full transient head profile matrix along x
reference: linearized 1D unconfined response to a piecewise west-boundary series
Solver Coverage#
Default solver: MODFLOW-NWT
Available variants: MODFLOW-NWT, MODFLOW 6, MODFLOW 6 irregular triangles, Boussinesq

Metrics - Space-time RMSE: 0.0055 m - Space-time max abs error: 0.0318 m - Final-profile RMSE: 0.0022 m - Cross-row head spread: 4.67e-07 m
Config file:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_modflownwt.tomlTolerances:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances.tomlExpected output: 24 periods, spatial shape 5 x 50
python -m validation_cases.analytical.transient.linearized_unconfined_boundary_piecewise_1d.run_case --no-show --solver modflow_nwt

Metrics - Space-time RMSE: 0.0055 m - Space-time max abs error: 0.0319 m - Final-profile RMSE: 0.0022 m - Cross-row head spread: 1.69e-06 m
Config file:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_modflow6.tomlTolerances:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances_modflow6.tomlExpected output: 24 periods, spatial shape 5 x 50
python -m validation_cases.analytical.transient.linearized_unconfined_boundary_piecewise_1d.run_case --no-show --solver modflow6

Metrics - Space-time RMSE: 0.0055 m - Space-time max abs error: 0.0323 m - Final-profile RMSE: 0.0023 m - Cross-row head spread: 0.00e+00 m
Config file:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_modflow6_irregular_tri.tomlTolerances:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances_modflow6_irregular_tri.tomlExpected output: 24 periods, spatial shape 5 x 50
python -m validation_cases.analytical.transient.linearized_unconfined_boundary_piecewise_1d.run_case --no-show --solver modflow6_irregular_tri

Metrics - Space-time RMSE: 0.0014 m - Space-time max abs error: 0.0082 m - Final-profile RMSE: 0.0004 m - Cross-row head spread: 1.12e-15 m
Config file:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_boussinesq.tomlTolerances:
validation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances.tomlExpected output: 24 periods, spatial shape 5 x 50
python -m validation_cases.analytical.transient.linearized_unconfined_boundary_piecewise_1d.run_case --no-show --solver boussinesq
Reproduce#
Run the underlying example or validation case with:
python -m validation_cases.analytical.transient.linearized_unconfined_boundary_piecewise_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. |
100 |
|
|
Baseline hydraulic head around which the linearized reference is expressed. |
10 m |
|
|
Sequence of west-boundary head levels used by the transient boundary forcing. |
[10.05 m, 10.05 m, 10.05 m, 10.05 m, 10.15 m, 10.15 m, 10.15 m, 10.15 m, 9.95 m, 9.95 m, 9.95 m, 9.95 m, 10.25 m, 10.25 m, 10.25 m, 10.25 m, 10.05 m, 10.05 m, 10.05 m, 10.05 m, 10 m, 10 m, 10 m, 10 m] |
|
|
Hydraulic conductivity used by the analytical or benchmark reference. |
0.0001 m/s |
|
|
Reference saturated thickness used by the linearized formulation. |
10 m |
|
|
Specific yield used by the transient reference formulation. |
0.1 - |
|
|
Axis along which the validation profile is extracted. |
0 |
|
|
Number of terms retained in the analytical series expansion. |
400 |
|
|
Reference time step used by the analytical evaluator. |
43200 s |
|
Common Numerical Setup#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Case-specific configuration field workflow used by the validation benchmark. |
simulation |
|
|
Simulation start time used by the benchmark. |
2003-01-01 00:00:00 |
|
|
Simulation end time used by the benchmark. |
2003-01-12 23:00:00 |
|
|
Nominal time step used by the benchmark. |
12 hour |
|
|
Synthetic-domain length along the x axis. |
100.0 m |
|
|
Synthetic-domain length along the y axis. |
10.0 m |
|
|
Grid cell count along the x axis. |
50 |
|
|
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 |
|
|
Depth model used to build the aquifer support. |
constant_thickness |
|
|
Aquifer or support thickness used by the benchmark. |
20.0 m |
|
|
Steady or transient flow regime used by the benchmark. |
transient |
|
|
Sink and source families activated in the benchmark. |
[] |
|
|
Boundary-condition families activated in the benchmark. |
[west_side, east_side] |
|
|
Hydraulic parameter families explicitly configured by the benchmark. |
[K, Ss, Sy] |
|
|
Parameterization mode used for K. |
homogeneous |
|
|
Homogeneous K value used by the benchmark. |
1e-4 m/s |
|
|
Parameterization mode used for Ss. |
homogeneous |
|
|
Homogeneous Ss value used by the benchmark. |
1e-10 m-1 |
|
|
Parameterization mode used for Sy. |
homogeneous |
|
|
Homogeneous Sy value used by the benchmark. |
0.10 - |
|
|
Initial-condition policy used to start the benchmark. |
custom |
|
|
Scalar initial head or top offset used by the initial-condition policy. |
10.0 m |
|
|
Case-specific configuration field flow.bc.dirichlet.west_side.forcing.mode used by the validation benchmark. |
csv |
|
|
Case-specific configuration field flow.bc.dirichlet.west_side.forcing.path_file used by the validation benchmark. |
west_head_piecewise.csv |
|
|
Case-specific configuration field flow.bc.dirichlet.west_side.forcing.aggregate used by the validation benchmark. |
last |
|
|
Fixed head applied on the east side boundary. |
10.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. |
50 |
|
|
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 |
|
|
Whether the first time period is treated as steady by MODFLOW-NWT. |
false |
|
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. |
50 |
|
|
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 |
|
|
Whether the first time period is treated as steady by MODFLOW 6. |
false |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Solver-specific override applied to MODFLOW 6. |
false |
|
|
Linear-solver complexity preset used by MODFLOW 6. |
COMPLEX |
|
|
Vertical anisotropy ratio passed to MODFLOW 6. |
1 |
|
|
Number of vertical layers used by MODFLOW 6. |
1 |
|
|
Whether the first time period is treated as steady by MODFLOW 6. |
false |
|
|
Committed unstructured mesh file used by the irregular-mesh solver variant. |
../../../shared/mesh_bundles/linearized_unconfined_drainage_irregular_tri_100x10/mesh_2d.msh |
|
|
Committed mesh-bundle directory used to recover support metadata for the irregular-mesh solver variant. |
../../../shared/mesh_bundles/linearized_unconfined_drainage_irregular_tri_100x10 |
|
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 number of stored time periods checked by the validation helper. |
24 |
|
|
Expected spatial shape for each stored time step. |
[5, 50] |
|
Acceptance Criteria by Solver#
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected output: 24 periods, spatial shape 5 x 50 |
|
|
Maximum accepted root-mean-square error for space time. |
0.05 |
|
|
Maximum accepted absolute error for space time. |
0.12 |
|
|
Maximum accepted cross-row spread for space time. |
2e-05 |
|
|
Maximum accepted root-mean-square error for final profile. |
0.05 |
|
|
Maximum accepted absolute error for final profile. |
0.1 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected output: 24 periods, spatial shape 5 x 50 |
|
|
Maximum accepted root-mean-square error for space time. |
0.008 |
|
|
Maximum accepted absolute error for space time. |
0.04 |
|
|
Maximum accepted cross-row spread for space time. |
1e-05 |
|
|
Maximum accepted root-mean-square error for final profile. |
0.003 |
|
|
Maximum accepted absolute error for final profile. |
0.005 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected output: 24 periods, spatial shape 5 x 50 |
|
|
Maximum accepted root-mean-square error for space time. |
0.05 |
|
|
Maximum accepted absolute error for space time. |
0.12 |
|
|
Maximum accepted cross-row spread for space time. |
0.012 |
|
|
Maximum accepted root-mean-square error for final profile. |
0.05 |
|
|
Maximum accepted absolute error for final profile. |
0.1 |
|
Field |
Meaning |
Value |
Source |
|---|---|---|---|
|
Expected output shape or time-space layout checked for this solver. |
Expected output: 24 periods, spatial shape 5 x 50 |
|
|
Maximum accepted root-mean-square error for space time. |
0.05 |
|
|
Maximum accepted absolute error for space time. |
0.12 |
|
|
Maximum accepted cross-row spread for space time. |
2e-05 |
|
|
Maximum accepted root-mean-square error for final profile. |
0.05 |
|
|
Maximum accepted absolute error for final profile. |
0.1 |
|
Source Pointers#
validation_cases/README.mdvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/README.mdvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/reference.pyvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/comparison.pyvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/plotting.pyvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/run_case.pyvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/metadata.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_modflownwt.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/runtime_boussinesq.pyvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances_modflow6.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/tolerances_modflow6_irregular_tri.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_modflow6.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_modflow6_irregular_tri.tomlvalidation_cases/analytical/transient/linearized_unconfined_boundary_piecewise_1d/config_boussinesq.tomlvalidation_cases/shared/boussinesq_uniform_strip.pyvalidation_cases/analytical/transient/linearized_unconfined_1d.py
Artifacts#
docs/source/_static/capability_gallery/validation/linearized_unconfined_boundary_piecewise_1d__modflownwt.pngdocs/source/_static/capability_gallery/validation/linearized_unconfined_boundary_piecewise_1d__modflow6.pngdocs/source/_static/capability_gallery/validation/linearized_unconfined_boundary_piecewise_1d__modflow6_irregular_tri.pngdocs/source/_static/capability_gallery/validation/linearized_unconfined_boundary_piecewise_1d__boussinesq.pngdocs/source/_static/capability_gallery/validation/linearized_unconfined_boundary_piecewise_1d_summary.jsonstores the displayed metrics plus source hashes used bypython -m tools.doc_gallery --check.