# Analytical Model of Wave Loads and Motion Responses for a Floating Breakwater System with Attached Dual Porous Side Walls.

1. IntroductionThe ocean-oriented activities, such as oil and gas extraction, aquacultural production, and recreations, have been identified to generate tremendous values to humans. Protection of the shoreline structures or offshore production facilities by degradation of the wave-structure interaction is therefore a vital issue in the coastal and offshore engineering applications. Various artificial structures have been designed and optimized for this purpose, of which the floating breakwaters or barriers are considered as valuable alternatives to conventional fixed gravity-type breakwaters [1]. The performance optimization of a breakwater, especially the type of floating breakwater, has received much attention with the advantage of energy dissipation [2]. Additionally, the porous structures may provide an added alternative to improve the performance of floating breakwaters.

Hydrodynamic analysis is widely applied to examine the wave-reduction performance of a floating breakwater. Mei and Black [3] investigated the wave forces acting on the rectangular floating obstacles by the variational-formulation-based numerical approach. Later, the wave radiation induced by small oscillations of a floating body was also studied [4]. A simplified analytical model to estimate the hydrodynamic performance of a long rectangular floating breakwater was developed under the assumption of a small gap condition [5]. Also, the diffraction and radiation problem of a rectangular buoy were analyzed under the conditions of heave oscillation [6], as well as in three dimensional scale, the sway, heave, and rolling motions [7]. Damping mechanisms of a floating breakwater were investigated numerically and experimentally considering the cases of a regular pontoon, a regular pontoon with wing plates attached, and a regular pontoon with wing plates penetrated [8]. Additionally, the hydrodynamic performances of rectangular boxes of fixed single and fixed double subject to the heave motion were estimated under the actions of either regular or irregular water waves [9]. Later, the wave absorbing effectiveness of different types of floating breakwaters, such as single box, double boxes, or board net, was also investigated [10].

Due to the advantage of energy dissipation, the water-wave absorbing phenomena have been observed when waves pass through a porous medium; that is, a porous structure can effectively attenuate both reflected and transmitted waves. Following the porous wavemaker theory of Chwang [11], considerable studies related to the interaction between porous structures and water waves were performed theoretically and experimentally, especially the porous barriers or the porous structures attached to the floating breakwaters. The wave energy dissipation of thin permeable wall was verified by conducting sets of experiments [12]. Analytical approaches have been carried out by Manam and Sivanesan [13] and Chwang and Dong [14], as well as Zhu and Chwang [15], to study the reflection and transmission problem of water waves propagating past vertical porous barriers. The hydroelastic issues between gravity waves and a submerged horizontal flexible porous plate were investigated by Behera and Sahoo [16]. The studies of perforated wall breakwaters and dual submerged horizontal porous plates were given in [17, 18], respectively. In order to examine the effectiveness of wave absorption by porous plates, the interactions between different types of porous plates and oblique monochromatic waves were investigated by Cho and Kim [19]. The concept of porous barrier was extended by Wang and Ren [20, 21] to analytically study the interaction between water waves and a flexible porous breakwater or, in a three-dimensional domain, a concentric porous cylinder system. The diffraction problems of either a circular porous plate submerged horizontally or a system with dual circular porous plates were investigated based on the linear potential flow theory [22, 23]. Recently, the scattering of water waves by a rectangular floating body that attached with two vertically porous walls was investigated analytically and experimentally by Qiao et al. [24]. In practices, the types of porous structure are various depending on the requirements of engineering problem; however, single or several vertical porous plates (including barrier and circular) and porous cylinder are the most common in engineering field. The vertical porous plates were often employed to absorb waves to protect coastal [25, 26]; furthermore, a set of porous plates were placed at the bottom of a wave tank to absorb waves [27]. The porous cylinder was employed to protect structure in Trafalgar offshore windfarm [28]; almost for the same purpose, a concentric two-cylinder system, involving an interior impermeable cylinder and external porous cylinder, was designed for Ekofisk gravity offshore structure in the North Sea [29]; then, the concept of single-layered porous cylinder was extended to double-layered perforated cylinder by Xiao et al. [30], and the experiments were conducted to investigate the hydrodynamic performances.

To achieve the combination of the advantages of a floating body and porous wall features, a new breakwater system that consists of a partially submerged pontoon-like body and two attached porous side walls is proposed in this study to provide a protected region sheltering from wave action. The results of wave reflection/transmission and wave induced hydrodynamic loads under the condition that the breakwater system is in a fixed position were given in Qiao et al. [24]. The investigation of wave interactions with the breakwater system is further extended to include the effect of the induced motions of the floating breakwater system, the so-called radiation problem. Under the assumptions of linear incident waves and small translational and rotational motions of surge, heave, and pitch, this study focuses on the investigation of hydrodynamic coefficients, such as the added mass and damping coefficient, and the dynamic responses of the floating breakwater system through developing an analytical model with solutions of diffracted and radiated waves and the equation of motion which is formulated and solved. The solutions of the velocity potentials in the defined regions, including the diffraction and radiation velocity potentials, are obtained by solving the formulated governing equations as well as the boundary and matching conditions. To verify the analytical solutions derived from this study, some limited numerical results published in literatures were adopted for comparisons. Additionally, the effects of the incident wave parameters and porous walls conditions on the hydrodynamic coefficients and RAOs (Response Amplitude Operators) are examined and discussed.

