Add Outflow boundary condition for Richards model
An Outflow boundary condition is a Dirichlet boundary condition that only allows a net flux in direction of the boundary normal vector, i.e., to the outside. This is particularly handy for simulation evaporative potentials. Depending on the actual potential at the boundary, setting a Dirichlet condition can lead to a flux into or out of the domain. An outflow boundary condition would prevent the inward flux and simply apply a (Neumann) no-flow condition in this case.
This boundary condition has already been implemented in !134 and has proven well in previous simulations.
Add boundary condition type
Outflow to Richards model. The boundary condition stores a matric head value which is only applied if the normal flux on the boundary is positive (i.e., outwards).
OutflowBoundaryConditionto the permitted BC types of the Richards model.
- Check for
Outflowboundary condition in Richards local operators. If the normal flux is negative, prohibit contribution to the local residual.
- Update documentation.
How to test the implementation?
Add a system test where an
Outflow BC is applied and a comparable Dirichlet BC would lead to a flux. Then compare the result with the result of an ODE solver which applies the no-flow BC.