Fenics：内部边界条件

Question

我想把斯托克斯方程和泊松方程（达西）与Beaver-Jospeh界面条件结合起来。我在单位正方形中创建了两个子域，一个来自[0，0.5]，另一个来自[0.5，1]。 在内部边界处，我要使用Beaver-Jospeh条件：

du/dy(x,0+)=c(u(x,0+)-u_D)

式中，u？D表示多孔区域中的达西速度，u（x，0+）表示流体域中的斯托克斯速度。c是常数。在

现在我不知道如何创建内部边界条件。 我发现了这样的东西：

^{pr2}$

但我想我必须使用内部方面？ 有谁能帮我做这个吗？我不知道怎么开始。在

谢谢！在