2. Theoretical Formulations and Solutions

The problem statement as illustrated in Figure 1 shows a train of right-going, small-amplitude, Airy waves with the frequency u>, wave height H, and wave length A to interact normally with a floating breakwater system that consists of a partially submerged rectangular structure and two attached fully submerged side porous walls. Subject to the wave induced hydrodynamic loads, the breakwater system responds with two translational motions and one rotational motion, i.e., surge, heave, and pitching motions. The Cartesian coordinates system is employed to define this two-dimensional boundary value problem. The x-axis is along the still water surface, where x = 0 is set at the center of the floating breakwater system, and the z-axis points vertically upwards with an origin located at the still water surface. The rectangular structure with a length of 2a and a submergence of s is initially positioned in a domain with undisturbed constant water depth h. Both of the thin porous walls have a length of d-s.

In this study, the fluid is assumed to be incompressible and inviscid and the motion irrational. Subsequently, the potential flow theory, providing that the velocity potential satisfies the Laplace Equation, with the defined boundary conditions could be applied to solve the waves and breakwater interaction problem. As shown in Figure 1, the entire fluid domain is divided into three subdomains, namely, regions I, II, and III, respectively. Region I is the subdomain upstream of the breakwater system (x [less than or equal to] a, -h [less than or equal to] z [less than or equal to] 0); region II is the subdomain between the two porous walls (-a < x < a, -d < z [less than or equal to] -s); and region III is the subdomain downstream of the breakwater system (a [less than or equal to] x, -h [less than or equal to] z [less than or equal to] 0). In the fluid domain, the Laplace equation of a complete velocity potential serves as the basic equation to describe the fluid motion, namely,

[[nabla].sup.2][phi] = 0; (1)

The fluid velocity potential f includes three components, incident wave potential [[phi].sub.I], diffracted wave potential fD, and the radiated wave potential [[PHI].sub.R], which covers the induced motions of the floating breakwater system. In this study, the dynamic motions of surge, heave, and pitching are considered. Thus, the velocity potential [phi] can be expressed as

[phi] = [[phi].sub.I] + [[phi].sub.D] + [3.summation over (j=1)] [[PHI].sup.(L).sub.R] (2)

where the subscript L = 1,2,3 donates the surge, heave, and pitching motions of the breakwater system, respectively. In this study, the incident wave potential is given as

[mathematical expression not reproducible] (3)

where [k.sub.0] is the wave number, [omega] is the wave frequency, and g is the gravitational acceleration. Here, [k.sub.0] and [omega] satisfy the usual dispersion relation

[[omega].sup.2] = [gk.sub.0] tanh ([k.sub.0]h) (4)

2.1. Theoretical Formulations and Solutions of the Incident and Diffracted Velocity Potentials. To solve the wave radiation problem in a linear wave system, the wave induced hydrodynamic forces and moment according to the wave diffraction approach by assuming the body is in a fixed position are required to be determined priorly. The wave diffraction solutions for the proposed breakwater system with two attached porous walls have been obtained in a previous study given in Qiao et al. [24] and Yip and Chwang [31]. Here, only limited formulations and solutions are summarized below. Let us define the wave diffraction potentials as [[phi].sub.Di]; i = 1,2,3 stands for the subdomains I, II, and III, respectively. The Laplace equations of [[phi].sub.Di] can be solved with the uses of the linearized kinematic and dynamic free surface boundary condition, bottom boundary condition, no flow condition at the body surface, and matching conditions at the interfaces of porous wall locations, including the porous flow conditions expressed as (Chwang [11])

[mathematical expression not reproducible] (5)

[mathematical expression not reproducible] (6)

where [b.sub.1] and [b.sub.2], having the dimension of length, are, respectively, material constants for the porous walls A and B, [p.sub.Di] (i= 1, 2, 3) are the dynamic pressures induced by the incident and diffracted velocity potentials, and g is the dynamic viscosity of the fluid. The wave diffraction solutions are derived as

[mathematical expression not reproducible] (7)

[mathematical expression not reproducible] (8)

[mathematical expression not reproducible] (9)

where [F.sub.m], [C.sub.m], [E.sub.m], and Tm (m = 0, 1, 2, 3, ...) are unknown complex coefficients to be determined. [I.sub.m], [[??].sub.m], and [r.sub.m] (m = 0, 1, 2, 3, ...) are expressed separately as

[I.sub.m] = {-I gH/2[omega], 0, 0, 0, ... ...} (10)

