Next Article in Journal
Assessing Coastal Road Flood Risk in Arctic Alaska, a Case Study from Hooper Bay
Next Article in Special Issue
Mean Drift Wave Forces on Arrays of Bodies Surrounded by Thin Porous Surfaces
Previous Article in Journal
Tailored Leaching Tests as a Tool for Environmental Management of Mine Tailings Disposal at Sea
Previous Article in Special Issue
A Semi-Analytical Model for Studying Hydroelastic Behaviour of a Cylindrical Net Cage under Wave Action
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrodynamics of a Moored Permeable Vertical Cylindrical Body

by
Dimitrios N. Konispoliatis
*,
Ioannis K. Chatjigeorgiou
and
Spyros A. Mavrakos
Laboratory for Floating Structures and Mooring Systems, Division of Marine Structures, School of Naval Architecture and Marine Engineering, National Technical University of Athens, 9 Heroon Polytechniou Avenue, GR 157-73 Athens, Greece
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(3), 403; https://doi.org/10.3390/jmse10030403
Submission received: 20 February 2022 / Revised: 7 March 2022 / Accepted: 8 March 2022 / Published: 10 March 2022
(This article belongs to the Special Issue Hydrodynamics of Fish Cages and Floating Permeable Structures)

Abstract

:
In this study, the problems of diffraction and radiation of water waves by a permeable vertical cylindrical body are formulated within the realm of the linear potential theory. The body, which is floating in constant water depth, is moored with a catenary mooring line system. The method of matched eigenfunction expansions for the prediction of the velocity potential in the fluid domain surrounding the body is applied. Furthermore, the static and dynamic characteristics of the mooring system are combined with the hydrodynamics of the body, to set up the coupled motion equations of the dynamical model, i.e., floater and mooring system, in the frequency domain. Numerical results obtained through the developed solution are presented. The results revealed that porosity plays a key role in reducing/controlling the exciting wave loads. As far as the mooring system is concerned, its quasi-static and dynamic characteristics, by employing several motion directions on the fairlead in accordance to varying environmental conditions, are examined, highlighting their effect on the body’s motions.

1. Introduction

Permeable floating structures have been widely applied by the marine sector to reduce the effect of incoming waves and to protect marine structures against the wave action, as they use their porous surface to decrease the transmission and reflection of wave heights. Hence, they become preferable to impermeable structures, due to their porosity, for applications such as harbor and shore protection [1,2,3,4]. Subsequently, several studies have followed concerning porous breakwaters and their capability in dissipating the wave energy, while minimizing the environmental impact [5,6,7,8,9,10,11].
Permeable floating structures are also related to aquaculture, which is gradually replacing ocean fishing. The shrinking availability of coastal sites, as well as the increased environmental impact of aquaculture, is forcing the latter into offshore areas, where the main challenge is to build a structure capable of withstanding the offshore severe environmental conditions, while being financially competitive with near-shore concepts. Cage farming has been practiced at an early phase for hundreds of years, initially in fresh water and later in seawater, whereas the development of modern cage systems has taken place in the last 30 years. Kawakami [12] was the first to evaluate the resistance of fish nets to currents using a simple analytical formula. Aarsnes et al. [13] calculated the current forces on cage systems and the deformation of nets on the basis of net-panel discretization with line finite elements in the plane of symmetry. Herein, an empirical formula for the drag coefficient of plane nets in a steady current was established. Continually, in [14], the stability and maneuverability problems of fishing gear were examined through the development of a dynamic study on submerged flexible reticulated surfaces. A multi-domain boundary element method based on Darcy’s fine pore model for the determination of the wave interaction with a horizontal porous flexible membrane was presented in [15]. Additionally, in [16] the scattering problem of a flexible fishnet, modeled as a porous, freely flexible, barrier displaced solely by hydrodynamic force and deformed similarly to a catenary, was studied. A nonlinear finite element method simulating the tension variation of a fishing net under a uniform current was presented in [17], whereas, in [18], a numerical calculation method for simulating the three-dimensional dynamic behavior of fishing nets was developed. Here, the fishing nets were modeled as a group of lumped point masses interconnected with weightless springs. In [19], a numerical model for the determination of the flow field around a 3D porous structure was presented. The model also calculated the responses of the structure taking into consideration the fluid velocity reduction due to the net. Bao et al. [20] applied a complicated dispersion equation for the evaluation of the eigenvalues of the fluid domain inside a semi-submerged porous circular cylinder. They derived an additional damping term created by the porosity.
In the recent decade, a wide range of studies on modeling the hydrodynamic problem of wave interaction with fish cages were presented, and several computational methods were developed. Indicatively, Faltinsen [21] presented a slender-body theory to simulate the wave effects on a floating fish cage with circular collar in the framework of potential flow theory without, however, considering current effects. Continually, in [22], a numerical model was developed to calculate viscous wave and current loads on aquaculture nets. Here, the net shape was determined by solving a time-stepping procedure which evaluates the unknown tensions at each timestep. The model was further developed in [23,24] to consider the induced wave and current mooring loads on a net cage. Furthermore, a cubic-shaped, elastic, fish cage was examined in [25] validating its mooring characteristics and behavior with experimental results. In [26,27], a hydroelastic analysis of surface gravity waves with floating flexible cage systems was presented considering the deflection of the porous and flexible structure, whereas, in [28], a semi-analytical model was developed describing the wave field around a submerged cylindrical fish net cage. Herein, the theory was based on linear approximation; hence, nonlinearities due to wave action or quadratic porous flow models were omitted. In addition, the wave–current effect on a moored flexible fish cage was studied in [29] through a semi-analytical model in the realm of the linearized wave theory, which was validated with numerical finite element model results. Recently, a semi-analytical model was proposed for wave interactions with a submerged flexible net cage governed by the membrane vibration equation [30].
Supplementary to aquaculture applications, permeable structures, considered as rigid bodies, were also examined as means of reducing/controlling the exciting wave loads. Specifically, in [31], a theoretical formulation was presented dealing with the wave interactions with arrays of bottom-seated, surface-piercing, circular cylinders. It was concluded that porosity significantly reduces the hydrodynamic loading, as well as the associated wave run up. The theoretical formulation was further developed in [32] to account for a freely floating circular cylinder, solving the corresponding radiation problem and concluding that the size of the porous region affects the hydrodynamic coefficients of the permeable body. In [33], the wave excitation forces on an array of freely floating, porous, circular cylinders were determined for various porous coefficients. In addition, a theoretical model, also validated experimentally, was developed in [34] to investigate the wave interactions with a truncated porous cylinder. Herein, the phenomenon of the sloshing mode inside the porous structure was observed. An array of porous circular cylinders with impermeable bottom, supplemented by horizontal, porous plates inside the bodies, was theoretically and experimentally examined in [35]. Furthermore, in [36], a 3D numerical model was developed describing the wave interactions with circular cylinders featuring a partially porous surface, leading to an optimal ratio of the porous portion to the impermeable portion, which minimizes the wave hydrodynamic impact. Furthermore, in [37,38], the effect of linear and quadratic resistance laws on porous bodies of arbitrary shape was examined. A porous bottom-seated vertical cylinder placed in front of a rigid vertical wall was studied in [39], in which a linear potential theory was assumed. It was concluded that the amplified wave effects due to the presence of the breakwater can be apparently reduced due to the presence of the porous cylinder. Recently, a CFD method was developed and experimentally validated in [40] concerning bottom-fixed porous cylinders. The method was compared with a boundary element model, applying a quadratic pressure drop condition, concluding that both methods can accurately replicate wave forces on porous cylinders.
The motivation of the present work is to evaluate the hydrodynamics of a moored, permeable body in waves on the basis of linear potential theory. The primary purpose of introducing this type of bodies is to dissipate the wave energy, effectively reducing the transmitted and reflected wave heights. Furthermore, permeable bodies are also related to aquaculture cages, the performance of which, in waves and currents, is strongly dependent on the mooring characteristics of the cage system. Therefore, in the present analysis, the problems of diffraction and radiation around a permeable vertical cylindrical body moored with a catenary mooring line system are formulated. The characteristics of the mooring system are combined with the hydrodynamics of the body, to set up the coupled motion equations of the dynamical model, i.e., floater and mooring system, in the frequency domain.
The organization of this paper is as follows: Section 2 formulates the solution of the corresponding semi-analytical diffraction and radiation problems. In Section 3, the coupled motion equations are formulated, whereas Section 4 is dedicated to the results and discussion. Lastly, the conclusions are presented in Section 5.

2. Material and Methods

2.1. Hydrodynamic Formulation

