Drag force: ShirgaonkarIB
Hi every one,
I tried to run a settling simulation (similar to the settlingTestMPI) with the IB-Solver and found that the particle doesnt settle as it should be (Stokes as comparison). So I had a look at the calculation of the drag force and I think there is a mistake in the equations.
(1) volVectorField h = nufField*fvc::laplacian(U_)-fvc::grad(p_);
...
(2) drag += h[cellI]*h.mesh().V()[cellI];
I think both terms in equation (1) (nufField*fvc::laplacian(U_) and fvc::grad(p_)) should be multiplied with the density of the fluid.
Forums:
- Read more about Drag force: ShirgaonkarIB
- 1 comment
- Log in or register to post comments