Mitigation of Hydroelastic Responses in a Very Large Floating Structure by a Connected Vertical Porous Flexible Barrier

: The hydroelastic response of an elastic thin plate combined with a vertical porous ﬂexible plate ﬂoating on a single- or a two-layer ﬂuid is analyzed in the two-dimensional Cartesian coordinate system. The vertical and the horizontal plates are placed in an inverted-L shape and rigidly connected together. The problem is studied with the aid of the method of matched eigenfunction expansions within the framework of linear potential ﬂow theory. The ﬂuid is assumed to be inviscid and incompressible, and the motion is assumed to be irrotational. Time–harmonic incident waves of the traveling mode with a given angular frequency are considered. Then, the least-squares approximation method and the inner product are used to obtain the expansion coefﬁcients of the velocity potentials. Graphical results show the interaction between the water waves and the structure. The effects of several physical parameters, including the length and the complex porous-effect parameter of the vertical plate, on the wave reﬂection and transmission are discussed. The results show that a vertical plate can effectively eliminate the hydroelastic response of the very large ﬂoating structure. The longer a vertical plate is, the more waves are reﬂected by the vertical plate. With the increase in the porous-effect parameter, the deﬂection of vertical plate decreases. Besides the effects of the ﬂexural rigidity, the lateral stress, the mooring line angle, the ﬂuid density ratio, and the position of interface on the wave reﬂection and transmission are discussed. Numerical results show the signiﬁcant mitigation effect due to the presence of the additional vertical plate.