The examined vertical cylindrical body is assumed to consist of a thin impermeable bottom and a permeable sidewall of radius α . The distance of the permeable surface bottom to the seabed is signified as h 1 , whereas the distance of the bottom of the structure to the seabed is signified as h . The water depth is set to be constant, i.e., d , as shown in Figure 1. The permeable body is exposed to the action of regular linear waves of amplitude H / 2 , angular frequency ω , and wave number k , propagating toward the positive direction. A cylindrical coordinate system ( r , θ , z ) is adopted to describe the problem. The origin is located at the seabed, coinciding with the cylinder’s vertical symmetry axis. The fluid’s and body’s motions are assumed to be small, so that linearized diffraction and radiation problems can be considered. The water domain is subdivided into three regions: (a) fluid region I : r a ,   0 z d ; (b) fluid region I I : 0 r a ,   h 1 z d ; (c) fluid region I I I : 0 r a ,   0 z h . Moreover, it is assumed that the fluid is inviscid and incompressible, while the flow is irrotational and can be represented in each fluid domain in terms of velocity potential, Φ y = Re ( φ y e i ω t ) ,   y = I , I I , I I I .
The velocity potential φ y ( r , θ , z ) can be decomposed into three terms associated with the corresponding diffraction and radiation problems. It follows that [41]
φ y = φ 0 + φ 7 y + j = 1 , 3 , 5 x ˙ j 0 φ j y = φ D y + j = 1 , 3 , 5 x ˙ j 0 φ j y .
In Equation (1), the velocity potential of the undisturbed incident harmonic wave is denoted as φ 0 , whereas the scattered wave potential for the permeable cylinder when it is assumed restrained to the wave impact is φ 7 y ,   y = I , I I , I I I . The radiation potential resulting from the forced oscillations of the body in the j -th direction, j = 1, 3, 5, with unit velocity amplitude is denoted as φ j y ,   y = I , I I , I I I . Furthermore, the body’s complex velocity amplitude in the j -th direction is x ˙ j 0 . It should be noted that the sum of the undisturbed incident harmonic wave potential with the scattered wave potential is equal to the diffraction velocity potential, i.e., φ D y = φ 0 + φ 7 y ,   y = I , I I , I I I .
The velocity potential φ j y ,   j = 1 , 3 , 5 , 7 ; y = I , I I , I I I satisfies the Laplace equation as its governing equation. In the framework of linear wave theory, φ j y satisfies a homogeneous free surface boundary condition at z = d and an impermeable boundary condition at z = 0 . In addition, the radiation and scattering potentials, j = 1 , 3 , 5 , 7 , have to be satisfy the Sommerfeld radiation condition at the far field [42].
Since the cylinder’s sidewall is permeable, and the Reynolds numbers of the flow through the permeable surface are low, Darcy’s law can be employed [43]. It is stated that the normal flow velocity is continuous through the porous boundary and proportional to the pressure difference through the porous boundary [31]; hence, the boundary condition on the sidewall forms
φ j I I r = n j + i k G [ φ j I I φ j I ]   on   r = a , h 1 z d ,   j = 1 , 3 , 5 , D .
In Equation (2), / r denotes the derivative with respect to r , k stands for the wave number, and n j is the generalized normal vector defined as ( n 1 , n 2 , n 3 ) = n ;   ( n 4 , n 5 , n 6 ) = r × n . Here, n is the unit normal vector pointing outward, r is the position vector regarding the coordinate system origin, and × the cross-product symbol.
Furthermore, G denotes the complex dimensionless porous effect parameter. The parameter can be composed by G = G r e + i G i m , where G r e represents the linearized drag effect of the permeable sidewall, and G i m is the inertial effect. Hence, for a real G , the resistance effects dominate over the inertia effects, whereas G attains complex values when the inertia effects dominate over the resistance ones [43]. The parameter G is also a measure of the sidewall porous effect, i.e., for G = 0 , the sidewall is total impermeable, whereas, as G approaches infinity, i.e., G 0 , the sidewall is completely permeable to fluid (i.e., no sidewall exists) [44]. Following [34], the porous effect parameter can be connected to the opening rate τ of the sidewall material, as well as the waveslope ε = k H / 2 , through Equation (3).
G = ( 17.8 ε + 143.2 ) τ 2 2 π ( 1 + 1.06 τ ) .
In Equation (3), the opening rate τ is equal to the ratio between the area of the opening holes and the total area.
Regarding the boundary conditions that have to be fulfilled by the velocity potentials φ j y ,   j = 1 , 3 , 5 , 7 ; y = I , I I , I I I on the cylinder’s impermeable wetted surface S (i.e., S is the surface of the impermeable bottom of the cylindrical body), they are formulated as follows [41]:
φ 7 I n = φ 0 n   on   S ,
φ j I n = n j   on   S ,   j = 1 ,   3 ,   5 .
The n j term in Equation (5) was defined in Equation (2).
In addition, the matching conditions at the interface between the fluid regions (see Figure 1) should be satisfied, which are given by
φ j I = φ j I I I ; r = a ; 0 z d ; j = 1 ,   3 ,   5 , D ,
φ j I r = φ j I I r ; r = a ;     h 1 z d ;   j = 1 ,   3 ,   5 , D ,
φ j I r = φ j I I I r ; r = a ;   0 z h ;   j = 1 ,   3 ,   5 , D .
The velocity potential of the undisturbed incident wave propagating along the positive x-axis can be expressed as follows [41]:
φ 0 ( r , θ , z ) = i ω H 2 cosh ( k z ) k sinh ( k d ) m = 0 ε m i m J m ( k r ) cos ( m θ ) .
In Equation (9), ε m denotes the Neumann’s symbol defined as ε 0 = 1 and ε m = 2 for m > 0 , and J m is the m -th order Bessel function of first kind.
Similar to Equation (9), the diffraction velocity potential around the permeable body can be obtained as
φ D y ( r , θ , z ) = i ω H 2 m = 0 ε m i m Ψ D m y ( r , z ) cos ( m θ ) ,   y = I , I I , I I I .
Furthermore, the radiation potentials φ j y ,   j = 1 , 3 , 5 can be expressed as follows [41]:
φ 1 y ( r , θ , z ) = Ψ 11 y ( r , z ) cos ( θ ) ,   φ 3 y ( r , θ , z ) = Ψ 30 y ( r , z ) ,   φ 5 y ( r , θ , z ) = Ψ 51 y ( r , z ) cos ( θ ) .
The functions Ψ D m y , Ψ 11 y , Ψ 30 y , Ψ 51 y in Equations (10) and (11) denote the principal unknowns of the corresponding diffraction and radiation problems. Here, the superscript indicates the fluid domain, y = I , I I , I I I . Moreover, the first subscript stands for the respective boundary value problem, D ,   1 ,   3 ,   5 , and the second stands for the number of m modes, which are applied in the solution procedure.
Equations (2) and (4)–(8) provide sufficient information for the treatment of the hydrodynamic problems (i.e., diffraction and radiation problems) in each fluid domain. Applying the method of separation of variables, the Laplace differential equation can be solved, and appropriate representations of the Ψ D m y , Ψ 11 y , Ψ 30 y , Ψ 51 y functions in each fluid domain can be established. The complete solution is obtained by applying the kinematic condition on the impermeable wetted surface, the porous boundary condition on the sidewall, and the matching relations (see Equations (6)–(8)) on the common cylindrical boundaries of the discrete fluid regions.
According to the presented formulation and similarly to [42], the following expressions for the Ψ l m y ,   l = D , 1 , 3 , 5 terms can be derived for the description of the induced flow field around the permeable cylindrical body.
Infinite fluid region, I : r a ,   0 z d
  1 δ l Ψ l m I = f l m I + i = 0 F l m , i I K m ( a i I r ) K m ( a i I a ) Z i ( z ) .
In Equation (12), K m is the m -th order modified Bessel function of second kind, whereas δ l = d , for l = D , 1 , 3 , and δ 5 = d 2 . The term f l m I is equal to
f 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 ( z ) ,   f 11 I = f 30 I = f 51 I = 0 ,
where H m is the m -th order Hankel function of the first kind.
In Equations (12) and (13), the terms Z i ( z ) are orthonormal functions in [ 0 , d ] defined as follows [41]:
Z 0 ( z ) = N 0 , I 1 / 2 cosh ( k I z ) ,
Z i ( z ) = N i , I 1 / 2 cos ( a i I z ) ,
where
N 0 , I = 1 2 [ 1 + sinh ( 2 k I d ) 2 k I d ] ,
N i , I = 1 2 [ 1 + sin ( 2 a i I d ) 2 a i I d ] .
The wave number k I k is related to the wave frequency through the dispersion equation, whereas a i I ,   i > 0 are the positive real roots of
ω 2 g + a i I tan ( a i I d ) = 0 ,
where the superscript denotes the infinite fluid region I .
Hence, it is convenient to write a 0 I = i k , resulting in the following [45]:
K 0 ( i k r ) = π 2 i H 0 ( k r ) .
Fluid region, I I : 0 r a ,   h 1 z d
  1 δ l Ψ l m I I = f l m I I + i = 0 F l m , i I I I m ( a i I I r ) I m ( a i I I a ) Z i ( z h 1 ) .
In Equation (20), I m is the m -th order modified Bessel function of first kind. In addition, the term Z i ( z h 1 ) is derived as follows [41]:
Z 0 ( z h 1 ) = N 0 , I I 1 / 2 cosh ( k I I ( z h 1 ) ) ,
Z i ( z h 1 ) = N i , I I 1 / 2 cos ( a i I I ( z h 1 ) ) ,
where:
N 0 , I I = 1 2 [ 1 + sinh ( 2 k I I ( d h 1 ) ) 2 k I I ( d h 1 ) ] ,
N i , I I = 1 2 [ 1 + sin ( 2 a i I I ( d h 1 ) ) 2 a i I I ( d h 1 ) ] .
Furthermore, it holds that
f D m I I = f 11 I I = 0 ,     f 30 I I = z d 1 + g d ω 2 ,   f 51 I I = r d 2 ( ( z d ) + g ω 2 ) .
The terms a i I I are the roots of
ω 2 g + a i I I tan ( a i I I ( d h 1 ) ) = 0 ,
with the imaginary a i I I = i k I I considered as first.
Fluid region, I I I : 0 r a ,   0 z h
  1 δ l Ψ l m I I I = f l m I I I + n = 0 ε m F l m , n I I I I m ( n π h r ) I m ( n π h a ) cos ( n π h z ) .
