How to apply boundary conditions inside the domain?

7 weeks ago by
Equation and domain are shown in the attached picture. If it is possible to apply the first and second types of boundary conditions inside the domain (points or lines as shown in the picture.)? If yes, what are the ways to do it?

Thank you very much!

Community: FEniCS Project
It should be possible to really just define the inner subdomains the way you would define a boundary subdomain. And then it's DirichletBC() for the first or some ds integral term in your weak form for the second condition. Have you already tried it out this way?
written 7 weeks ago by klunkean  
Thank you very much for your help! This way works only when we apply to a line (interface); however, I have no idea how to apply to a point.
written 7 weeks ago by Teeratorn Kadeethum  
For DirichletBCs you can add the keyword argument method = "pointwise".
Defining a flux at a point does not make immediate sense in my opinion, as in the context of weak forms in n dimensions a flux is always integrated over an n-1 dimensional object, e.g. surfaces in a 3D problem.
written 7 weeks ago by klunkean  
Thank you! I agree with your comments.
written 7 weeks ago by Teeratorn Kadeethum  
Please login to add an answer/comment or follow this question.

Similar posts:
Search »