Introduction
Very large floating structures (VLFSs) have broad research prospects in the utilization of ocean space to deal with insufficient land. VLFSs can modularly be constructed for airports, storage facility, military base, etc. Thus, a VLFS is potentially able to alleviate the land scarcity caused by urban development. The large floating ice sheets and floes also play a similar role to build polar research bases and ice airports. According to the arguments from [1][2][3][4], VLFS have following advantages: long design life, low environmental impact, and ease of expansion and mobility. Due to the characteristics that a VLFS can follow the motion of a wave, which may cause structural fatigue and failure, it is necessary to install additional structures for suppressing the response of VLFS.
According to [5], excessive deformation of the VLFS would disrupt its feasibility as a floating runway or a floating storage facility. Moreover, it is necessary to eliminate low-frequency vibration of residential VLFS to ensure the comfort of residence. In order to suppress the hydroelastic responses of a VLFS, several methods have been introduced by Wang et al. [5] such as different kinds of breakwaters and the submerged plate antimotion devices. A simple and effective anti-motion device is a submerged vertical plate attached with a VLFS, and it has been widely studied to suppress the hydroelastic response of a VLFS. Takagi et al. [6] concluded that, when a box-shaped body is attached to an edge of a VLFS, it reduces not only the deformation but also the shearing force and the moment of the platform, also attaining a good anti-motion performance. Masanobu et al. [7] proposed the additional wall with slits and the inverted-L type structures to mitigate the hydroelastic response of the VLFSs. They found that the drift force and the vertical displacement can be suppressed by providing slits for the anti-motion device. A composite grid method is adopted by Lee et al. [8] to solve the hydrodynamic forces and the nonlinear behaviors of fluid motion around a submerged plate. The added mass and damping forces on a VLFS were increased due to the generated vortex by the submerged plate, hence the structural responses were decreased. Taking a circle plate, for example, Pham et al. [9] utilized the modal expansion method to analyze the hydroelastic response of a pontoon-type circular VLFS attaching a horizontal submerged annular plate. The hydroelastic response reduction efficiency of VLFS edged with dual inclined plates has been investigated by Cheng et al. [10] with the aid of the method of matched eigenfunction expansions (MMEE). The finite element-boundary element method is utilized in research of a Wave Energy Converter (WEC)-type attachment to reduce hydroelastic responses of the VLFS [11,12]. Feng et al. [13] investigated the hydroelastic responses of a submerged horizontal solid/porous plate connected to a VLFS. On basis of understanding the above research, we think it is worth studying to suppress the deformation of VLFS with an additional vertical flexible plate.
Due to the feature that the horizontal size of the VLFS is much larger than its thickness, the flexible deformation of VLFS should be fully considered, while the rigid-body motion can be neglected. The pontoon-type VLFS was regarded by Ohmatsu [14] as a huge flat plate floating over the surface of water. Squire [15] claimed that there is a conspicuous overlap between the research in VLFS and in ice sheets. Thus, the VLFS model is also available for the analysis of floating ice sheets by using parameters of the ice. The MMEE was utilized to study the interaction between surface waves and an ice sheet by Fox and Squire [16], in which the unknowns coefficients were obtained by the error function method.
The key of MMEE, which is often used in hydroelastic problems, is to obtain the expansion coefficients. Another method to find the expansion coefficients is to derive an orthogonal inner product in the plate-covered region [17]. Furthermore, when the physical model has an open region, the orthogonality of the eigenfunctions in the open region can be utilized. Xu and Lu [18] directly dealt with the problem by using the inner product which involves the orthogonal eigenfunctions in open region. In a two-layer fluid model, the dispersion relation have two pairs of real roots which represent the surface and interfacial traveling wave modes, and the interaction between surface waves and interface waves is important. Xu and Lu [19] proposed a new inner product which utilizes the orthogonality of the eigenfunctions in the open region on a two-layer fluid. Meng and Lu [20] investigated the hydroelastic response of an elastic plate with lateral internal forces in a small amplitude waves on a three-layer fluid.
The porous structures which can dissipate a part of wave energy are studied to suppress the waves. Yu and Chwang [21] suggested a newly derived boundary condition with a complex porous-effect parameter for the porous breakwaters in a semi-circular harbor. Yip et al. [22] investigated the trapping of surface waves by submerged vertical porous and flexible barriers near the end of a semi-infinitely long channel of finite depth. The full reflection always occurs when the distance between the channel end-wall and the barrier is an integer multiple of half-wavelength. Kumar and Sahoo [23] studied the performance of a flexible porous plate breakwater in a two-layer fluid via the least-squares approximation method which is widely used in research of vertical barrier. The highest wave reflection and the lowest wave transmission amplitude in the surface and interface modes can be observed for the case which has zero porosity. The wave reflection and transmission in a two-layer fluid was found to be strongly dependent on the location of interface and the fluid density ratio. Mandal et al. [24] studied oblique wave scattering by multiple porous, flexible barriers in a two-layer fluid, and the condition for the Bragg resonance was derived in the case of multiple barriers. Singla et al. [25] analyzed the effectiveness of vertical solid permeable barriers located at a finite distance from a VLFS. The porosity of the structure can reduce 15-20% of the wave energy. However, the vertical flexible porous plate to mitigate the response of VLFS has not been considered in previous studies. Thus, a porous structure is also included in our research because it can dissipate the wave energy.
In the present paper, we investigate the hydroelastic responses of a semi-infinite floating plate combined with a vertical porous flexible plate under the wave action on a single-layer or a two-layer fluid of finite depth with the aid of the MMEE. The floating plate is much larger than the vertical plate. We neglect the rigid body motion of a VLFS. We assume that there is no gap between the horizonal elastic plate and the surface of water. The pressure from the elastic plate is included in the boundary conditions. Our focus is on how the vertical plates affects the hydroelastic response of VLFS. The inner product and the least square approximation are used to obtain a set of simultaneous equations for different boundary conditions, and the edge conditions are included as a part of the equation system. The conservation of energy is used to validate method of solution. In Section 2, we study the interaction between waves with semi-infinite or finite plates on a single-layer fluid. In Section 3, we extend the method to the case of a two-layer fluid. Numerical calculations for the analytical expressions and graphical results are performed in Section 4. Finally, conclusions are given in Section 5. We consider the hydroelastic interaction between incident gravity waves and a semiinfinite floating elastic thin plate combined with a vertical finite porous flexible plate in a single-layer fluid bounded by a flat rigid seabead, as shown in Figure 1. A twodimensional Cartesian coordinate system oxz is chosen in such a way that the x-axis points the horizontally rightward and the z-axis point the vertically upward. The fluid with the constant density ρ occupies the region −H < z < 0 with z = 0 being the undisturbed upper surface and z = −H a flat bottom. The horizontal semi-infinite elastic plate which floats on the fluid with zero draft covers the region 0 ≤ x < ∞. The vertical elastic plate is fixed with the horizontal plate at z = 0. The response of vertical plate is analyzed by assuming that the plate behaves as a one-dimensional beam. A clamped-mooring edge condition is used to describe the connection with the floating plate.  The fluid is assumed to be inviscid and incompressible and the fluid motion is considered irrotational. In the case of the simply time-harmonic motion with a frequency ω, we can decompose the time factor and write the velocity potential Re φ(x, z)e −iωt and the surface elevation Re ζ(x)e −iωt , where φ(x, z) is the spatial velocity potential, ζ(x) the spatial elevation on the surface and t the time variable. The whole flow domain is divided into two parts: the open water region (−∞ < x < 0) and the plate-covered region (0 < x < ∞), for which the corresponding spatial velocity potentials are denoted by φ L (x, z) and φ R (x, z), respectively. Thus, the governing equation is

