I have developed a parallel matrix-free monolithic FSI solver in the ALE frame of reference. Exemplary results
[1,
2], more details in my PhD thesis
[3].
Highlights:
-> The proposed monolithic predictor-corrector time integration scheme requires solving the generalized Stokes problem with discontinuous variable coefficients. This has to be done one or two times per time step, depending on the variant of the scheme.
-> I have introduced a new multilevel preconditioner, for a 3D case the time required to solve a system consisting of 421 million unknowns using 96 cores is 178 seconds.
-> As demonstrated by the movie[1], I am able to solve 3D problems at moderate Reynolds numbers (time step size of 0.01 did not result in any stability issues with Re=2k).
->The preconditioner from the thesis can still be greatly improved. I have a proof of concept demonstrating the basic properties of a new idea of smoother. The results look very promising and I'm curious how it would work. Moreover, the algorithm is suitable for GPU that could result in further speedup.
[1] 3D "Turek benchmark" at Re = 2000 (30M DoF)
https://youtu.be/7TaC1dTfAQY[2] Water hammer in elastic tube (12M DoF)
https://youtu.be/TRnlt_eC32QThose results were obtained as a part of my PhD thesis. I have finally managed to submit it in January, the defence will be held in 3 weeks (expect reviews next week).