Hi everybody, and developers of CFDEM,
Since the difference in hydro-static pressure ( p = ρgh) on the particle surface causes the buoyancy force, and the fluid force acting on the surface of immersed particle is computed in a fully resolved manner by integrating the pressure and viscous stress on the fluid-solid interface:
F = integral_S { -p*I + Tau } dS (on solid boundary)
= integral_V { -grad(p) + mu * laplacian(U) } dv (in solid region)
for in-compressible fluid.
So I wonder, is the buoyancy or Archimedes force already included in the pressure gradient term " -grad(p) "? Or it is just because the cfdemSolverIB is based on the openFOAM solver pisoFoam, in which the gravity force is ignored.
best,
David