1. Introduction
Wave generations in the laboratory, either wave channels or plane basins, are important related to performing experiments in studying wave characteristics or wave and structure interactions. Recent literature related to various types of wave generations, i.e., plunger or wedge type wavemakers, or wave characteristics regarding deep or intermediate waters, can be found in [
1,
2,
3]. Wavemaker theories are associated with wave generations; however, they are mostly steady and periodic solutions as far as analytical solutions are concerned. Transient wave generations concentrate on the entire wave development or decaying processes, but analytical solutions due to the inherited mathematical complexities are not so much available [
4]. For fully nonlinear water wave problems, inevitably, numerical methods are eventually good approaches to solve the problems. Jang and Sung [
5] used the smoothed particle hydrodynamics (SPH) method to simulate the piston-type wavemaker problem. In the numerical scheme, the fluid phase was discretized into a distribution of small fluid particles, whose dynamics were influenced by neighboring particles that lay within a support domain. For a given particle, the effect of the surrounding fluid was determined by a weighted sum over all the remaining particles in the system. Wei and Kirby [
6] on the other hand used an extended Boussinesq equation to govern the water wave motions, and a numerical code was given to simulate the time-dependent problems. De Padova et al. [
7] utilized the SPH numerical method to simulate the wave domain, while the buoyant jet was coexisting in the fluid domain. The characteristics of the buoyant jets affecting by the wave motion were investigated. Ozbulut et al. [
8] utilized the Navier–Stokes and continuity equations for governing the fluid motion through the Weakly Compressible SPH (WCSPH) approach to simulate wave generations in the wave tank. Gholamipoor and Ghiasi [
9] developed a numerical wave tank (NWT) by employing the local radial point interpolation collocation method, mixed Eulerian–Lagrangian approach, and the fourth order Runge–Kutta method. The potential theory was used for the mathematical description of the wave propagation problem. Vivanco et al. [
10] used the SPH simulation and introduced a framework for wave generation through the action of the pedal-wavemaking method, which generates long waves without effects of evanescent waves. The other category was the wide use of the open-source OpenFOAM CFD libraries, together with methods that were commonly used to deal with wave generation and absorption [
11,
12,
13].
For analytical solutions to unsteady wave generation problems, Lee et al. [
14] proposed a transient wavemaker theory. The time function was resolved by using the Laplace transform, the x function was expressed by the Fourier cosine function, and the z function was obtained by a direct integration method. The solution could be obtained without any difficulty, but the inversed Laplace transform was complicated and could involve numerical calculations for some mathematical functions. Reviews of some hydrodynamic characteristics of transient waves were given in references therein. Lee [
15] used the same solution methodology to solve for a piston-type wavemaker problem. Comparisons of unsteady wave characteristics generated by flap-type and piston-type wavemaker were discussed. Joo et al. [
16] also proposed an analytical solution to the initial and boundary-value problem of the wavemaker problems. A boundary condition on the wavemaker was presented. The unique part of the solution was a particular solution satisfying the time-dependent wavemaker boundary condition. Chang et al. [
17] solved the piston-type wavemaker problem using a similar solution approach, and Lin et al. [
18] solved the wavemaker problems generating irregular waves. Note that in this approach it was necessary to have a particular solution to satisfy the boundary condition that was in some cases difficult to derive.
The phenomena of undular bores at water surface also have transient characteristics while propagating along the water channel. The associated analytical solutions could possibly be applied for the transient wavemaker problems. Berry [
19] applied the quantum and Hawking effect analogies to describe analytically a model for undular tidal bore profile. EI et al. [
20] presented an asymptotic analytical description for undular bores, in which the Whitham modulation theory and the Su–Gardner system were used together with the Boussinesq governing equation. Landrini and Tyvand [
21] analytically and numerically studied the free surface flow induced by impulsive bottom movements. Marchant [
22] adopted the modified Korteweg–de Vries equation and presented analytical solutions for undular bores. Bestehorn and Tyvand [
23] solved nonlinear, two-dimensional Laplace equation and studied interactions between two undular bores in a long rectangular channel. Ali and Kalisch [
24] studied the dispersive undular bore model and concluded that there was no energy loss in the dispersive model. Chen et al. [
25] presented a mathematical model to simulate undular bores generated by blocking steady water currents in a flume, and an analytical solution was presented. Tong and Liu [
26] presented an analytical solution for studying transient water wave-induced responses inside an unsaturated poroelastic seabed of finite thickness. The transient waves, including the linear periodic wave, a solitary wave, and a bore, were considered.
In this study, a new analytical solution for the transient wavemaker problems in a wave channel is presented. The proposed method follows Chen et al. [
25], which combined the applicability approach of Lee et al. [
14] and Joo et al. [
16]. Our analytical solution is applied to simulate wave evolutions generated by the wavemaker motions. The generated steady state wave forms are compared with the steady state analytical solution given by Dean and Dalrymple [
27], and characteristics of evanescent waves are identified. Finally, we demonstrate capabilities of the present analytical solution by showing the generated wave forms along the wave channel at an integer number of periods of the wavemaker motion, and transient characteristics of wave generation regarding leading wave lengths and wave heights are investigated.
2. Problem Description and Solution
The problem of wave generation by a wavemaker in a wave channel was sketched, as shown in
Figure 1. The length of the channel was
, the constant water depth was
. A piston type wavemaker was specified at the left-hand side, and the end of the channel was a vertical wall. With an initially still water surface, and given that the wavemaker started with periodic motion, it was expected to have wave motions in front of the wavemaker, which also propagate toward the end of the channel. A Cartesian coordinate system with the origin located at the still water level was chosen, as indicated in
Figure 1. The positive
x axis pointed to the right and the positive
z axis pointed upward.
In this study, the wave motions in the channel are described by using the linear potential wave theory. The relation of the flow velocity
and the potential function
is given by
where
is the Laplacian operator. The boundary-value problem can be expressed as [
27].
The governing equation is the Laplace equation:
The free surface boundary condition and the bottom condition can be written, respectively, as
where
is the gravity constant. The wavemaker boundary condition at
x = 0 is
where
is the displacement function of the wavemaker and the upper dot indicates time derivative. At the right-hand end vertical wall, the boundary condition is
The initial conditions considered in this paper were the still water in the channel and no disturbances at the free surface. Equations (2)–(6), together with the initial conditions, represent the initial and boundary-value problem shown in
Figure 1. In the following, an analytical solution to the problem will be described.
The potential function
was a function of spatial coordinates
and time
. The first step of the present solution was to solve the
function. Equations (2)–(6) were reformulated using the Fourier cosine transform. The
function was expressed in terms of the finite Fourier cosine series, and the potential function could be expressed as [
28]
where
and
,
are defined, respectively, as
The use of the cosine functions can facilitate implements of boundary conditions at two ends of the channel, as can be observed later in the solution formulation.
The solution step then turned to solve the
function included in
and
. In the following, only the solution of
is given, as the solution to
is only a special case of
. The governing equation, the free surface, and bottom conditions, Equations (2)–(4), were rewritten, following the Fourier cosine transform of Equation (9), as
Equation (10) was a non-homogeneous second-order differential equation, and the solution can be obtained by using the method of variation of parameters [
28] and written as [
28]
in which
and
are the homogeneous solution and the particular solution, respectively. The solution to Equation (13) was
where
and
could be determined by utilizing the free surface and bottom boundary conditions, Equations (11) and (12).
The substitutions of Equations (13)–(15) into the free surface and bottom boundary conditions yield
where
. Equation (16) is a second-order differential equation in time solving for
, and Equation (17) is a relation between the coefficients
and
. The solution of Equation 16 can be written as
where
and
are the integration constants, whereas
The unknown constants and shown in Equation (18) could be determined by applying the initial conditions.
Consider that, initially in the wave channel, water is still and there are no disturbances at the free surface. The initial conditions can be written as
Equations (20) and (21) should be rewritten according to the Fourier cosine transform, Equation (9), and provide two equations solving for the two constants. Thus, we obtained undetermined constants, shown in Equation (18):
So far, the potential function
is completely determined. Similar to obtaining
, we also obtained
as
The initial boundary value problem described by Equations (2)–(6), (20), and (21) was solved.
Thus, the initial boundary value problem was solved. The wave potential function could be expressed explicitly as
For a piston type wavemaker, the displacement function of the wavemaker can be expressed as
where
is the stroke of the wavemaker motion,
, and
is the period of the wavemaker motion. Following the general analytical solution, Equation (25), the associated potential function for the generated wave can be written as
The associated free surface elevation can be calculated from Equation (27) by using the Bernoulli’s equation and can be expressed as
Note that, for the piston type wavemaker problem, an analytical transient solution was presented by Lee [
15] using a different solution methodology. In solving the initial and boundary-value problems, the Laplace transform was used to transform the time function, then solved the corresponding spatial functions, along with the boundary conditions and the final inverse Laplace transform to obtain the analytical solution. While we tried to utilize Lee [
15]’s approach to solve other wavemaker problems, we found that the inverse Laplace transform may not be able to work; therefore, we devised new and different solution methods to solve the problem. After we obtained the present analytical solution, we thought that there could be only one solution for the same initial and boundary-value problem, and we proceeded to give proof.
We show that the real part expression of the free surface elevation, Equation (28), in the present solution can be mathematically manipulated to show be identical with Equation (18), shown in Lee [
15].
The exponential function shown in Equation (28) can be expressed by
With the substitution of Equation (29) into Equation (28), we obtain
The real part of the free surface elevation, Equation (30), could be expressed as
Next, we adopted the variable definitions used in Lee [
15]. A cross reference of symbols between Lee [
15] and present formulation is shown in
Table 1.
After some mathematical rearrangements, Equation (31) can be rewritten as
Equation (32) is precisely identical to the expression shown in Lee [
15], which confirms the present analytical solution. We like to emphasize that, although the free surface elevation obtained in this study, Equation (28), was proved identical with the expression in Lee [
15], the derivation to obtain the analytical solution is entirely different from Lee [
15]. Furthermore, without having the inverse Laplace transform, the present analytical method could be applied to solve other problems more easily.
In the next section, convergence of the series expression shown in the present analytical solution is be discussed first. Then, the wave forms generated by the transient solution are shown to match the steady wave theory in the fully developed state. Additionally, the transient wave length is compared with the steady wave length to show transient characteristics of wave length variation. In particular, the first wave length in front of the wavemaker is shown to include both propagating and evanescent waves. Secondly, the present analytical solution is compared with numerical solution and experimental results to validate wave form distribution along the wave flume. Finally, transient characteristics of the generated leading waves are investigated.
3. Results and Discussion
The present analytical solution for the wavemaker problem contains the cosine series expression. To calculate the generated free surface elevation, a finite number of terms should be determined in the computation. A result of convergence test is shown in
Figure 2. The vertical axis is the dimensionless surface elevation, the subscript
indicates the number of terms used of the series expression, and
is the amplitude of the generated steady wave form. The conditions used in the calculation were
= 0.2 m,
= 0.035 m,
= 0.98 s. The relative water depth,
= 1.06, falls in the range of intermediate water depth. The generated steady fully developed wave height and wave length are 3.67 m and 1.18 m, respectively, as calculated from the steady wave theory [
27].
Figure 2 shows convergence of the series solution; where an error criterion 0.1% is adopted, the number of terms
= 178 should be used, as indicated by the dash vertical line in the figure. The relative error criterion is defined as
Time evolution of the generated surface waves at a fixed location of the wave channel is shown in
Figure 3, where
= 0.3 m,
= 0.052 m,
= 0.82 s, and
= 1.88. The observatory location is one wave length away from the wavemaker. The gain function,
, calculated from the steady wavemaker theory, is also plotted to show the envelopes of the developed wave amplitudes with time as the upper and lower dash lines.
Figure 3 shows that the fully developed wave heights follow the envelopes of the steady state solution. Furthermore, the wave forms take four waves to reach the fully developed state, and the developing rates of the wave form are illustrated.
The distribution of surface wave along the channel at a fixed time after waves are generated is plotted in
Figure 4. The horizontal axis is a dimensionless distance away from the wavemaker, and
is the steady wave length. The vertical axis is the dimensionless surface elevation. The conditions used are
= 0.3 m,
= 0.037 m, and
= 1.08 s, and 15 wave periods are calculated. The steady wave form calculated using the steady solution is plotted in phase with the present theory for comparison.
Figure 4 shows that the generated fifth wave almost has the same wave length as the steady result. The leading first wave has approximately 4.5 times the steady wave length, and the following wave lengths decrease gradually and reduce to the steady result.
In the wavemaker theory, motions of the wavemaker generate propagating waves and accompanied evanescent waves that decay while moving away from the wavemaker. In front of the wavemaker, the fully developed wave form calculated by the present analytical solution, the steady theory, and the propagating waves are plotted in
Figure 5.
Figure 5 shows that after the waves are fully developed, the present time domain solution (line) precisely matches the wave length and amplitude as the steady theory (circle). On the other hand, the difference between steady wave form and the propagating wave form (red dash line) indicates the evanescent wave in front of the wavemaker.
The generated surface waves by the wavemaker at 8.0 s, calculated using Equation (31), was validated with the numerical boundary element solution, as shown in
Figure 6. The result is cited from Figure 13 in Lee [
15]. The conditions used were water depth 3 m, length of the channel 30 m, stroke, and period of the wavemaker at 1.0 m and 2.0 s, respectively. The comparisons of the unsteady wave profiles were in reasonable agreements. The present analytical solution and the numerical solution matched very well for the first wave length. For the rest of the wave forms, the numerical solutions had slightly higher wave peaks and shallower wave troughs, with maximum 5% differences. The present analytical solution was further used to simulate experimental results presented by Gao [
29]. In the experiments, the water depth was 0.28 m and the stroke and period of the wavemaker motion were 0.02 m and 1 s, respectively. The length of the wave channel was 8.75 m. Variation of the free surface elevation with time at a location 0.55 m in front of the wavemaker was used for comparison. Calculated results using present analytical solution (line) and experimental results (discrete star) copied from Gao [
29] are plotted in
Figure 7. The comparison shows favorable agreements. Note that the wave amplitude calculated by using the steady wavemaker theory was 0.0125 m. At wave peaks, the two results match the steady theory; however, at wave troughs, the present analytical solution matches the steady theory, but the experimental results are shallower, and the differences are approximately 12%.
The present analytical solution can be used to simulate transient development of waves generated by the wavemaker in a channel. Time variations of surface waves in front of the wavemaker along the wave channel at one-period time, two-periods, six-periods, and 10-periods of the wavemaker motion are shown in
Figure 8. At time one-period, a developing wave form is induced, where a smaller wave peak is ahead of a lower trough propagating away from the wavemaker. At the two-period time, two wave forms are generated along the channel. A smaller wave is propagating ahead of a bigger wave, while waves are propagating away from the wavemaker. The wave forms are gradually developing, and it is expected that as more waves are continuously generated, the wave forms will reach fully developed conditions. At the six-period time, there are six unsteady waves in front of the wavemaker; among them there are four developed waves behind the two smaller leading waves. As for the 10-period time, there are four fully developed waves behind the four leading transient waves. We can notice that between fully developed waves and the leading waves, there is a highest wave. The simulated results indicated that there is a relatively bigger wave between developing waves and fully developed steady waves, which is about 15% bigger than the steady wave heights. As more time elapsed, there will be more fully developed waves generated following the leading waves.
Using the present analytical solution, we calculated the generated waves to investigate unsteady characteristics of the wave heights and wave lengths of the leading waves. The conditions used are water depth 0.35 m, the stroke of the wavemaker motion 0.019 m, and the wavemaker periods 1.14 s, 1.52 s, 1.95 s, and 2.56 s.
Figure 9 shows variations of dimensionless wave height versus dimensionless time, where
is the steady wave height and
is the wave height of the leading waves. Overall, the wave heights of the leading waves decrease as they propagate away from the wavemaker. For shorter waves, the wave heights decrease rapidly in the beginning five periods of time; on the other hand, for longer wave periods, the wave heights decrease rapidly in the beginning 15 periods of time. For 25 to 35 periods of time, the wave heights decrease and become gradual, showing a tendency of reaching constants. The shorter waves reach 0.1 and the longer waves reach 0.2.
Figure 10 shows variations of dimensionless wave lengths of the leading waves versus dimensionless time, where
is the steady wave length and
is the wave length of the leading waves. The wave length of the leading wave is defined as the distance starting from the raised water front to the first following zero-up crossing. In general, the wave lengths of the leading waves increase with time as they propagate away from the wavemaker. The calculated results indicate that for longer waves, tendency of wave length increase will reach 3.5; on the other hand, for shorter waves, the tendency will reach 7.