Hello
I'm solving a fluid dynamics problem for two phase flows. I'm working in c++.
I have to deal with the different domains separately, so it is not enough to define discontinuous coefficients. I should integrate the equations over the different subdomains using dx(1) and dx(2) in the FFC input file where I define forms. This should work, but I also need to impose boundary conditions on the boundary of the subdomains, probably I will also need to do this in weak form using Lagrange multipliers.
My question is if there is a way to define an internal boundary and integrate over it while defining forms, or even integrate over an arbitrary path in the domain to recover for example stresses.
Thank you in advance for any suggestion.
hello
Stefano
