Hydrodynamics of a Moored Permeable Vertical Cylindrical Body

: 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 ﬂoating 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 ﬂuid 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


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

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. (1) 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 + ϕ y 7 + ∑ j = 1,3,5 .
x j0 ϕ y j . (1) 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 ϕ y 7 , 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 ϕ y j , y = I, I I, I I I. Furthermore, the body's complex velocity amplitude in the j-th direction is .
x j0 . 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., ϕ y D = ϕ 0 + ϕ y 7 , y = I, I I, I I I. The velocity potential ϕ y j , 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, ϕ y j 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 ∂ϕ I I j ∂r = n j + ikG ϕ I I j − ϕ I j 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 re + iG im , where G re represents the linearized drag effect of the permeable sidewall, and G im 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 ε = kH/2, through Equation (3).
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 ϕ y j , 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]: ∂ϕ I j ∂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 The velocity potential of the undisturbed incident wave propagating along the positive x-axis can be expressed as follows [41]: 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 Furthermore, the radiation potentials ϕ y j , j = 1, 3, 5 can be expressed as follows [41]: The functions Ψ  (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 Ψ y Dm , Ψ y 11 , Ψ y 30 , Ψ y 51 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 Ψ y lm , 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, 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 I lm is equal to 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]: where 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 where the superscript denotes the infinite fluid region I. Hence, it is convenient to write a I 0 = −ik, resulting in the following [45]: Fluid region, I I: 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]: where: Furthermore, it holds that The terms a I I i are the roots of with the imaginary a I I i = −ik I I considered as first. Fluid region, I I I: 0 ≤ r ≤ a, 0 ≤ z ≤ h In Equation (27), the term f I I I lm is equal to The functions Ψ y lm , l = D, 1, 3, 5 have the advantage of being expressed by simple Fourier series representations, F I lm,i , F I I lm,i , F I I I lm,n , 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.

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 Substituting Equations (12) and (20) into Equation (29), we get 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, Substituting Equations (20) and (27) into Equation (31), it is derived that 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 Substituting Equations (12), (20), and (27) into Equations (33) and (34), it is derived that The terms L I 0 , L I i , L I I 0 , L I I i are presented in Appendix A. Similarly, the corresponding hydrodynamic reaction forces and pitching moment, F ij , 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 where S stands for the wetted surface, while n i was defined in Equation (2). In addition, Equation (37) can be written as In Equation (38), a ij and b ij 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 55 The terms M s and M b are defined by The terms L I 0 , L I i , L I I 0 , L I I i are presented in Appendix A.

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 i m , y i m , z i m ), 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 i m during the counterclockwise rotation of x in order to coincide with the positive axis of x i m .
The terms 0 , , 0 , are presented in Appendix A.

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 , diameter , and elasticity modulus . The geometry of a typical mooring line in 2D is shown in Figure 2a. A global mooring cable coordinate system ( , , ) is defined, located at the intersection of the body's vertical axis of symmetry, with the undisturbed free surface, ⃗, being the position vector of the i-th line's fairlead location with respect to the origin of the coordinate system and , = 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 ( , , ), = 1,2,3,4 are located at the connection points between the mooring lines and the body. It should be noted that the angle is formed by the positive axis of and during the counterclockwise rotation of in order to coincide with the positive axis of .
(a) (b) The mooring forces, , on the permeable cylindrical body in the -th direction, expressed in the global mooring cable coordinate system ( , , ), can be written in the frequency domain as The mooring forces, F m l , 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 In Equation (44), C m lj denotes the total restoring mooring stiffness to be applied on the permeable structure, whereas the term B m lj represents the corresponding mooring line damping coefficients. Both C m lj and B m lj 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 m lj , l, j = 1, . . . , 6, with respect to the global coordinate system, can be expressed through the restoring mooring stiffness matrix, c m,i lj , l, j = 1, . . . , 6, i = 1, 2, 3, 4, of each line, defined in the local mooring coordinate system. Hence, for the coefficients C m lj , l, j = 1, 2, 3, it holds that In Equation (45), the term C m lj is a 3 × 3 square matrix whose elements are C m lj , whereas λ i , c m,i lj are 3 × 3 square matrices defined by The matrix λ i T is the transpose of λ i .
The remaining terms C m lj of the restoring mooring stiffness matrix are presented in Appendix B.
The mooring restoring coefficients, c m,i lj , of each line can be evaluated using well-known quasi-static equations of a single mooring line, i.e., 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 d,i lj , 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]: where C i m,lj and B i m,lj , both frequency-and excitation-amplitude-dependent, stand for the real and imaginary parts of F d,i lj /A j expressed in the local mooring coordinate sys- In order to derive the total frequency-dependent mooring system restoring stiffnesses, C T m,lj , along with the corresponding total mooring line damping components, B T m,lj , 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 T m,lj and B T m,lj can be expressed as the sum of C i m,lj , B i m,lj , i = 1, . . . , 4 when the latter are formulated in the global coordinate system. Summarizing, the total mooring restoring coefficients, C m lj , are equal to C m lj *see Equation (45) and Appendix B (when only quasi-static considerations are taken into account), whereas they are equal to C T m,lj , when dynamics of mooring lines are included. In this last case, the additional term B m lj in Equation (44) represents the mooring line damping coefficients which are equal to B T m,lj /ω, when expressed in the global coordinate system.

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: In Equation (50), the term m lj denotes the generalized masses of the floater, a lj and b lj are the hydrodynamic added mass and damping coefficients, respectively (see Equation (38)), c lj represents the hydrostatic coefficients, and F l denotes the exciting forces and moments (see Equations (29), (31), (33), and (34)). C m lj represents the mooring restoring coefficients due to the mooring lines, and B m lj 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).
where k denotes the wave number. Equation (50) is solved through an iterative procedure in order for the amplitudedependent 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 environmentalgenerated 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 m,i lj , of each line and, consequently, the restoring mooring stiffness coefficients C m lj , 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 d,i lj , at the top of each line i = 1, 2, 3, 4, and subsequently the mooring line damping, B m lj , and the mooring restoring coefficients C m lj . 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 m lj and B m lj 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 m lj and B m lj are applied to the hydrodynamic formulation, and new values of RAO l (i.e., denoted by RAO 2 l ) are determined. The iterative procedure continues until where N stands for the iteration-cycle number, while the value of ε depends on the accuracy of the applied solution.

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. Here, the dimensionless porous effect parameter is equal to = 1.432, corresponding to an opening ratio of τ = 0.14 and wave steepness ε = 0.04633 (see Equation (3) (34)). In addition, the dimensionless hydrodynamic coefficients (added masses , and hydrodynamic damping coefficients , ), i.e., A , = , , B , = , (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 . 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.

Numerical Results
In the sequel, a permeable cylindrical body of radius = 35 m is considered moored with a symmetrically arranged four-catenary-mooring-line system at a water depth = 100 m. The draught of the permeable sidewall is equal to − ℎ = 17.5 m, whereas the thickness, ℎ − ℎ, 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) [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 .

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 unstretched length of each line is = 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 ( ) = 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.     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 (2) x = 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 ka at these wave frequencies are equal to 1.84, 5.33, and 8.53. These ka values are in the neighborhood of the wave numbers which zero the derivative of the Bessel function of first kind, i.e., J 1 (ka). 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 . maxima are shifted to higher values as 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 increases reversely to . On the other hand, for > 0.4, the moment F increases as increases, tending to the values of the submerged disc case (i.e., for ≫ 1). Moreover, for high values of ω, the porous effect parameter does not seem to affect F . Here, the cases of a truncated cylinder ( = 0 case) and of a submerged cylindrical plate of infinite thickness ( ≫ 1 case) are also considered. Regarding the added mass in surge direction, which is depicted in Figure 6a, it can be seen that A 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 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 is depicted. It can be seen that, for ω < 0.7, the added mass behaves proportionally with , since it increases as also increases tending to the values of the submerged plate case ( ≫ 1). On the other hand, the effect of on the values of A can be considered negligible for ω > 0.7. The pitch A and surge pitch A added masses, depicted in Figure 6c,d, respectively, follow a similar rational to A . Specifically, the tense oscillatory behavior at ≈ 0.7; 1.2; 1.55, etc. is notable, which smoothens as increases. In addition, the values of A , A tend to the corresponding values of the submerged plate case as increases. It should be noted that, due to symmetry, A = A , which also holds true for the permeable cylindrical bodies. 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 , B , B , B are presented versus the wave frequency. It can be seen from Figure 7a that B behaves reverse proportionally with for ω > 0.5. On the other hand, for small values of ω (i.e., ω < 0.5), the presented local maxima of B are shifted to lower values of ω as increases. Furthermore, the tense oscillatory behavior of B for small values of is also notable here, as in A (see Figure 6a). The variation of B is depicted in Figure  7b. The same conclusions can be drawn as in the A variation, concerning the decrease in B as increases, for small values of wave frequencies (ω < 0.25), as well as the negligible effect of on B for high values of ω (i.e., ω > 1). In Figure 7c,d, B and B are presented, respectively. It can be seen that the latter damping coefficients attain a similar 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 . tendency as A and A (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 = B . The effect of the porous parameter on the motions of the permeable cylindrical body is shown in Figure 8, where RAO , = 1,3,5, (see Equation (51)) is presented for various values of . 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 , 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. 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 variations behave reverse proportionally with 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 is notable, minimizing the surge motion regardless the value of at ≈ 0.7; 1.2; 1.55, etc. Concerning the maximization of RAO , 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 starts its variation from unity for ω tending to zero. Contrarily, as increases, the body's heave displacement decreases, whereas, for ω > 1, the porous parameter attains a negligible effect on RAO . In Figure 8c, the body's pitch motion is presented. In this figure, the mooring resonance at ω≈ 0.1 regardless of the values of should be noted, as well as the negligible effect of the sloshing phenomena on RAO .
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 = , , and vertical, Fd x = , , components 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 d, 1 11 , and vertical, Fd z x = F d,1 31 , 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 d, 1 33 , 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.  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 = 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 , = 1,3,5, (see Equation (51)  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. 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 , (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 follows in general the pattern of the quasistatic 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 from the quasi-static analysis attains similar results to the dynamic case. The variation of RAO 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 similar to that of RAO . Regarding the effect of 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 is negligible for ω > 0.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. 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 Furthermore, the vertical projection of the suspended mooring line length can be written as 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 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 where L stands for the total length of the mooring line.