Thank you very much! This is a piece of parallel code…
Please find attached the background mesh Th, the solid mesh Ths, together with the code.
I just found that the problem disappeared if only using one processor. However, it is still strange why the first 9 steps are OK for two or more processors, but then suddenly breaks down.
By the way, the method is similar to Prof. Daniele Boffi’s FSI method (two meshes plus Lagrangian multiplier) you implemented. However, I uses FEM interpolation rather than Lagrangian Multiplier. Prof. Olivier Pironneau is familiar with my work, who was my PhD thesis examiner. Of course, I am quite familiar with your monolithic FSI paper (2017) with Olivier. Olivier shared most of his FSI code, and I have implement your monolithic FSI method in parallel these days – but I came across stability issues, will ask you about that later.
Th_cylinder.mesh (639.7 KB)
Ths_cylinder.mesh (373.3 KB)
fsi_3d.edp (5.8 KB)