Python Fenics：内边界条件

我想结合斯托克斯方程和泊松方程，达西方程和比弗-约瑟夫界面条件。我在Unit Square中工作，创建了两个子域，一个来自[0,0.5]，另一个来自[0.5,1]。 在内部边界处，我想使用Beaver-Jospeh条件：

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

其中u_D表示多孔区域中的达西速度，ux，0+流体区域中的斯托克斯速度。c是一个常数

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

facet_domains = FacetFunction('size_t', mesh)
facet_domains.set_all(0)
for f in facets(mesh):
if any(ff.exterior() for ff in facets(f)):
facet_domains[f] = 1

u_in = Constant(-1.0)
bc_in = DirichletBC(V, u_in, facet_domains, 0)

但我想我必须使用内部面？ 有人能帮我怎么做吗？我不知道该怎么开始

谢谢