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.
We present a high-order nodal spectral element method for the two-dimensional simulation of nonlinear water waves. The model is based on the mixed Eulerian–Lagrangian (MEL) method. Wave interaction with fixed truncated structures is handled using unstructured meshes consisting of high-order iso-parametric quadrilateral/triangular elements to represent the body surfaces as well as the free surface elevation. A numerical eigenvalue analysis highlights that using a thin top layer of quadrilateral elements circumvents the general instability problem associated with the use of asymmetric mesh topology.We demonstrate how to obtain a robust MEL scheme for highly nonlinear waves using an efficient combination of (i) global L2 projection without quadrature errors, (ii) mild modal filtering and (iii) a combination of local and global re-meshing techniques. Numerical experiments for strongly nonlinear waves are presented. The experiments demonstrate that the spectral element model provides excellent accuracy in prediction of nonlinear and dispersive wave propagation. The model is also shown to accurately capture the interaction between solitary waves and fixed submerged and surface-piercing bodies. The wave motion and the wave-induced loads compare well to experimental and computational results from the literature.
An adaptive spectral/hp discontinuous Galerkin method for the two-dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non-conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p-1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h-type refinement, the parent element is subdivided into four similar sibling elements. The time-stepping is performed using a third-order Runge-Kutta scheme. The performance of the hp-adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p-adaptivity is more efficient than h-adaptivity with respect to degrees of freedom and computational time.