Hello,
I have a domain that is split in two (using the code by Guosheng:
ngsolve.org/forum/ngspy-forum/114-interface-integrals#491
)
In fluid 1 I define a FacetFESpace Q1 and in fluid 2 I define a FacetFESpace Q2.
Let p1 and q1 be test and trial functions in Q1.
Let p2 and q2 be test and trial functions in Q2.
On the interface between the fluids I compute the integral
100*p1*q1 + 100*p2*q2.
I then print the matrix and I see that I get different matrices depending on how I define Q2. If I define
Q2 = FacetFESpace(mesh, order=1)
I will get, for example, that A(54,54) = 16.6667, but if I define
Q2 = FacetFESpace(mesh, order=1, definedon=mesh.Materials("fluid2"))
I will get that A(54,54) = 0
I would prefer to use the second definition of Q2, but this gives the wrong matrix. Any suggestions how to fix this?
I have attached the code with Q2 defined in lines 52 and 53.
Thanks,
Sander