Natural laminar convective fluid flow has been simulated inside inclined rectangular cavities with and without internal heat generation for different aspect ratios and inclination angles. The most important basic dimensionless parameters for this problem are the external Rayleigh number (RaE) and the internal Rayleigh number (RaI), where RaE refers to the effects of the differential heating of the side walls and RaI refers to the amount of heat produced internally. Results were obtained for 4 cases with 192 tests: case (1), RaI = 0 without internal source generation, and cases (2, 3, and 4) with internal source generation for RaI = RaE, 10 RaE, and 100 RaE, respectively. In all cases, the parameters of study changed as 103 ≤ RaE ≤106, 0 ≤ RaI ≤ 107, inclination angle from 0 to 60 deg., and aspect ratios of the enclosure from 0.5 to 2. Results were represented graphically for flow and thermal fields as a streamline, isothermal contours, and Nusselt number. The computed results show that the strength of convection currents is measured by the internal energy. Finally, it is illustrated that by using a few grid points and a shorter CPU time for calculation, the present method can produce accurate numerical results. Also, increase in RaI leads to increasing heat transfer rate and its direction out from the cavity at both hot and cold walls. For lower values of RaI, heat transfer diffusion is more prominent, while for higher values of RaI, convection outweighs diffusion.
 HIGHLIGHTS
 
 Natural laminar convective fluid flow inside inclined rectangular cavities with and without internal heat generation for different aspect ratios and inclination angles has been simulated
 The most important basic dimensionless parameters, the external Rayleigh number (RaE) and the internal Rayleigh number (RaI) are studied
 DQ method performance was excellent
 The obtained computational results indicate that the strength of the convection currents depends on the internal energy
 Accurate numerical results can be obtained by the present method using a few grid points and shorter CPU time for calculation
 
 GRAPHICAL ABSTRACT