Hydroelastic Response on a Single-Layer
The boundary condition at the flat bottom is Within the framework of linear theory for small-amplitude waves, the combined kinematic and dynamic condition on the undisturbed surface (z = 0) for the open water region is given by where g is the gravitational acceleration.
Under the assumption of no gap between the fluid and the horizontal plate, the combined kinematic and dynamic conditions for the plate-covered region on z = 0 can be written as where M e = ρ e d, D = E d 3 /[12(1 − ν 2 )], and E, d, ν, and ρ e are the flexural rigidity, the effective Young modulus, the constant thickness, Poisson's ratio, and the density of the horizontal plate, respectively. As the vertical barrier is assumed to be a thin elastic plate of uniform mass subject to uniform rigidity, the equation of motion for the vertical plate acted upon by the fluid pressure is given by where ξ = ξ(z) is the horizontal deflection of barrier, M e = ρ e d, D = Ed 3 /[12(1 − ν 2 )], and E, d, ν, L, and ρ e are the flexural rigidity, the effective Young modulus, the constant thickness, Poisson's ratio, the length, and the density of the vertical plate, respectively. Q is related to the lateral stress of the plate (with compression at Q > 0 and stretch at Q < 0). With the aid of the Darcy law for describing the flow past a porous structure, the boundary condition on the vertical plate is given by where is the complex porous-effect parameter defined by Yu and Chwang [21], k 0 the wavenumber of surface progressive wave mode, ε the porosity, f the linearized resistance coefficient, and S the inertial force coefficient. It is assumed that the porous plate has negligible storage capacity.
The matching relations for the continuities of pressure and velocity and along x = 0 are given by The vanishing of deflection and slope at the clamped edge of the horizontal plate are given by With the requirements of the clamped-mooring edge conditions in Mandal et al. [24], zero deflection and zero slope at the fixed edge require the conditions as follows: At the edge near the mooring line, the bending moment is zero and the mooring tension relating to the elastic restoring force is added to the shearing force, which yields where θ is the mooring line angle and K m the stiffness of the mooring line.

