Transitional free surface flow profiles with a critical point occur with weak vorticity and viscosity effects and thus can be modeled with an irrotational flow approach. Important examples in hydraulic engineering include flow over low obstacles, transition structures in canals, and flow over high spillways. While solving Laplace’s equation is relatively simple, the determination of the unknown free surface and energy head of the flow is challenging. Both hydraulic quantities need to be iterated before solving the Laplace equation. Former models iterated the energy head on a trial-and-error basis, assuming that the linked free surface profile is smooth, i.e., free of waves. The iteration of the free surface for a given head is frequently accomplished using the Newton–Rapshon method, which is difficult for the challenging case of spillway flow, giving no solution in some cases. An alternative method of computing irrotational flow profiles in transitional flows involving a critical point is proposed in this work. The model contains three elements: mapping of Laplace’s equation to directly track the streamlines, determination of the critical point and unknown energy head using a critical flow condition for irrotational flows, and determination of the water surface position using an exact analytical solution. The proposed model is favorably compared with experimental data from different sources and CFD results, indicating a reasonable agreement.