A question about IBM used in Xcompact3D

[yzl429]

Dear developer of Xcompact3d,
It seems that the subroutine momentum-rhseq (dux1, duy1, duz1, rho1, ux1, uy1, uz1, ep1, phi1, divu3) of transeq.f90 does not take into account the forces exerted by solids on the flow field when IBM is turned on. I have a question, when using cubic spline interpolation, is it necessary to include the force of the solid on the fluid in the source term f? Or maybe this force is written elsewhere in src.
Looking forward to receiving your reply.

[mathrack]

Hello,

Each time a derivative is computed, the field is extended inside the solid. The extension will be compatible with the boundary condition at the fluid-solid interface. The resulting convection / diffusion should include the effect of the solid on the fluid. The pressure solver also takes into account the existence of the solid.

See https://github.com/xcompact3d/Incompact3d/blob/master/src/derive.f90

For more information, see https://theses.fr/2021POIT2302 and the references therein.

Please provide direct links to the source code instead of snapshots.

Best regards,
Mathrack