Next Article in Journal
Simulation of Electric Vehicle Charging Points Based on Efficient Use of Chargers and Using Recuperated Braking Energy from Trains
Previous Article in Journal
Examining the Relationship between Energy Consumption and Unfavorable CO2 Emissions on Sustainable Development by Going through Various Violated Factors and Stochastic Disturbance–Based on a Three-Stage SBM-DEA Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrodynamics of Moonpool-Type Floaters: A Theoretical and a CFD Formulation

by
Anargyros S. Mavrakos
1,
Dimitrios N. Konispoliatis
1,*,
Dimitrios G. Ntouras
2,
George P. Papadakis
2 and
Spyros A. Mavrakos
1
1
Laboratory for Floating Structures and Mooring Systems, School of Naval Architecture and Marine Engineering, National Technical University of Athens, GR 157-73 Athens, Greece
2
Laboratory for Naval and Marine Hydrodynamics, School of Naval Architecture and Marine Engineering, National Technical University of Athens, GR 157-73 Athens, Greece
*
Author to whom correspondence should be addressed.
Energies 2022, 15(2), 570; https://doi.org/10.3390/en15020570
Submission received: 29 November 2021 / Revised: 7 January 2022 / Accepted: 10 January 2022 / Published: 13 January 2022

Abstract

:
Moonpool-type floaters were initially proposed for applications such as artificial islands or as protecting barriers around a small area enabling work at the inner surface to be carried out in relatively calm water. In recent years, a growing interest on such structures has been noted, especially in relation to their use as heaving wave energy converters or as oscillating water column (OWC) devices for the extraction of energy from waves. Furthermore, in the offshore marine industry, several types of vessels are frequently constructed with moonpools. The present paper deals with the hydrodynamics of bottomless cylindrical bodies having vertical symmetry axis and floating in a water of finite depth. Two computation methods were implemented and compared: a theoretical approach solving analytically the corresponding diffraction problem around the moonpool floater and a computational fluid dynamics (CFD) solver, which considers the viscous effects near the sharp edges of the body (vortex shedding) as non-negligible. Two different moonpool-type configurations were examined, and some interesting phenomena were discussed concerning the viscous effects and irregularities caused by the resonance of the confined fluid.

1. Introduction

Moonpool-type floaters were initially proposed to be applied either as artificial islands or as structures enabling activities in their inner area, as they used their solid boundaries as protecting walls to the wave action. Nevertheless, in recent years, they have been associated with their use as wave energy converter (WEC) devices suitable for wave energy absorption, as well as floating pontoons for semi-submersible designs.
Several types of bottomless cylindrical hulls (known as moonpools) have been examined in the literature. The most common is a vertical cylinder with a vertical opening from the deck to keel through the hull [1], whereas square floaters with square moonpools have been also studied [2,3]. Additional types concern toroidal configurations consisting of a circular (core) section centered upon a large circle (ring) with an infinitesimally small [4] and finite [5] slenderness core-to-ring-radius ratio, as well as two coaxial surface piercing truncated circular cylinders [6]. The latter geometry was encompassed by an exterior partially immersed toroidal structure of finite volume and an interior coaxial free surface-piercing truncated cylinder. Thus, an internal annular free surface was developed, which was totally enclosed between the solid boundaries of the cylinders and open to the exterior fluid domain beneath the concentric bodies.
Hydrodynamics of the isolated truncated hollow cylindrical bodies have been investigated in the past. Garrett [7] presented an analytical solution of the wave scattering problem by a vertical bottomless partially immersed circular cylinder with infinitesimal wall thickness using matched axisymmetric eigenfunction expansions. Analytical computations of the exciting wave forces and moments on the hollow cylindrical shell structure, together with sway and pitch hydrodynamic coefficients, were derived by Miloh [8]. He also derived analytic expressions for the amplification ratios of heave (pumping) and pitch (seiche) modes of water motions inside the pond. The corresponding hydrodynamic characteristics of an open-ended vertically floating circular cylinder, as well as the induced water motion in the interior basin, accounting for the effect of the finite wall thickness, were presented by [1,9]. Furthermore, Liu et al. [10] investigated the second-order sum frequency wave loads on a hollow cylindrical shell structure in long crested irregular seas. Additionally, in [11,12] the solutions of the diffraction and radiation problems of a moonpool cylindrical body of infinite thickness and of a partially bottom-opened moonpool were presented. A moonpool with a restricted entrance was also studied. In [13], it was shown that the motion of the water inside the body was amplified as the size of the moonpool entrance was reduced. In addition, Liu et al. [14] measured the viscous effect of a moonpool cylindrical structure with a restricted entrance by introducing a quadratic dissipation assuming an additional dissipative disk at the moonpool entrance.
Regarding the two coaxial surface piercing truncated circular cylinder configuration, the solutions of the linearized diffraction and radiation problems around two independently moving concentric truncated cylinders were presented in [6,15]. Furthermore, the mean- and the time-dependent second-order wave loads on such a type of floating structures were evaluated in [16,17]. Recently, Chao and Kim [18] investigated the hydrodynamic performance of a two-body WEC by applying analytical solutions, whereas Kong et al. [19] examined a moonpool platform-wave energy buoy composed of an inner cylindrical buoy and an outer toroidal-cylinder for wave power extraction by the relative heave motions between the inner and outer buoys.
Apart from wave energy applications, concentric surface-piecing circular cylindrical bodies have also been investigated, combined with the trapped wave modes that they can develop, i.e., the localized oscillations of the enclosed formed free surface area, whereas the motion in the exterior fluid region decayed to zero. Abramson [20] presented a comprehensive review of liquid sloshing problems in moving containers, while Shipway and Evans [21] examined the wave trapping phenomena by two concentric circular cylindrical shells with a vertical axis and a zero thickness. As an extension of [21]’s work, McIver and Newman [22] studied the case of vertically axisymmetric trapping structures, of finite volumes, composed by two concentric interior free surfaces.
All the aforementioned research is based on potential flow methodologies concerning the hydrodynamic solutions of several types of bottomless cylindrical hulls. Nevertheless, over the last decade, there has been a significant interest in computational fluid dynamics (CFD) modelling due to its detailed results. CFD modelling solves the fluid flow problem using Navier–Stokes equations, where the hydrodynamics of random-shape floating structures can be calculated in detail, providing representative results of the flow physics. Several authors have studied the importance of including viscosity in moonpool problems using CFD modelling. Lu et al. [23] investigated the possibility of finding an artificial, empirically-based damping coefficient to the free surface condition inside a moonpool device applied to potential flow models. The latter exhibit overestimated magnitudes of wave forces as the fluid resonance takes place; hence, by introducing an artificial damping term, the potential flow model may work as well as the viscous fluid model. Lu et al. [24] extended the results from [23] on the wave elevation inside a moonpool body; whereas, in [25,26], the piston-mode in a 2D moonpool under forced heave oscillations was evaluated using a linear potential solver coupled to a viscous Navier–Stokes modelling. In addition, Fredriksen et al. [27] investigated, numerically and experimentally, how a low forward/incoming current speed influences the resonant piston-mode resonance in a moonpool. They concluded that the moonpool behavior was marginally reduced with a low forward velocity, whereas this reduction was dependent on the heave forcing amplitude. Recently, in [28] a viscous damping model for energy losses of a 2D fluid response in a moonpool due to wall friction and flow separation was proposed. Furthermore, the limitations of boundary element method solvers when viscous effects near floater’s sharp edges were taken into account and were exploited in [29] in the case of a moonpool type floater.
In the present article, two different geometric configurations of a moonpool-type floater were investigated; a floater encompassed by two coaxial, free-surface piercing toroidal cylinders with vertical symmetry axes (Configuration 1, see Figure 1a) and a floater consisting of an external toroidal cylinder and an internal coaxial surface piercing truncated circular cylinder (Configuration 2, see Figure 1b). In the annular and the cylindrical fluid areas formed between the cylinders’ vertical walls, oscillations of the enclosed water columns were developed. A semi-analytical model was used to solve the wave diffraction problem in the context of linear potential theory [30]. The method of matched axisymmetric eigenfunction expansions was applied to solve the relevant hydrodynamic problem, according to which the flow field around and inside the examined body was subdivided in coaxial ring-shaped fluid regions. Hence, appropriate series representations of the velocity potential can be established in each region, which were matched at the boundaries of adjacent fluid domains by enforcing the continuity conditions of the hydrodynamic pressure and the radial velocity in order to determine the unknown coefficients in the expressions of the diffraction potential. The exciting wave forces/moments and the wave elevation on both case studies were calculated and compared focusing on the resonant wave frequencies. Furthermore, the accuracy of the semi-analytical formulation was assessed in comparison to high fidelity CFD simulations.
The present work is structured as follows: in Section 2 the solution of the diffraction problem is formulated, along with the evaluation of the hydrodynamic forces and wave elevation around/inside the examined case studies. Section 3 is dedicated to the presentation of the used CFD solver, whereas, in Section 4, the outcomes of the two methods (i.e., semi-analytical and CFD solvers) are compared for the two moonpool arrangements. Finally, conclusions are drawn in Section 5.

2. Material and Methods

2.1. Velocity Potential Representation

