Hi, all
When DG methods are employed in the Stokes equation, we tend to use some codes like:
mesh = Mesh(geo.GenerateMesh(maxh=0.1, quad_dominated=False))
V = VectorL2(mesh, order=3, dgjumps=True)
n = specialcf.normal(mesh.dim)
u, v = V.TnT()
mean_grad_u = 0.5 * (grad(u)+grad(u.Other()))
jump_v = v-v.Other()
a = BilinearForm(V)
a += InnerProduct(mean_grad_u.trans*n, jump_v) * dx(skeleton=True)
My question is: for a fixed edge or facet, how do we know its inner element or outer element because we use jump_v = v-v.Other()
, and further the direction of its unit normal? Are the both consistent?
Best,
Di Yang