Abstract and Applied AnalysisVolume 2012 (2012), Article ID 214897, 12 pagesdoi:10.1155/2012/214897

Research Article

Nonlinear Numerical Investigation on Higher Harmonics at Lee Side of a Submerged Bar

D. Ning,1,2 X. Zhuo,¹ L. Chen,¹ and B. Teng¹

¹State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China²Waterway and Sediment Engineering Key Laboratory of Ministry of Transport, Nanjing Hydraulic Research Insititute, Nanjing 210029, China

Received 1 January 2012; Accepted 25 January 2012

Academic Editor: Muhammad Aslam Noor

Copyright © 2012 D. Ning et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The decomposition of a monochromatic wave over a submerged object is investigated numerically in a flume, based on a fully nonlinear HOBEM (higher-order boundary element method) model. Bound and free higher-harmonic waves propagating downstream the structure are discriminated by means of a two-point method. The developed numerical model is verified very well by comparison with the available data. Further numerical experiments are carried out to study the relations between free higher harmonics and wave nonlinearity. It is found that the nth-harmonic wave amplitude is growing proportional to the nth power of the incoming wave amplitude for weakly nonlinear wave condition, but higher-harmonic free wave amplitudes tend to a constant value for strong nonlinear wave condition.

1. Introduction

Submerged bars are common marine structures widely used in coastal engineering for different applications. Besides wave diffraction, higher-harmonic waves are generated due to wave shoaling effect, when the incoming waves propagate over a submerged structure. The generation of the superharmonic waves changes the swell spectrum because a significant part of the incoming wave energy may be transferred to higher frequencies. The free higher-harmonic waves are prominent at the lee side of the submerged bars and generally show strong nonlinear effect, whose wavelengths are much shorter than that of incident wave, wave amplitudes grow with increasing incoming wave amplitude, even to comparable to the incoming wave amplitude. Such free higher harmonics may debase sailing conditions and do harm to coastal structures. Thus, it is quite necessary to take account of their nonlinear characteristics.

Numerous researching works have been performed on the problem of submerged bar, experimentally and numerically. For example, Beji and Battjes [1], Luth et al. [2] made physical experiments of wave transformation over a submerged bar and accordingly stated the generation mechanism of higher-frequency wave energy. Jeng et al. [3] experimentally studied the wave propagation over a submerged breakwater in front of a vertical wall, and the corresponding wave surfaces are recorded and analyses. Cho et al. [4] reported their experimental study of wave reflection over submerged breakwaters in tandem, and concluded that reflection coefficients of impermeable submerged bars are bigger than those of permeable ones and trapezoid bars are advised to use in the view of applicability. Rambabu and Mani [5] numerically studied the effects of incident wave factors, water depth over submerged bar, and section width of submerged bar on the transmission coefficient, and they accordingly obtained the optimal section width of submerged bar. Madsen et al. [6], Madsen and Sφrensen [7], Engsig-Karup et al. [8], Fang and Zou [9] adopted various Boussinesq-type equations to model wave propagation and transformation over submerged bars or bottom topography. However, the analysis of the higher harmonic waves scatting by submerged bar is relatively scare. Grue [10] experimentally researched the decomposition of higher harmonic waves due to submerged rectangular step and horizontal cylinder in deep water and got the conclusion that superharmonic wave amplitudes grow up to a saturation values with the increasing of incoming wave amplitude. Brossard and Chagdali [11], Liu et al. [12] experimentally and numerically studied the harmonic generation by waves over a submerged plate, respectively.

Among various numerical studies, the Boussinesq-type equation has been commonly used in the analyses of nonlinear wave evolutions in shallow water. However, due to assumptions of both weak wave nonlinearity and weak wave dispersivity, the Boussinesq equation may not be valid for predictions of free higher harmonics in deep-water region and of highly nonlinear interaction between bound and free waves in the shallow-water region. On the other hand, applications of the higher-order boundary element method (HOBEM) to obtain a fully nonlinear numerical solution have become extremely successful in the past two decades. This method is achieved by directly solving the boundary integral equation with fully nonlinear free surface boundary conditions. The essential advantages include reducing the original -dimensional partial differential equation to an ()-dimensional boundary integral equation, which allows the numbers of introduced unknowns to be reduced significantly, and no limitation to water depth and wave height if there do not exist wave breaking.

In this paper, the free higher harmonics due to a monochromatic wave propagation over a submerged bar are addressed by using the authors’ developed fully nonlinear HOBEM [13]. This paper is arranged as follows. A brief description of the numerical wave tank model is given in Section 2, in which a two-point decomposition technique of higher harmonics from a periodic wave train is also introduced. In Section 3, the computed wave profiles at various locations and the higher harmonic wave amplitudes are compared to experimental observations in order to confirm the validity of the numerical model. The characteristics of free higher harmonics for weakly and strong nonlinear wave conditions are further addressed. Finally, the conclusions are presented in Section 4.