In Equation (27), the term f l m I I I is equal to
f D m I I I = f 11 I I I = 0 ,     f 30 I I I = z 2 0.5 r 2 2 h d ,   f 51 I I I = r ( z 2 0.25 r 2 ) 2 h d 2 .
The functions Ψ l m y ,   l = D , 1 , 3 , 5 have the advantage of being expressed by simple Fourier series representations, F l m , i I ,   F l m , i I I ,   F l m , n I I I , of all the types of ring regions. The system of equations for the unknown Fourier coefficients is derived by fulfilling the kinematic conditions at the vertical walls (i.e., permeable, and impermeable surfaces), as well as by the requirement for continuity of the potential and its first derivative. The formulation was described thoroughly in [42]. Hence, it is not further elaborated here.

2.2. Hydrodynamic Forces

The various hydrodynamic forces on the permeable cylindrical body are calculated by the pressure distribution given by the linearized Bernoulli’s equation. Thus, it can be written that
F 1 = 2 π i ρ ω 2 d H 2 a [ h d 1 d Ψ D 1 I ( a , z ) d z h 1 d 1 d Ψ D 1 I I ( a , z ) d z ] .
Substituting Equations (12) and (20) into Equation (29), we get
F 1 B = 2 i k d tan h ( k d ) [ N 0 , I 1 2 F D 1 , 0 I sinh ( k d ) sinh ( k h ) k a + i = 1 N i , I 1 2 F D 1 , i I sin ( a i I d ) sin ( a i I h ) a i I a N 0 , I I 1 2 F D 1 , 0 I I sinh ( k I I ( d h 1 ) ) k I I a i = 1 N i , I I 1 2 F D 1 , i I I sin ( a i I I ( d h 1 ) ) a i I I a ] ,
where B = π ρ g α 2 ( H 2 ) , where ρ is the water density, and g is the acceleration due to gravity. Similar, the vertical forces on the permeable cylindrical body are equal to the sum of the forces on the upper and lower surfaces, i.e., 0 r a ,   z = h 1 and z = h , respectively. Thus,
F 3 = 2 π ρ ω 2 d H 2 [ 0 a 1 d Ψ D 0 I I I ( r , h ) r d r 0 a 1 d Ψ D 0 I I ( r , h 1 ) r d r ] .
Substituting Equations (20) and (27) into Equation (31), it is derived that
F 3 B = 2 k d tan h ( k d ) [ 1 2 F D 0 , 0 I I I + 2 n = 1 ( 1 ) n F D 0 , i I I I 1 n π a h I 1 ( n π h a ) I 0 ( n π h a ) N 0 , I I 1 2 F D 0 , 0 I I 1 k I I a J 1 ( k I I a ) J 0 ( k I I a ) i = 1 N i , I I 1 2 F D 0 , i I I 1 a i I I a I 1 ( a i I I a ) I 0 ( a i I I a ) ] .
The moment on the permeable cylindrical body about a horizontal axis at an arbitrary distance z = e from the seabed can be decomposed into M S and M b terms arising from the pressure distribution on the body’s wetted surfaces (i.e., permeable and impermeable). It holds that
M S = 2 π i ρ ω 2 d H 2 a [ h d 1 d Ψ D 1 I ( a , z ) ( z e ) d z h 1 d 1 d Ψ D 1 I I ( a , z ) ( z e ) d z ] ,
M b = 2 π ρ ω 2 d H 2 [ 0 a 1 d Ψ D 1 I I I ( r , h ) r 2 d r 0 a 1 d Ψ D 1 I I ( r , h 1 ) r 2 d r ] .
Substituting Equations (12), (20), and (27) into Equations (33) and (34), it is derived that
M s B a = 2 i k d tan h ( k d ) [ N 0 , I 1 2 F D 1 , 0 I L 0 I + i = 1 N i , I 1 2 F D 1 , i I L i I N 0 , I I 1 2 F D 1 , 0 I I L 0 I I i = 1 N i , I I 1 2 F D 1 , i I I L i I I ] ,
M b B a = 2 k d t a n h ( k d ) [ 1 4 F D 1 , 0 I I I + 2 n = 1 ( 1 ) n F D 1 , i I I I 1 n π a h I 2 ( n π h a ) I 1 ( n π h a ) N 0 , I I 1 2 F D 1 , 0 I I 1 k I I a J 2 ( k I I a ) J 1 ( k I I a ) i = 1 N i , I I 1 2 F D 1 , i I I 1 a i I I a I 2 ( a i I I a ) I 1 ( a i I I a ) ] .
The terms L 0 I , L i I , L 0 I I , L i I I are presented in Appendix A.
Similarly, the corresponding hydrodynamic reaction forces and pitching moment, F i j , on the permeable cylindrical body in the i -th direction due to its sinusoidal motion with frequency ω and unit amplitude in the j -th direction are equal to
F i j = ρ ω 2 S Ψ j 0 y ( r , z ) cos ( m θ ) n i d S ,   y = I , I I , I I I ; i , j = 1 , 3 , 5 ,
where S stands for the wetted surface, while n i was defined in Equation (2).
In addition, Equation (37) can be written as
F i j = ω 2 a i j + i ω b i j .
In Equation (38), a i j and b i j denote the hydrodynamic added mass and damping coefficients (both real and dependent of ω) in the i -th direction due to the body’s unit sinusoidal motion in the j -th direction.
Substituting Equations (12), (20), and (27) into Equation (37), the following relations for the nondimensional hydrodynamic coefficients can be obtained:
a 11 ρ a 3 + i b 11 ρ ω a 3 = π d a [ N 0 , I 1 2 F 11 , 0 I sinh ( k d ) sinh ( k h ) k a + i = 1 N i , I 1 2 F 11 , i I sin ( a i I d ) sin ( a i I h ) a i I a N 0 , I I 1 2 F 11 , 0 I I sinh ( k I I ( d h 1 ) ) k I I a i = 1 N i , I I 1 2 F 11 , i I I sin ( a i I I ( d h 1 ) ) a i I I a ] ,
a 33 ρ a 3 + i b 33 ρ ω a 3 = 2 π d a [ ( h 4 d a 2 16 h d ) + 1 2 F 30 , 0 I I I + 2 n = 1 ( 1 ) n F 30 , i I I I 1 n π a h I 1 ( n π h a ) I 0 ( n π h a ) 1 2 ( h 1 d 1 + g ω 2 d ) N 0 , I I 1 2 F 30 , 0 I I 1 k I I a J 1 ( k I I a ) J 0 ( k I I a ) i = 1 N i , I I 1 2 F 30 , i I I 1 a i I I a I 1 ( a i I I a ) I 0 ( a i I I a ) ] ,
a 55 ρ a 5 + i b 55 ρ ω a 5 = π d 2 a 2 ( M s + M b ) .
The terms M s and M b are defined by
M s = [ N 0 , I 1 2 F 51 , 0 I L 0 I + i = 1 N i , I 1 2 F 51 , i I L i I N 0 , I I 1 2 F 51 , 0 I I L 0 I I i = 1 N i , I I 1 2 F 51 , i I I L i I I ] ,
M b = [ ( 1 8 a h d 2 1 48 a 3 d 2 h ) + 1 4 F 51 , 0 I I I + 2 n = 1 ( 1 ) n F 51 , i I I I 1 n π a h I 2 ( n π h a ) I 1 ( n π h a ) ( a g 4 d 2 ω 2 a ( d h 1 ) 4 d 2 ) N 0 , I I 1 2 F 51 , 0 I I 1 k I I a J 2 ( k I I a ) J 1 ( k I I a ) i = 1 N i , I I 1 2 F 51 , i I I 1 a i I I a I 2 ( a i I I a ) I 1 ( a i I I a ) ] .
The terms L 0 I , L i I , L 0 I I , L i I I are presented in Appendix A.

3. Motion Equations

3.1. Mooring Line Characteristics

