Abstract

Effects of underlying uniform current on the nonlinear hydroelastic waves generated due to an infinite floating plate are studied analytically, under the hypotheses that the fluid is homogeneous, incompressible, and inviscid. For the case of irrotational motion, the Laplace equation is the governing equation, with the boundary conditions expressing a balance among the hydrodynamics, the uniform current, and elastic force. It is found that the convergent series solutions, obtained by the homotopy analysis method (HAM), consist of the nonlinear hydroelastic wave profile and the velocity potential. The impacts of important physical parameters are discussed in detail. With the increment of the following current intensity, we find that the amplitudes of the hydroelastic waves decrease very slightly, while the opposing current produces the opposite effect on the hydroelastic waves. Furthermore, the amplitudes of waves increase very obviously for higher opposing current speed but reduce very slightly for higher following current speed. A larger amplitude of the incident wave increases the hydroelastic wave deflections for both opposing and following current, while for Young’s modulus of the plate there is the opposite effect.

1. Introduction

A floating elastic plate may be an appropriate physical model for a very large floating structure (VLFS) in the ocean or an ice sheet in the polar region due to its huge horizontal dimension compared with its vertical thickness. The nonlinear hydroelastic responses, without assuming that the waves are necessarily of small amplitude, become an indispensable factor during designing a VLFS or considering an ice sheet as a storage facility, a mobile offshore base or an aircraft runway, and so on. Forbes (1986) [1] studied preliminarily the nonlinear periodic waves beneath an elastic ice sheet in infinite water depths by the perturbation technology. Forbes (1988) [2] further optimized their previous method in the frame of perturbation technology by introducing Newton-Raphson techniques to approximate the Fourier coefficients. Vanden-Broeck and Părău (2011) [3] tried to extend the study of Forbes (1988) [2] to the large-amplitude periodic waves, while pointing out objectively that the perturbation method is not applicable to very steep waves for which it depended strongly on small parameters. The analytic solutions for the hydroelastic solitary waves propagating under an infinite elastic plate were derived by Milewski et al. (2011) [4], who used the asymptotic and numerical methods. Wang and Lu (2013) [5] investigated the nonlinear hydroelastic progressive waves traveling in an infinite elastic plate in deep water by the homotopy analysis method (HAM), which is completely independent of any small physical parameter [6]. Using the HAM in a similar way, the two-layer fluid was considered by Wang and Lu (2016) [7] to study the effects of the stratification of the fluid on the nonlinear hydroelastic waves traveling in an infinite elastic plate. It is found that a larger density ratio or depth ratio of the fluid layers can reduce the elastic plate deflection.

All the investigations mentioned above were under the hypothesis that there is no current in the oceanic environments. As we know, there are various reasons such as wind, thermal, earth’s rotation, and tidal effects that frequently generate ocean currents. Some scholars analytically investigated nonlinear hydroelastic waves with a current. Employing the linear potential-flow theory, Schulkes et al. (1987) [8] first considered the effect of uniform flow beneath a floating ice plate. Bhattacharjee and Sahoo (2007) [9] used dispersion relation to analyze detailed characteristics of the flexural gravity waves due to a floating elastic plate in the presence of a following current or an opposing current. Asymptotic depressions for the transient flexural gravity resulting from initial disturbances at a point were derived subsequently by Bhattacharjee and Sahoo (2008) [10]. Mohanty et al. (2014) [11] studied a combined effect of current and compressive force on time-dependent flexural gravity wave motion in both the cases of single and two-layer fluids and obtained the asymptotic results for Green’s function and the deflection of ice sheet using the method of stationary phase. Recently, Lu and Yeung (2015) [12] considered unsteady hydroelastic waves caused by the interaction of fixed concentrated line loads and the underlying current and found the flexural gravity wave motion depended on the ratio of current speed to phase or group speeds.

In this paper, we consider nonlinear hydroelastic waves generated due to a floating elastic plate in an underlying uniform current for the two-dimensional case. Two underlying currents are studied with the aid of the HAM. One is following current and the other is opposing current. Convergent analytical series solutions are deduced for the nonlinear hydroelastic response.

2. Mathematical Description