2. Mathematical Formulation

2.1. Governing Equation and Boundary Conditions

For a monochromic wave propagating over a submerged bar as shown in Figure 1, the numerical analysis are performed using a numerical wave tank model previously developed by the authors [13]. In the model, the incident wave is generated through a nonreflective wave generator which combines vertically distributed wave-making sources on , introduced by Brorsen and Larsen [14]. Two damping layers are located at both ends of the computational domain to dissipate the outgoing waves.

Figure 1: Sketch of the numerical wave tank model.

Fluid is assumed to be ideal, so that a velocity potential exists and satisfies the following Poisson equation inside the fluid domain : and the boundary conditions are expressed as in which is corresponding to the instantaneous free surface, the tank bottom, the surface of submerged body, the tank lateral surfaces. is the pulsating volume flux density of the internal source distribution and is prescribed second-order Stokes wave velocity. The damping function is adopted as the following: where is a reference value specifying the at-rest position of the fluid particle, α is damping coefficient, is length of damping zone. In this paper, , is one and half times incident wave length.

2.2. HOBEM

By using the second Green’s theorem, the prescribed boundary value problem can be transformed the following boundary integral equation: where and are source and field points, and is the solid angle. For the Green function , an image Green function is utilized so that the integration surface only includes the incident boundary and the free surface boundary and submerged bar . The Green function can be obtained by the superposition of the image of the Rankine source about the sea bed and the infinite images about the two lateral walls and written as (see [15]) where

is the tank width, and is the static water depth.

Then the boundary surface is discretized with a number of elements. The geometry of each element is represented by the shape functions, thus the entire curved boundary can be approximated by a number of higher-order elements. Within the boundary elements, physical variables are also interpolated by the same shape functions, that is, the elements are isoparametric. In the integration process, the solid angle, the single-layer and double-layer integration are directly resolved (see [16]). Since the discretized integral equation is always variant in time, all the boundary surfaces are remeshed and updated at each time step using the mixed Eulerian-Lagrangian scheme and 4th-order Runga-Kutta approach. Once (2.4) is solved, we can obtain the time history of surface evaluation at any position as described by Ning and Teng [17].

2.3. Separation of Free Waves from Phase-Locked Waves

After waves pass the submerged bar, higher harmonic waves generated by nonlinear wave-wave interactions in the shallow water over the bar will leave the bar leeward as free waves. So the surface elevation at at the lee side of the bar can be written as where is the amplitudes of the free transmitted waves with frequencies of integer times of the incident wave frequency; are the amplitudes of the nth-harmonic phase-locked waves, are the initial phase angles of the fundamental waves, and are the nth-harmonic free waves; k and are the wave number of the fundamental waves and the nth-harmonic free waves, and they satisfy the following dispersion relations: respectively. The fundamental wave amplitude, as well as the higher-order free and locked wave amplitudes, is obtained from the time histories of the surface elevation. The Fourier transform is introduced as follows: where and are the corresponding real and imaginary components, respectively.

Substituting (2.7) into (2.9) and making use of orthogonality of trigonometric function, we may obtain

Then, we can obtain the amplitudes of the nth-harmonic locked waves and free waves after applying two fixed points’ ( apart) surface elevations into (2.11) through an easy iterative solving process.

3. Numerical Results and Discussions

As a validation of the present model, a case for monochromic wave propagating a submerged trapezoid bar as shown in Figure 2 was considered firstly. The same problem was ever experimentally studied by Beji and Battjes [1] and Luth et al. [2]. In the present numerical simulation, some parameters for wave period , wave amplitude , and water depth are used. The corresponding computational domain is taken as ( denotes wave length), time step defined as T/80, and meshed with cells in , and directions after numerical convergent tests, in which the first and last are used as the damping layers.

Figure 2: Definition sketch of wave propagating over a trapezoid bar.

Figure 3 presents the distribution of wave elevations along the plane of symmetry at and . From the figure, it can be seen that there is a good agreement between these two times, which states that the numerical convergence has been obtained in the proposed model. It is also found that the wave profiles at the weather side of the submerge object are still regular and periodic, but they are quite skewed at the lee side. It means that a number of wave components are introduced the incident flow at the lee side.

Figure 3: Wave profile along the plane of symmetry at and 20T.