The mooring system constitutes a fundamental part for the motion response analysis of the permeable cylindrical body. In the present study, a four-point catenary mooring line system is considered. The lines are assumed as steel wires of unstretched length L , diameter d m , and elasticity modulus E . The geometry of a typical mooring line in 2D is shown in Figure 2a. A global mooring cable coordinate system ( x , y , z ) is defined, located at the intersection of the body’s vertical axis of symmetry, with the undisturbed free surface, r , being the position vector of the i-th line’s fairlead location with respect to the origin of the coordinate system and a i ,   i = 1 , 2 , 3 , 4 being the i-th line’s orientation angle in the horizontal plane (see Figure 2b). In addition, four local mooring coordinate systems ( x m i , y m i , z m i ), i = 1 , 2 , 3 , 4 are located at the connection points between the mooring lines and the body. It should be noted that the angle a i is formed by the positive axis of x and x m i during the counterclockwise rotation of x in order to coincide with the positive axis of x m i .
The mooring forces,   F l m , on the permeable cylindrical body in the l -th direction, expressed in the global mooring cable coordinate system ( x , y , z ), can be written in the frequency domain as
F l m = ( C l j m + i ω B l j m ) x j 0 .
In Equation (44), C l j m denotes the total restoring mooring stiffness to be applied on the permeable structure, whereas the term B l j m represents the corresponding mooring line damping coefficients. Both C l j m and B l j m terms are dependent on the exciting wave frequency. In order to evaluate these coefficients, the well-known quasi-static approximation, which is based on the static analysis of each line and allows prediction of the mooring stiffness imposed on the floating structure at zero excitation frequency, has to be properly extended by accounting of the line’s dynamic behavior.
Specifically, as far as the quasi-static approximation is concerned, the restoring mooring stiffness matrix (i.e., 6 × 6 square matrix), [ C ¯ l j m ] , l , j = 1 , , 6 , with respect to the global coordinate system, can be expressed through the restoring mooring stiffness matrix, c l j m , i , l , j = 1 , , 6 ,   i = 1 , 2 , 3 , 4 , of each line, defined in the local mooring coordinate system. Hence, for the coefficients C ¯ l j m , l , j = 1 , 2 , 3 , it holds that
[ C ¯ l j m ] = i = 1 4 [ λ i ] T [ c l j m , i ] [ λ i ] .
In Equation (45), the term [ C ¯ l j m ] is a 3 × 3 square matrix whose elements are C ¯ l j m , whereas [ λ i ] ,   [ c l j m , i ] are 3 × 3 square matrices defined by
[ λ i ] = [ cos a i sin a i 0 sin a i cos a i 0 0 0 1 ] ,   i = 1 , 2 , 3 , 4 ,
[ c l j m , i ] = [ c 11 m , i 0 c 13 m , i 0 c 22 m , i 0 c 31 m , i 0 c 33 m , i ] ,   i = 1 , 2 , 3 , 4 ; l , j = 1 , 2 , 3 .
The matrix [ λ i ] T is the transpose of [ λ i ] .
The remaining terms C ¯ l j m of the restoring mooring stiffness matrix are presented in Appendix B.
The mooring restoring coefficients, c l j m , i , of each line can be evaluated using well-known quasi-static equations of a single mooring line, i.e.,
c 11 m , i = δ P x δ x ¯ B ,     c 22 m , i = P x x B ,   c 33 m , i = δ P z δ z B ,   c 31 m , i = δ P z δ x ¯ B ,
where P x and P z denote the horizontal and vertical component of the tension force at the top of the line, x B is the projection of the suspended mooring line length in the horizontal direction, x ¯ B is the horizontal distance between the anchor and the fairlead, and z B is the vertical projection of the suspended mooring line length (see Figure 2a). Additional information on the terms of Equation (48) is presented in Appendix B.
The contribution of the mooring lines to the total damping of a moored permeable structure is a very important element for the evaluation of the body’s responses. Due to the line motions in the fluid domain, the phenomenon of energy dissipation appears, which offers to the moored body an additional amount of damping (i.e., mooring damping) originating from the drag and viscous forces on the mooring lines. The dynamic tension, F l j d , i , at the top of the line i , i = 1 , 2 , 3 , 4 , for sinusoidal motions of the upper end with amplitude A j can be written as follows [46]:
F l j d , i A j = C m , l j i + i B m , l j i , i = 1 , , 4 ; l , j = 1 , , 6 .
where C m , l j i and B m , l j i , both frequency- and excitation-amplitude-dependent, stand for the real and imaginary parts of F l j d , i / A j expressed in the local mooring coordinate system ( x m i , y m i , z m i ), i = 1 , 2 , 3 , 4 . In order to derive the total frequency-dependent mooring system restoring stiffnesses, C m , l j T , along with the corresponding total mooring line damping components, B m , l j T , by accounting for the mooring line’s dynamics, with respect to the global coordinate system, is subjected to the transformation expressed through Equation (45). Specifically, the coefficients C m , l j T and B m , l j T can be expressed as the sum of C m , l j i ,   B m , l j i ,   i = 1 , , 4 when the latter are formulated in the global coordinate system.
Summarizing, the total mooring restoring coefficients, C l j m , are equal to C ¯ l j m *see Equation (45) and Appendix B (when only quasi-static considerations are taken into account), whereas they are equal to C m , l j T , when dynamics of mooring lines are included. In this last case, the additional term B l j m in Equation (44) represents the mooring line damping coefficients which are equal to B m , l j T / ω , when expressed in the global coordinate system.

3.2. Equations of Motion

In the presented linear hydrodynamic analysis, the permeable cylindrical body is assumed to undergo small motions in its six degrees of freedom. Thus, the translational/rotational motions of the examined body are calculated by the following system of equations:
j = 1 6 [ ω 2 ( m l j + a l j ) i ω ( b l j + B l j m ) + ( c l j + C l j m ) ] x j 0 = F l .
In Equation (50), the term m l j denotes the generalized masses of the floater, a l j and b l j are the hydrodynamic added mass and damping coefficients, respectively (see Equation (38)), c l j represents the hydrostatic coefficients, and F l denotes the exciting forces and moments (see Equations (29), (31), (33), and (34)). C l j m represents the mooring restoring coefficients due to the mooring lines, and B l j m  represents the mooring line damping coefficients (see Equation (44)).
Under the assumption of symmetrical mass distribution and mooring arrangement, the examined body performs three degrees of freedom motions under the action of a regular wave train, i.e., two translations (i.e., surge and heave) and one rotation (i.e., pitch) in the wave propagation plane. Hence, the motions of the cylindrical body can be expressed in terms of the response amplitude operator (RAO).
RAO l = x j 0 H / 2 ,   l = 1 , 3 ;   RAO 5 = x j 0 k H / 2 ,
where k  denotes the wave number.
Equation (50) is solved through an iterative procedure in order for the amplitude-dependent mooring restoring coefficients and the mooring line damping to be determined. Regarding the quasi-static model, the mooring system is initially considered to undergo only pretention loads without any external excitation forces. Continually, the mooring system is displaced from its initial equilibrium position under the action of environmental-generated forces. Hence, the mooring characteristics of the system (i.e., tension forces, suspended mooring line length, horizontal distance between the anchor and the fairlead, and vertical projection of the suspended mooring line length) in its new displaced position are calculated. In addition, the mooring restoring coefficients, c l j m , i , of each line and, consequently, the restoring mooring stiffness coefficients C ¯ l j m , l , j = 1 , , 6 , in the global mooring coordinate system are defined (the followed procedure is presented in Appendix B). Regarding the dynamic modeling, it enables the evaluation of the dynamic tension, F l j d , i , at the top of each line i = 1 , 2 , 3 , 4 , and subsequently the mooring line damping, B l j m , and the mooring restoring coefficients C l j m . A detailed description of the evaluation of dynamic tensions was presented in [47,48], whereas, in [49], the motion responses of a CALM buoy under wave-current interactions were calculated. As far as the hydrodynamic calculations are concerned, the presented formulation, as described in Section 2, was applied and the hydrodynamic loads were fed into the motion Equation (50).
The coupling procedure of the mooring models with the hydrodynamic formulation was presented in [50] concerning a floating breakwater. Here, the procedure is extended for a permeable cylindrical body. Initially, the RAO l of the body is evaluated as if it were floating without mooring constraints, i.e., for zero C l j m and B l j m terms, and fed into the mooring dynamic model. Thus, the dynamic tensions are calculated for the specific values of the body’s motions. Subsequently, the corresponding values of C l j m and B l j m are applied to the hydrodynamic formulation, and new values of RAO l (i.e., denoted by RAO l 2 ) are determined. The iterative procedure continues until
RAO l N RAO l N 1 < ε ,
where N stands for the iteration-cycle number, while the value of ε  depends on the accuracy of the applied solution.

4. Results and Discussion

4.1. Methodology Validation

The numerical methods described above allow for the analysis of a great number of geometrical configurations of a permeable floating cylindrical body. For this purpose, the HAMVAB code [42,51] in FORTRAN was applied in the presented results.
Considering the available literature, the theoretical results from the presented methodology are compared with the results from [34]. In the latter study, comprehensive comparisons between theoretical and experimental results were performed. Specifically, a 1:2 scaled down model of a permeable cylinder, made of steel, was tested in various wave slope conditions and wave periods, as well as in an impact test case restrained to the wave and one in forced heave oscillations.
The examined permeable cylindrical body of radius a = 0.15 m and draught d h 1 = 0.3 m is subjected to an incident wave train at a water depth d = 5 m (see Figure 1). Here, the dimensionless porous effect parameter is equal to G = 1.432 , corresponding to an opening ratio of τ = 0.14 and wave steepness ε = 0.04633 (see Equation (3)). The validations of the results are made in terms of the dimensionless quantities of the surge, heave, and pitch wave exciting forces on the permeable body, i.e., F l = | F l | ρ g α 2 H / 2 ,   l = 1 , 3 and F 5 = | F 5 | ρ g α 3 H / 2 (see Equations (29), (31), (33), and (34)). In addition, the dimensionless hydrodynamic coefficients (added masses a 11 , a 33 and hydrodynamic damping coefficients b 11 , b 33 ), i.e., A 11 , 33 = a 11 , 33 ρ α 3 ,   B 11 , 33 = b 11 , 33 ω ρ α 3  (see Equations (39) and (40)), are also compared.
Figure 3 depicts the exciting force components in surge, heave, and pitch for the selected value of G . Moreover, the added mass and the damping coefficients for the specific porous effect parameter are presented. It can be noted an excellent correlation between the results of the present methodology and the results from [32]. Consequently, the presented theoretical formulation can effectively simulate the effect of the permeable sidewall of a cylindrical body.

4.2. Numerical Results

