hydromodpy.solver.boussinesq.runtimes.petsc_ts_vi_obstacle

Experimental PETSc TS + SNESVI runtime for a head-only obstacle problem.

The PETSc state vector contains only hydraulic head h. PETSc SNESVI bounds that state with z_bottom <= h <= z_top by default and PETSc TS performs fixed Backward-Euler steps inside the HydroModPy stress period. Surface and bottom rates are reconstructed after convergence from the remaining implicit residual; q_ex and q_dry are not primary unknowns and no Fischer-Burmeister equations are assembled here.

When an explicit positive Cauchy drainage conductance is supplied, the upper obstacle is relaxed so the drainage flux term carries the top boundary exchange.

Functions

solve_steady_problem(inputs)	The first TS VI prototype is transient-only.
solve_transient_step(inputs)	Solve one stress period using PETSc TS Backward Euler and SNESVI.