{[??].sub.m] = -i[k.sub.0], [k.sub.2], -[k.sub.2], -[k.sub.3], ... ..., -[k.sub.m] ... ...} (11)

[r.sub.m] = m[pi]/h - s (12)

Here, [k.sub.m] can be calculated from

[[omega].sup.2] = -[gk.sub.m] tanh ([k.sub.m]h) (13)

Substituting (7), (8), and (9) into the following derived matching equations

[mathematical expression not reproducible] (14)

[mathematical expression not reproducible] (15)

[mathematical expression not reproducible] (16)

[mathematical expression not reproducible] (17)

and utilizing the orthogonality properties of solution eigenfunctions, the unknown coefficients [C.sub.n], [D.sub.n], [R.sub.n], [T.sub.n] (n = 0, 1, 2, 3, ..., M) for the diffraction potentials can be determined by solving the 4(M+1) system of algebraic equations. Here, the terms with summation of infinite series are approximated by a large but finite number of M terms. In (16) and (17), [mathematical expression not reproducible] are given as

[mathematical expression not reproducible] (18)

[mathematical expression not reproducible] (19)

where [G.sub.01] = [rho][omega][b.sub.1]/([mu][k.sub.0]) and [G.sub.02] = [rho][omega][b.sub.2]/([mu][k.sub.0]) are dimensionless parameters defined according to Chwang [11] to represent the porous effect parameters of porous walls A and B, respectively.

2.2. Theoretical Formulations of the Radiated Velocity Potentials. For the radiation problem, the displacements of the floating breakwater system in motion are given as

[X.sub.j] = [X.sub.0j] [e.sup.-i[omega]t] j = 1, 2, 3 (20)

where [X.sub.oj] (j = 1,2, 3) are the amplitudes of surge, heave, and pitching motions, respectively. Let [[PHI].sup.(L).sub.Ri] = 1, 2, 3; i = 1, 2, 3) be the spatial radiated potentials describing the small-amplitude oscillations, where Lth indices denote, respectively, the induced surge, heave, and pitching motions while the ith indices represent the subdomains of I, II, and III, respectively. Then, the total radiated potentials in any of the flow domains can be expressed as

[mathematical expression not reproducible] (21)

As the velocity potentials [[phi].sup.(L).sub.R] satisfy the Laplace equation, we have

[[nabla].sup.2] [[phi].sup.(L).sub.Ri] = 0 (L = 1, 2, 3; I = 1, 2, 3) (22)

The boundary conditions, including the free surface boundary conditions, bottom boundary conditions, and the body surface conditions, for the surge, heave, and pitching oscillations, are listed as follows, respectively.

[mathematical expression not reproducible] (23)

[mathematical expression not reproducible] (24)

[mathematical expression not reproducible] (25)

on body surface

where n = ([n.sub.x], [n.sub.z]) is the unit normal vector and ([[xi].sub.x], [[xi].sub.z]) is the rotational center of the pitching oscillation. With the assumption that the rotational center is along the z axis, we can get [[xi].sub.x] = 0.

2.3. Analytical Solutions of the Radiated Velocity Potentials. Similar to the derivation of the solutions of the diffracted velocity potentials, the application of the method of separation variables to the governing equation (22) under the boundary conditions (23), (24), and (25) yields the radiated velocity potentials as follows:

[mathematical expression not reproducible] (26)

[mathematical expression not reproducible] (27)

[mathematical expression not reproducible] (28)

[mathematical expression not reproducible] (29)

[mathematical expression not reproducible] (30)

where

[[gamma].sub.m] = m[pi]/2a = (m = 1, 2, 3, ...) (31)

Combining the floating body surface condition at z = -s, -a [less than or equal to] x [less than or equal to] a and the orthogonality of cos[[[gamma].sub.m](x - a)] over the integral interval of -a [less than or equal to] x [less than or equal to] a, the expression of e can be derived as

[epsilon] = cos 2a[[gamma].sub.m]/sinh [[[gamma].sub.m](h - s)] a[[gamma].sub.m.sup.3] (32)

2.4. Determination of the Unknown Coefficients That Appeared in the Solutions of Radiated Velocity Potentials. The solutions of the unknown coefficients that appeared in the radiation potentials can be again obtained by utilizing the matching conditions at x = [+ or -]a under the principle of continuous pressure and normal velocity. We have the following.

(a) Matching Conditions for the Surge-Induced Potential

[mathematical expression not reproducible] (33a)

[mathematical expression not reproducible] (33b)

[mathematical expression not reproducible] (33c)

(b) Matching Conditions for the Heave-Induced Potential

[mathematical expression not reproducible] (34a)

[mathematical expression not reproducible] (34b)

[mathematical expression not reproducible] (34c)

(c) Matching Conditions for the Pitching-Induced Potential

[mathematical expression not reproducible] (35a)

[mathematical expression not reproducible] (35b)

[mathematical expression not reproducible] (35b)