In the sequel, a permeable cylindrical body of radius a = 35 m is considered moored with a symmetrically arranged four-catenary-mooring-line system at a water depth d = 100 m. The draught of the permeable sidewall is equal to d h 1 = 17.5 m, whereas the thickness, h 1 h , of the impermeable bottom is assumed infinite (see Figure 1). To investigate the effect of porosity, five different sidewalls surfaces with different opening ratios are considered, i.e., τ = 0.05, 0.13, 0.2, 0.40, and 0.60. It is reminded that τ is defined as the ratio of the opened area to the total sidewall area. In addition, wave steepness ε = 0.05 is assumed. Hence, from Equation (3) the dimensionless porous effect parameters are equal to G = 0.18, 1.22, 2.62, 8.92, and 17.48. The mass of the cylindrical body is assumed constant regardless the porosity of the sidewall surface, i.e., M = 19.72 t, the center of gravity (CoG) is located at the body’s vertical axis, at 17.497 m below the water free surface, and the mass moment of inertia relative to water free surface equals to I z = 12,080 t· m 2 .
The examined permeable cylindrical body is anchored to the seabed through a symmetric mooring system, composed of four identical mooring lines (see Figure 4). The line’s mass is equal to 48.76 kg/m, whereas the submerged weight per meter of each line is w = 415.25 N/m. The lines are assumed steel wires of diameter d m = 100 mm and of elasticity modulus E = 164.65 GPa. The unstretched length of each line is L = 700 m, and the minimum breaking tension is equal to 9147 kN. Furthermore, the formed angle by each line on the XY plane, with respect to the X-axis, is π/2.
The mooring system is initially considered to undergo only pretention loads. Continually, it is displaced from its initial equilibrium position under the action of horizontal environment forces F x ( 2 ) =  300 kN. Table 1 summarizes the mooring properties of each line, whereas, in Table 2, the locations of the attachment points of the lines on the cylindrical body in relation to the global coordinate system are presented.
In Figure 5, the exciting forces on the permeable cylindrical body are presented for various values of porous effect parameter, i.e., G = 0.18, 1.22, 2.62, 8.92, and 17.48. Here, G = 0 and G 1 are also considered. It should be reminded that, for G = 0 , the sidewall is assumed impermeable, whereas, for G 1 , the sidewall is considered fully permeable. The exciting forces are normalized by the terms ( π ρ g α 2 H / 2 ) and ( π ρ g α 3 H / 2 ), i.e., F l = | F l | π ρ g α 2 H / 2 ,   l = 1 , 3 , and F 5 = | F 5 | π ρ g α 3 H / 2 . It can be seen from the figure that, as the porous effect parameter increases, the horizontal exciting forces decrease (see Figure 5a). Specifically, for a zero porous parameter (i.e., impermeable sidewall case) the horizontal exciting forces attain generally higher values than those for G 0 , since, in the latter case, the wave energy is absorbed by the permeable surface. It should also be noted a peculiar behavior of F 1 at ω 0.7 ; 1.2 ; 1.55 , etc. This effect is notable for G = 0.18, 1.22, 2.62, 8.92, and 17.48, whereas it is not depicted for G = 0 and G 1 . In the vicinity of the corresponding wave frequencies, the F 1 values attain a sharp decrease, which is more profound for lower values of G . The values of k a at these wave frequencies are equal to 1.84, 5.33, and 8.53. These k a values are in the neighborhood of the wave numbers which zero the derivative of the Bessel function of first kind, i.e., J 1 ( k a ) . Hence, it can be concluded that sloshing phenomena do occur in permeable cylindrical bodies.
As far as the heave exciting forces presented in Figure 5b are concerned, it can be seen that, for G = 0 , the variation pattern of F 3 differs from that for G > 0 . Specifically, for G = 0 , F 3 tends to unity for ω tending to zero, whereas, for G > 0 , the vertical exciting force zeros when ω   0. Furthermore, the wave frequencies in which F 3 attains local maxima are shifted to higher values as G increases. Nevertheless, as ω increases (ω > 1), the effect of the porous effect parameter is no longer significant. Concerning the overturning moments (see Figure 5c), it can be seen that, at small wave frequencies (i.e., ω < 0.4 ), F 5 increases reversely to G . On the other hand, for ω > 0.4 , the moment F 5 increases as G increases, tending to the values of the submerged disc case (i.e., for G 1 ). Moreover, for high values of ω, the porous effect parameter does not seem to affect F 5 .
Figure 6 and Figure 7 depict the hydrodynamic coefficients of the examined permeable cylindrical body. The coefficients are normalized as A 11 , 33 = a 11 , 33 ρ α 3 ,     B 11 , 33 = b 11 , 33 ω ρ α 3 , A 15 = a 15 ρ α 4 ,   B 15 = b 15 ω ρ α 4 , and A 55 = a 55 ρ α 5 ,   B 55 = b 55 ω ρ α 5 . Here, the cases of a truncated cylinder ( G = 0 case) and of a submerged cylindrical plate of infinite thickness ( G 1 case) are also considered. Regarding the added mass in surge direction, which is depicted in Figure 6a, it can be seen that A 11 decreases as the porous effect parameter increases, tending to zero values for the submerged disc case. In addition, a sharp oscillation pattern can be seen at ω 0.7 ; 1.2 ; 1.55 , etc., which is more profound for lower values of G corresponding to sloshing phenomena occurring in the case of partially porous bodies (see also discussion of Figure 5a). In Figure 6b, the added mass in heave A 33 is depicted. It can be seen that, for ω < 0.7, the added mass behaves proportionally with G , since it increases as G also increases tending to the values of the submerged plate case ( G 1 ). On the other hand, the effect of G on the values of A 33 can be considered negligible for ω > 0.7. The pitch A 55 and surge pitch A 15 added masses, depicted in Figure 6c,d, respectively, follow a similar rational to A 11 . Specifically, the tense oscillatory behavior at ω 0.7 ; 1.2 ; 1.55 , etc. is notable, which smoothens as G increases. In addition, the values of A 55 ,   A 15 tend to the corresponding values of the submerged plate case as G increases. It should be noted that, due to symmetry,   A 15 = A 51 , which also holds true for the permeable cylindrical bodies.
The effect of the porous parameter on the hydrodynamic damping coefficients of the examined permeable cylindrical body is shown in Figure 7, where the variations of B 11 ,   B 33 ,   B 55 ,   B 15 are presented versus the wave frequency. It can be seen from Figure 7a that B 11 behaves reverse proportionally with G for ω > 0.5. On the other hand, for small values of ω (i.e., ω < 0.5), the presented local maxima of B 11 are shifted to lower values of ω as G increases. Furthermore, the tense oscillatory behavior of B 11 for small values of G is also notable here, as in A 11 (see Figure 6a). The variation of B 33 is depicted in Figure 7b. The same conclusions can be drawn as in the A 33 variation, concerning the decrease in B 33 as G increases, for small values of wave frequencies (ω < 0.25), as well as the negligible effect of G on B 33 for high values of ω (i.e., ω > 1). In Figure 7c,d, B 55 and B 15 are presented, respectively. It can be seen that the latter damping coefficients attain a similar tendency as A 55 and A 15 (see Figure 6c,d), regarding the decrease in the damping coefficients as the porous effect parameter increases, tending to the results of the submerged plate case. It also holds true that B 15 = B 51 .
The effect of the porous parameter on the motions of the permeable cylindrical body is shown in Figure 8, where RAO l ,   l = 1 , 3 , 5 , (see Equation (51)) is presented for various values of G . Here, a mooring quasi-static formulation is considered. The considered mooring characteristics are presented in Table 1. The mass of the body is equal to 19.72 t, regardless of the value of G , since the porous sidewall is assumed to be infinitesimally thin and have negligible mass. In the presented figures, the impermeable truncated cylinder and the submerged plate cases are not considered since they attain completely different mass and hydrostatic coefficients. Consequently, comparison of their motions with the corresponding ones of the examined permeable body seems meaningless.
Starting with the surge motions (Figure 8a), it can be seen that RAO 1 variations behave reverse proportionally with G for ω tending to zero (i.e., ω < 0.3). On the other hand, for ω > 0.3, the porous effect parameter seems to have a small effect on the body’s surge motions. Furthermore, the effect of the sloshing phenomena inside the porous sidewall on RAO 1 is notable, minimizing the surge motion regardless the value of G at ω 0.7 ; 1.2 ; 1.55 , etc. Concerning the maximization of RAO 1 , this can be traced back to the mooring restoring stiffness, which imposes a resonance location in the surge motion at ω     0.1 . The heave motions of the permeable cylindrical body are depicted in Figure 8b. It can be seen that RAO 3 starts its variation from unity for ω tending to zero. Contrarily, as G increases, the body’s heave displacement decreases, whereas, for ω > 1, the porous parameter attains a negligible effect on RAO 3 . In Figure 8c, the body’s pitch motion is presented. In this figure, the mooring resonance at ω     0.1 regardless of the values of G should be noted, as well as the negligible effect of the sloshing phenomena on RAO 5 .
Next, the static analysis is extended by accounting for the lines’ dynamic behavior. Hence, the dynamic mooring damping and dynamic mooring restoring coefficients were evaluated and included in Equation (50). The followed procedure was described in Section 3.2. In Figure 9a,b the horizontal, Fd x x = F 11 d , 1 , and vertical, Fd z x = F 31 d , 1 , components of the dynamic tensions (see Equation (49)) at the top of the mooring line for horizontal sinusoidal motions of its upper end with amplitude A 1 are presented. Similarly, in Figure 9c the corresponding vertical components of the dynamic tensions, Fd z z = F 33 d , 1 , at the top of the mooring line for vertical sinusoidal motions of amplitude A 3 of its upper end are depicted. Here, A j = 0.5 ,   1.0 ,   1.5 ,   2.0   m ,   j = 1 , 3 . The considered mooring characteristics are presented in Table 1. The results of Figure 9 demonstrate clearly that the motion amplitude of the mooring upper end (i.e., connection point of each mooring line with the permeable cylindrical body) affects the dynamic tensions on the latter location. Specifically, the dynamic tensions for the particular inertia, the geometric mooring line characteristics, and the examined wave frequency range seem to behave proportionally with A j ,   j = 1 , 3 .
The effect of the dynamic mooring damping and dynamic mooring restoring coefficients on the permeable body’s motions for several porous effect parameters is shown in Figure 10, where the variations of RAO l ,   l = 1 , 3 , 5 , (see Equation (51)) for G = 0.18, 1.22, 2.62, 8.92, and 17.48 as a function of ω are depicted. Here, the dynamic mooring restoring coefficients and the dynamic dumping coefficient terms are considered in the solution of Equation (50) for the determination of the permeable body’s motions. In the followed iterative procedure, the horizontal and vertical components of the dynamic tension at the line’s top end, depicted in Figure 9 (see also Equation (49)), are coupled with the body’s motion equations (see Equation (50)) with a convergence coefficient ε = 1   ×   10−10. In the case of surge response, RAO 1 , (Figure 10a) the dynamic mooring lines stiffness, as well as the mooring damping consideration, lead to a small decrease in the values compared to the nondynamic case. This decrease is more profound for ω < 1.25, whereas this effect decreases for ω > 1.25. Nevertheless, RAO 1 follows in general the pattern of the quasi-static case (see Figure 9a). Figure 10b depicts the variation of heave response. It can be seen that the dynamic characteristics of the mooring lines do not seem to affect the body’s vertical displacements, since RAO 3 from the quasi-static analysis attains similar results to the dynamic case. The variation of RAO 5 is shown in Figure 10c. A small decrease in pitch motions is also observed due to strong coupling between surge and pitch, which dictates a pattern of RAO 5 similar to that of RAO 1 . Regarding the effect of G on the pitch motions, it can be seen that the body’s rotations increase proportionally to the porous effect parameter. This is more profound for ω < 0.5. On the other hand, the effect of G is negligible for ω > 0.5.