As shown in Figure 1, Cartesian coordinates fixed on the earth are used, in which the -axis coincides with the undisturbed plate-water interface, while the -axis points vertically upward. The floating plate extends to the infinity along the -axis without draft, and an underlying following current is assumed to be moving from the left to right (an opposing current is from the right to left). Under the assumptions that the fluid is homogeneous, incompressible, and inviscid and the flow is irrotational, the governing equation iswhere is the wave velocity potential. is the hydroelastic wave elevation of the wave-plate interface and also is the vertical hydroelastic response of the floating plate. is the time from the start.

The bottom condition at deep water reads

The flexural rigidity of the plate is denoted by , where , , and are Young’s modulus, the constant thickness, and Poisson’s ratio of the plate, respectively. Let be the intensity of the uniform current, and let be the mass of the plate in a unit length, where is the uniform densities of elastic plate. is the uniform densities of the fluid. The nonlinear kinematic and dynamic boundary conditions on the unknown plate-water interface are given by

As for the progressive waves, we introduced an independent variable as follows: where and are the wave number and the angular frequency of the incident wave, respectively. So we can perform the velocity potential function and the hydroelastic wave profile .

To simplify the formulation, the physical variables are nondimensionalized as

By transformation (4) and nondimensionalization (5), (1)–(3) are changed into (after omitting the ) respectively, where A partial combination of (8) and (9) gives the boundary conditions on as follows:

For a given angular frequency or the current speed , the unknown velocity potential and the hydroelastic waves are obtained from (6), (7), (9), and (11) by using the HAM in the subsequent sections.

3. Analytical Approach Based on the HAM

3.1. Solution Expressions

We first predict a set of solution expressions of the unknown nonlinear hydroelastic waves. From physical points of view, unknown velocity potential and progressive hydroelastic wave deflections generated due to a thin elastic plate can be expressed by [5]with the base functions and , respectively, where and are all unknown coefficients to be derived.

According to the solution expressions (12), we can construct the initial estimation of the potential function aswhere is an unknown coefficient. We choose as the initial estimation of the hydroelastic wave deflection to simplify the computation [5, 6]. Considering the nondimensionalized linear dispersion relation [9] for deep water with a current, we think that the initial estimation of the nonlinear dispersion relation satisfies where is the initial estimation of angular frequency and is the initial estimation of the uniform current intensity .

3.2. Zeroth-Order Deformation Equations

Instead of solving these nonlinear equations (6), (7), (9), and (11) directly, we develop four topology homotopies , , , and . As the embedding parameter increases from to and deforms continuously from its initial estimation to the exact solution , varies continuously from to , is from to , and is from to . The homotopies and are governed by the zeroth-order deformation equations for (6), (7), (9), and (11) as follows: where is a nonzero convergence-control parameter. Based on the nonlinear boundary conditions (11) and (9), respectively, and are the nonlinear operators defined by respectively, where We choose the linear operators of the wave profile function in the nonlinear operator as the auxiliary linear operator : where .

Expanding these unknown variables , , , and into the Taylor series about at : where

We choose the linear operators of in as the auxiliary linear operator : where . In particular, we have great freedom to choose the auxiliary linear operator in the HAM, and based on the linear wave theory, the angular frequency and the current speed can be approximated as and in the operator , respectively.

Once is so properly chosen that the series in (23) converge at , we would obtain their convergent homotopy-series solutions.

3.3. High-Order Deformation Equations

Differentiating the zero-order deformation equations (16)–(19) times about , then dividing them by , and setting , we obtain the linear PDES for the unknown variables , , , and , namely, the th-order deformation equations: where . where and is the Heviside step function. The elaborated expressions for , , , and are given in Appendix.

In the first-order calculation, we still can not get the coefficient , so an additional equation is introduced to relate the hydroelastic wave height and the incident wave height as where is an even integer, an odd integer, and the dimensionless amplitude of first-order incident wave.

3.4. Optimal Convergence-Control Parameters

According to Liao [6, 13], it is the convergence-control parameter that essentially differs the HAM from all other analytic techniques. We can obtain the convergent approximate series solutions by choosing the optimal . For given the number of the discrete point, the optimal corresponds to the fastest decrease of the total squared residual defined bywherewhere and are the residual square errors of boundary conditions (18) and (19) at , respectively. . We choose here, and the optimal convergence-control parameter can be chosen by the minimum of in residual plot as shown in Figure 2.