Method of Solution
In the framework of the MMEE, the vertical eigenfunction V(k, z) for a single layer fluid reads Thus the spatial velocity potentials can be written as where ζ 0 is the amplitudes of surface incident wave; the reflection coefficients R 0 and R i and the transmission coefficients T 0 and T j are complex numbers to be determined, where (i = 1, 2, 3, . . . ; j = I, II, 1, 2, 3, . . .). Z 0 , Z i ,Z 0 , andZ j the vertical eigenfunctions. The symbols ∑ i and ∑ j are hereinafter defined as the summation signs for i = 1, 2, 3, . . . and j = I, II, 1, 2, 3, . . ., respectively. Let positive real numbers k 0 andk 0 be the wave numbers of progressive wave modes. The purely imaginary numbers k i andk j (i, j = 1, 2, · · · ) are the wave numbers of evanescent modes in the open water and the plate-covered regions, respectively. Furthermore,k I andk II correspond to two decaying progressive waves in the plate-covered regions. For a given frequency ω, the wave numbers k 0 , k i , andk j in Equations (14) and (15) satisfy where Γ = D/ρg and σ = M e /ρ. Equations (18) and (19) are called the dispersion relations for the waves in the open water and the plate-covered regions, respectively. By substituting the velocity potentials in Equations (14) and (15) into Equation (5), a general solution for the governing equation of the vertical plate is of the following form where with the unknown coefficients C n and S n (n = 1, 2) to be determined by the boundary conditions of the vertical plate from Equations (11) and (12). After substituting of the spatial velocity potentials Equations (14) and (15) into the matching relations Equations (8) and (9), it yields: In order to deduce the linear algebraic equations for the unknown expansion coefficients, we employ the inner product of the vertical eigenfunctions for a single layer fluid, which was defined by Xu and Lu [18] as follows, One can easily validate that Applying the inner product on both sides of Equation (25) with the vertical eigenfunction Z i (z) in the open water region, we have where Substituting of the spatial velocity potentials Equations (14) and (15) into the boundary condition on the vertical porous barrier Equation (6), we obtain Applying the least-squares method from [23], it is derived that Minimizing the above integral in Equation (33) with respect to R i , which are the unknowns in the open region, we obtain where the superscript " * " denotes the complex conjugate. We truncate the system such that the numbers of evanescent modes in the open water and plate-covered regions are N, and Equation (34) provides N + 2 linear equations. Substituting for the velocity potentials from Equations (14) and (15) to Equation (10), we obtain two other equations Thus, we derive a system of 2N + 8 equations from Equations (11), (12), (29), (30), (34), and (35) for 2N + 8 unknown coefficient R i , T j , C n , and S n (i = 0, 1, 2, 3, . . . , N; j = 0, I, II, 1, 2, 3, N, . . . ; n = 1, 2). The unknowns in the velocity potentials can be obtained by solving numerically this closed system.

Wave Interaction with a Finite Elastic Plate
In real situations, VLFS can also be modeled as a finite elastic thin plate. In this section, we will investigate the interaction of water waves with a finite floating plate combined with an attached vertical plate. The length of horizontal elastic plate isL, and the horizontal elastic plate has no draft floats on the fluid and covers the region (0 ≤ x ≤L). The whole flow domain is divided into three parts: the left open water region (−∞ < x < 0), the plate-covered region (0 ≤ x ≤L), and the right open water region (L < x < ∞). The corresponding spatial velocity potentials are denoted by φ L (x, z), φ M (x, z), and φ R (x, z), respectively. The governing equation, linearized boundary conditions, and the matching relations for the continuities of pressure and velocity for this case are the same as those given in Section 2.1.
The matching relations for the continuities of pressure and velocity at x = L are given by The clamped-free edge conditions of the horizontal plate read According to the method of eigenfunction expansion, the spatial velocity potentials φ L (x, z), φ M (x, z), and φ R (x, z) can be written as where R L i , R M j and T M j , T R j (i = 0, 1, 2, 3, . . . ; j = 0, I, II, 1, 2, 3, . . .) are the reflection and transmission coefficients to be determined, respectively. The derivation and calculation for the unknown coefficients R L i , R M j and T M j , T R j can be found in Appendix A.

Wave Interaction on a Two-Layer Fluid
Compared with the case of a single-layer fluid, the two-layer fluid has an interface which divides the upper and lower layers at z = −h. The densities of the upper and lower fluids are denoted by ρ 1 and ρ 2 , respectively. Next, we choose the case of the semi-infinite plate as an example, as shown in Figure 2. The matching relations of velocity and pressure at the interface (−∞ < x < ∞, z = −h 1 ) are given by where γ = ρ 1 /ρ 2 and K = ω 2 /g with 0 < γ < 1.
According to the method of eigenfunction expansion, φ L (x, z) and φ R (x, z) can be written as where The terms R 0 m , R i and T 0 m T j (i = 1, 2, · · · ; j = I, II, 1, 2, · · · ) are, respectively, the reflection and transmission coefficients to be determined, whereas ζ 0 and η 0 are the amplitudes of incident waves on the surface and interface, respectively. The vertical eigenfunction V(k, z) in Equation (52) for a two-layer fluid is defined by Xu and Lu [19] as follows The inner product on the two-layer fluid is from [19] The derivation and calculation for the unknown coefficients When a VLFS is modeled as a finite elastic thin plate, φ L (x, z), φ M (x, z) and φ R (x, z), can be written as We can also use the least-squares method as Equations (A19) and (A22) in which R i is replaced with R L i , for m = 1, 2 (1 + e ik 0nL )ωã 0 n + e ik 0nL Gk 0 1 − e ik 0nLk 0 n R M

