Gatin, Inno

Jasak, Hrvoje

Vukčević, Vuko

FAMENA Repository

CFD Validation of a Container Ship in Calm Waterand Head SeasInno Gatin1, Vuko Vukcˇevic´1 and Hrvoje Jasak11University of Zagreb, Zagreb/Croatia, inno.gatin@stud.fsb.hr, {vuko.vukcevic, hrvoje.jasak}@fsb.hr ,1 IntroductionAdded resistance of ships in waves is one of increasingly important problems in naval engineer-ing due to energy efficiency regulations. In this work, validation of the Naval Hydro pack inOpenFOAM is performed by conducting simulations of a KRISO container ship (KCS) in calmwater and head waves. Steady resistance and dynamic sinkage and trim at different Froudenumbers are compared to experimental data. Mesh refinement study has been carried out fordesign Froude number in order to asses numerical uncertainty. Seakeeping simulations of theship in head waves are carried out for a number of different wave parameters at design Froudenumber. Added resistance is compared with experimental results. All simulations are performedwith fully non–linear, turbulent, two–phase CFD solver. Wave modelling is performed usingSpectral Wave Explicit Navier–Stokes Equations (SWENSE) [1] with implicit relaxation zones[5] that are used to prevent wave reflection. In addition to the validation runs, 3-D freak wavesimulation encountering a KCS has been performed using the SWENSE solver coupled with adirectional Higher Order Spectrum (HOS) [2] method for nonlinear wave propagation.2 Mathematical modelIn this section mathematical model of incompressible two-phase flow is presented. Two–phasemodel is based on a single set of mixture equations taking into account jump conditions at thefree surface implicitly [4]. Volume of Fluid (VOF) method will be used for interface capturing[9] in calm water simulations, while the Level Set (LS) method will be used for head wavessimulations, where SWENSE procedure is employed.Three sets of governing equations are used to model the two phase flow: continuity equation,momentum equation and interface capturing equation. The continuity equation, used to for-mulate dynamic pressure equation, states:∇•u = 0. (1)Mixture momentum equation reads:∂u∂t+∇•(uu)−∇• (νe∇u) = −1ρ∇pd , (2)where νe stands for effective dynamic viscosity, allowing for the general turbulence models [11].Interface capturing equations used in the calculations can be found in Rusche (2002) [7] forVOF method, and in Sun & Beckermann (2007) [8] for LS method.Apart from the governing equations, special attention is given to the jump conditions at the freesurface. In the present model, normal stress balance is modelled exactly, while the tangentialstress is taken into account approximately. Surface tension is neglected. This approach isjustified for large scale two–phase problems encountered in naval hydromechanics. Normalstress balance yields two free surface jump conditions [4]. First condition prescribes the dynamicpressure jump across the free surface:p−d − p+d = −(ρ− + ρ+)g•x , (3)where superscript ”+” denotes the heavier fluid (water), while superscript ”−” denotes thelighter fluid (air). Second condition prescribes the jump of dynamic pressure gradient overdensity to be zero:1ρ−∇p−d −1ρ+∇p+d = 0. (4)Jump conditions are used to discretize dynamic pressure terms (pressure gradient and laplacian)at the free surface.3 Numerical modelNumerical model is based on second–order accurate polyhedral FV method used in foam–extend.Rigid body motion is solved using quaternion based six degrees of freedom (6DOF) equations.Mesh is modelled as a rigid body with special boundary conditions. Coupling of pressure,velocity, free surface and 6DOF equations is performed in a segregated manner using PIMPLEalgorithm.4 Numerical simulationsIn this section simulation results regarding KCS model are compared to experimental data.Steady state resistance with dynamic sinkage and trim is simulated for a range of Froudenumbers. Seakeeping simulations are performed in regular head waves with different incidentwave parameters at design Froude number. Finally, a simulation of a freak wave encounteringa full scale KSC is shown.4.1 Steady resistance, simulations with dynamic sinkage and trimCFD results computed on 4.6 million cells unstructured mesh are compared with experimentaldata for following Froude numbers: 0.1516, 0.1949, 0.2274, 0.2599 and 0.2816. KCS modelparticulars and experimental settings details are available at Tokyo Workshop On CFD inShip Hydrodynamics website [6]. Symmetry boundary condition is used at the central verticallongitudinal plane, and only half of the domain is simulated. Average mesh non–orthogonalityis 6.6 degrees while the maximum is 76.1 degrees. Results comparison for total resistancecoefficient is shown in Figure 1. Dynamic sinkage and trim comparison is shown in Figure 2and Figure 3, respectively.Grid convergence study is performed for 0.282 Froude number. Table 1 shows the to-tal resistance coefficient, sinkage and trim obtained with different grid sizes. Relative er-ror compared to experimental result is also presented. Relative error is calculated as  =(CtEFD − CtCFD)/CtEFD. Calculated order of spatial accuracy for total resistance [10] isp = 3.98 for three finest grids, while grid uncertainty is 0.029× 10−3.4.2 Seakeeping simulations in regular head wavesAdded resistance is compared with experimental data for various incident wave parameters.Simulations are performed for design Froude number and head waves. Comparison is performedin Fourier space, comparing individual harmonic amplitudes to experimental data [6]. KCSmodel particulars can be found on Tokyo Workshop On CFD in Ship Hydrodynamics website[6]. Figure 4 shows comparison of mean, first and second order of total resistance coefficient fordifferent incident wave parameters (λ/Lpp).0.1 0.15 0.2 0.25 0.3Fr3.43.63.844.24.4C T·103CT, EFDCT, CFDFigure 1: Comparison of total resistance coefficient for a range of Froude numbers.0.1 0.15 0.2 0.25 0.3Fr-2-1.5-1-0.50σ, cmσ, EFDσ, CFDFigure 2: Comparison of dynamic sinkage for a range of Froude numbers.Table 1: Grid convergence study results.Grid No. 1 2 3 4 EFDCt × 103 4.837 4.736 4.608 4.560 4.501Ct ,% -7.476 -5.234 -2.381 -1.313 /σ, cm -1.665 -1.650 -1.655 -1.653 -1.702σ,% 2.194 3.067 2.756 2.867 /τ,o -0.145 -0.147 -0.147 -0.146 -0.159τ ,% 8.813 7.264 7.847 7.974 /It can be seen in Figure 4 that the trend of all the curves is well captured. Mean valuesobtained with CFD generally overshoot the experimental data, which is consistent with the0.1 0.15 0.2 0.25 0.3Fr-0.2-0.15-0.1-0.050τoτ, EFDτ, CFDFigure 3: Comparison of dynamic trim for a range of Froude numbers.0.6 0.8 1 1.2 1.4 1.6 1.8 2λ/LPP0510152025C T· 103CFD - meanEFD - meanCFD - 1st orderEFD - 1st orderCFD - 2nd orderEFD - 2nd orderFigure 4: Comparison of total resistance coefficient for various incident wave parameters.steady state resistance simulations. First order resistance is undershot for all cases.4.3 Freak wave impact on a full scale floating KCSA three-dimensional, realistic freak wave is initialized using long time evolution of wave spectrumwith Higher Order Spectrum (HOS) potential nonlinear wave theory [2] implemented in NavalHydro pack. Initial condition used for HOS simulation is JONSWAP directional spectrum.Wave spectrum parameters are: Hs = 10.5 m, Tp = 9.5 s, and directionality constant n = 8 [3].Freak wave height is H = 21.8 m. Simulation is performed for full size KCS, Table 2. Figure 5sequentially shows the encounter of the freak wave on the KCS.Table 2: Full scale KCS characteristics.L, m B, m D, m T, m xG, m yG, m KG, m rx, m ry, m rz, m230 32.2 19 10.8 111.6 0 14.32 12.9 57.5 57.5(a) (b)(c) (d)Figure 5: Three–dimensional freak wave encountering a full scale KCS.5 ConclusionSimulations of KCS in calm water and head waves are carried out and the results are com-pared with experimental data. Steady state simulations show good agreement for resistanceand dynamic sinkage and trim. Larger discrepancies of dynamic sinkage and trim for lowerFroude numbers are caused by coarse mesh considering the absolute values of sinkage and trim.Preliminary results of added resistance in head waves are compared to experimental data, be-ing an important design criteria. Comparison showed that added resistance can be predictedfairly accurately. Finally, an example of three–dimensional freak wave simulation is performedto show the capability of the HOS wave theory coupled with CFD simulation using SWENSEmethod.References[1] Ferrant, P. and Gentaz, L. and Le Touze´, D. A new ranse/potential approach for waterwave diffraction. Proc. Numerical Towing Tank Symposium, 2002.[2] I. Gatin. Master Thesis; Implementation and validation of the Higher Order Spectrum(HOS) method for Simulation of Surface Waves. Faculty of Mechanical Engineering andNaval Architecture, Zagreb, 2015.[3] L. H. Holthuijsen. Waves in Oceanic and Costal Waters. Delft University of Techonlogyand UNESCO-IHE, 2007.[4] J. Huang, P. M. Carrica, and F. Stern. Coupled ghost fluid/two–phase level set methodfor curvilinear body–fitted grids. Int. J. Numer. Meth. Fluids, 44:867–897, 2007.[5] H. Jasak, V. Vukcˇevic´, and I. Gatin. Numerical Simulation of Wave Loads on StaticOffshore Structures. In CFD for Wind and Tidal Offshore Turbines, pages 95–105. SpringerTracts in Mechanical Engineering, 2015.[6] National Maritime Reasearch Inistitute (NMRI). Tokyo 2015: A Workshop on CFD in ShipHydrodynamics. http://www.t2015.nmri.go.jp/, 2015. [Online; accessed 20 August2015].[7] H. Rusche. Computational Fluid Dynamics of Dispersed Two - Phase Flows at High PhaseFractions. PhD thesis, Imperial College of Science, Technology & Medicine, London, 2002.[8] Y. Sun and C. Beckermann. Sharp interface tracking using the phase–field equation. Jour-nal of Computational Physics, 220:626–653, 2007.[9] O. Ubbink and R. I. Issa. A method for capturing sharp fluid interfaces on arbitrarymeshes. Journal of Computational Physics, 153:26–50, 1999.[10] H. K. Versteeg and W. Malalasekera. An introduction to Computational Fluid Dynamics:The finite volume method. Pearson Education Limited, 1995.[11] D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, 1993.

CFD Validation of a Container Ship in Calm Water and Head Seas

Abstract

Similar works

Full text

Available Versions

FAMENA Repository