Incorporating particle-level correlations for heat/mass transfer and thermal buoyancy based on Boussenisq approximation into two-way coupled direct numerical simulations (DNS) (an advanced CFD technique) to study the fluid and particle phase dynamics.