For the porous walls A and B at the boundaries x = [+ or -]a, the condition of continuity of normal velocity for flows passing through the porous walls A and B leads to

[mathematical expression not reproducible] (36a)

[mathematical expression not reproducible] (36b)

where the porous flow velocities [W.sub.A](z, t) and [W.sub.B](z, t) are assumed to obey Darcy's law; that is, they are linearly proportional to the pressure difference induced by the radiated motions across their corresponding porous wall (Chwang [11]). With the application of the linearized Bernoulli equation, we have

[mathematical expression not reproducible] (37a)

[mathematical expression not reproducible] (37b)

Utilizing the orthogonality property of eigenfunctions, cos [[??].sub.n](z + h), for [[phi].sup.(L).sub.R1] and [[phi].sup.(L).sub.R3] (L = 1, 2, 3) over the integral interval of -h [less than or equal to] z [less than or equal to] 0, the following equations can be formulated by the matching conditions of (33b), (33c), (34b), (34c), (35b), and (35c) as

[mathematical expression not reproducible] (38)

[mathematical expression not reproducible] (39)

[mathematical expression not reproducible] (40)

[mathematical expression not reproducible] (41)

[mathematical expression not reproducible] (42)

[mathematical expression not reproducible] (43)

Additionally, by combining (36a), (36b), (37a), and (37b) with the matching conditions given in (33a), (34a), and (35a), respectively, and applying the orthogonality property of eigenfunctions, cos [r.sub.n](z + h), for [[phi].sup.(L).sub.R2] (L = 1, 2, 3), over the integral interval of -h [less than or equal to] z [less than or equal to] -s, we have

[mathematical expression not reproducible] (44)

[mathematical expression not reproducible] (45)

where

[mathematical expression not reproducible] (46)

[mathematical expression not reproducible] (47)

[mathematical expression not reproducible] (48)

[mathematical expression not reproducible] (49)

[mathematical expression not reproducible] (50)

[mathematical expression not reproducible] (51)

Again, the summation of terms with the infinite series is approximated by a large but finite numbers of terms, e.g., a total of M terms. Then, the unknown coefficients [A.sup.(1).sub.m], [C.sup.(1).sub.m], [E.sup.(1).sub.m], [B.sup.(1).sub.m] (m = 0, 1, 2, M) [4(M+1) unknowns totally] can be determined by solving the 4(M+1) system of algebraic equations formulated in (38), (39), (44), and (45) (for L=1). Similarly, the unknown coefficients [A.sup.(2).sub.m], [C.sup.(2).sub.m], [E.sup.(2).sub.m], [B.sup.(2).sub.m] (m = 0, 1,2, ..., M) can be calculated from the 4(M+1) system of algebraic equations expressed in (40), (41), (44), and (45) (for L=2) while the coefficients [A.sup.(3).sub.m], [C.sup.(3).sub.m], [E.sup.(3).sub.m], [B.sup.(3).sub.m], (m = 0, 1, 2, ..., M) can be evaluated by solving the 4(M+1) system of algebraic equations formulated in (42), (43), (44), and (45) (for L=3). After the determination of the above unknown coefficients, a series of engineering properties important to practical applications could be calculated, for example, the hydrodynamic coefficients and the motion induced dynamic responses of the breakwater system.

3. Exciting Forces and Hydrodynamic Coefficients

3.1. Exciting Forces. To determine the hydrodynamic coefficients related to the induced motions of the floating breakwater system, the exciting forces and moment of the combined incident and diffracted waves on the breakwater need to be computed as inputs for solving the equation of motion. By integrating the dynamic pressures with the use of the velocity potentials of incident and diffracted waves (see (7)(9)) over the wet surface, the hydrodynamic forces along the x direction ([F.sub.D1]) and along the z direction ([F.sub.D2]) can be formulated, respectively, as (Qiao et al. [24])

[mathematical expression not reproducible] (52)

[mathematical expression not reproducible] (53)

Also, the wave exciting moment (FD3) that is referenced to a rotational center (0, [[xi].sub.z]) can be derived as

[mathematical expression not reproducible] (54)

3.2. Hydrodynamic Coefficients. The radiation forces induced by the translational and rotational motions of the floating breakwater system can be calculated from the derived radiation potentials utilizing the linearized Bernoulli equation as

[mathematical expression not reproducible] (55)

where [m.sub.L,k] is the radiation-induced added mass and the [d.sub.L,k] denotes the radiation damping. The expressions are, respectively,

[mathematical expression not reproducible] (56)

[mathematical expression not reproducible] (57)

[[??].sub.Lk] can be computed with formulations given below

[mathematical expression not reproducible] (58)

[mathematical expression not reproducible] (59)

[mathematical expression not reproducible] (60)

[mathematical expression not reproducible] (61)

3.3. Equation of Motion and RAOs. The derived solutions for the diffraction and radiation problems could now be integrated with the equation of motion of the floating breakwater system to analyze the dynamic responses of the structural system. We have

