Search

International audience; This paper studies if adaptive grid refinement combined with finite-volume simulation of the incompressible RANS equations can be used to obtain grid-independent solutions of realistic flow problems. It is shown that grid adaptation based on metric tensors can generate series of meshes for grid convergence studies in a straightforward way. For a two-dimensional airfoil and the flow around a tanker ship, the grid convergence of the observed forces is sufficiently smooth for numerical uncertainty estimation. Grid refinement captures the details of the local flow in the wake, which is shown to be grid converged on reasonably-sized meshes. Thus, grid convergence studies using automatic refinement are suitable for high-Reynolds incompressible flows.

CFD simulation results for self propelled full scale ship are compared to sea trial measurements in this work. Two–phase RANS based CFD numerical model used in this work is based on the Ghost Fluid Method for numerically robust treatment of discontinuities at the free surface and the algebraic Volume–of–Fluid method for interface capturing. The propeller is modelled as a pressure–jump based actuator disc, allowing CPU time efficient simulations while preserving the accuracy of integral results. The numerical model is implemented in foam–extend, a community driven fork of the OpenFOAM software. The comparison with sea trials includes achieved forward speed, thrust and torque for given shaft speed (in RPM) for a general cargo carrier.

Split Trailing Suction Hopper Dredgers (Split TSHD) have longitudinally-split hull, which symmetrically opens when executing gravity-driven unloading of the cargo, while being ex- posed to various environmental conditions. Even though they have variable hull geometry, their hydrostatic and stability characteristics are usually calculated for initial and unchanged loading conditions only, which is a requirement imposed by classification society stability regulations for TSHD ships [2, 3, 4]. In order to investigate the significance of the discharge process dynamics on actual ship stability, unsteady numerical simulations were performed with the Discrete Element Method (DEM) for symmetrical hopper opening during cargo discharge procedure, without the hull opening failure modes examined. The ship hydrostatic properties, which are pre-calculated analytically using Radial Basis Functions (RBF) for all possible states [11], are used in combin- ation with the solver in order to compute the righting moment and the righting arm, which are affected by the dynamics of the cargo and the loss of displacement. The dynamics of the cargo discharge process was simulated with a DEM solver implemented for Graphics Processing Units (GPUs), Blaze-DEMGPU [8]. Spherical shapes of particulate elements were employed to model the soil cargo, with both cohesion and buoyancy effects included for wetted elements. The simu- lations of the discharging were performed for various loading conditions. Numerical simulations indicate that the dynamics of the cargo during its discharging should not be ignored due to its effect on the transverse stability of the ship. Therefore, an incoming wave and other environ- mental loads in combination with a hull opening failure during the discharge could lead to inapt unstable situation of the ship. Non-symmetrical Split TSHD ship openings will be examined in future work, with an investigation of its influence on ship stability and safety of cargo discharge procedures in failure modes.

In this paper a monolithic algorithm for coupling rigid body motion equations and pressure field within a Finite Volume framework is presented. Monolithic coupling enables fewer number of pressure–velocity loops per time–step, thus reducing the required computational time. The presented method is compared to conventional partitioned coupling approach in terms of computational efficiency and accuracy. The results are compared to experimental data as well.

An important consideration in the catamaran design is the distance between the hulls. Arrangement of the hulls in catamaran configuration can have strong influence on the wave making resistance and thus on the total resistance of a catamaran. The hydrodynamic interaction between hulls becomes significant when spacing between hulls is sufficiently small. In this paper, numerical simulations of viscous flow around monohull and catamaran model are performed utilizing commercial software package STAR-CCM+, in order to investigate the influence of spacing between hulls on the interference resistance. A mathematical model based on Reynolds Avaraged Navier-Stokes (RANS) equations, k-ε turbulence model and Volume of Fluid (VOF) method for describing the motion of two-phase media are briefly described. Numerical simulations are performed for Series 60 monohull and two catamaran configurations with CB=0.6 for different values of Froude number. Results of performed numerical simulations are compared with experimental results available in the literature and satisfactory agreement has been achieved. It has been shown that CFD is a very useful tool in preliminary catamaran design.

To realise the full potential of Computational Fluid Dynamics (CFD) within ma- rine science and engineering, there is a need for continuous maturing as well as verification and validation of the numerical methods used for free surface and interfacial flows. One of the distinguishing features here is the existence of a water surface undergoing large deformations and topological changes during transient simulations e.g. of a breaking wave hitting an off- shore structure. To date, the most successful method for advecting the water surface in marine applications is the Volume-of-Fluid (VOF) method. While VOF methods have become quite advanced and accurate on structured meshes, there is still room for improvement when it comes to unstructured meshes of the type needed to simulate flows in and around complex geometric structures. We have recently developed a new geometric VOF algorithm called isoAdvector for general meshes and implemented it in the OpenFOAM interfacial flow solver called interFoam. We have previously shown the advantages of isoAdvector for simple pure advection test cases on various mesh types. Here we test the effect of replacing the existing interface advection method in interFoam, based on MULES limited interface compression, with the new isoAd- vector method. Our test case is a steady 2D stream function wave propagating in a periodic domain. Based on a series of simulations with different numerical settings, we conclude that the introduction of isoAdvector has a significant effect on wave propagation with interFoam. There are several criteria of success: Preservation of water volume, of interface sharpness and shape, of crest kinematics and celerity, not to mention computational efficiency. We demonstrate how isoAdvector can improve on many of these parameters, but also that the success depends on the solver setup. Thus, we cautiously conclude that isoAdvector is a viable alternative to MULES when set up correctly, especially when interface sharpness, interface smoothness and calcula- tion times are important. There is, however, still potential for improvement in the coupling of isoAdvector with interFoam’s PISO based pressure-velocity solution algorithm.

Validation and verification of a KCS model in head seas is presented in this paper for five wave conditions. Validation is performed by comparing added resistance, heave and pitch motions with experimental data. Numerical model is verified via grid refinement studies, assessing grid uncertainty for each test case. Periodic uncertainty by means of moving window FFT is introduced and assessed. The decomposition model is based on SWENSE (Spectral Wave Explicit Navier-Stokes Equations) method and implicit relaxation zones. Jump conditions at the interface are used to obtain interface-corrected interpolation schemes for density and pressure, embedding the free surface into two single-phase equations for water and air. The method fully resolves issues with parasitic air velocities and ensures one-cell-sharp jump of dynamic pressure and density across the free surface. Implicitly redistanced Level Set equation is used for interface capturing, while the k-omega SST model is used for turbulence. The method is implemented in the Naval Hydro pack based on foam-extend, a fork of the open source software OpenFOAM.

This paper presents detailed validation and verification study for the KCS in calm water conditions. Validation is performed by comparing the resistance, sinkage and trim with experimental data for six Froude numbers, while verification is performed for each case by assessing grid and iterative uncertainties. A single set of governing equations for water and air is used where jump conditions at the free surface are used to discretise discontinuous density and dynamic pressure fields. Volume of Fluid (VOF) method is used for interface capturing, and the k-omega SST model is used for turbulence modelling. Ship motion is solved with quaternion based six degrees-of-freedom (6 DOF) rigid body motion equations. The method is implemented in the Naval Hydro pack based on foam-extend, a community driven fork of OpenFOAM software.