Rate of Energy Flux
We take the case on the two-layer fluid as an example. For the sake of clarity, the total depth of the two-layer fluid H, the gravitational acceleration g and the density of the upper fluid ρ 1 are chosen as the characteristic quantities to non-dimensionalize relative quantities. The nondimensional variables and parameters arê The symbol "ˆ" over the nondimensional variables and parameters will be dropped hereinafter for clarity. The parameters used for the computation are d = d = 0.01, ρ e = ρ e = 0.9, D = D = 0.05, K m = 0.01, and θ = 45 • . We assume that the incident waves predominantly consisting of surface and interfacial modes have the same frequency ω.
According to the kinematic boundary conditions, the amplitudes of surface and interfacial elevations, A S (x) and A I (x), can be given by The amplitudes of the bending moment M(x) and the shear force Q(x) on the horizontal plate can be given by The amplitudes of the bending moment M(z) and the shear force Q(z) on the vertical plate can be given by To validate the above-mentioned method of solution, numerical examples are presented here. The conservation of the energy for the two-layer fluid under consideration can be written as where Φ m (x, y, z, t) is the velocity potential for the upper (m = 1) and lower (m = 2) fluids, and the superscripts "O" and "P" denote the open water and plate-covered regions, respectively.

Response on a Single-Layer Fluid
In the present subsection, we consider the case of the semi-infinite plate model on a single-layer fluid, and set ζ 0 = 0.01, ω = 1.25, G = 0, L = 0.9, and N = 15 for calculations. The relative error of the energy is Figure 3 illustrates the variation of the amplitudes of surface versus x and the horizontal deflection of the vertical plate versus z for various values of N, where N is the number of terms for the evanescent modes. It is obvious in Figure 3a that the changes of the curves are tiny enough to neglect. Figure 3b displays the horizontal deflection of the vertical plate which decreases slowly with the increment of N. It is depicted that the curves keep almost the same for N ≥ 15 which reveals that the results of surface amplitudes converge quickly and the higher-order terms have a tiny contribution to the solution. Therefore we choose the parameter N = 15 in the following cases.          Figure 8 displays the variations of the bending moments and the shear forces along the vertical and the horizontal plates. It is observed that the maximal amplitudes of the bending moments and the shear forces have a small variation between two plate models in Figure 8a,b. For the finite plate model, the bending moments and the shear forces is zero at the free edges. The reflected wave from the right edge caused the periodic changes in the bending moments and the shear forces. The additional vertical plate significantly decreases the bending moments and the shear forces on the floating plate, and it can reduce the load on the horizontal plate. Due to the clamped-mooring edge conditions, the bending moments and the shear forces decrease from top to bottom of the vertical plate in Figure 8c,d. The bending moments measure zero at the bottom of the vertical plate. Figure 8d shows a small non-zero shear force caused by the mooring line at the bottom of the vertical plate.
The porosity of a submerged plate has a significant effect in mitigating the wave loads on the structure. Figure 9 shows the variations of the amplitudes of surface and the horizontal deflection of the vertical plate with different values of the porous-effect parameter G of the vertical porous flexible plate. The highest wave reflection and lowest wave transmission are observed for the case with G = 0. With the increment of the porouseffect parameter, the amplitudes in the plate-covered region and the right open water region increase, as shown in Figure 9a. Thus, the effect of the vertical plate on reducing the hydroelastic response is weakened. The reason for the phenomenon is that the higher porous-effect parameter G means higher porosity, which indicates that more water can flow through the vertical plate. Figure 9b illustrates that the deflection of the vertical plate near the surface decreases with the porous effect parameter G = 0. With the introduction of porous structure, the wave loads on the vertical board has been reduced, and a part of the total wave energy is dissipated by the porous structure. Thus, in the design of marine structures, the reduction in wave loads and the effect of wave elimination need to be chosen according to the specific situation.     Figure 10a shows that, as the flexural rigidity increases, the amplitude ratio of surface in the region behind the wall decreases, and an opposite trend can be observed on the left open water region. It can be observed that the horizontal deflection and rotating angle of vertical plate decreases while the flexural rigidity D increases in Figure 10b. Thus, the increment of D is mainly affected in the region sheltered by vertical plates.  Figure 11 demonstrates the variations of the amplitudes of surface versus x and the deflection of vertical plate versus z with different values Q for the lateral stress of the vertical plate (with compression at Q > 0 and stretch at Q < 0). In this case, we have Figure 11a displays that the fluctuation of the surface elevations decreases slightly on the plate-covered region when Q becomes bigger. With the increment of Q, the horizontal deflection of vertical plate in the interface increases in Figure 11b. In general, the lateral stress has a little effect to reduce the deformation of the wall, unless the rigidity of the anchor chain is greatly increased.    Figure 12b shows, with the increment of θ, the horizontal deflection of vertical plate decreases. In fact, the effect of the mooring line angle is tiny enough to ignore.