{-[[omega].sup.2] [M + m] - i[omega] [d] + [s]} [[X.sub.0]] = [[F.sub.0]] (62)

where [[[X.sub.0]] = [[X.sub.01], [X.sub.02], [X.sub.03]].sup.T] and are vectors indicating the motion amplitudes and the wave exciting forces, respectively. [m] and [d] are added mass and damping matrices, [M] is the mass matrix, and [S] denotes the hydrostatic restoring force matrix; they can be expressed as

[mathematical expression not reproducible] (63)

[mathematical expression not reproducible] (64)

Here [m.sub.0] is the mass per unit length of the floating breakwater and [z.sub.G] is the z-coordinate of the center of gravity. The mass moment of inertia about y-axis I can be computed by I = [m.sub.0][r.sub.G.sup.2], in which [r.sub.G] is the radius of gyration in pitching oscillation and [bar.GM] is the metacentric height of the floating breakwater system. Once [X.sub.0j] (j=l, 2, 3) are determined by solving (62), the RAOs (Response Amplitude Operators) of the breakwater system can be calculated as

[RAO.sub.j] = [X.sub.0j]/[H.sub.j] (j = 1, 2, 3) (65)

4. Results and Discussion

The validation of the wave diffraction solutions in terms of reflection/transmission coefficients and hydrodynamic loads for a train of monochromatic wave propagating through the proposed floating breakwater system that consists of a partially submerged rectangular body and with or without the attached two porous side walls has been carried out in a previous study (Qiao et al. [24]). Here, additional comparisons of wave diffraction results with other published numerical solutions for the limiting case of neglecting the effect of two porous walls are again presented. For the wave radiation problem, the motion induced hydrodynamic coefficients, such as added mass and damping coefficient, are also investigated for their variations with changes of wave parameters. As a validation, those results for the cases without porous walls are compared with limited numerical solutions published in literatures. Subsequently, the effects of various dimensionless parameters, for example, porous properties, incident wave conditions, and the geometry of floating body, to wave-induced hydrodynamic forces, hydrodynamic coefficients, and RAOs were also examined. In our computations, the first 40 terms in the infinite series of the diffraction potentials and radiation potentials are taken, the same as [7]; namely, M=40.

As an example showing the results comparisons for the diffraction problem, Figure 2 illustrates the variations of the dimensionless exciting forces in horizontal and vertical directions versus dimensionless wavelength, [k.sub.0]s, for the case of h/s = 2.0, a/s = 1.0, [G.sub.01] = [G.sub.02] [right arrow] [infinity]. The exciting moments with reference to the rotational center of (0,0) are also presented in Figure 2. Here, [G.sub.01] = [G.sub.02] [right arrow] [infinity] represents the case of a floating breakwater system without dual porous side walls. The numerical results computed by Black et al. [4] are also plotted in Figure 2. The comparisons indicate that the exciting forces and moments are well predicted by the present analytical model. Both the analytical solutions and the numerical results reveal that the horizontal exciting forces and the exciting moments show rapidly increasing trend initially until a peak value is reached and, then, are followed by a gradual decrease of the values with the further increase of [k.sub.0]s. However, the vertical forces as shown in Figure 2(b) decrease continuously with an increase in [k.sub.0]s. According to the results given in the previous study by Qiao et al. [24], it indicated that the effectiveness of the porous walls is reflected with the reduction of the vertical hydrodynamic forces and when the breakwater system is rotated with reference to the bottom center of the body, the induced moment can also be reduced. However, due to the added additional length of the porous walls, the horizontal forces are slightly greater than those without porous walls.

Figure 3 presents the variations of the dimensionless added mass and damping coefficients in surge only motion versus [k.sub.0]h for the setting of s/h = 1/3, d/h = 1/2, and a/h = 1/6 under the conditions of various porous-effect parameters. When [G.sub.01] = [G.sub.02] = [infinity], it represents the case without considering the effect of dual porous walls. The present analytical solutions for this limiting case can be adopted to compare to the limited numerical results available in the study of [7]. The comparison plot is shown in Figure 4 where the nearly identical matches between the present analytical solutions and the numerical results of [7] can be observed. The results shown in Figure 4 suggest that the addition of porous walls tends to enhance the added mass and damping coefficients. In general, with a decrease of porous effect parameter, the added mass and damping coefficients increase. For the variation of added mass coefficient, with porous walls, it can be seen that a peak value occurred at a long wave (small [k.sub.0]h value) condition. However, for the radiation damping coefficient, it is shown to increase rapidly to its maximum value with an increase of [k.sub.0]h; then, a gradual decrease trend is followed by a further increase of [k.sub.0]h.