5. Conclusions

In the present paper, a moored permeable cylindrical body under the action of regular wave trains was investigated. Frequency analysis formulation was applied for the investigation of the effect of the porous parameter on the body’s hydrodynamics. Furthermore, the quasi-static and dynamic behavior of wire mooring lines on the body’s motions was studied. Special attention is given to the evaluation of dynamic stiffness and damping of the mooring lines through an iterative procedure. The presented numerical schemes account for the diffraction exciting forces and moments, the body’s hydrodynamic characteristics, and its translational and rotational motions. It is shown that the porous effect parameter plays a key role in reducing the wave loads on the permeable body by dissipating the wave energy. In addition, a small but significant decrease in the body’s surge and pitch motions is observed through the insertion of dynamic mooring characteristics on the system’s motion equation. On the other hand, the effect of the mooring dynamic characteristics on the body’s heave motion can be considered negligible.
The present theoretical formulation will be further developed in order to study the hydrodynamics of a moored cylindrical fish cage, similar to the examined permeable body, under the consideration of the sidewall’s flexibility in regular waves.

Author Contributions

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

Funding

This research received financial support by Greek national funds through the Operational Program “Competitiveness, Entrepreneurship, and Innovation” of the National Strategic Reference Framework (NSRF)—Research Funding Program: “MATISSE: Study of the appropriateness and the adequacy of modern materials for offshore fish cage—numerical and experimental investigation in realistic loading conditions”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Values of integrals L 0 I , L i I , L 0 I I , L i I I  defined in Equation (35):
L 0 I = k ( d e ) sinh ( k d ) k ( h e ) sinh ( k h ) cosh ( k d ) + cosh ( k h ) ( k a ) 2 .
L i I = a i I ( d e ) sin ( a i I d ) a i I ( h e ) sin ( a i I h ) + cos ( a i I d ) cos ( a i I h ) ( a i I a ) 2 .
L 0 I I = 1 k I I ( d e ) sin h ( k I I ( d h 1 ) ) + cosh ( k I I ( d h 1 ) ) ( k I I a ) 2 .
L i I I = 1 + a i I I ( d e ) sin ( a i I I ( d h 1 ) ) + cos ( a i I I ( d h 1 ) ) ( a i I I a ) 2 .

Appendix B

The C ¯ l j m  terms of the restoring mooring stiffness matrix are determined by
C ¯ 11 m = i = 1 4 c 11 m , i ( cos a i ) 2 + c 22 m , i ( sin a i ) 2 ,   C ¯ 22 m = i = 1 4 c 11 m , i ( sin a i ) 2 + c 22 m , i ( cos a i ) 2 ,   C ¯ 33 m = c 33 m , i ,
C ¯ 12 m = i = 1 4 c 11 m , i ( cos a i ) ( sin a i ) c 22 m , i ( cos a i ) ( sin a i ) ,
C ¯ 13 m = i = 1 4 c 13 m , i cos a i ,   C ¯ 23 m = i = 1 4 c 13 m , i sin a i ,
C ¯ 21 m = C ¯ 12 m ,     C ¯ 31 m = C ¯ 13 m   ,     C ¯ 32 m = C ¯ 23 m .
The remaining terms of C ¯ l j m  can be calculated by Equations (A5)–(A8).
C ¯ 41 m = C ¯ 31 m y 1 C ¯ 21 m z 1 ,   C ¯ 51 m = C ¯ 11 m z 1 C ¯ 31 m x 1 ,   C ¯ 61 m = C ¯ 21 m x 1 C ¯ 11 m y 1 .
C ¯ 42 m = C ¯ 32 m y 1 C ¯ 22 m z 1 ,   C ¯ 52 m = C ¯ 21 m z 1 C ¯ 32 m x 1 ,   C ¯ 62 m = C ¯ 22 m x 1 C ¯ 21 m y 1 .
C ¯ 43 m = C ¯ 33 m y 1 C ¯ 32 m z 1 ,   C ¯ 53 m = C ¯ 31 m z 1 C ¯ 33 m x 1 ,   C ¯ 63 m = C ¯ 32 m x 1 C ¯ 31 m y 1 .
C ¯ 44 m = C ¯ 33 m y 1 2 2 C ¯ 32 m y 1 z 1 + C ¯ 22 m z 1 2 .
C ¯ 54 m = C ¯ 31 m y 1 z 1 C ¯ 21 m z 1 2 C ¯ 33 m y 1 x 1 + C ¯ 32 m x 1 z 1 .
C ¯ 64 m = C ¯ 32 m x 1 y 1 C ¯ 31 m y 1 2 C ¯ 22 m x 1 z 1 + C ¯ 21 m y 1 z 1 .
C ¯ 55 m = C ¯ 11 m z 1 2 2 C ¯ 31 m x 1 z 1 + C ¯ 33 m x 1 2 .
C ¯ 65 m = C ¯ 21 m x 1 y 1 C ¯ 32 m x 1 2 C ¯ 11 m y 1 z 1 + C ¯ 31 m x 1 y 1 .
C ¯ 66 m = C ¯ 22 m x 1 2 2 C ¯ 21 m x 1 y 1 + C ¯ 11 m y 1 2 .
In the above Equations (A9)–(A17), ( x 1 , y 1 , z 1 ) are the coordinates of the mooring line attaching points on the permeable structure with respect to the global coordinate system located at the body’s vertical axis and at the undisturbed free surface.
The projection of the suspended mooring line length on the horizontal direction, x B , presented in Equation (48) is equal to
x B = P x w { sin h 1 ( P z P x ) sin h 1 ( P z w l s P x ) } + P x l s E A .
Furthermore, the vertical projection of the suspended mooring line length can be written as
z B = P x w { 1 + ( P z P x ) 2 1 + ( P z w l s P x ) 2 } + P z l s E A ,
where l s is the suspended mooring length, P x and P z are the horizontal and vertical components of the tension force at the top of the line, respectively, E is the elasticity modulus, w is the mooring line weight per meter in water, and A  is the line’s cross-sectional area. That is, for steel wires, it holds that
A = π 4 d m 2 ,
where d m  is the diameter of the steel wire.
The horizontal distance between the anchor and the top mooring-line attaching point x ¯ B  is equals to
x ¯ B = L l s + x B ,
where L  stands for the total length of the mooring line.