Two coaxial, free-surface piercing toroidal cylinders with vertical symmetry axes are considered floating in constant water depth d. The cylinders are exposed to the action of a regular wave-train of wave frequency ω and amplitude H / 2 . The draughts of the moonpool bodies are denoted by h 1 and h 2 , respectively. Cylindrical co-ordinates ( r , θ , z ) are introduced with the z axis positive upwards and with origin on the seabed (see Figure 1). In the present formulation, viscous effects were neglected, whereas the fluid was assumed incompressible. Furthermore, the fluid’s motions were assumed to be small, so that the linearized diffracted problem could be considered.
The fluid flow is described by the velocity potential expressed as:
Φ ( r , θ , z ; t ) = R e [ φ ( r , θ , z ) e i ω t ]
with
φ ( r , θ , z ) = φ D ( r , θ , z ) = φ 0 ( r , θ , z ) + φ 7 ( r , θ , z )
Here, φ 0 is the velocity potential of the undisturbed incident harmonic wave and φ 7 is the diffraction potential for the restrained body in the wave impact. The diffraction potential is described by the velocity potential φ D .
The undisturbed velocity potential of an incident wave train propagating along the positive axis can be expressed in cylindrical coordinates as:
φ 0 ( r , θ , z ) = i ω H 2 Z 0 ( z ) Z 0 ( d ) m = 0 ε m i m J m ( k r ) cos ( m θ )
In Equation (3) J m stands for the mth order Bessel function of first kind, ε m is the Neumann’s symbol, i.e., ε 0 = 1 ,   ε m = 2 for m 1 , and Z 0 is derived by:
Z 0 ( z ) = [ 1 2 [ 1 + sin h ( 2 k d ) 2 k d ] ] 1 2 cos h ( k z )
with Z 0 ( d ) being its derivative at z = d. The wave number k and the wave frequency ω are related by the dispersion equation i.e.,: ω 2 = k g tan h ( k d ) .
In accordance with Equation (3) the total diffracted potential of the flow field around the examined moonpool can be written in the form:
φ D ( r , θ , z ) = i ω H 2 m = 0 ε m i m Ψ D , m ( r , z ) cos ( m θ )
The Ψ D , m denote the principal unknowns of the problem. The subscript D stands for the diffraction problem, whereas the subscript m indicates the modes that should be taken into consideration within the solution of the corresponding problem.
The potential φ D must satisfy:
Δ φ D = 2 φ D r 2 + 1 r φ D r + 1 r 2 2 φ D θ 2 + 2 φ D z 2 = 0
with the linearized boundary conditions:
on the free surface z = d:
ω 2 φ D g φ D z = 0
on the seabed z = 0:
φ D z = 0
on the body’s bottom ( z = d h 1 , b r a and z = d h 2 ,   l r c ):
φ D z = 0
on the body’s side walls ( d h 1 z d , r = a , b   and   d h 2 z d ,   r = c , l ):
φ D r = 0
It is also required that the scattered potential satisfies an appropriate radiation condition as r which has the form [1,31]:
lim r r ( r i k ) φ 7 = 0
In addition, both the velocity potential and its derivatives φ D r must be continuous at the boundaries r = a , b , c , l of the neighboring ring elements shown in Figure 1a. This results in: for 0 < z < d h 1 , and r = a :
Ψ D , m I ( a , z ) = Ψ D , m I I ( a , z )
Ψ D , m I ( r , z ) r | r = a = Ψ D , m I I ( r , z ) r | r = a
for 0 < z < d h 1 , and r = b :
Ψ D , m I I ( b , z ) = Ψ D , m I I I ( b , z )
Ψ D , m I I ( r , z ) r | r = b = Ψ D , m I I I ( r z ) r | r = b
for 0 < z < d h 2 , and r = c :
Ψ D , m I I I ( c , z ) = Ψ D , m I V ( c , z )
Ψ D , m I I I ( r , z ) r | r = c = Ψ D , m I V ( r , z ) r | r = c
for 0 < z < d h 2 , and r = l :
Ψ D , m I I I ( c , z ) = Ψ D , m I V ( c , z )
Ψ D , m I I I ( r , z ) r | r = c = Ψ D , m I V ( r , z ) r | r = c
The superscripts I, II, III, IV, V imply quantities corresponding to respective types of ring elements.
The procedure for the determination of the velocity potential around a moonpool body that is exposed to the action of surface gravity waves has been presented thoroughly in previous studies: [1,5,6,9,15], to name a few. This procedure is extended here to the case of two coaxial surface piercing toroidal vertical cylinders. The following equations present proper series representations of the functions Ψ D , m i ,   i = I , I I , I I I , I V , V in each fluid region around the body:
(a)
Infinite fluid domain, Type I ( a r ,   0 z d ):
1 d Ψ D , m I = { J m ( k r ) J m ( k a ) H m ( k a ) H m ( k r ) } Z 0 ( z ) d Z 0 ( d ) + j = 0 F m j I K m ( a j r ) K m ( a j a ) Z j ( z )
where H m ,   K m is the m-th order Hankel function of the first kind and the modified Bessel function of the second type, respectively. In addition, Z j ( z ) are orthonormal functions in [ 0 , d ] defined by Equation (4) for j = 0 and for j > 0:
Z j ( z ) = [ 1 2 [ 1 + sin ( 2 a j d ) 2 a j d ] ] 1 2 cos ( a j z ) ,   j > 0
The terms a j are the roots of the transcendental equation: ω 2 + a j g tan ( a j d ) = 0 , which proposes one imaginary a 0 = i k ,   k > 0 and infinite number of real roots.
(b)
Second fluid domain, Type II ( b r a ,   0 z d h 1 ):
1 d Ψ D , m I I = n = 0 ε n [ R m n I I ( r ) F m n I I + R m n * I I ( r ) F m n * I I ] cos ( n π z d h 1 )
where ε n stands for the Neumann’s symbol and the functions R m n I I ( r ) ,   R m n * I I ( r ) are presented in Appendix A.
(c)
Third fluid domain, Type III ( c r b ,   0 z d ):
1 d Ψ D , m I I I = i = 0 [ R m i I I I ( r ) F m i I I I + R m i * I I I ( r ) F m i * I I I ] Z i ( z )
where Z i are orthonormal functions in [ 0 , d ] defined by Equations (4) and (21) for i = 0 ,   i > 0 , respectively. The functions R m i I I I ( r ) , R m i * I I I ( r ) are presented in Appendix A.
(d)
Fourth fluid domain, Type IV ( l r c ,   0 z d h 2 ):
1 d Ψ D , m I V = s = 0 ε s [ R m s I V ( r ) F m s I V + R m s * I V ( r ) F m s * I V ] cos ( s π z d h 2 )
Here ε s stands for the Neumann’s symbol and the functions R m s I V ( r ) ,   R m s * I V ( r ) are defined in Appendix A.
(e)
Fifth fluid domain, Type V ( 0 r l ,   0 z d ):
1 d Ψ D , m V = q = 0 F m q V I m ( a q r ) I m ( a q l ) Z q ( z )
where Z q are orthonormal functions in [ 0 , d ] defined by Equations (4) and (21) for q = 0 ,   q > 0 , respectively, whereas the term I m is the mth order modified Bessel function of the first kind. We note that the modified Bessel function of second type K m ( a q r ) does not appear in the velocity potential representation of this fluid domain type, as the function K m has a singular behavior at r 0 . The latter has namely a contribution proportional to ( l / r ) m for m 0 and log ( r ) for m = 0 at r 0 [32].
Based on the derived velocity potentials in each fluid domain, the continuity requirements of the hydrodynamic pressure and radial velocity at the boundaries of neighboring fluid regions, expressed through Equations (12)–(19), along with the kinematic boundary conditions at the vertical walls of the body, are then to be fulfilled (see Equation (10)). Hence, the linear system of equations for the determination of the unknown Fourier coefficients F m j l ,   l = I , I I , I I I , I V , V , and F m j * s ,   s = I I , I I I , I V , in each fluid domain, was derived. The solution procedure is presented in the following Section 2.2.

2.2. Solution Procedure

In this subsection, the solution procedure for determining the unknown Fourier coefficients in each fluid domain is presented. Towards this goal, we first expressed the velocity potentials of the fluid domains at their vertical boundaries. To this end, by accounting for the values of the functions R m n I I ( r ) ,   R m n * I I ( r ) ,   R m i I I I ( r ) ,   R m i * I I I ( r ) ,   R m s I V ( r ) ,   R m s * I V ( r ) (see Appendix A) at r = a ,   b ,   c   and   r = l , the potential functions Ψ D , m i ,   i = I , I I , I I I , I V , V , see Equations (20), (22)–(25), can be reduced to the following simple series expressions at r = a , b , c , l :
1 d Ψ D , m I ( a , z ) = j = 0 F m j I Z j ( z )
1 d Ψ D , m I I ( a , z ) = n = 0 ε n F m n I I cos ( n π z d h 1 )
1 d Ψ D , m I I ( b , z ) = n = 0 ε n F m n * I I cos ( n π z d h 1 )
1 d Ψ D , m I I I ( b , z ) = i = 0 F m i I I I Z i ( z )
1 d Ψ D , m I I I ( c , z ) = i = 0 F m i * I I I Z i ( z )  
1 d Ψ D , m I V ( c , z ) = s = 0 ε s F m s I V cos ( s π z d h 2 )
1 d Ψ D , m I V ( l , z ) = s = 0 ε s F m s * I V cos ( s π z d h 2 )
1 d Ψ D , m V ( l . z ) = q = 0 F m q V Z q ( z )
Next, by taking into consideration Equations (26)–(33), which represent the velocity potential at the vertical boundaries of the various fluid regions, as well as Equations (12), (14), (16), (18), which express the conditions for continuity of the potential function at r = a , b , c , l , we multiply both sides of Equations (12) and (14) by 1 / ( d h 1 ) cos ( v π z / ( d h 1 ) ) , v = 0 , 1 , 2 , , n . , and of Equations (16) and (18) by 1 / ( d h 2 ) cos ( v π z / ( d h 2 ) ) , v = 0 , 1 , 2 , , s . , and then we integrate the resulting expressions over their regions of validity. The resulting expressions relating the Fourier coefficients of adjacent fluid regions are:
F m s * I V = q = 0 L s q F m q V
F m s I V = i = 0 L s i F m i * I I I
F m n * I I = i = 0 L n i F m i I I I
F m n I I = j = 0 L n j F m j I
The L s q , L s i , L n i , L n j terms are defined in Appendix B.
Furthermore, the conditions for the continuity of the first derivatives of the potential at r = a , b , c , l (see Equations (13), (15), (17) and (19)), as well as the kinematical boundary conditions at the body’s vertical boundaries (see Equation 10), must also be fulfilled. Hence, multiplying both sides of Equations (10), (13), (15), (17) and (19) by ( 1 / d ) Z v ( z ) , v = 0 , 1 , 2 , , q , . , integrating over the regions of their validity and adding the resulting expressions, the flowing set of equations can be obtained:
F m j I Λ j I = d h 1 d n = 0 ε n L n j [ Λ m n I I F m n I I + Λ m n * I I F m n * I I ] ,
  for   r = a ,   0 z d Λ m i I I I F m i I I I + Λ m i * I I I F m i * I I I = d h 1 d n = 0 ε n L n i [ D m n I I F m n I I + D m n * I I F m n * I I ] ,   for   r = b ,   0 z d  