The comparisons of wave elevation time series at specific points marked in Figure 2 among the present results, experimental data [1], and the solutions of the extended Boussinesq model of Nwogu [9] are given in Figure 4. From the figures, it can be seen that there are good agreements for the three methods at points , and . But for the points at the lee side such as , the Boussinesq model cannot give accurate description because the wave nonlinearity is stronger behind the shoal. Therefore, the fully nonlinear model, such as the present scheme, is advised to analyze such strong nonlinear wave problem.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: Comparisons of free surface displacement among the experiment, Boussinesq model, and present method.

The results in Figures 2 and 3 show that much higher-harmonic waves are involved at the lee side of the submerged bar. So it is necessary to analyze these higher harmonics. Equation (2.7) is used to separate various waves at point by truncating to 7. The downstream neighboring point is also used here for solving (2.11). Figure 5 gives the comparison of free surface at point between the numerical solution and the fitting curve for weakly nonlinear condition (, , ) and strong nonlinear condition (, ), respectively. The other parameters are the same with the above. It proves the validity of the present two-point separation method from good agreement. Figure 5(a) shows a regular wave profile similar to incident wave. But Figure 5(b) presents a strong deformed wave profile, even secondary wave peaks can be apparently found.

(a) Weakly nonlinear condition (, )

(b) Strong nonlinear condition (, )

Figure 5: Comparison of free surface at between the numerical solution and the fitting curve.

Once the parameters , and initial phase angles are solved, the various order wave records can be obtained as shown in Figures 6 and 7. From the figure, it can be seen that the periods of second and third harmonic waves are corresponding to 1/2 and 1/3 times wave period, and phase-locked waves and free waves keep same periods but different phases. Locked waves are released to the corresponding free wave in a large degree after the wave propagating over the submerged bar, especially for strong nonlinear condition shown in Figure 7. With the increasing of wave nonlinearity, the higher harmonic free waves are found to be more prominent at the obstacle’s lee side. The second-harmonic free wave amplitude is comparable to the incoming wave amplitude as shown in Figure 7.

(a) 1st-harmonic free wave

(b) 2nd-harmonic waves

(c) 3rd-harmonic waves

Figure 6: Time histories of various order harmonic waves for weakly nonlinear condition (, ).

(a) 1st-harmonic free wave

(b) 2nd -harmonic waves

(c) 3rd-harmonic waves

Figure 7: Time histories of various order harmonic waves for strong nonlinear condition (, ).

From above results, it shows that the incoming wave is mainly broken up into superharmonic free waves at the lee side of submerged structure due to nonlinear free surface effect. Thus, the higher harmonic free waves are especially researched in the following section. Figure 8 gives the distribution of first-, second-, and third-harmonic free wave amplitudes with the increasing of incident wave amplitudes for static water depth . In the figure, the real lines are the linear, quadratic, and cubic fitting curves obtained from the computed data. It means that the nth-harmonic wave amplitude is growing proportional to the nth power of the incoming wave amplitude for weakly nonlinear wave condition. Figure 9 shows the distribution of first-, second-, and third-harmonic free wave amplitudes with the increasing of incident wave amplitudes for static water depth . With the water depth decreased, the wave nonlinearity is increased. The relation between harmonic free wave amplitudes and incoming wave amplitudes does not obey the rules shown in Figure 8 again. It is observed that the generation of higher harmonic free waves is more powerful than that in Figure 8. The various harmonic free wave amplitudes tend to a constant value as the incoming wave amplitude . Surprisingly, the portion of the second-harmonic free wave amplitude is even larger than that of first-harmonic wave amplitude since .

(a) Fundamental free wave

(b) 2nd- and 3rd-harmonic free waves

Figure 8: First-, second-, and third-harmonic free wave amplitudes at the lee side of submerged bar for .

(a) Fundamental wave

(b) 2nd- and 3rd-harmonic free waves

Figure 9: First-, second-, and third-harmonic free wave amplitudes at the lee side of submerged bar for .

To further testify the conclusion obtained in Figure 9, another case for monochromic wave propagating a submerged rectangular step is considered. The rectangular shelf has a cross-section 0.5 m long and 0.41 m high. The distance between the uppermost point of the geometry and the mean free surface is chosen as 0.0375 m and 0.05 m, respectively. The incident wave frequency is defined as 5.969 and 6.597 rad/s, respectively. Figure 10 shows the distribution of higher free harmonic wave amplitudes with different incident wave amplitudes, in which the experimental data from Grue [10] are also given and the longitudinal coordinates are dimensionless divided by the incident wave amplitude . From the figure, it can be seen that there are better agreements between numerical solutions and experimental data. The similar conclusions with those in Figure 9 are also drawn.

(a) (, )

(b) (, )

Figure 10: Distribution of higher free harmonic wave amplitudes with incident wave amplitude.

4. Conclusions