References

  1. Sollitt, C.; Cross, R. Wave transmission through permeable breakwaters. In Proceedings of the 13th Coastal Engineering Conference, Vancouver, BC, Canada, 10–14 July 1972; pp. 1827–1846. [Google Scholar]
  2. Madsen, O. Wave transmission through porous structures. J. Waterw. Ports Coast. Ocean. Eng. Div. 1974, 100, 169–188. [Google Scholar] [CrossRef]
  3. Solitt, C.; Cross, R. Wave Reflection and Transmission at Permeable Breakwaters; Tech. Paper 76-8; US Army Corps of Engineers, Coastal Engineering Research Center: Fort Belvoir, VA, USA, 1976. [Google Scholar]
  4. Sulisz, W. Wave reflection and transmission at permeable breakwaters of arbitrary cross section. Coast. Eng. 1985, 9, 371–386. [Google Scholar] [CrossRef]
  5. Pengzhi, L.; Karunarathna, S. Numerical study of solitary wave interaction with porous breakwaters. J. Waterw. Port Coast. Ocean. Eng. 2007, 133, 352–363. [Google Scholar]
  6. Huang, Z.; Li, Y.; Liu, Y. Hydraulic performance and wave loadings of perforated/slotted coastal structures: A review. Ocean Eng. 2011, 38, 1031–1053. [Google Scholar] [CrossRef]
  7. Lan, Y.; Hsu, T.; Lai, J.; Chwang, C.; Ting, C. Bragg scattering of waves propagating over a series of poro-elastic submerged breakwaters. Wave Motion 2011, 48, 1–12. [Google Scholar] [CrossRef]
  8. Liu, Y.; Li, H.J. Wave reflection and transmission by porous breakwaters: A new analytical solution. Coast. Eng. 2013, 78, 46–52. [Google Scholar] [CrossRef]
  9. Pereira, E.; The, H.; Manoharan, L.; Lim, C. Design optimization of porous box-type breakwater subjected to regular waves. MATEC Web Conf. 2018, 203, 01018. [Google Scholar] [CrossRef] [Green Version]
  10. Dai, J.; Wang, C.M.; Utsunomiya, T.; Duan, W. Review of recent research and developments on floating breakwaters. Ocean Eng. 2018, 158, 132–151. [Google Scholar] [CrossRef]
  11. Han, M.M.; Wang, C.M. Hydrodynamics study on rectangular porous breakwater with horizontal internal water channels. J. Ocean Eng. Mar. Energy 2020, 6, 377–398. [Google Scholar] [CrossRef]
  12. Kawakami, T. The theory of designing and testing fishing nets in model. In Modern Fishing Gear of the World (2); Fishing News (Books) Ltd.: London, UK, 2006; pp. 471–481. [Google Scholar]
  13. Aarsnes, J.V.; Rudi, H.; Loland, G. Current forces on cage, net deflection. Engineering for offshore fish farming. In Proceedings of the a Conference Organized by the Institution of Civil Engineers, Glasgow, UK, 17–18 October 1990. [Google Scholar]
  14. Bessonneau, J.S.; Marichal, D. Study of the dynamics of submerged supple nets (applications to trawls). Ocean. Eng. 1998, 25, 563–583. [Google Scholar] [CrossRef]
  15. Cho, I.H.; Kim, M.H. Interactions of horizontal porous flexible membrane with waves. J. Waterw. Port Coast. Ocean Eng. 2000, 126, 245–253. [Google Scholar] [CrossRef]
  16. Chan, A.; Lee, S.W.C. Wave characteristics past a flexible fishnet. Ocean Eng. 2001, 28, 1517–1529. [Google Scholar] [CrossRef]
  17. Wan, R.; Hu, F.; Tokai, T. A static analysis of the tension and configuration of submerged plane nets. Fish. Sci. 2002, 68, 815–823. [Google Scholar] [CrossRef]
  18. Tagaki, T.; Suzuki, K.; Hiraishi, T. Development of the numerical simulation method of dynamic fishing net shape. Nippon Suisan Gakkaishi 2002, 68, 320–326. [Google Scholar] [CrossRef] [Green Version]
  19. Fredheim, A.; Faltinsen, O.M. Hydroelastic analysis of a fishing net in steady inflow conditions. In Proceedings of the 3rd International Conference on Hydroelasticity in Marine Technology 2003, Oxford, UK, 15–17 September 2003. [Google Scholar]
  20. Bao, W.; Kinoshita, T.; Zhao, F. Wave forces acting on a semi-submerged porous circular cylinder. Proc. Inst. Mech. Eng. Part M J. Eng. Marit. Environ. 2009, 223, 349–360. [Google Scholar] [CrossRef]
  21. Faltinsen, O.M. Hydrodynamic aspects of a floating fish farm with circular collar. In Proceedings of the 26th International Workshop on Water Waves and Floating Bodies (IWWWFB 2011), Athens, Greece, 17–20 April 2011. [Google Scholar]
  22. Kristiansen, T.; Faltinsen, O.M. Modelling of current loads on aquaculture net cages. J. Fluids Struct. 2012, 34, 218–235. [Google Scholar] [CrossRef]
  23. Kristiansen, T.; Faltinsen, O.M. Experimental and numerical study of an aquaculture net cage with floater in waves and current. J. Fluids Struct. 2015, 54, 1–26. [Google Scholar] [CrossRef] [Green Version]
  24. Shen, Y.G.; Greco, M.; Faltinsen, O.M.; Nygaard, I. Numerical and experimental investigations on mooring loads of a marine fish farm in waves and current. J. Fluids Struct. 2018, 79, 115–136. [Google Scholar] [CrossRef]
  25. Ito, S.; Kinoshita, T.; Bao, W. Hydrodynamic behaviors of an elastic net structure. Ocean Eng. 2014, 92, 188–197. [Google Scholar] [CrossRef]
  26. Mandal, S.; Sahoo, T. Wave interaction with floating flexible circular cage system. In Proceedings of the 11th International Conference on Hydrodynamics (ICHD 2014), Singapore, 19–24 October 2014. [Google Scholar]
  27. Su, W.; Zhan, J.M.; Huang, H. Wave interactions with a porous and flexible cylindrical fish cage. Procedia Eng. 2015, 126, 254–259. [Google Scholar] [CrossRef] [Green Version]
  28. Ma, M.; Zhang, H.; Jeng, D.S.; Wang, C.M. A semi-analytical model for studying hydroelastic behaviour of a cylindrical net cage under wave action. J. Mar. Sci. Eng. 2021, 9, 1445. [Google Scholar] [CrossRef]
  29. Liu, Z.; Mohapatra, S.C.; Guedes Soares, C. Finite element analysis of the effect of currents on the dynamics of a moored flexible cylindrical net cage. J. Mar. Sci. Eng. 2021, 9, 159. [Google Scholar] [CrossRef]
  30. Ma, M.; Zhang, H.; Jeng, D.S.; Wang, C.M. Analytical solutions of hydroelastic interactions between waves and submerged open-net fish cage modeled as a porous cylindrical thin shell. Phys. Fluids 2022, 34, 017104. [Google Scholar] [CrossRef]
  31. Williams, A.N.; Li, W. Water wave interaction with an array of bottom-mounted surface-piercing porous cylinders. Ocean Eng. 2000, 27, 841–866. [Google Scholar] [CrossRef]
  32. Williams, A.N.; Li, W.; Wang, K.H. Water wave interaction with a floating porous cylinder. Ocean Eng. 2000, 27, 1–28. [Google Scholar] [CrossRef]
  33. Park, M.S.; Koo, W.; Choi, Y. Hydrodynamic interaction with an array of porous circular cylinders. Int. J. Nav. Archit. Ocean Eng. 2010, 2, 146–154. [Google Scholar] [CrossRef] [Green Version]
  34. Zhao, F.; Bao, W.; Kinoshita, T.; Itakura, H. Theoretical and experimental study on a porous cylinder floating in waves. ASME J. Offshore Mech. Arct. Eng. 2011, 133, 011301. [Google Scholar] [CrossRef]
  35. Zhao, F.; Kinoshita, T.; Bao, W.G.; Huang, L.Y.; Liang, Z.; Wan, R. Interaction between waves and an array of floating porous circular cylinders. China Ocean Eng. 2012, 26, 397–412. [Google Scholar] [CrossRef]
  36. Park, M.S.; Koo, W. Mathematical modeling of partial-porous circular cylinders with water waves. Hindawio Publishing Corporation. Math. Probl. Eng. 2015, 2015, 903748. [Google Scholar] [CrossRef]
  37. Dokken, J.; Grue, J.; Karstensen, P. Wave analysis of porous geometry with linear resistance law. J. Mar. Sci. Appl. 2017, 16, 480–489. [Google Scholar] [CrossRef] [Green Version]
  38. Dokken, J.; Grue, J.; Karstensen, P. Wave forces on porous geometries with linear and quadratic pressure-velocity relations. In Proceedings of the 32nd International Workshop on Water Waves and Floating Bodies, IWWWFB, Dalian, China, 23–26 April 2017. [Google Scholar]
  39. Cong, P.; Bai, W.; Teng, B. Analytical modeling of water wave interaction with bottom-mounted surface-piercing porous cylinder in front of a vertical wall. J. Fluids Struct. 2019, 88, 292–314. [Google Scholar] [CrossRef]
  40. Qiao, D.; Mackay, E.; Yan, J.; Feng, C.; Li, B.; Feichtner, A.; Ning, D.; Johanning, L. Numerical simulation with a macroscopic CFD method and experimental analysis of wave interaction with fixed porous cylinder structures. Mar. Struct. 2021, 80, 103096. [Google Scholar] [CrossRef]
  41. Kokkinowrachos, K.; Mavrakos, S.A.; Asorakos, S. Behavior of vertical bodies of revolution in waves. Ocean Eng. 1986, 13, 505–538. [Google Scholar] [CrossRef]
  42. Konispoliatis, D.; Chatjigeorgiou, I.; Mavrakos, S. Theoretical hydrodynamic analysis of a surface-piercing porous cylindrical body. Fluids 2021, 6, 320. [Google Scholar] [CrossRef]
  43. Sankar, A.; Bora, S.N. Hydrodynamic coefficients for a floating semi-porous compound cylinder in finite ocean depth. Mar. Syst. Ocean Technol. 2020, 15, 270–285. [Google Scholar] [CrossRef]
  44. Sankar, A.; Bora, S.N. Hydrodynamic forces and moments due to interaction of linear water waves with truncated partial-porous cylinders in finite depth. J. Fluids Struct. 2020, 94, 120898. [Google Scholar]
  45. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions, 9th ed.; Dover Publication: Washington, DC, USA, 1970. [Google Scholar]
  46. Mavrakos, S.A.; Chatjigeorgiou, I.K. Mooring-induced damping on floating structures. In Proceedings of the 1st International Conference on Marine Industry, Varna, Bulgaria, 2–7 June 1996; Volume II, pp. 365–378. [Google Scholar]
  47. Triantafyllou, M.S. Preliminary design of mooring systems. J. Ship Res. 1982, 26, 25–35. [Google Scholar] [CrossRef]
  48. Triantafyllou, M.S.; Bliek, A.; Shin, H. Static and Fatigue Analysis of Multi-Leg Mooring Systems; Technical report; MIT Press: Cambridge, MA, USA, 1986. [Google Scholar]
  49. Amaechi, C.V.; Wang, F.; Ye, J. Numerical studies on CALM buoy motion responses and the effect of buoy geometry cum skirt dimensions with its hydrodynamic waves—Current interactions. Ocean. Eng. 2022, 244, 110378. [Google Scholar] [CrossRef]
  50. Loukogeorgaki, E.; Angelides, D. Stiffness of mooring lines and performance of floating breakwater in three dimensions. Appl. Ocean Res. 2005, 27, 187–208. [Google Scholar] [CrossRef]
  51. Mavrakos, S.A. User’s manual for the computer code HAMVAB. In Laboratory for Floating Bodies and Mooring Systems; National Technical University of Athens: Athens, Greece, 2001. [Google Scholar]
