The current study explored the capability of a discrete particle method known as dissipative particle dynamics with energy conservation (eDPD) to simulate combined convection heat transfer in a vertical lid driven cavity. The study investigated two cases of aiding and opposing buoyancy mechanisms in the lid driven cavity. The eDPD results were compared against the finite volume solutions for the range of Richardson number, 10 −2 ≤Ri≤10 2 . The method showed good comparison for the range of Richardson number 10 −2 ≤Ri≤10 1 . However, the eDPD method showed deviation from the FV solutions for a high value of Richardson number, Ri=10 2 , and this deviation is attributed to the compressibility of eDPD system experienced at such high value of Richardson number. Parametric study on the influence of the Richardson number (Ri) on the eDPD compressibility was conducted and presented via temperature isotherms, streamlines, velocity contours, velocity vectors, temperature and velocity profiles.