The phenomenon of wave propagation over a submerged bar is examined by a powerful numerical wave tank method. The transmitted waves are separated into the free waves and the locked waves with different harmonic frequencies by a two-point method. Numerical results show that, at the leeside of the bar, the portions of the free waves are much larger than those of the locked waves with the same frequency. For weakly nonlinear wave conditions, the nth-harmonic free wave amplitude is growing proportional to the th power of the incoming wave amplitude and the fundamental free wave plays an important role. However, for strong nonlinear wave condition, the fundamental free wave amplitude is decreased and second- and third-harmonic free wave amplitudes are increased with the increasing of incoming wave amplitude, and all tend to a constant value, that is, a kind of saturation statues. Even the phenomena of higher-harmonic free wave amplitude exceeds the lower-harmonic free wave amplitude occur.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Basic Research Program of China (973 Program, Grant no. 2011CB013703), the National Natural Science Foundation of China (Grant nos. 51179028, 11072052, and 50921001), and the Open Fund of Waterway and Sediment Engineering Key Laboratory of Ministry of Transport of China.

References

S. Beji and J. A. Battjes, “Experimental investigation of wave propagation over a bar,” Coastal Engineering, vol. 19, no. 1-2, pp. 151–162, 1993.

H. R. Luth, R. Klopman, and N. Kitou, “Kinematics of waves breaking partially on an offshore bar; LDV measurements of waves with and without a net onshore current,” Tech. Rep. H-1573, Delft Hydraulics, 1994.

D. S. Jeng, C. Schacht, and C. Lemckert, “Experimental study on ocean waves propagating over a submerged breakwater in front of a vertical seawall,” Ocean Engineering, vol. 32, no. 17-18, pp. 2231–2240, 2005.

Y. S. Cho, J. I. Lee, and Y. T. Kim, “Experimental study of strong reflection of regular water waves over submerged breakwaters in tandem,” Ocean Engineering, vol. 31, no. 10, pp. 1325–1335, 2004.

C. Rambabu and J. S. Mani, “Numerical prediction of performance of submerged breakwaters,” Ocean Engineering, vol. 32, no. 10, pp. 1235–1246, 2005.

P. A. Madsen, R. Murray, and O. R. Sφrensen, “A new form of Boussinesq equaitons with improved linear dispersion characteristics,” Coastal Engineering, vol. 15, no. 4, pp. 371–388, 1991.

P. A. Madsen and O. R. Sφrensen, “Bound waves and triad interactions in shallow water,” Ocean Engineering, vol. 20, no. 4, pp. 359–388, 1993.

A. P. Engsig-Karup, J. S. Hesthaven, H. B. Bingham, and T. Warburton, “DG-FEM solution for nonlinear wave-structure interaction using Boussinesq-type equations,” Coastal Engineering, vol. 55, no. 3, pp. 197–208, 2008.

K. Z. Fang and Z. L. Zou, “Boussinesq-type equations for nonlinear evolution of wave trains,” Wave Motion, vol. 47, no. 1, pp. 12–32, 2010.

J. Grue, “Nonlinear water waves at a submerged obstacle or bottom topography,” Journal of Fluid Mechanics, vol. 244, pp. 455–476, 1992.

J. Brossard and M. Chagdali, “Experimental investigation of the harmonic generation by waves over a submerged plate,” Coastal Engineering, vol. 42, no. 4, pp. 277–290, 2001.

C. R. Liu, Z. H. Huang, and S. K. Tan, “Nonlinear scattering of non-breaking waves by a submerged horizontal plate: experiments and simulations,” Ocean Engineering, vol. 36, no. 17-18, pp. 1332–1345, 2009.

D. Z. Ning, B. Teng, R. Eatock Taylor, and J. Zang, “Nonlinear numerical simulation of regular and focused waves in an infinite water depth,” Ocean Engineering, vol. 35, no. 8-9, pp. 887–899, 2008.

M. Brorsen and J. Larsen, “Source generation of nonlinear gravity waves with the boundary integral equation method,” Coastal Engineering, vol. 11, no. 2, pp. 93–113, 1987.

J. N. Newman, “The approximation of free-surface Green functions,” in Wave Asymptotics, P. A. Martin and G. R. Whickham, Eds., pp. 107–142, Cambridge University Press, Cambridge, UK, 1992.

B. Teng, Y. Gou, and D. Z. Ning, “A higher order BEM for wave-current action on structures—direct computation of free-term coefficient and CPV integrals,” China Ocean Engineering, vol. 20, no. 3, pp. 395–410, 2006.

D. Z. Ning and B. Teng, “Numerical simulation of fully nonlinear irregular wave tank in three dimension,” International Journal for Numerical Methods in Fluids, vol. 53, no. 12, pp. 1847–1862, 2007.