Skip to content

GitLab

The TS-GitLab will have to shut down towards the end of the year — please think about migrating your projects to GitLab.com or GitHub.
(This is still a very early message, meant to keep you informed. There will be more urgent ones in the future.)

Closed
Created by Lukas Riedel@lriedelOwner

Add Outflow boundary condition for Richards model

Description

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.

Proposal

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).

  • Add OutflowBoundaryCondition to the permitted BC types of the Richards model.
  • Check for Outflow boundary 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.

Related issues

See #

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information