Accurate numerical prediction of boiling phenomena is critical for many thermal applications due to its high heat transfer capacity and constant temperature nature. Particularly, in nucleate boiling, vapor bubble formation significantly improves heat dissipation. The surface conditions of the boiling site strongly influence nucleate boiling. Therefore, to enhance heat transfer, the formation of microscale cavities on the boiling surface is a viable approach. The present study numerically investigates the effect of cavity formation on nucleate boiling heat transfer, with a central focus on a rectangular microscale cavity geometry. The Finite Volume Method (FVM) is used to model nucleate boiling with and without rectangular cavities. Rohsenow’s boiling correlation was implemented with the Volume of Fluid (VOF) Boiling Method. The geometries are modeled in 2D, and the working liquid is selected to be liquid nitrogen. The liquid temperature is set to the evaporation temperature of nitrogen in atmospheric conditions and the wall temperature is kept constant. The nucleate boiling characteristic is examined with transient simulations and the results are compared based on key parameters such as the heat transfer rate and the mass flow rate of the boiling gas.