Particle Drag Force Output in MFiX-DEM

Hello everyone,

I am new to MFIX, and I am trying to use the DEM method to calculate the drag force acting on each particle.

Based on some relevant posts on the forum, I added the following lines at the appropriate location in the `drag_gs_des1.f` file and rebuilt the solver:

DES_USR_VAR(NP,1) = D_FORCE(1)

DES_USR_VAR(NP,2) = D_FORCE(2)

DES_USR_VAR(NP,3) = D_FORCE(3)

However, my fluidized bed contains a large number of particles, so I need to use the DMP mode for parallel computation. Before running the DMP simulation, I have two questions:

1. Are any additional modifications required for DMP parallel computation? For example, do I need to modify any other files related to parallel communication or data output?

2. I noticed that the spatial distribution of the particle drag force exhibits an obvious regional or block-like pattern. It does not appear to be as smooth and continuous as the drag coefficient field. I am wondering whether this is caused by the numerical algorithm, the interpolation method, the mesh resolution, or the way the drag force is calculated and output in MFIX.

Thank you very much for your help.