Two methods for generating nonlinear solitary waves on constant water depth in order to study their further runup over a steep sloping beach are presented. The first method is based on a prescribed wave paddle motion inducing a desired free surface profile. The second generation procedure is based on an appropriate mass source term added to the mass conservation equation in an internal source line. These two wave generation methods are implemented in the industrial Computational Fluid Dynamics code (CFD): PHOENICS (Parabolic Hyperbolic Or Elliptic Numerical Integration Code Series). The generated solitary wave's profile and their non-breaking runup over an impermeable steep beach are in good agreement with experimental data of Lin, Chang and Liu (1999). No appreciable discrepancies are observed between experimental and simulated velocities during the runup and rundown processes except for the horizontal velocity component at the final stage of the second runup. The moving paddle and mass source methods performed equally well on reproducing non-breaking runup despite the difference in their wave generation procedures.