How to simulate a spouted bed with an orifice using diffusion + Gaussian kernel method in MFIX?

Hi, everyone,
I want to simulate a spouted bed with an orifice. The MFIX version I used is 19.2.2. The dp = 3 mm, the length of orfice is 5 mm. The orifice is divided into 2 grids, and then the dx/dp = 0.84. The diffusion equation with a Gaussian kernel distribution was used in my case. However, the case runs well on my computer with a single CPU while it crushes with parallel computation (4 CPU or 8 CPU). Is it the bug of the implementation of diffusion equation + Gaussian kernel in MFIX under parallelization simulation?

The bug occurs when the case is run under parallelization using 8 processors:

The file is also attached.
SFB_test.zip (39.0 KB)