Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
D
dorie
  • Project overview
    • Project overview
    • Details
    • Activity
    • Releases
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 31
    • Issues 31
    • List
    • Boards
    • Labels
    • Service Desk
    • Milestones
  • Merge Requests 9
    • Merge Requests 9
  • Operations
    • Operations
    • Incidents
    • Environments
  • Analytics
    • Analytics
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Members
    • Members
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • dorie
  • dorie
  • Issues
  • #182

Closed
Open
Opened Mar 24, 2020 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
Assignee
Assign to
None
Milestone
None
Assign milestone
Time tracking
None
Due date
None
Reference: dorie/dorie#182