Hello, everyone!
If we want to output the drag force on spherical particles, call drag_gs_des1.f file by default. This file calls the des_drag_gp.f file subroutines of drag_gp_des.f file, and this subroutine calls the drag_gs.f file subroutines to obtain the drag coefficient F_GP of the particle set for the drag model, then calculate the drag force D_FORCE.
But in the SQP method, is the output drag force the same? Output D_FORCE in drag_gs_des1.f file. The drag coefficient of non-spherical particles is calculated in sq_drag.f file instead of drag_gs.f file. So I want to ask software developers, is the traction value output by this method correct? If not, I would like to consult software developers on how to correctly output the drag force acting on non-spherical particles.
Interestingly, the collision force acting on particles is much greater than the drag force. The following compressed package contains examples and code. Thank you all for your help！
0.212.zip (27.1 MB)
Even though the drag coefficient for non-spherical particles is different (different file) from that of spherical particles, this does not affect the calculation of drag force (same file).
1 Like
Thank you for your reply!