Simulating the free decay motion and wave radiation from a heaving semi-submerged sphere poses significant computational challenges due to its three-dimensional complexity. By leveraging axisymmetry, we reduce the problem to a two-dimensional simulation, significantly decreasing computational demands while maintaining accuracy. In this paper, we exploit axisymmetry to perform a large ensemble of Computational Fluid Dynamics (CFDs) simulations, aiming to evaluate and maximize both accuracy and efficiency, using the Reynolds Averaged Navier–Stokes (RANS) solver interFOAM, in the opensource finite volume CFD software OpenFOAM. Validated against highly accurate experimental data, extensive parametric studies are conducted, previously limited by computational constraints, which facilitate the refinement of simulation setups. More than 50 iterations of the same heaving sphere simulation are performed, informing efficient trade-offs between computational cost and accuracy across various simulation parameters and mesh configurations. Ultimately, by employing axisymmetry, this research contributes to the development of more accurate and efficient numerical modeling in ocean engineering.