Response on the Two-Layer Fluid
In this section, we analyze the hydroelastic response on a two-layer fluid. We calculate the wave scattering due a semi-infinite plate and set ζ 0 = 0.01, η 0 = 0.0001, γ = 0.9, h = 0.5, G = 0, L = 0.9, ω = 1.25, L = 12, and N = 15. Figure 13 illustrates the variations of the amplitudes of surface and interfacial elevations versus x and the horizontal deflection of vertical plate versus z with different values of the porous-effect parameter G of the vertical porous flexible plate. It is similar to the situation in the single-layer fluid. The highest wave reflection and the lowest wave transmission are observed for the case with G = 0. With the increment of the porous-effect parameter, the amplitude in the plate-covered and right open water region increase. With the introduction of porosity, the deflection of the wall decreases with the increment of porous-effect parameter G from Figure 13c. Thus, the main function of porosity is to reduce the load on the wall, and a medium porosity is a suitable choice to keep the balance between the reduction in the load and the suppression of the waves. Figure 14 shows the variations of the amplitudes of surface and interfacial elevations versus x and the horizontal deflection of vertical plate versus z for different depth h of the interface. It is observed that, as h increases, the interfacial elevations decrease. The reason for this situation is that the wave energy decreases with a decreasing depth, and the wave energy in the interface is correspondingly reduced. The change of the surface amplitude with the increment of h is tiny enough to ignore. The deflection of the vertical plate changes slightly with the increment of h in Figure 14c. In general, the position of the interface has a small effect on the results. Figure 15 illustrates the variations of the amplitudes of surface and interfacial elevations versus x and the horizontal deflection of vertical plate versus z for various values of L in the case of the incident waves of interfacial wave mode, where L is the length of the vertical plates. It reveals, in Figure 15a, that the interfacial incident wave has a small effects on the surface. Figure 15b shows the interfacial waves is mostly reflected by the vertical plate for L = 0.7 when the vertical barrier pierces through the interface. The interfacial waves are almost unhindered for L = 0.3.     Figure 16b shows that, with decrease in γ, the interfacial elevations decrease in the open region and increase in the plate-cover region. The interfacial incident wave has a small effect on the surface for γ = 0.9 in Figure 16a. The amplitudes of the surface significantly increases with the decrement of γ. The deflection of the vertical plate also increases with the decrement of γ, as shown in Figure 16c. Figure 17 demonstrates the variations of the amplitudes of surface and interfacial elevations versus x and the horizontal deflection of vertical plate versus z for various values of γ. Figure 17b shows that with the increase in γ, the interfacial elevations decrease, and it reveals that the density ratio significantly affects the wave amplitude in the interface. The amplitudes of surface remain almost the same for γ = 0.7 and γ = 0.9 in Figure 17a. In the case of γ = 0.5, it is observed that the amplitudes of surface vary periodically due to the interfacial incident wave. The reason for the situation is that the energy of interfacial waves, according to the wave number displayed on the graph, transfers to surface and propagates in surface wave mode. The deflection of the vertical plate has little change due to the increment of γ, as shown in Figure 17c.