The dimensionless heave-heave added mass ([m.sub.22]/2[rho]as) and radiation damping coefficient ([d.sub.22]/2[rho][omega]as) versus [k.sub.0]h under various [G.sub.01] and [G.sub.02] values for the case of s/h = 1/3, d/h = 1/2, a/h = 1/6 are presented in Figure 4. Under the limiting condition of [G.sub.01] = [G.sub.02] = [infinity] (without porous walls), the numerical results produced by the study of [7] are also plotted for comparison. Again, good agreements between the present analytical solutions and the published numerical results can be noticed. The results in Figure 4(a) suggest that the variations of added mass for the shorter wavelength waves (large [k.sub.0]h values) are less sensitive to the change of porous parameters. The presence of dual porous walls increases the heave-heave added mass, which tends to reach its maximum value at the long wave extreme. Also, a decrease in porous parameters results in an increase in heave induced added mass. As shown in Figure 4(b), the heave-heave damping coefficients decrease rapidly when [k.sub.0]h increases initially. Then, a gradual changing trend is followed by a further increase in [k.sub.0]h. Additionally, the presence of dual porous walls tends to reduce the values of heave-heave damping coefficients.

For the pitching induced added mass and damping coefficients, the results are plotted against [k.sub.0]h in Figure 5 for various porous-effect parameters, [G.sub.01] and [G.sub.02], and under the settings of s/h = 1/3, d/h = 1/2, a/h = 1/6. When the rotational center is set at (0,0), as shown in Figure 5, the present analytical solutions match well with the published numerical results [7] for the limiting case of no porous wall (i.e., [G.sub.01] = [G.sub.02] = [infinity]). It can be observed that the pitching oscillation induced added mass and damping coefficients for the floating breakwater system with dual porous side walls are greater than those without porous walls. With decrease of the porous effect parameters, the added mass and damping coefficients increase. For the added mass coefficients, it is found that the variations trends are more complex for waves with relatively smaller [k.sub.0]h under the conditions of different porous-effect parameters. Generally, in cases of weaker porous-effect (for example, [G.sub.01] = [G.sub.02] [less than or equal to] 2 as shown in Figure 5(a)), the value of added mass coefficient tends to approach 1.4 under the case considered. For the variations of damping coefficient as shown in Figure 5(b), the varying trends are similar to those of damping coefficient in surge motion.

The results presented in Figure 6 are the variations of the dimensionless surge-pitch added mass and damping coefficients versus [k.sub.0]h for varying [G.sub.01] and [G.sub.02]. The case is for s/h = 1/3, d/h = 1/2, a/h = 1/6 with reference to the rotational center of (0,0). Clearly, the added mass and damping coefficients are shown to increase with an increase in the porous-effect parameters of both porous walls A and B, while the damping coefficients reach their absolute maximum values when [k.sub.0]h is within the range of 1 and 2. The magnitude of surge-pitch hydrodynamic coefficients, especially the damping, is greatly influenced by the permeability of the dual porous walls compared to the situation with surge or pitching motion only.

The analytical solutions for the Response Amplitude Operators (RAOs) of the floating breakwater system proposed in this study are illustrated in Figure 7 as a function of dimensionless wavelength [k.sub.0]h for various porous-effect parameters with the conditions of s/h = 0.3, d/h = 0.5, a/h = 0.25, [r.sub.G]/a = 0.9, and GM/a = 0.1, where the rotational center is located at (0,0). It is found that the varying trends of surge RAOs are similar to those of pitching RAOs, although the peak values of surge RAOs are smaller than those of pitching RAOs. When the porous-effect parameters decrease, the main peak values of either surge or pitching RAOs that occurred are shifted to the wave conditions with shorter wavelength. Furthermore, for the heave RAO, it experiences a typical damped response for the case of single degree of freedom system and its maximum values of RAO would move to the lower frequency wave region when the porous-effect parameters decrease.

5. Summary and Conclusions

The hydrodynamic performance investigations, including the hydrodynamic loads, radiation induced added mass and damping coefficients, and dynamic responses related RAOs, are critical to the study that is associated with the problems of wave and structure interactions. In this study, the hydrodynamic problems considering a floating breakwater system that consists of a rectangular body and two vertical attached porous walls are investigated analytically based on the two-dimensional potential flow theory. The expressions of the diffraction and radiation based velocity potentials and the associated unknown coefficients for each subdomain are derived by the MEEM approach. Subsequently, the exciting hydrodynamic forces, added mass and radiation damping coefficients in terms of surge, heave and pitching motions, and the equation of motion for the floating breakwater system are formulated and calculated for various case scenarios by varying the wave conditions and porous effect parameters. Considering the limiting case without porous walls, good agreements between the present analytical solutions and published numerical results are achieved. This suggests that the present analytical model is demonstrated to be a reliable engineering tool to provide reasonable estimations for the hydrodynamic performances of a floating breakwater system with dual side porous walls. It is worthy to notice that the wave loads on the floating structures, especially the moments when referenced to the rotational center of bottom center, would be reduced with an increase in the porous-effect parameters of the porous walls. Additionally, the existence of the dual porous walls could move the peaks of surge and pitching RAOs to the waves of higher frequency when the porous effect parameters decrease. Overall, it can be concluded that a properly adopted porous-wall system when attached to a floating body can effectively improve its hydrodynamic performance as a breakwater. The analytical approach developed in this study can also be further extended to assist the theoretical analysis for other similar but improved floating breakwater systems.

https://doi.org/10.1155/2018/1295986

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by "Fundamental Research Funds for the Central Universities" (Grand No. 017180224).

References

[1] I.-H. Cho, "Transmission coefficients of a floating rectangular breakwater with porous side plates," International Journal of Naval Architecture and Ocean Engineering, vol. 8, no. 1, pp. 53-65, 2016.

[2] H. Hu, K.-H. Wang, and A. N. Williams, "Wave motion over a breakwater system of a horizontal plate and a vertical porous wall," Ocean Engineering, vol. 29, no. 4, pp. 373-386, 2001.

[3] C. C. Mei and J. L. Black, "Scattering of surface waves by rectangular obstacles in water of finite depth," Journal of Fluid Mechanics, vol. 38, no. 3, pp. 499-511, 1969.

[4] J. L. Black, C. C. Mei, and M. C. G. Bray, "Radiation and scattering of water waves by rigid bodies," Journal of Fluid Mechanics, vol. 46, no. 1, pp. 151-164, 1971.

[5] N. Drimer, Y. Agnon, and M. Stiassnie, "A simplified analytical model for a floating breakwater in water of finite depth," Applied Ocean Research, vol. 14, no. 1, pp. 33-41, 1992.

[6] J.-F. Lee, "On the heave radiation of a rectangular structure," Ocean Engineering, vol. 22, no. 1, pp. 19-34, 1995.

[7] Y. H. Zheng, Y. G. You, and Y. M. Shen, "On the radiation and diffraction of water waves by a rectangular buoy," Ocean Engineering, vol. 31, no. 8-9, pp. 1063-1082, 2004.

[8] E. D. Christensen, H. B. Bingham, A. P. Skou Friis, A. K. Larsen, and K. L. Jensen, "An experimental and numerical study of floating breakwaters," Coastal Engineering Journal, vol. 137, pp. 43-58, 2018.

[9] E. Koutandos, P Prinos, and X. Gironella, "Floating breakwaters under regular and irregular wave forcing: Reflection and transmission characteristics," Journal of Hydraulic Research, vol. 43, no. 2, pp. 174-188, 2005.

[10] G. H. Dong, Y. N. Zheng, Y. C. Li, B. Teng, C. T. Guan, and D. F. Lin, "Experiments on wave transmission coefficients of floating breakwaters," Ocean Engineering, vol. 35, no. 8-9, pp. 931-938, 2008.

[11] A. T. Chwang, "A porous-wavemaker theory," Journal of Fluid Mechanics, vol. 132, pp. 395-406, 1983.

[12] Y. C. Li, Y. Liu, and B. Teng, "Porous effect parameter of thin permeable plates," Coastal Engineering Journal, vol. 48, no. 4, pp. 309-336, 2006.

[13] S. R. Manam and M. Sivanesan, "Scattering of water waves by vertical porous barriers: an analytical approach," Wave Motion. An International Journal Reporting Research on Wave Phenomena, vol. 67, pp. 89-101, 2016.

[14] A. T. Chwang and Z. Dong, "Wave-trapping due to a porous plate," in Proceedings of 15 ONR symposium on Naval Hydrodynamics, pp. 407-414, 1985.

[15] S. Zhu and A. T. Chwang, "Analytical study of porous wave absorber," Journal of Engineering Mechanics, vol. 127, no. 4, pp. 326-332, 2001.

[16] H. Behera and T. Sahoo, "Hydroelastic analysis of gravity wave interaction with submerged horizontal flexible porous plate," Journal of Fluids and Structures, vol. 54, pp. 643-660, 2015.

[17] T. L. Yip and A. T. Chwang, "Perforated wall breakwater with internal horizontal plate," Journal of Engineering Mechanics, vol. 126, no. 5, pp. 533-538, 2000.

[18] Y. Liu and H.-J. Li, "Wave scattering by dual submerged horizontal porous plates: Further results," Ocean Engineering, vol. 81, no. 5, pp. 158-163, 2014.

[19] I. H. Cho and M. H. Kim, "Wave absorbing system using inclined perforated plates," Journal of Fluid Mechanics, vol. 608, pp. 1-20, 2008.

[20] K.-H. Wang and X. Ren, "Water waves on flexible and porous breakwaters," Journal of Engineering Mechanics, vol. 119, no. 5, pp. 1025-1047, 1993.

[21] K. H. Wang and X. Ren, "Wave interaction with a concentric porous cylinder system," Ocean Engineering, vol. 21, no. 4, pp. 343-360, 1994.

[22] F. Zhao, T. Zhang, R. Wan, L. Huang, X. Wang, and W. Bao, "Hydrodynamic loads acting on a circular porous plate horizontally submerged in waves," Ocean Engineering, vol. 136, pp. 168-177, 2017.

