Fluctuations embedded force in two-phase flow

[Nathaniel]

Hi everyone!

I'm currently trying to combine moving embedded boundaries with CLSVOF (or VOF). First of all, a huge thanks to Ghigo for his work, it's been incredibly helpful!

My current code runs to some extent, but when I check the forces on the embedded boundary, I'm seeing abnormal high-frequency fluctuations, especially after the solid exits the water.

I suspect this might be because the vof interface isn't being properly truncated by the solid boundary. I've tried using tavares's 'contact-embed.h'. For CLSVOF, my workaround was basically adding some simple cs checks in the redistance and related files (though this really isn't my area of expertise...). As shown in the attached figures, they didn't seem to help much.

I'm wondering about the best direction: should I keep focusing on resolving the three-phase truncation issue, or could this actually be a problem with how the embedded force is being calculated?

I would really appreciate any advice or discussion on this!