Conclusions
We have analytically studied the hydroelastic response of semi-infinite and finite floating elastic plate combined with a vertical elastic porous plate in a single-layer and a two-layer fluids with the aid of the method of matched eigenfunction expansions. The vertical plate is fixed to an edge of the floating plate and moored at its lower edge. The wave reflection and transmission by the structure was found to be strongly dependent on the length L and the porous effect parameter G of the vertical plate. Due to the low wave energy near the seabed, the changes of the curves between L = 0.7 and L = 0.9 is small, hence a further extension of the vertical plate has little effect on the solution. The vertical plate can decrease the bending moments and the shear forces on the floating plate.
The highest wave reflection and lowest wave transmission are observed for the case with the zero porous-effect parameter (G = 0). For a higher porosity of the wall, more water will flow through the vertical plate. Thus, the amplitudes of surface and interfacial elevations increase as the value of G increases. The structural rigidity, the lateral stress and the mooring line angle of the vertical plate mainly affect the response of the vertical plate, and have a modest effect on the response of the floating plate. For a two-layer fluid system, the fluid density ratio and the position of interface mainly effect on the interfacial waves. As the lower position of interface descends, the interfacial amplitude decreases. When the vertical barrier pierces through the interface, the interfacial waves is mostly reflected by the vertical plate. With the decrease in γ, the interfacial incident waves have more effects on the surface. In particular, in the case of γ = 0.5, the amplitudes of surface on the left open water region vary periodically due to the fact that the wave energy transfers from interfacial modes to surface ones and propagates in surface wave mode. The rate of energy flux in the results is approximatively conserved, and the method has high rate of convergence.
These observations are of significant importance in the design for eliminating the vibration of VLFS. Based on the discussion of the results presented here, we can conclude that, in order to keep a balance between the reduction in the wave loads on the wall and the suppression of the deflection on the VLFS, the addition of a submerged plate with the length of L = 0.7 and a medium porous effect parameter is a suitable choice in design.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Derivation Process of the Unknown Coefficients for a Finite Plate Model in a Single-Layer Fluid
The derivation process of the unknown coefficients in the potential function expansion in Equations (41)-(43) is shown in this appendix. Substituting of the velocity potentials from Equations (41)-(43) into Equation (5), a general solution for Equation (5)  (A1) Substituting for the velocity potentials from Equations (41)-(43) into the matching relations Equations (8), (9), and (38), and performing the inner product on both sides of to obtain a set of simultaneous equations, we obtain −k l R L l P ll =k 0 e ik 0L T M T R l P ll = e ik 0L T M l Q 0l + i ∑ j T M j e ik jL Q jl + R M 0 Q 0l + i ∑ j R M j Q jl , (l = 0, 1, 2, · · · ). (A5) Applying the least-squares method, we obtain where the superscript " * " denotes the complex conjugate, (1 + e ik jL )ωã j + e ik jL Gk 0 − e ik jLk j R M jZj + ωW. (A8)

Appendix B. Derivation Process of the Unknown Coefficients in the Two-Layer Fluid
In this appendix, we take the case with a semi-finite plate to show how to dertermine the unknown coefficients in a two-layer fluid. Substituting of Equations (48) and (49) into Equation (8), and performing the inner product on both sides of to obtain a set of simultaneous equations, we have When the vertical barrier pierces through the interface (namely L > h), the deflections of the vertical plate above (−h < z < 0) and beneath (−L < z < −h) the interface are denoted by ξ 1 (z) and ξ 2 (z) as follows The vertical plate is continuous at z = −h, so we have [C mn cosh(κ n z) + S mn sinh(κ n z)]. (A16) Here, C mn and S mn (m = 1, 2; n = 1, 2) are constants determined by the boundary conditions of the vertical plate, and a im andã jm are defined by By applying the least-squares method, we obtain where