Setting external body force in Mfix-DEM

MFiX
,
1

Hello,

The simulation consists of generated particles falling due to gravity.
X-cyclic condition enabled with set outlet velocity.

The issue is:

  1. Any input values for particle-particle rolling friction is not giving the right velocity contour

  2. Outlet velocity should be controlled by body force. The body force is added to gradually reduce the deviation between outlet velocity and set outlet velocity of 0.06m/s. The file .f needs to be reviewed.

Formulation:

Fx,add = ms (g+(vx- v_out)/ delta_t) - Fx

where ms = mass of particle, g = gravity, vx = particle velocity in x direction, v_out = particle outlet velocity, t = time, Fx = force in x direction

I have attached the below .f file for setting up external body force.

The particle velocity must be updated according to the force on the particle at each time step.

Outlet velocity is set as 0.06m/s.

Kindly help me fix these issues also review the .f file

Thanks

pinfin_3mm_9.97_0.06_walls.zip (37.9 MB)

bodyforce_mod.f (2.6 KB)