A detailed mathematical model of a two-phase heat spreader with axial microchannels is developed in which the fluid flow is considered along with the heat and mass transfer processes during evaporation and condensation. The model is based on the equations for the mass, momentum and energy conservation, which are written for the evaporator, adiabatic, and condenser zones. The model, which permits to simulate several shapes of microchannels, can predict the maximum heat transfer capacity of the two-phase heat spreader, the optimal fluid mass, and the temperatures and pressure gradients along the microchannel. The effect of shear stresses at the free liquid surface in a microchannel due to the frictional liquid-vapor interaction on the liquid flow is taken into consideration. The heat transfer through the liquid films in both evaporator and condenser is accounted for in the model, which is described with respect to the disjoining pressure, interfacial thermal resistance, surface roughness, and curvature. The thermal resistances of the evaporator and condenser are determined by accounting for the longitudinal distribution of the meniscus curvature, which is dependent on heat load and heat spreader inclination. NOMENCLATURE A Section of liquid passage (m²), Hamaker constant (J) A' Dispersion constant (J) d Width of the microchannel (m) Dh Hydraulic diameter (m) ep Thickness of the wall (m) f Friction coefficient (-) g Gravity acceleration (m/s²) hs Microchannel height (m) hc Condenser heat transfer coefficient (W/m².K) hey Evaporator heat transfer coefficient (W/m².K) K Curvature of liquid-vapor interface (m 1 ) La Length of the adiabatic zone (m) Las Length of the dried zone (m) Lb Length of the blocked zone (m) Lc Length of the condenser (m) Le Length of the evaporator (m) Lt Total length of the heat spreader (m) M Mass (kg) e m Evaporated mass flow rate per unit surface (kg/s.m²) c m Condensed mass flow rate per unit surface (kg/s.m²)