Commit 4b3029a5 authored by Lukas Riedel's avatar Lukas Riedel

Add outflow BC to Richards FV local operator

parent 59d7b337
#ifndef DUNE_DORIE_RICHARDS_OPERATOR_FV_HH
#define DUNE_DORIE_RICHARDS_OPERATOR_FV_HH
#include <dune/common/float_cmp.hh>
#include <dune/geometry/referenceelements.hh>
#include <dune/pdelab/common/quadraturerules.hh>
......@@ -256,6 +258,11 @@ public:
// Water flux in normal direction w.r.t the intersection
const auto water_flux_n = - cond_i*dudn;
// No contribution for inward flux
if (BoundaryCondition::isOutflow(bc)
and Dune::FloatCmp::lt(double(water_flux_n), 0.0))
return;
// Contribution to residual from Dirichlet boundary
r_i.accumulate(lfsv_i, 0, water_flux_n*volume_f);
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment