Hello,
My question is how one can put a boundary condition in the equation-based environment on an internal boundary. To be more presize: assume there are two domains with a common boundary. The dependent variable is u. Let us assume for simplicity that I am making a simulation of a Laplace equation in the both domains:
Laplace(u)=0
There may be several types of boundary conditions on the internal boundary such as u1=u2 or a*u1=u2
where u1 and u2 are the boundary values on the internal boundary in the domains 1 and 2 correspondingly. a is a constant.
The second type is related to the flux: n1*grad(u1)=n2*grad(u2) or a*n1*grad(u1)=n2*grad(u2)
where n1 and n2 are the unit normal vector to the boundary in the domains and grad stays for gradient.
There should be a trick to implement such types of the boundary conditions.
Please let me know, if there is one.
Alexei.
My question is how one can put a boundary condition in the equation-based environment on an internal boundary. To be more presize: assume there are two domains with a common boundary. The dependent variable is u. Let us assume for simplicity that I am making a simulation of a Laplace equation in the both domains:
Laplace(u)=0
There may be several types of boundary conditions on the internal boundary such as u1=u2 or a*u1=u2
where u1 and u2 are the boundary values on the internal boundary in the domains 1 and 2 correspondingly. a is a constant.
The second type is related to the flux: n1*grad(u1)=n2*grad(u2) or a*n1*grad(u1)=n2*grad(u2)
where n1 and n2 are the unit normal vector to the boundary in the domains and grad stays for gradient.
There should be a trick to implement such types of the boundary conditions.
Please let me know, if there is one.
Alexei.