D m i I I I F m i I I I + D m i * I I I F m i * I I I = d h 2 d s = 0 ε s L s i [ Λ m s I V F m s I V + Λ m s * I V F m s * I V ] , for   r = c ,   0 z d  
F m q V Λ q V = d h 2 d s = 0 ε s L s q [ D m s I V F m s I V + D m s * I V F m s * I V ] ,
for   r = l ,   0 z d
The terms Λ j I , Λ q V ,   Λ m n I I , Λ m n * I I ,   D m n I I , D m n * I I , Λ m i I I I , Λ m i * I I I ,   D m i I I I , D m i * I I I , Λ m s I V , Λ m s * I V ,   D m s I V , D m s * I V in Equations (38)–(41) are described in Appendix B.
The solution of the infinite system described by the linear Equations (34)–(37) and (38)–(41) will provide the unknown Fourier coefficients, F m j l ,   l = I , I I , I I I , I V , V , and F m j * s ,   s = I I , I I I , I V , in each fluid domain. It should be noted that for the numerical implementation of the method, infinite series (20), (22)–(25), are truncated after J, N, I, S, Q terms, respectively. A similar procedure has also been used in [1,5,6].

2.3. Hydrodynamic Exciting Forces/Moment and Wave Elevation

Having evaluated the diffraction velocity potential, the hydrodynamic pressure in each fluid region could be estimated by the linearized Bernoulli equation. Hence, the exciting forces and moments on the moonpool cylindrical body could be determined by integrating the hydrodynamic pressure, p , over the body’s wetted surface, as follows:
F k ( t ) = S 0 p n k d S = i ω ρ e i ω t S 0 φ D n k d S = ω 2 ρ H 2 e i ω t m = 0 ε m i m S 0 Ψ D , m ( r , z ) cos ( m θ ) n k d S , for   k = 1 , 3
and
F 5 ( t ) = M 1 ( t ) + M 3 ( t )
where:
M k ( t ) = S 0 p ( r × n k ) d S = i ω ρ e i ω t S 0 φ D ( r × n k ) d S = ω 2 ρ H 2 e i ω t m = 0 ε m i m S 0 Ψ D , m ( r , z ) cos ( m θ ) ( r × n k ) d S , for   k = 1 , 3
Here S 0 is the mean wetted surface of the body; ρ is the water density; n k are the generalized normal components defined by n = ( n 1 , n 2 , n 3 ) ;   r ×   n = ( n 4 , n 5 ,   n 6 ) , r being the position vector of a point on the wetted surface with respect to the reference co-ordinate system of the body. Analytical representations of the exciting forces and moments on the body are presented in Appendix C.
The free surface elevation outside and in the annulus and moonpool areas inside the body is given by the linearized boundary condition on the free surface (see Equation (7)). It holds:
Z ( r , θ , z ; t ) = R e { ζ ( r , θ , d ) e i ω t } = R e { i ω g φ D ( r , θ , d ) ( cos ( ω t ) i sin ( ω t ) ) }

3. Computation Fluid Dynamics—MaPFlow Solver

In this section, a brief description of the CFD modelling applied to the present bodies’ configurations is presented. For a more detailed review of the solver’s methodology, the reader may refer to [29,33]. In the previous mentioned publications, a series of test cases have been examined, including a moonpool configuration excited by incident waves that prove the liability of the solver in this type of simulations.
MaPFlow treats free surface flows as flows of two immiscible and incompressible fluids. The presence of the two fluids is described using the Volume of Fluid (VoF) approach [34], while the requirement of the incompressibility is satisfied using the artificial compressibility method (ACM) [35]. The governing equations are presented in Equation (46). The system of equations, in 3 dimensions, consists of 5 scalar equations. The equations are augmented by the pseudo-time derivatives of the variables. The aim of the numerical procedure is to drive these derivatives to zero; thus the original unsteady system of equations will be retrieved. The coupling of the equations is performed during the pseudo-time, where a relation between the density and the pressure field is assumed. The coupling is controlled through the relation ρ p | τ = 1 β , where β is a free parameter that in typical free surface flows ranges from 5 to 15 [33,36].
Γ Ω Q τ d Ω + Γ e Ω Q t d Ω + Ω ( F c F v ) d S = Ω S q d Ω
In the above integral equation, Ω is a control volume with boundary Ω , Q = [ p , u , a ] T is the vector of the unknown variables (pressure p , velocity u and volume fraction a ), vector S q contains the various source terms of the equations such as the gravitational field, while t and τ denote the real and fictitious time, respectively. The original ACM when applied to two-phase flows scales locally with the density field and this results to a stiff numerical behavior. In order to remedy this dependency, the preconditioner Γ of Kunz [37] is used; it is presented in Equation (47). Furthermore, to recover the conservative form of the equations, the unsteady terms of Equation (46) are supplemented with the matrix Γ e , which is also presented in Equation (47).
Γ = [ 1 β ρ m 0 0 0 ρ m I u Δ ρ α β ρ m 0 1 ] ,   Γ e = [ 0 0 0 0 ρ m I u Δ ρ 0 0 1 ]
Finally, F c is the vector of the convective fluxes and F v the vector of the viscous fluxes. The two vectors are given by Equation (48).
F c = [ V n ρ m u Δ V + p n x ρ m v Δ V + p n y ρ m w Δ V + p n z a V n ] , F v = [ 0 τ x x n x + τ x y n y + τ x z n z τ y x n x + τ y y n y + τ y z n z τ z x n x + τ z y n y + τ z z n z 0 ]
In Equation (48), Δ V is the velocity difference between the contravariant velocity V n = u · n and the grid face velocity due to the mesh motion V g = u v o l · n . In the present study, because the structures were considered fully restrained, the mesh was stationary; thus, V g = 0 . Furthermore, the viscous fluxes are computed through the viscous stresses τ i j , which in turn were expressed using the Boussinesq approximation, by Equation (49).
τ i j = ( μ t + μ   m ) ( u i x j + u j x i ) 2 3 ρ m k δ i j
In the above equation, μ   m is the viscosity of the mixture, μ t is the turbulent viscosity, k the turbulent kinetic energy and δ i j the Kronecker’s symbol. Moreover, since here only small amplitude linear waves were investigated, the turbulence modelling was switched off ( μ t = 0 ,   k = 0 ).
The generation and damping of the free surface waves were performed by adding source terms to the momentum equations in user-specified regions near the boundaries of the computational domain. The form of this source term is described by Equation (50). The source terms scale according to the coefficient C n w t , which is controlled though the factor α n w t [ 1 s ] and the function f ( x n w t ) defined inside the wave generation or damping zone. f ( x n w t ) has typically an exponential or polynomial form [38]; α n w t is a free parameter that depends on the flow itself (e.g., the waves entering the damping zone) and the discretization process. In [33], a parametric analysis of the effects of the coefficient C n w t is presented. The source term S n w t drives the solution to the desired velocity field u t a r . The target velocity in the case of generation zone was the velocity field provided by an analytical solution of free surface waves (linear, Stokes nth order waves, etc.); in the case of damping zone, the usual strategy was to set the normal velocity to the boundary equals to zero.
S n w t = C n w t ρ m   ( u u t a r     ) ,                 C n w t = α n w t f ( x n w t )
The above system of equations was discretized using the finite volume method. In each cell of the computational mesh, a control volume Ω I was defined, where the unknown variables had a constant value and the control volume was represented by the geometric center of the computational mesh (cell-centered approximation). The surface integrals were decomposed into a sum of constant surface integral over the faces that make up the cell. The semi-discretized form of the equations is presented in Equation (51). The source terms of the equations are considered constant in each cell.
Γ e ( Q Ω I ) t + Γ Ω I Q τ = R Ω I
R Ω I is the spatial residual of the discretization process (see Equation (52)).
R Ω I j N f ( Ω I ) ( F c F v ) j Δ S j + Ω I S q
The convective fluxes are computed using the approximate Riemann solver of Roe [39,40], whereas a second order interpolation scheme supplemented by a directional derivative to account for the skewness of the mesh has been applied for the computation of the viscous fluxes.
For the time discretization, a second order backwards differentiation scheme was used [41] for the real time derivatives, while, for the pseudo-time derivatives, a first order explicit Euler scheme was utilized. In order to facilitate convergence, the pseudo-timestep scales locally in each cell based on a given Courant-Friedrichs-Lewy (CFL) number and the characteristics of the flow. For an analytic description of the time discretization process the reader can refer to [33].
Finally, regarding the boundary conditions, at the wall boundaries, the no-slip boundary condition was applied; the symmetry plane was treated as a slip boundary, while the far field boundaries were approximated using ghost cells. At every far field boundary face, an approximate Riemann problem was solved between the state inside the domain and the far field conditions imposed by the ghost cell. It should be noted that, in the case of wave generation, the ghost cell’s values were given by the current wave solution.
The discretization of the equations led to a non-linear system of equations. The linearization process took places in pseudo-time. The resulting linear system was solved using a Gauss–Seidel iterative method, while it was permutated using the Reverse Cuthill-Mckee reordering scheme to accelerate convergence.
The above methodology was implemented as part of the CFD solver MaPFlow [42,43]; MaPFlow was developed in NTUA, solved the URANS equations in unstructured meshes and is capable of running up to 1000 parallel processes utilizing the MPI protocol.

4. Results

