Objective: Incorporating particle-level correlations for heat/mass transfer and thermal buoyancy based on Boussenisq approximation into two-way coupled direct numerical simulations (DNS) (an advance CFD technique) for particle-laden flows to study the turbulence modulation and particle clustering.