implement virtual mass force

Submitted by gelinhan on Mon, 01/15/2018 - 01:56

Hello developers,

I am simulating a liquid-solid flow in which the virtual mass force is significant. The equation of the virtual mass force is like f = 0.5*rhof*Vp*(Duf/Dt-dup/dt). I looked into the virtual mass force model implemented in CFDEM and found that only the material derivative of fluid phase (Duf/Dt) is included. The acceleration of particle (dus/dt) is missing. Is there any simple method to calculate the dup/dt term?

Thank you.
Linhan

Forums:

paul | Wed, 01/17/2018 - 23:16

The present implementation is a bit clumsy in that regard, introducing inaccuracies even when particles cross processor boundaries.

May I suggest using the total force acting on a particle to calculate du/dt? You could have that communicated from LIG to CFDEM where it could be used by your new force model.

- Paul

gelinhan | Thu, 01/18/2018 - 06:11

Thanks for the idea. I am using the force communicated from LIGGGHTS now. It is working fine. I am doing some further test with this implementation.

Thanks again for the advice.

Linhan

carl | Wed, 03/09/2022 - 21:43

Hi Linhan,

I am interested in your code of the virtual mass force, is it available for sharing?

- Carl

Carl Henrik Dahmén