In accordance with the procedure presented in the previous sections, it is obvious that the accuracy of the method and the CPU time required was affected by the evaluation procedure of the Fourier coefficients in each fluid domain. In the present calculations, 80 terms were used for the series expansions of the velocity potential in the fluid domains of the I, III and V types, whilst 150 terms are retained for the velocity representation in the II, IV types.
Two different geometric configurations of a moonpool-type floater were investigated; a floater encompassed by two coaxial, free-surface piercing moonpool cylinders with vertical symmetry axes (Configuration 1, see Figure 1a) and a floater consisting of an external toroidal cylinder and an internal coaxial surface piercing truncated circular cylinder (Configuration 2, see Figure 1b) [30]. Specifically, for the Configuration 1 (i.e., coaxial hollow cylinders) it held: α = 13 m, b = 12 m, c = 9 m, l = 6.083 m, and h1 = 14 m, h2 = 5.5 m (see Figure 2a), whereas, for the Configuration 2 (i.e., coaxial truncated cylinder), it was assumed: α = 13 m, b = 12 m, c = 9 m, h1 = 14 m and h2 = 5.5 m (see Figure 2b). Both arrangements were considered restrained to the wave impact at a water depth of d = 70 m.
Regarding the numerical results that were produced by the in-house CFD solver, the two different configurations were excited by linear waves with wave frequencies ranging from 0.6 rad/s to 1.2 rad/s, with a constant wave height of 0.6 m. Based on the linear wave theory of Airy, the velocity field of the wave was imposed on the one end of the domain, while, at the other end, wave damping was used. The waves were propagated along the x-axis, while the free surface was placed perpendicular to the z-axis. In order to save computational resources, only the half domain was simulated by setting symmetry boundary conditions to the xz-plane. Due to the large deviation of the wave lengths, two different computational domains were used, for each case. For wave frequencies between 0.6 and 0.8 rad/s (which corresponds to wave lengths between 170 to 95 m, respectively) the domain had a total length of 600 m, while, for the wave frequencies between 0.85 and 1.20 rad/s (where the wave lengths range between 85 and 45 m), a domain of total length of 300 m is used.
The domain discretization is depicted in Figure 3. The discretization of the domain upwind of the configuration was based on the wave characteristics. Near the free surface, the domain was uniform along the z-axis, with a length that corresponded to 18 cells per wave height. The mesh along the wave propagation was also uniform, with a discretization that corresponded to at least 100 cells per wavelength in each domain for the studied wave frequencies. Furthermore, special care was taken near the sharp edges of the structure, as illustrated in Figure 3a, to account for the flow separation that was expected in these regions. Finally, the mesh coarsened as it reached the lateral and the outflow boundary, to make sure that no wave reflections occurred. In order to capture the viscous effects close to the wall, the first point was placed at 0.5 mm normal to the wall.
Figure 4 presents the grid independency study. The M2 grid has the aforementioned mesh characteristics and a total size of approximately 5 million cells. Additionally, three more grids were examined. The M1 grid consisted of 1.6 million cells. The main difference compared to the M2 grid is that it had 8 cells along the wave height and 50 cells at the smallest wavelength. Furthermore, the characteristic size of the mesh close to the wall is two times larger, compared to the previous mesh. The M2-w mesh had the same characteristics as the M2 grid; however, the first node from the wall was placed at 1 mm normal to the wall. Finally, the M3 had a total size of 8.5 million cells. The M3 grid differed from M2 in the characteristics of the mesh near the structure, as well as upwind of the structure in the direction of the wave propagation. The cells per wave height were the same as M2; the total cells per smallest wavelength were 200; the mesh near the structure had two times smaller characteristic size. From the independency study, the M2, M2-w and M3 produced similar results. The characteristics of the M2 grid were chosen for all simulations. The time discretization was chosen as 1000 timesteps for the corresponding wave period [33].
Figure 5 and Figure 6 depict the horizontal and vertical exciting forces on the moonpool floater of the 1st examined Configuration, derived by the presented potential flow formulation. The results are also compared with the corresponding exciting forces on the outer moonpool floater (i.e., Body 2) (see Figure 2a), without the presence of the inner moonpool floater (i.e., Body 1), as well as on Body 1, without the presence of Body 2 (i.e., isolated bodies). The results of the exciting forces are normalized by the term πρgα2(H/2).
The peculiar behavior of the horizontal exciting forces (see Figure 5) on the Body 1 and Body 2, in the neighborhood of ω ~ 1.74 rad/s and ω ~ 1.23 rad/s, respectively, was due to the resonance pitch oscillations of the interior water basin of each Body 1 and 2. The same behavior has been reported in [1] regarding moonpool floaters. It should be further noted that the peculiar behavior of the horizontal exciting force on the total Configuration 1 (i.e., floater of two coaxial, free-surface piercing toroidal cylinders) was moved at a lower wave frequency, i.e., ω ~ 1.13 rad/s. This transpose could be attributed to the annular water area that formed between the two toroidal cylinders, which resonated in a pitch at a lower wave frequency than the corresponding one of the independent moonpool floaters (i.e., Body 1 and Body 2). In addition, a secondary resonance of the horizontal exciting force on the total Configuration 1 was observed in the neighborhood of ω ~ 1.74 rad/s. In the latter wave frequency, the resonance pitch oscillations of the interior water basin inside the Body 1 occurred.
With regard to the vertical forces (see Figure 6) on Body 1 and Body 2, resonances appeared in the vicinity of ω ~ 1.05 rad/s and ω ~ 0.7 rad/s, respectively. These peaks were due to the pumping resonance of the fluid motion in the interior water area of the moonpools Body 1 and Body 2. Concerning the vertical exciting force on the total Configuration 1, it can be seen (Figure 6) that its pumping resonant wave frequency was located in between the corresponding pumping wave frequencies of Body 1 and Body 2 (i.e., ω ~ 0.8 rad/s).
Next, the horizontal and vertical wave exciting forces on the examined Configuration 2 (see Figure 2b), derived from the theoretical formulation, are presented. Specifically, in Figure 7 and Figure 8 the values of Fx/(πρgα2(H/2)) and Fz/(πρgα2(H/2) on the Configuration 2 are depicted, respectively. The results are also compared with the corresponding exciting forces on the isolated Body 1 and Body 2, when the latter were considered alone in the wave field (i.e., Body 1 was assumed alone in the wave field without the presence of Body 2, and vice versa for Body 2).
The resonance of the horizontal exciting force on the Configuration 2 (see Figure 7), which occurred in the neighborhood of ω ~ 1.125 rad/s due to the pitch fluid motion in the annular fluid area between the interior cylinder and the exterior toroidal body, was moved to ω ~ 1.24 rad/s for the Body 2 case, when the latter was assumed isolated in the wave field. As far as the forces on the isolated Body 1 are concerned, it can be seen that they followed a smooth variation pattern at the examined wave frequency range, without any sharp maximizations.
Regarding the vertical exciting forces on the Configuration 2 (see Figure 8), the pumping resonance of the fluid motion in the annular area was depicted at ω ~ 0.85 rad/s [20], whereas, for the isolated Body 2, this pumping resonance lied near ω ~ 0.7 rad/s (see discussion of Figure 5). It should also be noted that the smooth variation pattern of the vertical exciting force on the isolated cylinder follows (i.e., Body 1).
Comparing Figure 5 and Figure 7, it can be concluded that both configurations 1 and 2 attained, in the non-resonance regime, comparable results concerning the horizontal exciting wave forces. Furthermore, their resonance pitch oscillations occurred at the vicinity of ω ~ 1.13 rad/s. Hence, it could be derived that the wave motion inside the pond of Body 1 in the first configuration did not seem to affect the floater’s horizontal forces. On the other hand, the secondary resonance of the horizontal exciting wave force on the Configuration 1 (i.e., ω~1.74 rad/s) was not observed in the case of Configuration 2. Regarding the vertical exciting forces Fz/(πρgα2(H/2)), see Figure 6 and Figure 8; in the case of Configuration 2, due to the truncated inner cylindrical body, they attained higher values compared to the vertical loads on Configuration 1. Hence, assuming both Configurations as heaving point absorbers for wave energy exploitation, Configuration 2 seems to be more efficient compared to Configuration 1.
In the sequel, the results from the theoretical formulation were supplemented by numerical ones, which took into consideration the viscous effects near the sharp edges of the examined moonpool-type floaters. The numerical results presented follow the configuration terminology described in the previous section. A grid consisting of 5 million cells was employed with symmetry conditions (half of the moonpool is modelled). Since the moonpool was considered fixed and the waves were essentially in the linear region, the flow is assumed laminar [29,44].
Regarding Configuration 1 the comparison of the non-dimensional horizontal and vertical exciting forces could be seen in Figure 9 and Figure 10. The illustration of the horizontal forces (Figure 9) suggested that the comparison was relatively good for both methods. CFD predictions follow the analytical method’s trends with minor discrepancies. In the proximity of the peak excitation frequency (between 1.0–1.2 rad/s), the differences between the predictions became larger. More specifically, in the near peak frequency region, CFD results predicted earlier the peak of the Body 1 than the analytical method, leading to a rise in total force. At 1.13 rad/s, both methods predicted the maximum horizontal excitation force. However, it was evident that the amplitude of the CFD predictions was significantly smaller. This could be justified by both the phase difference of the peak excitation loads between the two bodies in the CFD predictions (discussed in the following) induced by the viscous effects that were taken into account in the CFD computations.
Regarding the vertical excitation force in Figure 10, the comparison between the methodologies was fair, with the CFD simulations predicting smaller peak loads and shifted towards the lower frequencies. This has also been noted in [29]. Nevertheless, the qualitative comparison for both the total force and the individual body forces was considered to be in good agreement.
The respective comparison for Configuration 2 (external toroidal cylinder with internal coaxial truncated cylinder) can be found in Figure 11 and Figure 12. The agreement between the numerical results is good in both the horizontal and vertical forces. As Figure 11 suggests, CFD predictions followed the analytical method trends; nevertheless, the peak excitation frequency was predicted earlier. As opposed to Configuration 1, the individual body forces peak excitation frequencies were predicted in the same frequency, suggesting that the forces on the two bodies were in phase. This was also the case for the analytical method’s predictions. Finally, regarding the vertical force in Figure 12, again the CFD predicted the peak excitation frequency earlier; however, as in Configuration 1, the predicted trends were the same.
In order to gain a better insight of this phase difference between each body’s horizontal force in Configuration 1, instantaneous free surface elevation contours at four snapshots of the wave period are presented in Figure 13 at ω = 1.1 rad/s.
As Figure 13 suggests, the sloshing mode was excited between the outer (Body 1) and the inner (Body 2) cylinder, which explains the phase difference between the horizontal force of each individual body.
Figure 14 presents the free surface elevations inside the moonpool (between Body 1 and 2) for both configurations with respect to the incident wave frequency. The space averaged free surface elevation inside the moonpool is presented, as well as in three circumferential positions at 0, 90 and 180 degrees (as shown in Figure 13a). In both configurations up to ω = 1 rad/s, it was evident that the free surface elevation did not greatly vary between the various spatial positions, suggesting that the piston mode dominated the free surface elevation inside the moonpool. As the frequency increased beyond 1 rad/s the gauge at 90 degrees followed the spaced averaged elevation; however, the other two locations at 0 and 180 degrees clearly deviated (they increase in amplitude), thus indicating the appearance of a sloshing mode inside the moonpool. As the radial frequency was further increased, this mode dampens out. Furthermore, it is worth noting that peak frequencies of the piston mode resulted in an increase of the vertical forces, while the sloshing mode resulted in the resonance of the horizontal excitation forces.
This is further illustrated in Figure 15, where the phase angle of the signal in the three wave gauges is presented. In smaller frequencies, the three signals had no phase difference; however, in larger frequencies, the phases started to deviate. In the case of ω = 1.1 rad/s, the signal in the upwind gauge (0 deg) and in the downwind gauge (180 deg) had a 90 degrees phase difference compared to the middle gauge (90 deg) and 180 degrees between them. In higher frequencies, there were still some phase differences; however, the amplitude of free surface was relatively small (see Figure 13).
Lastly, Figure 16 presents the y-vorticity field in the symmetry plane in the case of the Configuration 1 for different wave frequencies. Close to the resonance of the piston mode, excessive flow separation was noted near the sharp edges of the structure (Figure 16b). Due to the loss of energy in this region of recirculating flow, the viscous solver predicted smaller amplitude of the piston mode (and thus of the vertical forces) compared to the one of the analytical solution. Contrary to Figure 16b, the vorticity field in larger frequencies (where the amplitude of the piston mode is reduced significantly), the vortex shedding was not significant; thus, minor discrepancies were noted between the two approaches. It was also noted that the performance of inviscid potential models is influenced by the radiation effect of the inner water area, leading to discrepancies between the viscous solver and the analytical solution.

5. Conclusions

A semi-analytical method has been developed to solve the linearized diffraction problem of two types of moonpool floaters, namely: (a) two coaxial toroidal cylinders configuration and (b) a floater encompassed by an external toroidal cylinder and a coaxial internal truncated cylinder. The potential flow theoretical formulation was supplemented by a numerical modelling that considered the viscous effects near the sharp edges of the body as non-negligible.
The accuracy of the two presented methods, the theoretical formulation and the numerical approximation, have been validated in both examined Configurations 1 and 2 and in a wide range of wave frequencies, i.e., 0.6 rad/s to 1.2 rad/s. It can be concluded that, in the non-resonance regime, the comparisons between the two methods were excellent; hence, both models predicted accurate results. However, in the critical wave frequencies near resonance, the CFD modelling predicted the peak excitation frequency earlier, whereas the amplitude of the predictions was significant smaller than the corresponding ones from the theoretical analysis. This can be traced back to the viscous effects that were taken into account in the CFD computations, which altered both the phase difference of the peak excitation loads between the two bodies in the CFD predictions. Specifically, an excessive flow separation was noted in the resonance regime, near the sharp edges of the structure; hence, the viscous solver predicted smaller force amplitudes. On the other hand, in larger frequencies the vortex shedding was negligible; thus, minor discrepancies were noted between the two approaches.
The present research will be continued further by comparing the accuracy of both methodologies to arrays of floating bodies in which the exciting forces on the structures are affected not only by the shape of the floaters but also from their position in the array towards the incoming waves.

Author Contributions

Conceptualization, S.A.M. and G.P.P.; methodology, A.S.M. and D.G.N.; software, S.A.M., D.G.N. and G.P.P.; validation, D.G.N. and A.S.M.; investigation, A.S.M., D.G.N., G.P.P. and D.N.K.; writing—original draft preparation, D.N.K.; writing—review and editing, S.A.M., G.P.P., D.G.N. and A.S.M.; visualization, D.N.K. and D.G.N.; supervision, S.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Greek General Secretariat for Research and Innovation, Program Reward 2019, Title: Marine Structures for the exploitation of renewable energy sources and dissemination of results.

Acknowledgments

MaPFlow computations were supported by computational time granted from the Greek Research and Technology Network (GRNET) in the National HPC facility—ARIS—under project “SHIPFLOW” with ID pr010039_thin.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Expressions of the functions R m n I I ( r ) ,   R m n * I I ( r ) in Equation (22):
R m n I I = K m ( n π b d h 1 ) I m ( n π r d h 1 ) K m ( n π r d h 1 ) I m ( n π b d h 1 ) K m ( n π b d h 1 ) I m ( n π a d h 1 ) K m ( n π a d h 1 ) I m ( n π b d h 1 ) ,   n 0
R m n * I I = K m ( n π r d h 1 ) I m ( n π a d h 1 ) K m ( n π a d h 1 ) I m ( n π r d h 1 ) K m ( n π b d h 1 ) I m ( n π a d h 1 ) K m ( n π a d h 1 ) I m ( n π b d h 1 ) ,   n 0
whereas, for n = 0 the terms R m 0 I I ,   R m 0 * I I , equal to:
R m 0 I I = ( r / b ) m ( b / r ) m ( a / b ) m ( b / a ) m , R m 0 * I I = ( a / r ) m ( r / a ) m ( a / b ) m ( b / a ) m   ,   m 0 ,   n = 0
and for n = 0 ,   m = 0 it holds:
R 00 I I = ln ( r / b ) ln ( a / b ) , R 00 * I I = ln ( a / r ) ln ( a / b ) ,   m = 0 ,   n = 0
The functions R m i I I I ( r ) , R m i * I I I ( r ) in Equation (23), are given by:
R m i I I I = K m ( a i c ) I m ( a i r ) K m ( a i r ) I m ( a i c ) K m ( a i c ) I m ( a i b ) K m ( a i b ) I m ( a i c ) ,   i 0
R m i * I I I = K m ( a i r ) I m ( a i b ) K m ( a i b ) I m ( a i r ) K m ( a i c ) I m ( a i b ) K m ( a i b ) I m ( a i c ) , i 0
Here a i are the roots of the transcendental equation: ω 2 + a i g tan ( a i d ) = 0 .
Similar to Equations (A1) and (A2), the expressions of the functions R m s I V ( r ) ,   R m s * I V (r) in Equation (24) are given by:
R m s I V = K m ( s π l d h 2 ) I m ( s π r d h 2 ) K m ( s π r d h 2 ) I m ( s π l d h 2 ) K m ( s π l d h 2 ) I m ( s π c d h 2 ) K m ( s π l d h 2 ) I m ( s π c d h 2 ) ,   s 0
R m s * I V = K m ( s π r d h 2 ) I m ( s π c d h 2 ) K m ( s π c d h 2 ) I m ( s π r d h 2 ) K m ( s π l d h 2 ) I m ( s π c d h 2 ) K m ( s π l d h 2 ) I m ( s π c d h 2 ) ,   s 0
R m 0 I V = ( r / l ) m ( l / r ) m ( c / l ) m ( l / c ) m ,     R m 0 * I V = ( c / r ) m ( r / c ) m ( c / l ) m ( l / c ) m   ,   m 0 ,   s = 0
R 00 I V = ln ( r / l ) ln ( c / l ) ,     R 00 * I V = ln ( c / r ) ln ( c / l ) ,   m = 0 ,   s = 0

Appendix B

The terms L s q , L s i , L n i , L n j defined in Equations (34)–(37) equal to:
L s q = 1 d h 2 0 d h 2 Z q ( z ) cos ( s π z d h 2 ) d z = ( 1 ) s N q 1 / 2 a q ( d h 2 ) a q 2 ( d h 2 ) 2 s 2 π 2 sin ( a q ( d h 2 ) )   for   q 1
L s i = 1 d h 2 0 d h 2 Z i ( z ) cos ( s π z d h 2 ) d z = ( 1 ) s N i 1 / 2 a i ( d h 2 ) a i 2 ( d h 2 ) 2 s 2 π 2 sin ( a i ( d h 2 ) )   for   i 1
L s 0 = 1 d h 2 0 d h 2 Z 0 ( z ) cos ( s π z d h 2 ) d z = ( 1 ) s N 0 1 / 2 k ( d h 2 ) k 2 ( d h 2 ) 2 + s 2 π 2 sin h ( k ( d h 2 ) )   for   i , q = 0
L n i = 1 d h 1 0 d h 1 Z i ( z ) cos ( n π z d h 1 ) d z = ( 1 ) n N i 1 / 2 a i ( d h 1 ) a i 2 ( d h 1 ) 2 n 2 π 2 sin ( a i ( d h 1 ) )   for   i 1
L n j = 1 d h 1 0 d h 1 Z j ( z ) cos ( n π z d h 1 ) d z = ( 1 ) n N j 1 / 2 a j ( d h 1 ) a j 2 ( d h 1 ) 2 n 2 π 2 sin ( a j ( d h 1 ) )   for   j 1
L n 0 = 1 d h 1 0 d h 1 Z 0 ( z ) cos ( n π z d h 1 ) d z = ( 1 ) n N 0 1 2 k ( d h 1 ) k 2 ( d h 1 ) 2 + n 2 π 2 sin h ( k ( d h 1 ) )   for   j , i = 0
In Equations (A11)–(A16) the terms N l ,   N 0 equal to:
N l = 1 2 [ 1 + sin ( 2 a l d ) 2 a l d ] , N 0 = 1 2 [ 1 + sin h ( 2 k d ) 2 k d ]
Also, k is related to the wave frequency ω through the dispersion equation and a l ,   l > 0 are the positive real roots of ω 2 + a l g tan ( a j d ) = 0 .
The terms Λ j I ,   Λ m n I I , Λ m n * I I in Equation (38) equal to:
Λ j I = a K m ( a j r ) | r = a K m ( a j a ) ,               Λ m n I I = a R m n I I r | r = a , Λ m n * I I = a R m n * I I r | r = a
Moreover, the terms D m n I I , D m n * I I , Λ m i I I I , Λ m i * I I I in Equation (39) can be written as:
D m n I I = b R m n I I r | r = b ,   D m n * I I = b R m n * I I r | r = b ,   Λ m i I I I = b R m i I I I r | r = b ,   Λ m i * I I I = b R m i * I I I r | r = b
The terms D m i I I I , D m i * I I I , Λ m s I V , Λ m s * I V of Equation (40) equal to:
D m i I I I = c R m i I I I r | r = c ,   D m i * I I I = c R m i * I I I r | r = c , Λ m s I V = c R m s I V r | r = c ,   Λ m s * I V = c R m s * I V r | r = c
and the terms Λ q V ,   D m s I V , D m s * I V of Equation (41) can be written as:
D m s I V = l R m s I V r | r = l ,   D m s * I V = l R m s * I V r | r = l , Λ q V = l I m ( a q r ) | r = l I m ( a q l )
In (A18) and (A21) primes denote derivatives with respect to the radius

Appendix C

The first order horizontal exciting wave forces on the body, based on Equation (42) can be written as:
F 1 = ω 2 ρ H 2 a e i ω t m = 0 ε m i m d h 1 d 0 2 π Ψ D , m I ( r , z ) cos ( m θ ) cos θ d θ d z + + ω 2 ρ H 2 b e i ω t m = 0 ε m i m d h 1 d 0 2 π Ψ D , m I I I ( r , z ) cos ( m θ ) cos θ d θ d z ω 2 ρ H 2 c e i ω t m = 0 ε m i m d h 2 d 0 2 π Ψ D , m I I I ( r , z ) cos ( m θ ) cos θ d θ d z + + ω 2 ρ H 2 l e i ω t m = 0 ε m i m d h 2 d 0 2 π Ψ D , m V ( r , z ) cos ( m θ ) cos θ d θ d z = = 2 i k d ( π ρ g a 2 H 2 ) tan h ( k d ) { N 0 1 2 ( k a ) ( F 10 I b a F 10 I I I ) [ sin h ( k d ) sin h ( k ( d h 1 ) ) ] + + j = 1 J 1 N j 1 / 2 ( α j a ) F 1 j I [ s i n ( α j d ) s i n ( α j ( d h 1 ) ) ] b a i = 1 I 1 N i 1 / 2 ( α i a ) F 1 i I I I [ sin ( a i d ) sin ( a i ( d h 1 ) ) ] + ( c a ) 2 N 0 1 2 ( k c ) ( F 10 III l c F 10 V ) [ sin h ( k d ) sin h ( k ( d h 2 ) ) ] + + ( c a ) 2 i = 1 I 1 N i 1 / 2 ( α i c ) F 1 i I I I [ sin ( α i d ) sin ( α i ( d h 2 ) ) ] ( c a ) 2 l c q = 1 Q 1 N q 1 / 2 ( α q c ) F 1 q V [ sin ( α q d ) sin ( α q ( d h 2 ) ) ] }
The first three terms in the parenthesis express the horizontal exciting wave force on the external torus, while the rest three the corresponding force imposed on the internal toroidal cylinder.
Similarly, the vertical exciting wave forces on the body equal to:
F 3 = ω 2 ρ H 2 e i ω t m = 0 ε m i m b a 0 2 π Ψ D , m I I ( r , z ) cos ( m θ ) ( 1 ) r d r d θ ω 2 ρ H 2 e i ω t m = 0 ε m i m l c 0 2 π Ψ D , m I V ( r , z ) cos ( m θ ) ( 1 ) r d r d θ = = 2 k d ( π ρ g a 2 H 2   ) tan h ( k d )   { 1 2 F 00 I I ( 1 1 ( b / a ) 2 2 l n ( a / b ) ) 1 2 F 00 * I I ( ( b a ) 2 1 ( b / a ) 2 2 l n ( a / b ) ) + + 2 n = 1 N 1 ( 1 ) n ( d h 1 n π a ) 2 [ F 0 n I I ( Λ 0 n I I D 0 n I I ) + F 0 n * I I ( Λ 0 n * I I D 0 n * I I ) ] + 1 2 F 00 I V ( c a ) 2 ( 1 1 ( l / c ) 2 2 l n ( c / l ) ) 1 2 F 00 * I V ( c a ) 2 ( ( l c ) 2 1 ( l / c ) 2 2 l n ( c / l ) ) + + 2 ( c a ) 2 s = 1 S 1 ( 1 ) s ( d h 2 s π c ) 2 [ F 0 s I V ( Λ 0 s I V D 0 s I V ) + F 0 s * I I ( Λ 0 s * I V D 0 s * I V ) ]   }
The first three terms in the parenthesis express the vertical exciting wave force on the external torus, while the rest three the corresponding force imposed on the internal toroidal cylinder.
The overturning moment about a horizontal axis lying at an arbitrary distance z = e from the seabed arises from the pressure distribution on the body’s vertical walls and on its bottom surface, based on Equation (43), i.e.,
M 1 = ω 2 ρ H 2 a e i ω t m = 0 ε m i m d h 1 d 0 2 π Ψ D , m I ( r , z ) cos ( m θ ) cos θ ( z e ) d θ d z + + ω 2 ρ H 2 b e i ω t m = 0 ε m i m d h 1 d 0 2 π Ψ D , m I I I ( r , z ) cos ( m θ ) cos θ ( z e ) d θ d z ω 2 ρ H 2 c e i ω t m = 0 ε m i m d h 2 d 0 2 π Ψ D , m I I I ( r , z ) cos ( m θ ) cos θ ( z e ) d θ d z + + ω 2 ρ H 2 l e i ω t m = 0 ε m i m d h 2 d 0 2 π Ψ D , m V ( r , z ) cos ( m θ ) cos θ ( z e ) d θ d z = 2 i k d ( π ρ g a 3 H 2 ) tan h ( k d ) { N 0 1 2 ( k a ) 2 ( F 10 I b a F 10 I I I ) x [ k ( d e ) sin h ( k d ) k ( d h 1 e ) sin h ( k ( d h 1 ) ) cos h ( k d ) + cos h ( k ( d h 1 ) ) ] + + j = 1 J 1 N j 1 2 ( α j a ) 2 F 1 j I [ α j ( d e ) sin ( α j d ) α j ( d h 1 e ) sin ( α j ( d h 1 ) ) + + cos ( α j d ) cos ( α j ( d h 1 ) ) ] b a i = 1 I 1 N i 1 2 ( α i a ) 2 F 1 i I I I [ α i ( d e ) sin ( α i d ) α i ( d h 1 e ) sin ( α i ( d h 1 ) ) + cos ( α i d ) cos ( α i ( d h 1 ) ) ] + ( c a ) 3 N 0 1 2 ( k c ) 2 ( F 10 I I I l c F 10 V ) x [ k ( d e ) sin h ( k d ) k ( d h 2 e ) sin h ( k ( d h 2 ) ) cos h ( k d ) + cos h ( k ( d h 2 ) ) ] + + ( c a ) 3 i = 1 I 1 N i 1 2 ( α i c ) 2 F 1 i I I I [ α i ( d e ) sin ( α i d ) α i ( d h 2 e ) sin ( α i ( d h 2 ) ) + cos ( α i d ) cos ( α i ( d h 2 ) ) ] ( c a ) 3 l c q = 1 Q 1 N q 1 2 ( α q c ) 2 F 1 q V x [ α q ( d e ) sin ( α q d ) α q ( d h 2 e ) sin ( α q ( d h 2 ) ) + cos ( α q d ) cos ( α q ( d h 2 ) ) ] }
M 3 = ω 2 ρ H 2 e i ω t m = 0 ε m i m b a 0 2 π Ψ D , m I I ( r , z ) cos ( m θ ) ( 1 ) r 2 d r d θ ω 2 ρ H 2 e i ω t m = 0 ε m i m l c 0 2 π Ψ D , m I V ( r , z ) cos ( m θ ) ( 1 ) r 2 d r d θ = = 2 i k d ( π ρ g a 3 H 2   ) tan h ( k d )   { 1 4 F 10 I I ( 1 ( b a ) 2 ) + 1 4 F 10 * I I ( b a ) ( 1 ( b a ) 2 ) + + 2 n = 1 N 1 ( 1 ) n ( d h 1 n π a ) 2 . [ F 1 n I I ( Λ 1 n I I b a D 1 n I I 1 ) + F 1 n * I I ( Λ 1 n * I I + b a ( 1 D 1 n * I I ) ) ] + + 1 4 F 10 I V ( c a ) 3 ( 1 ( l c ) 2 ) + 1 4 F 10 * I V ( c a ) 3 ( l c ) ( 1 ( l c ) 2 ) + + 2 ( c a ) 3 s = 1 S 1 ( 1 ) s ( d h 2 s π c ) 2 . [ F 1 s I V ( Λ 1 s I V l c D 1 s I V 1 ) + F 1 s * I I ( Λ 1 s * I V + l c ( 1 D 1 s * I V ) ) ] }

References

  1. Mavrakos, S. Wave loads on a stationary floating bottomless cylindrical body with finite wall thickness. Appl. Ocean Res. 1985, 7, 213–224. [Google Scholar] [CrossRef]
  2. Molin, B. On the piston and sloshing modes in moonpools. J. Fluid Mech. 2001, 430, 27–50. [Google Scholar] [CrossRef]
  3. Zheng, S.; Gao, W.; Liu, Y.; Cheng, S. Wave diffraction from a truncated cylinder with a moonpool of arbitrary cross-section: A semi-analytical study. Ocean Eng. 2021, 237, 109573. [Google Scholar] [CrossRef]
  4. Wu, T.Y.; Yates, G.T. Finite-Amplitude Unsteady Slender-Body Flow Theory. In Proceedings of the 11th Symposium on Naval Hydrodynamics, London, UK, 28 March–2 April 1976. [Google Scholar]
  5. Mavrakos, S. Hydrodynamic characteristics of floating toroidal bodies. Ocean Eng. 1997, 24, 381–399. [Google Scholar] [CrossRef]
  6. Mavrakos, S.A. Hydrodynamic Characteristics of Two Concentric Surface-Piecing Floating Circular Cylinders. In Proceedings of the 11th International Congress of International Maritime Association of the Mediterranean, IMAM, Lisboa, Portugal, 25–30 September 2005. [Google Scholar]
  7. Garrett, C.J.R. Bottomless harbours. J. Fluid Mech. 1970, 43, 433–449. [Google Scholar] [CrossRef]
  8. Miloh, T. Wave Loads on a Floating Solar Pond. In Proceedings of the International Workshop on Ship and Platform Motions, Berkeley, CA, USA, 26–28 October 1983; pp. 110–131. [Google Scholar]
  9. Mavrakos, S. Hydrodynamic coefficients for a thick-walled bottomless cylindrical body floating in water of finite depth. Ocean Eng. 1988, 15, 213–229. [Google Scholar] [CrossRef]
  10. Liu, Y.; Yue, D.K.; Kim, M.-H. First- and second-order responses of a floating toroidal structure in long-crested irregular seas. Appl. Ocean Res. 1993, 15, 155–167. [Google Scholar] [CrossRef]
  11. Zhou, S.; Mitchell, L. Diffraction of ocean waves around a hollow cylindrical shell structure. Wave Motion 2009, 46, 78–88. [Google Scholar] [CrossRef] [Green Version]
  12. Zhou, H.; Zhang, H. Radiation and diffraction analysis of a cylindrical body with a moon pool. J. Hydrodyn. 2013, 25, 196–204. [Google Scholar] [CrossRef]
  13. Chen, X.B.; Liu, H.X.; Duan, W.Y. Semi-analytical solutions to wave diffraction of cylindrical structures with a moonpool with a restricted entrance. J. Eng. Math. 2014, 90, 51–66. [Google Scholar] [CrossRef]
  14. Liu, H.-X.; Chen, H.-L.; Zhang, L.; Zhang, W.-C.; Liu, M. Quadratic dissipation effect on the moonpool resonance. China Ocean Eng. 2017, 31, 665–673. [Google Scholar] [CrossRef]
  15. Mavrakos, S.A. Hydrodynamic coefficients in heave of two concentric surface-piercing truncated circular cylinder. Appl. Ocean. Res. 2004, 26, 84–97. [Google Scholar] [CrossRef]
  16. Mavrakos, S.A.; Chatjigeorgiou, I.K. Second-order hydrodynamic effects on an arrangement of two concentric truncated vertical cylinders. Mar. Struct. 2009, 22, 545–575. [Google Scholar] [CrossRef]
  17. Mavrakos, S.A.; Chatjigeorgiou, I.K. Second-Order Wave Induced Loads on Vertical Bodies of Revolution. In Marine Technology and Engineering; Carlos, G., Soares, Y., Garbatov, N., Fonseca, A.P., Teixeira, Eds.; CRC Press: Boca Raton, FL, USA; Taylor & Francis Group: London, UK, 2011; pp. 479–501. ISBN 978-0-415-69808-5. [Google Scholar]
  18. Cho, I.; Kim, M. Hydrodynamic performance evaluation of a wave energy converter with two concentric vertical cylinders by analytic solutions and model tests. Ocean Eng. 2016, 130, 498–509. [Google Scholar] [CrossRef]
  19. Kong, F.; Liu, H.; Su, W.; Ao, J.; Chen, H.; Jing, F. Analytical and Numerical Analysis of the Dynamics of a Moonpool Platform–Wave Energy Buoy (MP–WEB). Energies 2019, 12, 4083. [Google Scholar] [CrossRef] [Green Version]
  20. Abramson, H.N. The Dynamic Behavior of Liquids in Moving Containers; NASA SP-106; Scientific and Technical Information Division, National Aeronautics and Space Administration: Washington, DC, USA, 1996.
  21. Shipway, B.J.; Evans, D.V. Wave Trapping by Axisymmetric Concentric Cylinders. In Proceedings of the 21st International Conference on Offshore Mechanics and Artic Engineering (OMAE 2002), Oslo, Norway, 23–28 June 2002. [Google Scholar]
  22. McIVER, P.; Newman, J.N. Trapping structures in the three-dimensional water-wave problem. J. Fluid Mech. 2003, 484, 283–301. [Google Scholar] [CrossRef] [Green Version]
  23. Lu, L.; Teng, B.; Sun, L.; Chen, B. Modelling of multi-bodies in close proximity under water waves—Fluid forces on floating bodies. Ocean Eng. 2011, 38, 1403–1416. [Google Scholar] [CrossRef]
  24. Lu, L.; Teng, B.; Cheng, L.; Sun, L.; Chen, X. Modelling of multi-bodies in close proximity under water waves—Fluid resonance in narrow gaps. Sci. China Ser. G Phys. Mech. Astron. 2010, 54, 16–25. [Google Scholar] [CrossRef]
  25. Kristiansen, T.; Faltinsen, O.M. Gap Resonances Analyzed by a Domain-Decomposition Method. In Proceedings of the 26th International Workshop on Water Waves and Floating Bodies, Athens, Greece, 17–20 April 2011. [Google Scholar]
  26. Kristiansen, T.; Faltinsen, O.M. Gap resonances analyzed by a newdomain-decomposition method combining potential and viscous flow. Appl. Ocean. Res. 2012, 34, 198–208. [Google Scholar] [CrossRef]
  27. Fredriksen, A.G.; Kristiansen, T.; Faltinsen, O.M. Experimental and numerical investigation of wave resonance in moonpools at low forward speed. Appl. Ocean Res. 2014, 47, 28–46. [Google Scholar] [CrossRef]
  28. Tan, L.; Lu, L.; Tang, G.-Q.; Cheng, L.; Chen, X.-B. A viscous damping model for piston mode resonance. J. Fluid Mech. 2019, 871, 510–533. [Google Scholar] [CrossRef]
  29. Ntouras, D.; Manolas, D.; Papadakis, G.; Riziotis, V. Exploiting the limit of BEM solvers in moonpool type floaters. J. Phys. Conf. Ser. 2020, 1618, 082059. [Google Scholar] [CrossRef]
  30. Mavrakos, A.S. Hydrodynamic Parametric Investigation of the Chamber’s Geometry on the Efficiency of Floating OSCILLATING Water Column (OWC) Wave Energy Devices. Master′ Thesis, School of Naval Architecture and Marine Engineering, National Technical University of Athens, Athens, Greece, July 2019. [Google Scholar]
  31. Wehausen, J.V. The Motion of Floating Bodies. Annu. Rev. Fluid Mech. 1971, 3, 237–268. [Google Scholar] [CrossRef]
  32. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions, 9th ed.; Dover Publications: Mineola, NY, USA, 1970. [Google Scholar]
  33. Ntouras, D.; Papadakis, G. A Coupled Artificial Compressibility Method for Free Surface Flows. J. Mar. Sci. Eng. 2020, 8, 590. [Google Scholar] [CrossRef]
  34. Hirt, C.W.; Nichols, B.D. Volume of fluid (VOF) method for the dynamics of free boundaries. J. Comput. Phys. 1981, 39, 201–225. [Google Scholar] [CrossRef]
  35. Chorin, A.J. A Numerical Method for Solving Incompressible Viscous Flow Problems. J. Comput.-Tional Phys. 1967, 2, 12–26. [Google Scholar] [CrossRef]
  36. Nichols, D. Development of a Free Surface Method Utilizing an Incompressible Multi-Phase Algorithm to Study the Flow about Surface Ships and Underwater Vehicles. Ph.D. Thesis, Department of Aerospace Engineering, Mississippi State University, Starkville, MS, USA, 2002. [Google Scholar]
  37. Kunz, R.F.; Boger, D.A.; Stinebring, D.R.; Chyczewski, T.S.; Lindau, J.W.; Gibeling, H.J.; Venkateswaran, S.; Govindan, T. A preconditioned Navier–Stokes method for two-phase flows with application to cavitation prediction. Comput. Fluids 2000, 29, 849–875. [Google Scholar] [CrossRef]
  38. Peric, R. Minimizing Undesired Wave Reflection at the Domain Boundaries in Flow Simulations with Forcing Zones. Ph.D. Thesis, Technischen Universität Hamburg, Hamburg, Germany, 2019. [Google Scholar]
  39. Roe, P.L. Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys. 1981, 43, 357–372. [Google Scholar] [CrossRef]
  40. Kelecy, F.; Pletcher, R. The Development of a Free Surface Capturing Approach for Multidimensional Free Surface Flows in Closed Containers. J. Comput. Phys. 1997, 138, 939–980. [Google Scholar] [CrossRef]
  41. Biedron, R.; Vatsa, V.; Atkins, H. Simulation of Unsteady Flows Using an Unstructured Navier-Stokes Solver on Moving and Stationary Grids. In Proceedings of the 23rd AIAA Applied Aerodynamics Conference, Toronto, ON, USA, 6–9 June 2005; p. 5093. [Google Scholar]
  42. Papadakis, G. Development of a Hybrid Compressible Vortex Particle Method and Application to External Problems including Helicopter Flows. Ph.D. Thesis, Department of Mechanical Engineering, National Technical University of Athens, Athens, Greece, November 2014. [Google Scholar]
  43. Diakakis, K. Computational Analysis of Transitional and Massively Separated Flows with Application to Wind Turbines. Ph.D. Thesis, Department of Mechanical Engineering, National Technical University of Athens, Athens, Greece, 2019. [Google Scholar]
  44. Fredriksen, A.G.; Kristiansen, T.; Faltinsen, O.M. Wave-induced response of a floating two-dimensional body with a moonpool. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2015, 373, 20140109. [Google Scholar] [CrossRef]
Figure 1. Two examined moonpool case studies exposed to the action of regular wave trains. Definitions and discretization of the flow field around the bodies: (a) two coaxial surface piercing moonpool vertical toroidal cylinders; (b) moonpool with a coaxial surface piercing truncated circular cylinder.
Figure 1. Two examined moonpool case studies exposed to the action of regular wave trains. Definitions and discretization of the flow field around the bodies: (a) two coaxial surface piercing moonpool vertical toroidal cylinders; (b) moonpool with a coaxial surface piercing truncated circular cylinder.
Energies 15 00570 g001
Figure 2. Schematic representation of the two moonpool-type case studies: (a) floater of two coaxial, free-surface piercing moonpool cylinders with vertical symmetry axes; (b) floater of an external toroidal cylinder and an internal coaxial truncated cylinder.
Figure 2. Schematic representation of the two moonpool-type case studies: (a) floater of two coaxial, free-surface piercing moonpool cylinders with vertical symmetry axes; (b) floater of an external toroidal cylinder and an internal coaxial truncated cylinder.
Energies 15 00570 g002
Figure 3. Computational mesh used for the CFD simulations, in case of Configuration 1. The mesh of the Configuration 2 has the corresponding characteristics; (a) symmetry plane; (b) z plane cut; (c) surface mesh.
Figure 3. Computational mesh used for the CFD simulations, in case of Configuration 1. The mesh of the Configuration 2 has the corresponding characteristics; (a) symmetry plane; (b) z plane cut; (c) surface mesh.
Energies 15 00570 g003aEnergies 15 00570 g003b
Figure 4. Grid Independency study. Forces acting on Configuration 2 during the 35th wave period. The wave excitation frequency is ω = 0.80 rad/s; (a) horizontal excitation force; (b) vertical excitation force.
Figure 4. Grid Independency study. Forces acting on Configuration 2 during the 35th wave period. The wave excitation frequency is ω = 0.80 rad/s; (a) horizontal excitation force; (b) vertical excitation force.
Energies 15 00570 g004
Figure 5. Horizontal exciting forces on the Configuration 1. The results are also compared with the corresponding forces on the isolated toroidal Bodies 1 and 2.
Figure 5. Horizontal exciting forces on the Configuration 1. The results are also compared with the corresponding forces on the isolated toroidal Bodies 1 and 2.
Energies 15 00570 g005
Figure 6. Vertical exciting forces on the Configuration 1. The results are also compared with the corresponding forces on the isolated toroidal Bodies 1 and 2.
Figure 6. Vertical exciting forces on the Configuration 1. The results are also compared with the corresponding forces on the isolated toroidal Bodies 1 and 2.
Energies 15 00570 g006
Figure 7. Horizontal exciting forces on the Configuration 2. The results are also compared with the corresponding forces on the isolated Bodies 1 and 2.
Figure 7. Horizontal exciting forces on the Configuration 2. The results are also compared with the corresponding forces on the isolated Bodies 1 and 2.
Energies 15 00570 g007
Figure 8. Vertical exciting forces on the Configuration 2. The results are also compared with the corresponding forces on the isolated Bodies 1 and 2.
Figure 8. Vertical exciting forces on the Configuration 2. The results are also compared with the corresponding forces on the isolated Bodies 1 and 2.
Energies 15 00570 g008
Figure 9. Comparison between the theoretical and CFD results of the non-dimensional horizontal excitation force, in case of the Configuration 1; (a) wave frequencies between 0.05 and 2 rad/s; (b) wave frequencies between 0.8 and 1.2 rad/s (detailed view).
Figure 9. Comparison between the theoretical and CFD results of the non-dimensional horizontal excitation force, in case of the Configuration 1; (a) wave frequencies between 0.05 and 2 rad/s; (b) wave frequencies between 0.8 and 1.2 rad/s (detailed view).
Energies 15 00570 g009
Figure 10. Comparison between the theoretical and CFD results of the non-dimensional vertical excitation force, in case of the Configuration 1; (a) wave frequencies 0.05 and 2 rad/s; (b) wave frequencies between 0.6 and 1.0 rad/s (detailed view).
Figure 10. Comparison between the theoretical and CFD results of the non-dimensional vertical excitation force, in case of the Configuration 1; (a) wave frequencies 0.05 and 2 rad/s; (b) wave frequencies between 0.6 and 1.0 rad/s (detailed view).
Energies 15 00570 g010
Figure 11. Comparison between the theoretical and CFD results of the non-dimensional horizontal excitation force, in case of the Configuration 2; (a) wave frequencies 0.05 and 2 rad/s; (b) wave frequencies between 0.8 and 1.2 rad/s (detailed view).
Figure 11. Comparison between the theoretical and CFD results of the non-dimensional horizontal excitation force, in case of the Configuration 2; (a) wave frequencies 0.05 and 2 rad/s; (b) wave frequencies between 0.8 and 1.2 rad/s (detailed view).
Energies 15 00570 g011
Figure 12. Comparison between the theoretical and CFD results of the non-dimensional vertical excitation force, in case of the Configuration 2; (a) wave frequencies 0.05 and 2 rad/s; (b) wave frequencies between 0.6 and 1.0 rad/s (detailed view).
Figure 12. Comparison between the theoretical and CFD results of the non-dimensional vertical excitation force, in case of the Configuration 2; (a) wave frequencies 0.05 and 2 rad/s; (b) wave frequencies between 0.6 and 1.0 rad/s (detailed view).
Energies 15 00570 g012
Figure 13. Snapshots of the non-dimensional surface elevation (ζ/(H/2)) during a wave period for Configuration 1 and ω = 1.1 rad/s; (a) t = 0T/4; (b); t = T/4; (c) t = 2T/4; (d) t = 3T/4.
Figure 13. Snapshots of the non-dimensional surface elevation (ζ/(H/2)) during a wave period for Configuration 1 and ω = 1.1 rad/s; (a) t = 0T/4; (b); t = T/4; (c) t = 2T/4; (d) t = 3T/4.
Energies 15 00570 g013aEnergies 15 00570 g013b
Figure 14. Analysis of the instantaneous free surface elevations inside the moonpool of Configuration 1 (a) and Configuration 2 (b) The position of the wave gauges is shown in Figure 13a. The non-dimensional amplitude of the free surface elevation (ζ/(H/2)) is plotted for the various wave frequencies.
Figure 14. Analysis of the instantaneous free surface elevations inside the moonpool of Configuration 1 (a) and Configuration 2 (b) The position of the wave gauges is shown in Figure 13a. The non-dimensional amplitude of the free surface elevation (ζ/(H/2)) is plotted for the various wave frequencies.
Energies 15 00570 g014
Figure 15. Analysis of the phase angle of the downwind (0 deg), middle (90 deg) and the upwind (180 deg) wave gauges for the two configurations. The position of the gauges is shown in Figure 13a; (a) Configuration 1; (b) Configuration 2.
Figure 15. Analysis of the phase angle of the downwind (0 deg), middle (90 deg) and the upwind (180 deg) wave gauges for the two configurations. The position of the gauges is shown in Figure 13a; (a) Configuration 1; (b) Configuration 2.
Energies 15 00570 g015
Figure 16. Snapshot of the y-vorticity field in the symmetry plane for two different wave frequencies, in case of Configuration 1; (a) vorticity field at the wave frequency where a resonant behavior of the piston mode is predicted (ω = 0.75 rad/s); (b) vorticity field at the wave frequency where a resonant behavior of the pitch sloshing mode is predicted (ω = 1.1 rad/s).
Figure 16. Snapshot of the y-vorticity field in the symmetry plane for two different wave frequencies, in case of Configuration 1; (a) vorticity field at the wave frequency where a resonant behavior of the piston mode is predicted (ω = 0.75 rad/s); (b) vorticity field at the wave frequency where a resonant behavior of the pitch sloshing mode is predicted (ω = 1.1 rad/s).
Energies 15 00570 g016
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mavrakos, A.S.; Konispoliatis, D.N.; Ntouras, D.G.; Papadakis, G.P.; Mavrakos, S.A. Hydrodynamics of Moonpool-Type Floaters: A Theoretical and a CFD Formulation. Energies 2022, 15, 570. https://doi.org/10.3390/en15020570

AMA Style

Mavrakos AS, Konispoliatis DN, Ntouras DG, Papadakis GP, Mavrakos SA. Hydrodynamics of Moonpool-Type Floaters: A Theoretical and a CFD Formulation. Energies. 2022; 15(2):570. https://doi.org/10.3390/en15020570

Chicago/Turabian Style

Mavrakos, Anargyros S., Dimitrios N. Konispoliatis, Dimitrios G. Ntouras, George P. Papadakis, and Spyros A. Mavrakos. 2022. "Hydrodynamics of Moonpool-Type Floaters: A Theoretical and a CFD Formulation" Energies 15, no. 2: 570. https://doi.org/10.3390/en15020570

APA Style

Mavrakos, A. S., Konispoliatis, D. N., Ntouras, D. G., Papadakis, G. P., & Mavrakos, S. A. (2022). Hydrodynamics of Moonpool-Type Floaters: A Theoretical and a CFD Formulation. Energies, 15(2), 570. https://doi.org/10.3390/en15020570

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop