Aiming to uncover the propulsion mechanisms underlying a cruising crocodile, we conduct computational fluid dynamics (CFD) simulations on the flow around a simplified three-dimensional model of the Crocodylus siamensis. The locomotion of the crocodile model is realized through undulating its body and tail, mimicking a crocodile-like swimming pattern. At a cruising speed of U∞ = 0.5 m/s (corresponding to a Reynolds number Re = 9.95 × 105 based on U∞ and the body length L), the hydrodynamics of the crocodile model are investigated, taking into account effects of the undulation parameters (i.e., amplitude A and frequency f). The normalized undulation parameters cover broad ranges of 0.6 ≤ A* = A/W ≤ 1.0 and 0.25 ≤ f * = fW/U∞ ≤ 0.625, where W is the body width. The CFD simulations are conducted in ANSYS Fluent, with the SST k–ω turbulence model and user-defined functions for dynamic mesh being used. Numerical results reveal that A* and f * render profound effects on the hydrodynamic performance of the crocodile model. The time-mean axial force coefficient (CA¯) and power coefficient (C¯Power) exhibit rapid growth with increasing A* and/or f *, while the root mean square lateral force coefficient (Cy,rms) is more dependent on f * than on A*. It is further found that, irrespective of A*, CA¯ and C¯Power can be well scaled with Strouhal number St (= 2fA/U∞) or St2(1 − U∞/c). Furthermore, distinct flow patterns are observed in the wake of the crocodile model undulating at different St, corresponding to the drag, transition (or cruising), and thrust type swimming, respectively. Discussion is made on the wake flow structures and their connections to the generation of the hydrodynamic forces. The findings from this work contribute to the understanding of the propulsion mechanisms of the swimming crocodile, meaningful for the design of efficient biomimetic amphibious robots.