Force balance

Submitted by Yingzhi on Fri, 11/27/2020 - 01:20

Hi everyone

I try to use dump custom fx fy fz f_dragforce(1) f_dragforce(2) f_dragforce(3) to record the total forces and fluid drag forces on each particle.

But the output is so wired:ITEM:
ATOMS id type type x y z ix iy iz vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] omegax omegay omegaz radius
1 1 1 24.9037 89.5883 2300.29 0 0 0 -6.94822e-07 3.52545e-06 -8.67465e-05 -7.7927e-20 3.06414e-19 9.97466e-18 -1.14887e-06 -1.62293e-06 1.65673e-06 0 0 0 1.5
2 1 1 50.4235 200.241 2292.43 0 0 0 7.8989e-07 5.26336e-06 -8.79128e-05 7.68815e-20 5.86253e-19 -9.97466e-18 8.78281e-08 -1.62379e-07 1.07349e-06 0 0 0 1.5

I do not understand why there are so huge differences between f[] and f_dragforce. There is no contact force because the particle is Particles suspended in the fluid.

Could someone please help me solve or explain this discrepancy?
Thank you very much!