4. Results and Discussion

The total residual square error at several different orders is illustrated for a fluid with the uniform current versus to show the validity and accuracy of our analytic solutions in Figure 2 with the case of , , , , , and . It is found that, as the order increases gradually, reduces to a minimum in the interval and then the optimal is about . Further, we find from Table 1 that, as , the square error decreases quickly to at the 7th-order approximation.

We consider the effects of several important physical parameters, including the uniform current intensity and direction, the amplitude of incident wave, and Young’s modulus of the floating plate, on the nonlinear dynamics characteristic of the hydroelastic waves. Firstly, Figure 3 shows variation of the hydroelastic wave deflection for different current intensity. As the maximum current speed is about in ocean [14], we select four dimensionless current intensities (without currents), , , and in our computations. From the comparison between the following current and the opposing current, it is shown that the amplitudes of the hydroelastic waves decrease slightly with the increment of the following current intensity, while the opposing current has obvious opposite effects. Further, the hydroelastic wave is more greatly affected by the opposing current.

In Figures 4 and 5, three dimensionless amplitudes are considered to examine the variation of the hydroelastic waves by increasing from to for different current. We find that for both a following current and an opposing current a larger increases the hydroelastic waves profile. In the following current case, we can observe from Figure 4 that a larger current intensity reduces the hydroelastic wave profile very slightly, while for the opposing current case, we note from Figure 5 that amplitude of the hydroelastic waves increases obviously with the increment of the current intensity .

Finally, we consider the effects of Young’s modulus of the floating plate on the hydroelastic waves by increasing from to for different current. It is seen from Figures 6 and 7 that for both a following current and an opposing current a larger decays the waves’ amplitude. It is because the plate becomes stiffer so that a large reactive force can be generated to resist the deformation of the hydroelastic waves. Particularly in the following current case, Figure 6 indicates that for a larger the wave amplitude reduces very slightly, while for the opposing current case, we observe from Figure 7 that amplitude of the hydroelastic wave increases very obviously with the increment of current intensity .

5. Conclusions

We extend the study by Wang and Lu (2013) [5] to the nonlinear hydroelastic waves on an underlying uniform current, which describes wave-current interaction in the real ocean. In the frame of the HAM, we construct the zeroth-order deformation equations, in which we choose the auxiliary linear operator containing only the derivatives of in nonlinear operator and containing only the derivatives of in nonlinear operator . In particular, we choose the corresponding linear solutions as the initial estimations of the dispersion relation, the velocity potential , and the hydroelastic waves profile , respectively. The convergent solutions are analytically achieved by solving the th-order deformation equations, and numerical results indicate the convergence and accuracy of our HAM-based analytical expressions by choosing the optimal .

The influences of important physical parameters on the nonlinear hydroelastic waves in the underlying uniform current are discussed in detail. It is found that in the following current case a larger reduces the amplitude of the hydroelastic waves very slightly but increases it obviously in an opposing current. The hydroelastic waves become steeper for both following and opposing current as the incident wave amplitudes increase, but Young’s modulus has the opposite effect on the hydroelastic waves. Further, the hydroelastic wave is more greatly affected by the opposing current. Particularly, the hydroelastic wave is more greatly affected by the opposing current than by the following current.

Appendix

Detailed Derivation of , , , and in (27) and (28)

Let We have a Maclaurin series for any as follows: For , it follows from (A.1) and (A.2) that where Thus, for , we have where

Let

Substituting series expansions (A.1), (A.5), and (A.7) into conditions (18) and (19) and then equating the like-power of , we acquire two linear boundary conditions (27) and (28), respectively. Then the explicit expressions for , , , and in these two conditions are derived by where

Disclosure

The authors themselves programed the symbolic computation software named Mathematica 8.0 independently to gain the approximate analytical solutions of the PDEs considered here.

Conflicts of Interest

There are no conflicts of interest regarding this manuscript.

Acknowledgments

This research was financially supported by the National Natural Science Foundation of China under Grant no. 51374134 and Project of Shandong Province Higher Educational Science and Technology Program no. J14LJ53. Qingdao Postdoctoral Applied Research Project no. 020022034 is also acknowledged.