Abstract
The scattering of incident surface water waves due to small bottom undulation on the porous bed of a laterally unbounded ocean in the presence of surface tension at the free surface is investigated within the framework of two-dimensional linearized water wave theory. Perturbation analysis in conjunction with the Fourier transform technique is employed to derive the first-order reflection and transmission coefficients in terms of integrals involving the shape function representing the bottom undulation. One special type of bottom topography is considered as an example and the related coefficients are determined in detail. These coefficients are presented in graphical forms. The theoretical observations are validated computationally. The results for the problem involving scattering of water waves by bottom deformations on an impermeable ocean bed are obtained as a particular case.
1. Introduction
The interaction problems of surface waves with a preexisting (fixed) pattern of undulation on an otherwise flat bed are important for their possible applications in the areas of coastal and marine engineering, and as such these are being studied for a long time. The problem of reflection of surface waves by patches of bottom undulations has received an increasing amount of attention as its mechanism is important in the development of shore-parallel bars.
The problems of water wave diffraction by undulating bottom topography of a sea bed have been considered by many scientists and several papers have been published in this direction (see Davies and Heathershaw [1], Mandal and Basu [2], Martha and Bora [3], Miles [4], and many others).
All the above works focused only on the wave motion over the impermeable bottom topography of sea bed in absence of surface tension at free surface. Nowadays, due to many interesting applications in the theory of water waves scattering, many researchers have turned their attention to the problems related to porous bed rather than an impermeable one. Mase and Takeba [5], Zhu [6], and Silva et al. [7] considered water-wave reflection and/or transmission problems where a porous medium was assumed to lie on a sea bed of varying quiescent depth. Martha et al. [8] focused on the problem involving wave scattering by small undulation on a porous bottom topography without considering the surface tension at the free surface.
Here, we investigate the scattering problem of surface water waves by small undulation on an porous ocean-bed of finite depth in the presence of surface tension at the free surface. The motion of the fluid below the porous ocean bed is not analyzed here and it is assumed that the fluid motions are such that the resulting boundary condition on the ocean bed as considered here holds good and depends on a known parameter , called porosity parameter, in this analysis. Due to the small parameter (≪1) being present in the representation of the small undulation of the bottom, the perturbation analysis can be employed. By this analysis, the velocity potential, the reflection, and transmission coefficients appearing in the governing boundary value problem (BVP) can be expanded as a power series involving this parameter . Equating the coefficients of identical powers of the parameter from both sides of all equations appearing in the coupled BVP, a number of BVPs can be obtained. However, the BVP up to first order is considered here for its solutions as the BVP of higher order can be solved successively in principle. The solution of the zeroth order BVP is obvious. Using Fourier transform technique, the solution involving the velocity potential for the first order BVP is obtained. By comparing this solution with far-field behaviour of the velocity potential, the first-order reflection and transmission coefficients are derived in terms of integrals involving the shape function representing the bottom undulation. These coefficients are computed numerically for a special form of bottom topography, namely, sinusoidal bottom, and depicted graphically. The results for problem involving scattering of water waves by bottom deformations on an impermeable ocean bed are obtained as a particular case.
2. Problem Formulation
A right-handed rectangular Cartesian coordinate system is considered in which -axis is the position of the undisturbed free surface of the ocean and -axis is positive vertically downward from the undisturbed free surface (see Figure 1). The ocean bed with small undulation, composed of porous material of specific type, is described by , where is a bounded and continuous function, describing the undulation of the ocean bed and as so that the ocean is of uniform finite depth far away from the undulation on either side; the nondimensional number (≪1) is a small parameter giving a measure of the smallness of the undulation. It is also assumed that the fluid is incompressible and inviscid, and the motion is irrotational. The usual assumptions of linear water wave theory and the removal of the harmonic time independence lead to the following governing equations for the time-independent complex-valued velocity potential function : where and , is the angular frequency of the incoming water-wave train, is acceleration due to gravity, is the coefficient of surface tension at the free surface of the ocean, is the density of water, denotes the normal derivative at a point on the bottom, and is the porous effect parameter corresponding to the ocean bed under consideration.
It is assumed that a train of progressive wave is incident to the bottom undulation from the direction of , then it is partially reflected by and partially transmitted over the bottom undulation so that has the far-field behaviour given by where and represent, respectively, the unknown reflection and transmission coefficients which are to be determined here.
3. Methodology and Results
3.1. Perturbation Technique
Assuming to be very small, the bottom boundary condition (3) can be expressed as If there is no undulation at the bottom, the surface wave train propagates without any hindrance and there is total transmission. In view of this along with the approximate form of relation (5), can be expressed in terms of the small parameter as follows: Substituting the expansions (6a)–(6c) into (1)-(2) and (4)-(5), we find, after equating the coefficient of and from both the sides, that the functions and , respectively, satisfy the following BVPs.
BVP-I. The function satisfies
BVP-II. The function satisfies The solution of the BVP-I corresponds to the incident wave over finite depth and can be written as where , the wave number of the normal incident wave, is the unique positive root of the following dispersion relation: Note that, is the negative root of the dispersion relation (11).
3.2. Fourier Transform Technique
The solution of BVP-II is obtained by using Fourier transform technique. To solve this BVP, we now assume to have a small positive imaginary part so that decreases exponentially as . This ensures the existence of Fourier transform of with respect to and defined as follows: with the following inverse transform: Applying Fourier transform defined in relation (12), to (8a), (8b) and (8c) and solving them, we obtain as given by where is the fourier transform of .
Taking the inverse fourier transform, the solution for the first order velocity potential can be written as where Now the first order reflection and transmission coefficients, and , due to normally incident wave, are obtained by letting in the relation (15) and comparing with (9).
To calculate the first order reflection coefficient , we let in relation (15). As , the behaviour of can be obtained by rotating the path of the integral involving the term into a contour in the first quadrant so that we must include the residue term at the pole . The path of the integral involving the term in the relation (15) is rotated into a contour in the fourth quadrant so that the integral involving the term does not contribute as . Then comparing the resultant integral value with (9), we obtain the value of as where and is given in the relation (16) with denoting the derivative of with respect to .
Similarly, to calculate the first order transmission coefficient , we take in relation (15). As , the behaviour of can be obtained by rotating the path of the integral involving the term into a contour in the first quadrant so that we must include the residue term at the pole . The path of the integral involving the term in (15) is rotated into a contour in the fourth quadrant so that the integral involving the term does not contribute as . Then comparing the resultant integral value with (9), we obtain the value of as These resulting integrals and can be evaluated once the bottom profile is known. In the Section 5, we consider a special form for the shape function .
4. Limiting Case
In the absence of the surface tension and porosity (i.e., , ), the reflection and transmission coefficients in relations (17) and (19) agree with the ones obtained by Miles [4], Davies and Heathershaw [1], Martha and Bora [3], and Mandal and Basu [2] (for , impermeable bed).
In the absence of the surface tension, that is, when , the above results coincide exactly with the results of Martha et al. [8].
5. Special Form of the Bottom Profile
In 1982 Davies [9] found that an undulating bed has the ability to reflect incident wave energy which has an important implication in the application of coastal protection and in the case of possible ripple growth if the bed is erodible. Here we consider a special form of the shape function in the form of a patch of sinusoidal ripples as the bottom undulation. The shape function is given by where with as the amplitude of the sinusoidal ripples, the wave number of the sinusoidal ripples with an arbitrary phase angle , and , positive integers. For this case, the reflection coefficient and the transmission coefficient , respectively, are obtained as In the situation in which there is an integer number of ripple wavelengths in the patch such that and , we find and , respectively, as Relation (23) illustrates that for a given number of ripples, the first order reflection coefficient is an oscillatory in nature. Furthermore, when the bed wave number is approximately twice the surface wave number (i.e., when ), Bragg-resonance described earlier by Davies [9], Mei [10], and Davies and Heathershaw [1] takes place between the bed forms and the surface waves. The relation (23) becomes of the form when (at resonance), and in such a situation, we employ L'Hôpital's rule to compute the limiting value and find that the reflection coefficient is given by the following expression: Thus the reflection coefficient becomes a constant multiple of , the number of ripples in the patch. Hence, the reflection coefficient increases linearly with . This indicates that relatively few bottom undulation with its wave number equal to approximately twice the surface wave number may give rise to a very substantial reflected wave.
6. Numerical Results and Discussion
In this section, the first order reflection coefficient as given by (23) is computed numerically for various values of different dimensionless parameter and they are presented graphically.
In Figure 2, the first order reflection coefficient is depicted against wave number for different values of the porous effect parameter and for the dimensionless parameters , , , . From this figure, it is clear that is oscillatory in nature and the peak value of is attained when the wave number of the bottom undulation becomes approximately twice the surface wave number , which validate the theoretical observation. Also the peak value of the reflection coefficient increases as the value of the porous effect parameter increases.
In Figure 3, the first order reflection coefficient is plotted against the wave number for different number of ripples of the bed () and for , , , . From this figure, it is concluded that when the number of ripples of the bed increases, the peak value of also increases, which also validate the theoretical observation of (25). It is also clear that its oscillatory nature is more pronounced with the number of zeros of increased but the general feature of remains the same.
The effect of the surface tension on reflection coefficient is shown for , , , in Figure 4 and for the same set of values with in Figure 5. From these figures, it is clear that the value of the first order reflection coefficient decreases as the surface tension increases.
7. Conclusion
Scattering of water waves by small undulation of the porous bottom of an ocean in the presence of surface tension is investigated here. Perturbation analysis in conjunction with the Fourier transform technique is used to derive the first order velocity potential, the reflection, and transmission coefficients. The main advantage of this method is that we need to solve relatively easier ordinary differential equation to find the Fourier transform of the velocity potential. The derived results for reflection and transmission coefficients are exactly coincide with the known ones obtained earlier in the case when the bed has no porous effect and when the surface tension at the free surface is neglected. The results for the reflection and transmission coefficients for sinusoidal bottom topography are determined in detail and demonstrated graphically. From the computational results, it is clear that the peak value of reflection coefficient increases when the porosity of the uneven ocean bed increases. Another important conclusion is that the value of the reflection coefficient is found to be decreasing with the surface tension at the free surface. The analysis of the present work which is developed mainly for two-dimensional problems is also applicable to problems for the case of oblique incidence (three dimensions), with appropriate assumptions and modifications and this will be the subject matter of our future work.
Acknowledgment
S. Panda is grateful to the Council of Scientific and Industrial Research (CSIR), Government of India, for providing the Junior Research Fellowship for pursuing Ph. D. degree at the Indian Institute of Technology Ropar, India. S. C. Martha also thanks the Indian Institute of Technology Ropar, India, for the support through ISIRD grant and for providing all the necessary facilities. The authors thank Professor Aloknath Chakrabarti, Indian Institute of Science, Bangalore, India, for his valuable suggestions and comments during the preparation of the paper.