hydromodpy.solver.boussinesq.runtimes.petsc_mixed#

PETSc runtime for the Boussinesq backend on Linux.

This backend switches from the historical head-only regularized overflow law to one semi-explicit mixed double-obstacle formulation:

  • h remains the differential unknown,

  • q_ex is one algebraic unknown per cell for saturation excess,

  • q_dry is one algebraic unknown per cell for the lower obstacle.

Time integration still uses one backward-Euler step per stress period, but each step now solves the coupled nonlinear DAE residual with PETSc SNES.

Functions

solve_steady_problem(inputs)

Solve one steady mixed complementarity system with PETSc SNES.

solve_transient_step(inputs)

Solve one implicit transient DAE step with PETSc SNES.