[23] A. Mondal and R. Gayen, "Wave interaction with dual circular porous plates," Journal of Marine Science and Application, vol. 14, no. 4, pp. 366-375, 2015.

[24] W. L. Qiao, K.-H. Wang, and Y. Q. Sun, "Scattering of water waves by a floating body with two vertically attached porous walls," Journal of Engineering Mechanics, vol. 144, no. 2, Article ID 04017162. 2018.

[25] S. W. Twu and D. T. Lin, "On a highly effective wave absorber," Coastal Engineering Journal, vol. 15, no. 4, pp. 389-405, 1991.

[26] X. P. Yu, "Diffraction of water waves by porous breakwaters," Journal of Waterway, Port, Coastal, and Ocean Engineering, vol. 121, no. 6, pp. 275-282, 1995.

[27] I. J. Losada, M. A. Losada, and A. Baquerizo, "An analytical method to evaluate the efficiency of porous screens as wave dampers," Applied Ocean Research, vol. 15, no. 4, pp. 207-215, 1993.

[28] G. Besio and M. A. Losada, "Sediment transport patterns at Trafalgar offshore windfarm," Ocean Engineering, vol. 35, no. 7, pp. 653-665, 2008.

[29] H. Song and L. Tao, "Short-crested wave interaction with a concentric porous cylindrical structure," Applied Ocean Research, vol. 29, no. 4, pp. 199-209, 2007

[30] L. Xiao, Y. Kou, L. Tao, and L. Yang, "Comparative study of hydrodynamic performances of breakwaters with double-layered perforated walls attached to ring-shaped very large floating structures," Ocean Engineering, vol. 111, pp. 279-291, 2016.

[31] T. L. Yip and A. T. Chwang, "Water wave control by submerged pitching porous plate," Journal of Engineering Mechanics, vol. 124, no. 4, pp. 428-434, 1998.

Weiliang Qiao [ID], (1) Keh-Han Wang, (2) Wenqi Duan, (1) and Yuqing Sun (1)

(1) Marine Engineering College, Dalian Maritime University, Dalian 116026, China

(2) Dept. of Civil and Environmental Engineering, Univ. of Houston, Houston, TX 77024, USA

Correspondence should be addressed to Weiliang Qiao; qiao_wl@hotmail.com

Received 17 September 2018; Accepted 6 November 2018; Published 19 November 2018

Academic Editor: Roman Lewandowski

Caption: Figure 1: Schematic diagram of a floating body with attached two porous walls under the action of an incident wave and forced surge, heave, and pitching motions.

Caption: Figure 2: Comparison of the dimensionless wave exciting forces and moment versus [k.sub.0]s for h/s = 2.0, a/s = 3.0, [G.sub.01] = [G.sub.02] [right arrow] [infinity]: (a) dimensionless horizontal exciting force; (b) dimensionless vertical exciting force; (c) dimensionless exciting moment about the rotational center of (0,0).

Caption: Figure 3: Variation of the dimensionless surge-surge hydrodynamic coefficients ((a) for the added mass and (b) for the radiation damping) versus [k.sub.0]h under the conditions of different porous-effect parameters for s/h = 1/3, d/h = 1/2, a/h = 1/6.

Caption: Figure 4: Variation of the dimensionless heave-heave hydrodynamic coefficients ((a) for the added mass and (b) for the radiation damping) versus [k.sub.0]h under the conditions of different porous-effect parameters for s/h = 1/3, d/h = 1/2, a/h = 1/6.

Caption: Figure 5: Variation of the dimensionless pitch-pitch hydrodynamic coefficients ((a) for the added mass and (b) for the radiation damping) versus [k.sub.0]h under the conditions of different porous-effect parameters for s/h = 1/3, d/h = 1/2, a/h = 1/6 with the rotational center set at (0, 0).

Caption: Figure 6: Variation of the dimensionless surge-pitch hydrodynamic coefficients ((a) for the added mass and (b) for the radiation damping) versus [k.sub.0]h under the conditions of different porous-effect parameters for s/h =1/3, d/h = 1/2, a/h =1/6 with the rotational center set at (0,0).

Caption: Figure 7: Variation of the dimensionless RAO ((a) for surge RAO, (b) for heave RAO, and (c) for pitch RAO) versus [k.sub.0]h under the conditions of different porous-effect parameters for s/h = 0.3, d/h = 0.5, a/h = 0.25 [r.sub.G]/a = 0.9, GM/a = 0.1 with rotation center at (0,0).

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Qiao, Weiliang; Wang, Keh-Han; Duan, Wenqi; Sun, Yuqing |

Publication: | Mathematical Problems in Engineering |

Date: | Jan 1, 2018 |

Words: | 6662 |

Previous Article: | Stochastic Collocation Applications in Computational Electromagnetics. |

Next Article: | Contractual Efficiency of PPP Infrastructure Projects: An Incomplete Contract Model. |