A high fluidized bed is modeled using MFIX-DEM. In order to test the case, the lower zone is computed first, and it gets stable results. However, when the computational domain increases to higher zone, computation always goes to divergence. Is it a meshing problem or other? How to deal with it? Thanks. This is the mfx document. FB1.mfx (13.2 KB)