我想将斯托克斯和泊松方程（达西）与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是常数。
现在我不知道如何创建内部边界条件。 我找到了类似的东西：
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)
但我认为我必须使用内部方面？ 有人可以帮我怎么做？我不确切知道如何开始。
谢谢！