Infragravity waves, which have the frequency of 0.004–0.04 Hz, are generated indirectly by the wind and they obtain energy from sea-swell short waves through nonlinear energy transfer. The better understanding of infragravity wave dynamics can help to provide insights into coastal flooding and erosion. Accurate estimation of infragravity wave related processes requires the efficient and reliable numerical modeling approach. In the context of RANS-based simulations of large-scale wave processes, nonhydrostatic modeling is balanced in efficiency and accuracy in which the free surface is simplified with a single-valued function of the horizontal coordinate. The open-source NHWAVE model is developed to account for a second-order correction that suppresses the spurious harmonics when a first-order wave generation method is used. The model is validated with the GLOBEX dataset. The adopted method can generate bound infragravity wave at the incident boundary and can reproduce infragravity wave height evolution in the shoaling and surf zone.