I am trying so solve diffusion model using DG method. But my solution “u” between interfaces from sub domains are discontinuous. u(“+”) = function(u(“-”)) or u(“-”) = function(u(“+”)). For example u(“+”) = 10*(u(“-”)). I don't know how can I impose jump condition on the system. Is it possible?
I marked the line between domains by 3.
u(“+”) = 10*(u(“-”)) in interior line 3.
Follow DG formulation:
a = dot(D*S*grad(v), grad(u))*dx \
- dot(avg(grad(v)), jump(u, n))*dS \
- dot(jump(v, n), avg(grad(u)))*dS \
+ alpha/h_avg*dot(jump(v, n), jump(u, n))*dS \
- dot(D*S*grad(v), u*n)*ds(1) - dot(v*n, D*S*grad(u))*ds(1)\
- dot(D*S*grad(v), u*n)*ds(2) - dot(v*n, D*S*grad(u))*ds(2)\
+ (gamma/h)*v*u*ds(1) + (gamma/h)*v*u*ds(2)\
L = v*f*dx- u0*dot(D*S*grad(v), n)*ds(1) + (gamma/h)*u0*v*ds(1) \
- u1*dot(D*S*grad(v), n)*ds(2) + (gamma/h)*u1*v*ds(2) \
How can i impose this jump condition?
