Dear all,
Recently I have read the code of virtual mass force, but there is a question that:
the original force = 0.5*rho_f*V_p*(Du/Dt-du_s/dt)=0.5*rho_f*V_p*(du/dt+u*div(u)-du_s/dt)
is that means in our virtual mass force model we use:
if splitUrelCalculation, using (Du/Dt) instead of (Du/Dt-du_s/dt),
if not splitUreCalculation, using (Ur-Urold)/dt instead of (Du/Dt-du_s/dt).
Is my understanding of this right?
Thank you!
Detian
2015.7.18