I need to compute the harmonic extension u2 of a field u1 (defined in a domain Omega_1) on a domain Omega_2, sharing a common portion of the boundary. By calling u1.set_allow_extrapolation(True) and setting the BC as DirichletBC(W, u1, bdry()) it perfectly works (in serial).
However, if I run the code in parallel, even if the code end without errors, the results are bad: the Dirichlet condition is wrongly assigned on a portion of the boundary.
I’ve found old posts online saying that interpolation between different meshes doesn’t work in parallel. Has this been fixed in the latest version? Is there any workaround?
I guess that doing ALE for FSI the same problem shows up. How is it solved in that case?