Water electrolysis in neutral medium is a promising avenue to reduce the overall cost of industrial hydrogen production. However, the complex interplays among various processes, including electrochemical reaction, mass transport, bubble detachment, two-phase flow, and homogeneous reaction of the buffer, result in a high nonlinearity in solving the multiphysics of hydrogen evolution reaction (HER). In this paper, a two-phase flow and mass transport model is developed for HER under nearneutral pH condition. The mathematical model is solved to obtain the velocity and concentration profiles at various reaction rates, buffer concentrations, gravitational accelerations, and electrode properties. The results indicate that the homogeneous reaction of the buffer remarkably facilitates the transport of H + and OH − , thereby prompting the HER performance. The porosity plays a major role in facilitating ion transport inside/outside the electrode, compared to superaerophobic design and supergravity. In addition, the phosphate buffer exhibits a better feasibility and capacity than the bicarbonate buffer in near-neutral pH medium. This paper concludes that the phosphate/bicarbonate buffer and two-phase flow concurrently contribute to the improvement of ion transport and HER performance. The developed mathematical model provides a deeper understanding of the fundamental processes of HER in near-neutral pH medium.