This paper presents an analysis of forced convection heat transfer in microchannel heat sinks for electronic system cooling. In view of the small dimensions of the microstructures, the microchannel is modeled as a fluid-saturated porous medium. Numerical solutions are obtained based on the Forchheimer-Brinkman-extended Darcy equation for the fluid flow and the two-equation model for heat transfer between the solid and fluid phases. The velocity field in the microchannel is first solved by a finite-difference scheme, and then the energy equations governing the solid and fluid phases are solved simultaneously for the temperature distributions. Also, analytical expressions for the velocity and temperature profiles are presented for a simpler flow model, i.e., the Brinkman-extended Darcy model. This work attempts to perform a systematic study on the effects of major parameters on the flow and heat transfer characteristics of forced convection in the microchannel heat sink. The governing parameters of engineering importance include the channel aspect ratio (alpha(s)), inertial force parameter (Gamma), porosity (epsilon), and the effective thermal conductivity ratio (k(Gamma)). The velocity profiles of the fluid in the microchannel, the temperature distributions of the solid and fluid phases, and the overall Nusselt number are illustrated for various values of the problem parameters. It is found that the fluid inertia force alters noticeably the dimensionless velocity distribution and the fluid temperature distribution, while the solid temperature distribution is almost insensitive to the fluid inertia. Moreover, the overall Nusselt number increases with increasing the values of alpha(s) and epsilon, while it decreases with increasing k(Gamma). (c) 2006 Elsevier Ltd. All rights reserved.