In this study, a flat heat pipe with asymmetric boundary conditions is modeled numerically in two-dimensional space. Vapor and wick regions are incorporated in the analysis. The Navier-Stokes equations along with the energy equation are solved under steady state, subsonic, laminar, incompressible and Newtonian flow assumptions. Furthermore, Brinkman-Forchheimer extended Darcy model is employed for the liquid flow in the wick structure. SIMPLE algorithm is employed with staggered grid arrangement to overcome the complications caused by the pressure terms in the momentum equations. The model accounts for uniform injection and suction in the evaporator and condenser sections and it is assumed that evaporation and condensation only occur in evaporator and condenser sections, respectively. Wick structure is presumed to be isotropic and saturated with working fluid. All thermo physical properties are taken as constant and no body forces are included in the model. By considering all of the above information, an in-house code is developed by using the software, MATLAB in order to perform the analysis. The major results introduced are the distribution of velocity components along with their profiles at different locations, temperature distribution and pressure drops along the heat pipe as well as their distributions within the domain. Additionally, parametric analyses are performed to investigate the effects of heat input and heat pipe aspect ratio on the required capillary pressure.