Apply BC only on part of a sideset (spatial activation/deactivation)

[josebastiase]

Check these boxes if you have followed the posting rules.

• Q&A General is the most appropriate section for my question
• I have consulted the posting Guidelines on the Discussions front page
• I have searched the Discussions forum and my question has not been asked before
• I have searched the MOOSE website and the documentation does not answer my question
• I have formatted my post following the posting guidelines (screenshots as a last resort, triple back quotes around pasted text)

Question

Hi,

I’m working on an environmental problem where a boundary condition is defined on a sideset, but it should only act on certain spatial portions of that boundary. The active region can change in space and time and is not known a priori.

When the BC is not acting, it should be fully turned off (not just set to zero).

Splitting the boundary into many small predefined sidesets is not practical in this case.

What is the recommended way in MOOSE to activate/deactivate a BC only on part of a sideset based on spatial coordinates (and possibly time)?

Thanks!

[GiudGiud]

Hello

is this an integratedBC (creating an integral contribution) or a dirichlet BC?

[josebastiase]

dirichlet BC

[josebastiase]

well... actually is a linear sink

[GiudGiud]

right so it derives from PorousFlowSink which is an IntegratedBC not a dirichletBC.
So if is trying to enfore a dirichlet condition it is either doing it weakly or using a penalty method (which is weakly in a way).

But that is good news.
To not enforce an integratedBC, all you need to do is to have it contribute 0 to the integral of the BC.
So in the compteQpIntegral, have the desired logic for "should-contribute-or-not" and then return 0 if it should not be active

[josebastiase]

Hi @GiudGiud that worked! I used a penalty function to replace a dirichlet condition and works as expected. Thank you