Figure 1. A 3D representation of the examined permeable vertical cylindrical body.
Figure 1. A 3D representation of the examined permeable vertical cylindrical body.
Jmse 10 00403 g001
Figure 2. Representation of mooring cable and global coordinate system: (a) schematic representation of a 2D typical mooring line; (b) depiction of the mooring cable and global coordinate system orientation.
Figure 2. Representation of mooring cable and global coordinate system: (a) schematic representation of a 2D typical mooring line; (b) depiction of the mooring cable and global coordinate system orientation.
Jmse 10 00403 g002
Figure 3. Comparison of the present analytical results with Zhao et al. [32] in terms of (a) dimensionless exciting forces F 1 , F 3 ,   F 5 and (b) dimensionless hydrodynamic coefficients A 11 , B 11 ,   A 33 , B 33 . Here, G = 1.432 .
Figure 3. Comparison of the present analytical results with Zhao et al. [32] in terms of (a) dimensionless exciting forces F 1 , F 3 ,   F 5 and (b) dimensionless hydrodynamic coefficients A 11 , B 11 ,   A 33 , B 33 . Here, G = 1.432 .
Jmse 10 00403 g003
Figure 4. A 3D representation of the examined symmetric mooring system, composed by four identical mooring lines.
Figure 4. A 3D representation of the examined symmetric mooring system, composed by four identical mooring lines.
Jmse 10 00403 g004
Figure 5. Dimensionless exciting forces F 1 , F 3 ,   F 5 on the permeable cylindrical body against wave frequency ω for various values of porous effect parameter G = 0.18, 1.22, 2.62, 8.92, and 17.48: (a) F 1 ; (b) F 3 ; (c) F 5 . Here, G = 0 and G 1  cases are also presented.
Figure 5. Dimensionless exciting forces F 1 , F 3 ,   F 5 on the permeable cylindrical body against wave frequency ω for various values of porous effect parameter G = 0.18, 1.22, 2.62, 8.92, and 17.48: (a) F 1 ; (b) F 3 ; (c) F 5 . Here, G = 0 and G 1  cases are also presented.
Jmse 10 00403 g005aJmse 10 00403 g005b
Figure 6. Dimensionless hydrodynamic added mass A 11 , A 33 ,   A 55 ,   A 15 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter: (a) A 11 ; (b) A 33 ; (c) A 55 ; (d) A 15 . The results are also compared with the corresponding ones of a truncated cylinder ( G = 0   case) and of a submerged cylindrical plate ( G 1 case).
Figure 6. Dimensionless hydrodynamic added mass A 11 , A 33 ,   A 55 ,   A 15 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter: (a) A 11 ; (b) A 33 ; (c) A 55 ; (d) A 15 . The results are also compared with the corresponding ones of a truncated cylinder ( G = 0   case) and of a submerged cylindrical plate ( G 1 case).
Jmse 10 00403 g006
Figure 7. Dimensionless hydrodynamic damping coefficients B 11 , B 33 ,   B 55 ,   B 15 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter: (a) B 11 ; (b) B 33 ; (c) B 55 ; (d) B 15 . The results are also compared with the corresponding ones of a truncated cylinder ( G = 0 case) and of a submerged cylindrical plate ( G 1 case).
Figure 7. Dimensionless hydrodynamic damping coefficients B 11 , B 33 ,   B 55 ,   B 15 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter: (a) B 11 ; (b) B 33 ; (c) B 55 ; (d) B 15 . The results are also compared with the corresponding ones of a truncated cylinder ( G = 0 case) and of a submerged cylindrical plate ( G 1 case).
Jmse 10 00403 g007
Figure 8. Dimensionless motions RAO 1 , RAO 3 ,   RAO 5 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter G = 0.18, 1.22, 2.62, 8.92, and 17.48: (a) RAO 1 ; (b) RAO 3 ; (c) RAO 5 .
Figure 8. Dimensionless motions RAO 1 , RAO 3 ,   RAO 5 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter G = 0.18, 1.22, 2.62, 8.92, and 17.48: (a) RAO 1 ; (b) RAO 3 ; (c) RAO 5 .
Jmse 10 00403 g008
Figure 9. Horizontal and vertical dynamic tensions at the top of the mooring line against wave frequency ω for various values of horizontal motion amplitudes A j = 0.5 ,   1.0 ,   1.5 ,   and   2.0 m: (a) Fd x x ; (b) Fd z x ; (c) Fd z z .
Figure 9. Horizontal and vertical dynamic tensions at the top of the mooring line against wave frequency ω for various values of horizontal motion amplitudes A j = 0.5 ,   1.0 ,   1.5 ,   and   2.0 m: (a) Fd x x ; (b) Fd z x ; (c) Fd z z .
Jmse 10 00403 g009
Figure 10. Dimensionless motions RAO 1 , RAO 3 ,   RAO 5 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter G = 0.18, 1.22, 2.62, 8.92, and 17.48: (a) RAO 1 ; (b) RAO 3 ; (c) RAO 5 . Here, the dynamic mooring characteristics were taken into consideration.
Figure 10. Dimensionless motions RAO 1 , RAO 3 ,   RAO 5 of the permeable cylindrical body against wave frequency ω for various values of porous effect parameter G = 0.18, 1.22, 2.62, 8.92, and 17.48: (a) RAO 1 ; (b) RAO 3 ; (c) RAO 5 . Here, the dynamic mooring characteristics were taken into consideration.
Jmse 10 00403 g010
Table 1. Mooring characteristics.
Table 1. Mooring characteristics.
w Submerged weight per unit length415.25 N/m
d Line diameter0.1 m
L Mooring length700 m
E A Axial stiffness1.293 × 109 N
F x ( 2 ) Horizontal environmental force300 × 103 N
l B Bottom seated mooring length79.11 m
P x i Horizontal   force   at   the   cable s   top ,   i = 1 , , 4 780 × 103 N
P z i Vertical   force   at   the   cable s   top ,   i = 1 , , 4 258 × 103 N
x B Horizontal projection of the suspended mooring line length610.5 m
c 11 m Mooring stiffness in x-direction due to the unit translational motion of its attachment point in the x-direction135.9 kN/m
c 22 m Mooring stiffness in y-direction due to the unit translational motion of its attachment point in the y-direction1.27 kN/m
c 33 m Mooring stiffness in z-direction due to the unit translational motion of its attachment point in the z-direction4.84 kN/m
c 13 m Mooring stiffness in x-direction due to the unit translational motion of its attachment point in the z-direction21.87 kN/m
Table 2. Mooring attachment points.
Table 2. Mooring attachment points.
Line 1Line 2Line 3Line 4
a 0 π /2 π 3 π /2
x m i 350−350
z m i 0000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Konispoliatis, D.N.; Chatjigeorgiou, I.K.; Mavrakos, S.A. Hydrodynamics of a Moored Permeable Vertical Cylindrical Body. J. Mar. Sci. Eng. 2022, 10, 403. https://doi.org/10.3390/jmse10030403

AMA Style

Konispoliatis DN, Chatjigeorgiou IK, Mavrakos SA. Hydrodynamics of a Moored Permeable Vertical Cylindrical Body. Journal of Marine Science and Engineering. 2022; 10(3):403. https://doi.org/10.3390/jmse10030403

Chicago/Turabian Style

Konispoliatis, Dimitrios N., Ioannis K. Chatjigeorgiou, and Spyros A. Mavrakos. 2022. "Hydrodynamics of a Moored Permeable Vertical Cylindrical Body" Journal of Marine Science and Engineering 10, no. 3: 403. https://doi.org/10.3390/jmse10030403

APA Style

Konispoliatis, D. N., Chatjigeorgiou, I. K., & Mavrakos, S. A. (2022). Hydrodynamics of a Moored Permeable Vertical Cylindrical Body. Journal of Marine Science and Engineering, 10(3), 403. https://doi.org/10.3390/jmse10030403

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