This study reports the results of the numerical investigation of steady-state thermal developing forced convection between parallel plate channels. The Darcy-Forchheimer-Brinkman model is used for momentum equation. Based on the local thermal equilibrium taking into account the viscous dissipation and the heat generation in solid and fluid phases, temperature profiles of the solid and fluid are obtained numerically in a saturated porous medium, with walls being at constant temperature and constant heat flux. The temperature and the Nusselt number profiles are shown in figures for different dimensionless fluid parameters such as Prantdl, Darcy, Forchheimer, and Brinkman number and the effects of thermal parameters such as the dimensionless heat exchange coefficient, thermal conductivity ratio, fluid heat number and solid heat numbers on the length, and the heat transfer of the thermal developing region are revealed by parameter study. Furthermore, four terms contributing to the viscous dissipation are investigated and the comparisons between the effects of combinations of these terms with temperature boundary conditions are highlighted.