“…We note that (i) Kriebitzsch et al (2013b) reported the main difference between the results obtained from both DNS and CFD-DEM method via the bed expansion characteristics, we therefore reported the CFD results of H only; (ii) due to the restitution coefficient used is unity, which means there is no energy dissipation during particle-particle and particle-wall collisions, therefore, there is no bubbling structure formed (the results are not shown here for brevity), although the superficial gas velocity is roughly two times of the minimum fluidization velocity, this is in agreement with the conclusion made from both DNS and CFD-DEM method ; (iii) The value of H is calculated as in present study, where is the number of time step used to calculated the mean value ( , meaning 10 seconds have been used), is the total number of particles used in simulations, N is the number of computational cells, is the volume of cell i, is the volume of particle, is the height of the center of cell i and is the solid concentration. Further simulations have been carried out to compare with the DNS results of Third et al (Third et al, 2016), the relevant parameters are summarized in table 2. For a gas velocity of 0.3 m/s, the values of H calculated using equations (7) and (10) are 15.8 mm and 15.5 mm, respectively, both of which are in a good agreement with the DNS result (H=15.7 mm); For a gas velocity of 0.24 m/s, the values of H calculated using equations (7) and (10) are 14.6 mm and 13.2 mm, respectively, the former is very close to the DNS result (H=14.4 mm), whereas the latter is much less than the DNS result.…”