Commit 9e1b498f authored by Kilian Weishaupt's avatar Kilian Weishaupt

[coupling-ff-pm] Improve docu in ff problem

parent c2c04c9c
Pipeline #365 passed with stage
in 47 seconds
......@@ -135,8 +135,15 @@ public:
BoundaryTypes boundaryTypes(const Element &element, const SubControlVolumeFace &scvf) const
{
BoundaryTypes values;
// set Neumann BCs to all boundaries first
values.setAllNeumann();
// set the coupling boundary condition at the interface
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
values.setAllCouplingNeumann();
// set a Dirichlet boundary condition at the bottom
if (onLowerBoundary_(scvf.center()))
values.setAllDirichlet();
......@@ -153,6 +160,7 @@ public:
*/
PrimaryVariables dirichlet(const Element &element, const SubControlVolumeFace &scvf) const
{
// set p = 0 at the bottom
PrimaryVariables values(0.0);
values = initial(element);
......@@ -175,8 +183,10 @@ public:
const ElementVolumeVariables& elemVolVars,
const SubControlVolumeFace& scvf) const
{
// no-flow everywhere ...
NumEqVector values(0.0);
// ... except at the coupling interface
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
values[Indices::conti0EqIdx] = couplingManager().couplingData().massCouplingCondition(fvGeometry, elemVolVars, scvf);
......
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