This numerical work presents a geometrical investigation of a corrugated isothermal surface placed in a two-dimensional cavity subjected to unsteady, turbulent pool boiling flows. The main purposes are maximizing the heat transfer rate between the isothermal surface and the surrounding water flow, and the volume of vapor generated into the cavity. The geometric investigation followed the constructal design method, being the ratio Hi/Li (i = 1, 2 or 3) of the corrugations varied for three different numbers of corrugations: N = 1, 2, and 3, keeping constant the corrugations area. The volume of fluid (VOF) and Lee's evaporation-condensation models are used to estimate the volume fractions of water vapor/liquid and interfacial mass transfer. The unsteady Reynolds Averaged Navier Stokes (URANS) continuity, momentum and conservation of energy equations, and volume fraction transport equation, are solved using the finite volume method (FVM) available in software Ansys FLUENT. For closure of turbulence, the k - ɛ model is adopted. For validation of the model, the heat flux and convection heat transfer coefficient obtained for a pool boiling bared surface case are compared with Rohsenow's correlations, and differences lower than 7.0 % are reached. Results indicated a strong influence of the ratio Hi/Li and number of corrugations (N) in the heat transfer rate per unit depth (qs) and dimensionless volume of vapor (Vdim) generated into the cavity. The highest intrusion of the corrugations led to the generation of few large and many small scales, benefiting the thermal performance, regardless of the performance indicator employed. The optimal configuration, N = 3 and H3/L3 = 2.0 improved 49 % and 188 % the Vdim and qs compared with the worst corrugated case, showing the importance of the geometry of the corrugation in this problem.