# Asymptotic Formulas for the Reflection/Transmission of Long Water Waves Propagating in a Tapered and Slender Harbor

Authors: Eric-Gustavo Bautista; Federico Méndez; Oscar Bautista
Journal: Journal of Applied Mathematics
 (2015)

---

## Abstract

We obtain asymptotic formulas for the reflection/transmission coefficients of linear long water waves, propagating in a harbor composed of a tapered and slender region connected to uniform inlet and outlet regions. The region with variable character obeys a power-law. The governing equations are presented in dimensionless form. The reflection/transmission coefficients are obtained for the limit of the parameterκ2≪1, which corresponds to a wavelength shorter than the characteristic horizontal length of the harbor. The asymptotic formulas consider those cases when the geometry of the harbor can be variable in width and depth: linear or parabolic among other transitions or a combination of these geometries. For harbors with nonlinear transitions, the parabolic geometry is less reflective than the other cases. The reflection coefficient for linear transitions just presents an oscillatory behavior. We can infer that the deducted formulas provide as first approximation a practical reference to the analysis of wave reflection/transmission in harbors.


---

## Body

## 1. Introduction

Long water waves can propagate over large distances from the coast line along slender harbors. In some cases, the geometry of these structures can be strongly tapered, producing the amplification of the wave amplitude; see Li et al. [1]. A property of oscillations in harbors is that even relatively small vertical motions can be accompanied by large horizontal motions of the water, Wiegel [2]; under these circumstances, the oscillations could create unacceptable vessel motions and excessive mooring forces leading to breaking of mooring lines as it was detailed by Rabinovich [3].Inside the harbors can be presented different configurations and connections between two or more basins through a narrow passage [4, 5]; at least one of these is open to the sea, and it is identified as outer harbor and the others are identified as inner harbors. The outer harbor can act like a filter and reduces the energy of waves propagating towards the inner harbor; see, for instance, Bellotti [6].When the water wave is travelling from the outer harbor entrance, it is strongly affected by the geometry of the harbor. In harbors, a small portion of the water waves are reflected repeatedly by the interior boundaries. Some portions of the wave energy escape to the ocean, while others stay inside the harbor, Mei et al. [7]. The wave reflected is an important mechanism that may endanger harbor access, due to sea states at the entrance, Zanuttigh et al. [8]. In addition, if the harbor geometry converges in a narrow passage, we can appreciate the motion of waves with large amplitudes. The reflection and the amplification of the water waves by tapered and slender harbors are effects that should be taken into account in their design and are characterized by the reflection and transmission coefficients.Different researchers conducted fundamental studies related to the propagation of water waves in complex geometries with the aid of experimental and theoretical models [9, 10], made specially to obtain the reflection and transmission coefficients. Among others, one important case that has been widely studied is the propagation of water waves over an infinite step [11, 12] and the deformation of long water waves onto a shelf with abrupt and gradual changes in the depth; see, for example, Goring [13]. On the other hand, the linear long water wave approximation has been used to estimate the reflection and transmission coefficients for different submerged trapezoidal breakwaters, Chang and Liou [14]. In this direction, Jung et al. [15] obtained an approximated analytical solution for the mild-slope equation, applied to the analysis of wave reflection for an asymmetric trench. The aforementioned work was improved by Xie and Liu [16] by obtaining an exact analytical solution to the modified mild-slope equation. Bender and Dean [17] studied the reflection and transmission of normally incident waves by two-dimensional trenches and shoals. Wang et al. [18] using the linear shallow water approximation analyzed the longitudinal oscillations in a rectangular harbor with a hyperbolic-cosine squared bottom induced by incident perpendicular waves. Lin and Liu [19] obtained an analytical solution for linear long wave reflection by an obstacle of general trapezoidal shape. In the same direction, Xie et al. [20] explored the linear long wave reflection by a rectangular obstacle with two scour trenches of power function profile and later Liu et al. [21] proposed an exact analytical solution in the form of a Taylor series for linear long wave reflection by a submerged quasi-idealized breakwater or trench. In addition, there are studies applying perturbation techniques to predict the reflection and modulation of long water waves propagating over geometries with variable depth and in channels with constant depth and parabolic width transition; see, for example, [22, 23].In this work, we obtain asymptotic formulas for the reflection/transmission coefficients of linear long water waves, propagating in a harbor composed of a tapered and slender region connected to uniform inlet and outlet regions. The coefficients are derived in the limit for which the wavelength is smaller than the characteristic horizontal length of a harbor. In order to obtain the above formulas, we have adopted the analytical solution for the free surface elevation of linear long water waves propagating in a tapered and slender channel proposed by Bautista et al. [24], which is obtained with the WKB singular perturbation technique. We consider obviously those cases where the geometry of the harbor can be variable in width and depth: linear or parabolic among other transitions or a combination of these geometries. To verify the derived formula for the reflection coefficient, it is compared against the classical formula for the reflection of linear long water waves propagating on a depth with linear transition. The influence of the dimensions and shape of the harbors on the reflection and transmission coefficients is investigated. The easy application of the reflection/transmission coefficients allows us to clarify which geometrical transitions of the harbor are less or more reflective.


## 2. Formulation

Let us consider a linear long water wave propagating from right to left hand side in a tapered and slender harbor. We assume that the harbor is formed by three regionsR1, R2, and R3, where region R2 defines just the variable character of the cross section, while regions R1 and R3 assume uniform values for the width and depth of the harbor; see Figures 1 and 2. In the selected Cartesian coordinate system, the x-axis is positive to the right with origin placed in the intersection of regions R1 and R2, y-axis is pointing out in the transversal direction of the harbor, and the z-axis is pointing outwards in a normal direction to the mean sea water level. The interval where the transition region R2 occurs is 0≤x≤L2, while the regions R1 and R3 are defined by the intervals -L1≤x≤0 and L2≤x≤L3, respectively, where L1 and L2 are the lengths of the regions R1 and R2, respectively, while L3 corresponds to a length measure from the origin of the reference system x=0 to the inlet of the region R3. In the present analysis, we assume that the walls of the harbor are impermeable. However, the vertical boundaries on the right and left sides of the system are completely open to the flow of the water waves. In particular, we assume the sum of incident and reflected waves on the right and only a transmitted wave on the left. The prescribed widths and depths of the harbors are given by the following relationships:(1)bx=b1for-L1≤x≤0b2x=βxm+b1for0≤x≤L2b3forL2≤x≤L3hx=h1for-L1≤x≤0h2x=αxn+h1for0≤x≤L2h3forL2≤x≤L3.In the above equations, the slenderness ratios α and β are defined by(2)β=b3-b1L2m,α=h3-h1L2n.Figure 1
The top view diagram of the physical model studied. The water wave amplitudeAw for w=I, R, and T represents the incident (travelling from right to the left hand side), reflected, and transmitted water waves, respectively. The harbor is formed for three regions R1, R2, and R3.Figure 2
The side view diagram of the physical model studied. The water wave amplitudeAw for w=I, R, and T represents the incident (travelling from right to the left hand side), reflected, and transmitted water waves, respectively. The harbor is formed for three regions R1, R2, and R3.In addition, the geometry shapes of the width and the depth of the harbors are given by the values ofm≥1 and n≥1, respectively, with m,n~O1. The arguments, for which m and n can not be less than one, are explained in Section 3. The power-law functions for the width and the depth, given by (1), are taken as geometrical representations or approximations to tapered and slender harbors.Under some specific circumstances, in harbors, the water depth is very much smaller compared with the wavelength. Therefore, the wave oscillations can be approximated as shallow water waves, for which water flows are mostly in the horizontal planes, while the vertical velocity variations are weak and the pressure is hydrostatic as a first approximation. The fluid flow is governed by the unidimensional linear continuity and momentum equations, Mei et al. [7, p. 124]:(3)∂ζx,t∂t+1bx∂Ux,thxbx∂x=0,(4)∂Ux,t∂t+g∂ζx,t∂x=0,where ζx,t is the water surface elevation and U(x,t) is the mean water velocity in the x direction; for simplicity, hereafter h=hx, t is the time and g is the gravitational acceleration.The continuity and momentum equations (3) and (4) can be combined in only one. By this purpose, we can derive (3) with respect to the time and replace (4) in the resulting equation, obtaining the following equation:(5)∂2ζx,t∂t2-gbx∂∂xhxbx∂ζx,t∂x.For infinite regions of the wave propagation, typically the wavelength (but not necessarily) is selected as the characteristic horizontal length, Svendsen [25, p. 328]; however, for the present work, which has a well defined geometry, the horizontal characteristic scale over which changes the magnitude of the dependent variables is the length L2, the characteristic time scale of the wave is the period T, and the characteristic scale for the wave amplitude is the incident wave amplitude AI. For nondispersive linear long water waves the characteristic wave celerity is given by c=gh1=ω/k, where the wave frequency is ω=2π/T, the wave number is given by k=2π/λ, and the wavelength in region R1 is denoted by λ.Therefore, for regionR2, we can introduce the following dimensionless variables:(6)χ2=xL2,Δ2=ζx,tAI,τ=ωt.Substituting the previous dimensionless variables into (5), the dimensionless governing equation for linear long waves propagating in a tapered and slender harbor (b3/L2≪1) is given by(7)∂2Δ2∂τ2-κ2b¯χ2∂∂χ2Gχ2∂Δ2∂χ2=0,where Gχ2=b¯χ2h¯χ2 is the dimensionless area of the transversal section of the harbor in the transition region R2, the dimensionless width is written as b¯χ2=1-ɛwχ2m+ɛw, and the dimensionless depth is expressed as h¯χ2=1-ɛhχ2n+ɛh;ɛh=h1/h3 and ɛw=b1/b3. In the above equation, the dimensionless parameter κ2 is defined as (8)κ2=λ2πL22=gh3ωL22.The parameterκ2=λ/2πL22 represents the competition between the wavelength, λ, and the characteristic horizontal length of the harbor, L2. From a physical point of view, the limit of κ2≪1 corresponds to the horizontal projection L2, very large compared with λ.Along the three regions of the harbor, we assume that the waves are harmonic functions in time and can be expressed as(9)Δjχj,τ=Reδjχje-iτ,where j=1,2, and 3 indicates the regions R1, R2, and R3, respectively. The dimensionless amplitudes of the water waves at any cross section of the harbor are denoted by δjχj,i=-1 and Re “is the real part of.”For the regionR2 and from (9), we suggest the following initial conditions:(10)Δ2χ2,0=Reδ2χ2,∂Δ2χ2,0∂τ=Re-iδ2χ2.It should be noted that, in (10), the unknown function δ2χ2 must be determined as a part of the problem.In addition, (7) requires appropriated boundary conditions for Δ2(χ2,τ) at χ2=0 and χ2=1. These conditions are obtained from the continuity condition for the total head pressure at the junctions [7, 26].We define the following dimensionless variables:(11)χ1=-xL1,χ3=x-L2L3-L2,Z¯=zh3,P¯jχj,Z¯,τ=Pjx,z,tρgh3,Δ1χ1,τ=ζ1x,tAI,Δ3χ3,τ=ζ3x,tAI,where j=1,2, and 3 indicates the regions R1, R2, and R3, respectively, Pj(x,z,t) is the pressure, ρ is the water density, and ζ1(x,t) and ζ3(x,t) are the free surface elevations in the regions R1 and R3, respectively.The pressure under linear long water waves is still hydrostatic; that is,Pj(x,z,t)=ρg[ζj(x,t)-z], Dean and Dalrymple [27, p. 133]. Substituting the dimensionless variables (11) into the previous relationship, we obtain the dimensionless version of the pressure(12)P¯jχj,Z¯,τ=ɛΔjχj,τ-Z¯,where ɛ=AI/h3.Considering that the pressures at the connection between the regionsR1-R2 and R2-R3 are equals, we have that(13)P¯2χ2=0,Z¯,τ=P¯1χ1=0,Z¯,τ,(14)P¯2χ2=1,Z¯,τ=P¯3χ3=0,Z¯,τ.Substitution of (12) into (13) and (14) leads to(15)Δ2χ2=0,τ=Δ1χ1=0,τ,Δ2χ2=1,τ=Δ3χ3=0,τ.Introducing (9) into (7) and (15), the following boundary value problem is obtained:(16)κ2d2δ2χ2dχ22+κ2Gχ2dGχ2dχ2dδ2χ2dχ2+δ2χ2h¯χ2=0,with the following boundary conditions:(17)δ2χ2=0=δ1χ1=0,δ2χ2=1=δ3χ3=0.In the above boundary conditions (17), δ1χ1=0 and δ3χ3=0 are the unknown dimensionless amplitudes of the linear long water waves propagating in regions R1 and R3, respectively, which are presented in paragraphs below.The asymptotic solution of (16) for the limit of κ2≪1, obtained with the WKB perturbation technique and reported in detail by Bautista et al. [24], is adopted for the present work and it is given by(18)δ2χ2=1h¯1/4χ2b¯1/2χ2Ceiϑx+De-iϑx.Equation (18) is derived for a channel composed of a region with variables width and depth and connected to inlet and outlet regions which extend with infinite lengths, conversely to the equivalent regions analyzed in the present work, where we assume well defined lengths; however, the infinite lengths have no influence on the behavior of the waves in the region of the harbor with variable section, because the patching boundary conditions are not functions of the lengths of the regions with constant width and depth. Therefore, the above analytical solution is maintained valid for the present work.In (18), the constants C and D must be determined from the appropriated patching boundary conditions. On the other hand, ϑx=Fχ2,n/κ21/2 and Fχ2,n is defined as(19)Fχ2,n=∫1h¯χ2dχ2.Following a similar procedure as that developed for the regionR2 and taking into account relationships (11) and (9), the dimensionless surface elevation for regions R1 and R3 can be determined by solving the following differential equation:(20)d2δjχjdχj2+κj2δjχj=0,where j=1 and 3 indicates the regions R1 and R3, respectively. The analytical solutions of (20) for the water waves amplitude in regions R1 and R3 are trivial and are given by(21)δ1χ1=βTeiκ1χ1,(22)δ3χ3=βReiκ3χ3+e-iκ3χ3,with κ1=ωL1/gh11/2, κ3=ωL3-L2/gh31/2, βR=AR/AI, and βT=AT/AI.In the above equations,AR and AT are the unknown complex amplitudes of the reflected and transmitted waves, respectively. Evaluating (21) at χ1=0, we can obtain that δ1χ1=0=βT and, for χ3=0, δ3χ3=0=βR+1. These relationships are just the unknown boundary conditions in (17).


## 3. Reflection and Transmission Coefficients for the Limit ofκ2≪1

To obtain the values of the coefficientsC, D, βR, and βT, we require two additional boundary conditions, which can be obtained from the continuity of mass flux at the connections between the regions R1-R2 at χ2=0 and R2-R3 at χ2=1. Taking into account the above comments the corresponding patching boundary conditions are given by the following dimensionless relationships:(23)dδ1dχ1χ1=0=-Φdδ2dχ2χ2=0,dδ2dχ2χ2=1=Φ~dδ3dχ3χ3=0,where Φ=L1/L2 and Φ~=L2/L3-L2.From (17), (18), and (23), we can easily derive the following system of equations:(24)eiϑ1e-iϑ1-101+a2ieiϑ1-1-a2ie-iϑ1-10eiϑ2e-iϑ20-a01+a0a1ɛh1/2ieiϑ2-1-a0a1ɛh1/2ie-iϑ20a0CDβRβT=1-100with(25)ϑ1=Fχ2=1,nκ21/2,ϑ2=Fχ2=0,nκ21/2,a0=ɛh1/4ɛw1/2,a1=α1κ21/2,a2=α2κ21/2,α1=m1-ɛwqχ2=0,m2a0ɛw1/4+n1-ɛhgχ2=0,n4a0ɛh,α2=m1-ɛw2+n1-ɛh4.In the above set of equations,Fχ2=0,n and Fχ2=1,n are the solutions of integral (19), evaluated at χ2=0 and χ2=1, respectively, for different values of the exponent n with n≥1. In addition, in relationships (25) for α1, gχ2,n=χ2n-1, and qχ2,m=χ2m-1, these functions for values of m,n>1 and m,n~O1 must be evaluated at χ2=0 obtaining that gχ2=0,n=0 and qχ2=0,m=0 and, for harbors with linear transition, (m=1 and/or n=1), are reduced to gχ2=0,n=qχ2=0,m=1. It is very important to note that the exponents n and m can take different values, depending on the harbor shape geometry, that is, 1,1.5,2,2.5,3,5. In order to satisfy the condition of λ≪L2, represented by κ2≪1 and slender harbor, given by b3/L2≪1, the values of m and n must be restricted to be of order unity. Harbors with m,n≫1 degenerate in an abrupt transition as χ2→1, and therefore the slender assumption is no longer satisfied. It is possible to combine different shapes of the width and the depth. In this mathematical model, the values of m and n must be greater than or equal to unity, because the power-law distribution has the origin located at χ2=0, just in the intersection of regions R1 and R2 and, for the case of n<1 and m<1, the functions gχ2,n and qχ2,m are divergent and, therefore, (18) will be also undefined.The solution of the system of (24) is trivial and is given by(26)C=-i2a1a02ɛh1/2Ωeiϑ1,D=2a02+a1a0ɛh1/2iΩeiϑ1+2ϑ2,βR=1Ωa02a1ɛh1/2a2-2iei2ϑ1-a0a2a0a1ɛh1/2-2iei2ϑ2,βT=4Ωeiϑ1+ϑ2,with Ω=a0a2+2ia0a1ɛh1/2-2iei2ϑ2-a02a1a2ɛh1/2ei2ϑ1.To obtain the reflection and transmission coefficients, we make an energy flow balance: atχ2=0 only a transmitted wave is presented, while at χ2=1 there are incident and reflected waves; therefore, in a simplified form, the energy flux balance is given by(27)CR2+ɛwɛh1/2CT2=1.In the above equation, the reflection and transmission coefficientsCR=βR and CT=βT can be written as(28)CR=1Ωa02a1ɛh1/2a2-2iei2ϑ1-a0a2a0a1ɛh1/2-2iei2ϑ2,(29)CT=4Ω,where  represents the module of a complex number.Taking into account relationships (25), for harbors with width transition m>1 and ɛh=1, (28) and (29) can be written in a simplified form as(30)CR≈m1-ɛwκ21/24,(31)CT≈1ɛw1/21+α22κ2/4,respectively. In a similar form as before, for harbors with depth transition n>1 and ɛw=1, (28) and (29) are reduced to(32)CR≈n1-ɛhκ21/28,(33)CT≈1ɛh1/41+α22κ2/4,respectively.Equations (28)–(33) are developed for slow variations of the depth and width of the harbor, (h3/L2≪1 and b3/L2≪1), and for a characteristic length of the harbor bigger than the wavelength κ2≪1, and for sufficiently small amplitude AI/h3≪1.For the limit ofɛw→1 and ɛh→1, the coefficients take values of CR≈0 and CT≈1.Now, from the dimensionless generalized Green’s law, Lighthill [28, p. 122],(34)δ2χ2=0~1h¯1/4χ2=0b¯1/2χ2=0,and, from the definition of the transmission coefficient CT=βT~Oδ2χ2=0, we have that(35)CT~1h¯1/4χ2=0b¯1/2χ2=0,with h¯χ2=0=ɛh and b¯χ2=0=ɛw; therefore, for the limits of ɛh≪1 and ɛw≪1, the second term on the left side of (27) is ɛwɛh1/2CT2→1, which implies that, to preserve the energy flow, for the above limits, the reflection coefficient must be CR≪1, as is shown in Results and Discussion.From the above comments, in this work, for the limits ofɛw≪1 and ɛh≪1, the reflection coefficient can not be approximated to one and the transmission coefficient can not be approximated to zero, because the present mathematical model is derived for tapered and slender harbors, as was previously commented on. The geometrical conditions for which CR could be approximated to one are those cases of rapidly varying cross section or abrupt transitions, geometrical conditions that are out of the scope of the present work and are not treated here.Equations (28) and (29) allow obtaining values of the reflection and transmission coefficients for harbors with linear or parabolic transitions among other transitions or a combination of these geometries.For the limit ofκ2≪1, (28) and (29) can reproduce the classical case of linear long wave reflection by a depth with linear transition reported by Dean [29], as can be verified in Section 4.2 of Results and Discussion. Although the values of the coefficient CR, obtained for the limit of κ2≪1 and calculated from (28), are small, these allow to identifying those geometrical transitions that are less or more reflective.


## 4. Results and Discussion

In order to investigate the influence of the harbor geometry on the reflection and transmission of linear long water waves, as an example, we use the following typical physical data: the length and width of the harbor areL2=3500m and b3=300m, respectively, and a depth h3=25m, Zwamborn and Grieve [30]; however, because the present mathematical formulation is expressed in dimensionless form, many other representative values of the physical variables can be used. In addition, we propose two widths b1(=300 m and 75 m) and two depths h1(=25 m and 6.25 m); with these values we obtain the width and depth ratios ɛw=1,0.25 and ɛh=1,0.25, respectively. In order to satisfy the long water wave theory, the wavelength in region R3 is λ=695.4 m; therefore, the kinematical parameter takes a value of κ2=0.001. In this manner, the selected data satisfy the assumption of long water wave theory, which in this case is equivalent to h3/λ=0.036≤1/20; see Rahman [31, p. 161] and the fundamental assumption is fulfilled and satisfied.


### 4.1. Wave Reflection by a Harbor with Constant Depth and Variable Width

In Figure3, we show the results of the reflection coefficient CR as a function of the width ratio ɛw for different values of m=1,2,3.5,5, ɛh=1, and κ2=0.001. The results of the coefficient CR for a harbor with a linear transition width m=1 are obtained by using (28) and for those cases of harbors with m>1 and ɛh=1, we use (30). From (30) it can be deduced that the reflection coefficient for a specific geometry of the width transition m increases or decreases its magnitude depending on the assumed values of the parameter κ2=λ/2πL22. For values of the parameter κ2 such that L2→λ, the harbors are more reflective than the case with L2≫λ. On the other hand, for constant values of the parameters κ2 and ɛw, the reflection coefficient increases or decreases its magnitude proportionally in terms of the values of m, where always m>1 and m~O1.Figure 3
Reflection coefficient as a function of the dimensionless parameterɛw, for waves propagating in a harbor with width transitions m=1,2,3.5,5, with ɛh=1 and κ2=0.001.In general, from Figure3 we show that a harbor with linear width transition is less reflective than harbors with nonlinear transitions; however, as ɛw→0, the reflection coefficient for m=2,3.5,5 takes finite values of CR=0.01581,0.02766,0.03949, respectively. While the harbor with m=1, the reflection coefficient increases its values exponentially.In Figure4, the behaviors of the reflection CR and transmission CT coefficients as a function of the parameter κ2 are plotted for four values of m=1,2,2.5,5 with constant values of ɛw=0.25 and ɛh=1. For the cases of m>1 the results of reflection coefficients are calculated from (30) and for m=1 we use (28). The transmission coefficients CT for the values of m>1 and m=1 are obtained from (31) and (29), respectively.Figure 4
Reflection/transmission coefficients as a function of the kinematical parameterκ2 with ɛh=1 and ɛw=0.25 for four values of m=1,2,2.5,5.From (31) with κ2→0, which means that L2≫λ, the transmission coefficient is CT~Oɛw-1/2 and the harbor geometry has no influence on the wave transmission; however, for values of κ2 such that L2→λ the term α22κ2/4 must be retained in order to identify the influence of harbor geometry on the wave transmission.In general, Figure4 reveals that for large values of the parameter κ2 the coefficient CR increases and CT decreases. From a physical point of view, this condition is presented only for L2→λ; otherwise, for L2≫λ the harbors are less reflective. In the limit of CR→0, (27) is simplified to CT=ɛh-1/4ɛw-1/2; therefore, we can substitute the values of ɛw=0.25 and ɛh=1 to obtain that CT=2, which is the maximum value observed in Figure 4. Under these conditions, the maximum value of the transmission coefficient for this case is CT~Oɛw-1/2. It can be appreciated that, in the interval plotted for the parameter κ2 and with the exception of m=1, the other transition shapes do not have an oscillatory behavior, this characteristic also can be observed in (30) and (31), which do not include oscillating terms, and these are only presented in the generalized formulas used to obtain the coefficients CR and CT for the harbor with linear transition width, (28) and (29), respectively.


### 4.2. Wave Reflection by a Harbor with Constant Width and Variable Depth

In Figure5, we show the values of the reflection coefficient CR as a function of the depth ratio ɛh with constant values of κ2=0.001 and ɛw=1 and four values of n=1,2,3.5,5. From (32), we obtain similar results to those obtained in Section 4.1, for specific values of the parameters n, ɛh, and κ2. For constant values of the parameters κ2 and ɛh the reflection coefficient increases its magnitude for large values of n>1 and n~O1. Comparing the nonlinear transitions, just the shape geometry with n=2 is less reflective than those cases of n(=3.5 and 5); see Figure 5.Figure 5
Reflection coefficient as a function of the dimensionless parameterɛh, for waves propagating on a depth with transitions n=1,2,3.5,5 with ɛw=1 and κ2=0.001.For the particular case ofn=1, the reflection coefficient is calculated from (28) and it is compared against the analytical model proposed by Dean [29], which is given by(36)CR=γ0M1+γ1M0-2v0v1-8/π2E02ɛh1/2γ0M1+γ1M0-2v0v1+8/π2E02ɛh1/2,where(37)E0=4πL2λ11-ɛh,E1=4πL2λɛh1/21-ɛh,γ0=J02E0+J12E0,γ1=J02E1+J12E1,M0=Y02E0+Y12E0,M1=Y02E1+Y12E1,ν0=J0E0Y0E0+J1E0Y1E0,ν1=J0E1Y0E1+J1E1Y1E1.In the above equations,J0 and J1 are the Bessel functions of first kind of orders 0 and 1, respectively. Y0 and Y1 are Bessel functions of second kind of orders 0 and 1, respectively. We must note that the results are properly adjusted. In addition, Figure 5 reveals that, for ɛh→0, the reflection coefficient CR for the lineal transitions, n=1, grows exponentially.In Figure6, we show the reflection coefficient CR and transmission coefficient CT as a function of the parameter κ2 with constant values of ɛw=1 and ɛh=0.25.Figure 6
Reflection/transmission coefficients as a function of the kinematical parameterκ2 with ɛw=1 and ɛh=0.25 for four values of n=1,2,2.5,5.The depth transitions are given for four values ofn=1,2,2.5,5. It should be noted that the reflection and transmission coefficients for a linear bottom transition have an oscillatory behavior with strong oscillations as κ2→0. For those cases of n>1 the parabolic transition (n=2) is less reflective than the other cases. The complementary transmission coefficients are shown at the top of the same figure. For values of the parameter κ2→0.008, the bottom shape with n=5 has a transmission coefficient smaller than that of n=1,2,2.5. In the case of κ2=0.008 the transmission coefficient for n=5 is CT=1.41297 while that for n=2 is CT=1.41401, which represents a difference of 7.4%.


### 4.3. Wave Reflection by a Harbor with Variables Width and Depth

In Figure7, we show the reflection and transmission coefficients for those harbors with variable geometry distribution with m=n=1,2,2.5,5 as a function of the parameter κ2, with constant values of ɛh=ɛw=0.25. It is found that, for κ2=0.0005, the maximum value of the transmission coefficient is near CT=2.827 where always CR→0, for the four harbor geometries. In the above limit, from (27) the maximum value of the transmission coefficient is CT~Oɛw-1/2ɛh-1/4, which depends only on the width and depth ratios, ɛw and ɛh, respectively. On the other hand, for a bigger value than κ2=0.0005, that is, κ2=0.008, the coefficient CR for m=n=5 takes a value of CR=0.125 and the corresponding one for m=n=2 is CR=0.02, which represents that a harbor with geometry distribution m=n=5 is 6.25 times greater than the case of m=n=2.Figure 7
Reflection/transmission coefficients as a function of the kinematical parameterκ2 with ɛh=ɛw=0.25 for four values of m=n=1,2,2.5,5.


## 4.1. Wave Reflection by a Harbor with Constant Depth and Variable Width

In Figure3, we show the results of the reflection coefficient CR as a function of the width ratio ɛw for different values of m=1,2,3.5,5, ɛh=1, and κ2=0.001. The results of the coefficient CR for a harbor with a linear transition width m=1 are obtained by using (28) and for those cases of harbors with m>1 and ɛh=1, we use (30). From (30) it can be deduced that the reflection coefficient for a specific geometry of the width transition m increases or decreases its magnitude depending on the assumed values of the parameter κ2=λ/2πL22. For values of the parameter κ2 such that L2→λ, the harbors are more reflective than the case with L2≫λ. On the other hand, for constant values of the parameters κ2 and ɛw, the reflection coefficient increases or decreases its magnitude proportionally in terms of the values of m, where always m>1 and m~O1.Figure 3
Reflection coefficient as a function of the dimensionless parameterɛw, for waves propagating in a harbor with width transitions m=1,2,3.5,5, with ɛh=1 and κ2=0.001.In general, from Figure3 we show that a harbor with linear width transition is less reflective than harbors with nonlinear transitions; however, as ɛw→0, the reflection coefficient for m=2,3.5,5 takes finite values of CR=0.01581,0.02766,0.03949, respectively. While the harbor with m=1, the reflection coefficient increases its values exponentially.In Figure4, the behaviors of the reflection CR and transmission CT coefficients as a function of the parameter κ2 are plotted for four values of m=1,2,2.5,5 with constant values of ɛw=0.25 and ɛh=1. For the cases of m>1 the results of reflection coefficients are calculated from (30) and for m=1 we use (28). The transmission coefficients CT for the values of m>1 and m=1 are obtained from (31) and (29), respectively.Figure 4
Reflection/transmission coefficients as a function of the kinematical parameterκ2 with ɛh=1 and ɛw=0.25 for four values of m=1,2,2.5,5.From (31) with κ2→0, which means that L2≫λ, the transmission coefficient is CT~Oɛw-1/2 and the harbor geometry has no influence on the wave transmission; however, for values of κ2 such that L2→λ the term α22κ2/4 must be retained in order to identify the influence of harbor geometry on the wave transmission.In general, Figure4 reveals that for large values of the parameter κ2 the coefficient CR increases and CT decreases. From a physical point of view, this condition is presented only for L2→λ; otherwise, for L2≫λ the harbors are less reflective. In the limit of CR→0, (27) is simplified to CT=ɛh-1/4ɛw-1/2; therefore, we can substitute the values of ɛw=0.25 and ɛh=1 to obtain that CT=2, which is the maximum value observed in Figure 4. Under these conditions, the maximum value of the transmission coefficient for this case is CT~Oɛw-1/2. It can be appreciated that, in the interval plotted for the parameter κ2 and with the exception of m=1, the other transition shapes do not have an oscillatory behavior, this characteristic also can be observed in (30) and (31), which do not include oscillating terms, and these are only presented in the generalized formulas used to obtain the coefficients CR and CT for the harbor with linear transition width, (28) and (29), respectively.


## 4.2. Wave Reflection by a Harbor with Constant Width and Variable Depth

In Figure5, we show the values of the reflection coefficient CR as a function of the depth ratio ɛh with constant values of κ2=0.001 and ɛw=1 and four values of n=1,2,3.5,5. From (32), we obtain similar results to those obtained in Section 4.1, for specific values of the parameters n, ɛh, and κ2. For constant values of the parameters κ2 and ɛh the reflection coefficient increases its magnitude for large values of n>1 and n~O1. Comparing the nonlinear transitions, just the shape geometry with n=2 is less reflective than those cases of n(=3.5 and 5); see Figure 5.Figure 5
Reflection coefficient as a function of the dimensionless parameterɛh, for waves propagating on a depth with transitions n=1,2,3.5,5 with ɛw=1 and κ2=0.001.For the particular case ofn=1, the reflection coefficient is calculated from (28) and it is compared against the analytical model proposed by Dean [29], which is given by(36)CR=γ0M1+γ1M0-2v0v1-8/π2E02ɛh1/2γ0M1+γ1M0-2v0v1+8/π2E02ɛh1/2,where(37)E0=4πL2λ11-ɛh,E1=4πL2λɛh1/21-ɛh,γ0=J02E0+J12E0,γ1=J02E1+J12E1,M0=Y02E0+Y12E0,M1=Y02E1+Y12E1,ν0=J0E0Y0E0+J1E0Y1E0,ν1=J0E1Y0E1+J1E1Y1E1.In the above equations,J0 and J1 are the Bessel functions of first kind of orders 0 and 1, respectively. Y0 and Y1 are Bessel functions of second kind of orders 0 and 1, respectively. We must note that the results are properly adjusted. In addition, Figure 5 reveals that, for ɛh→0, the reflection coefficient CR for the lineal transitions, n=1, grows exponentially.In Figure6, we show the reflection coefficient CR and transmission coefficient CT as a function of the parameter κ2 with constant values of ɛw=1 and ɛh=0.25.Figure 6
Reflection/transmission coefficients as a function of the kinematical parameterκ2 with ɛw=1 and ɛh=0.25 for four values of n=1,2,2.5,5.The depth transitions are given for four values ofn=1,2,2.5,5. It should be noted that the reflection and transmission coefficients for a linear bottom transition have an oscillatory behavior with strong oscillations as κ2→0. For those cases of n>1 the parabolic transition (n=2) is less reflective than the other cases. The complementary transmission coefficients are shown at the top of the same figure. For values of the parameter κ2→0.008, the bottom shape with n=5 has a transmission coefficient smaller than that of n=1,2,2.5. In the case of κ2=0.008 the transmission coefficient for n=5 is CT=1.41297 while that for n=2 is CT=1.41401, which represents a difference of 7.4%.


## 4.3. Wave Reflection by a Harbor with Variables Width and Depth

In Figure7, we show the reflection and transmission coefficients for those harbors with variable geometry distribution with m=n=1,2,2.5,5 as a function of the parameter κ2, with constant values of ɛh=ɛw=0.25. It is found that, for κ2=0.0005, the maximum value of the transmission coefficient is near CT=2.827 where always CR→0, for the four harbor geometries. In the above limit, from (27) the maximum value of the transmission coefficient is CT~Oɛw-1/2ɛh-1/4, which depends only on the width and depth ratios, ɛw and ɛh, respectively. On the other hand, for a bigger value than κ2=0.0005, that is, κ2=0.008, the coefficient CR for m=n=5 takes a value of CR=0.125 and the corresponding one for m=n=2 is CR=0.02, which represents that a harbor with geometry distribution m=n=5 is 6.25 times greater than the case of m=n=2.Figure 7
Reflection/transmission coefficients as a function of the kinematical parameterκ2 with ɛh=ɛw=0.25 for four values of m=n=1,2,2.5,5.


## 5. Conclusions

In this work, for the limit ofκ2≪1, we obtain simple asymptotic formulas for the reflection and transmission coefficients of linear long water waves, propagating in a tapered and slender harbor. The present mathematical formulation can be used as a first approximation to identify which harbors are less or more reflective than other geometries. The new formulas for the reflection and transmission coefficients are general and are very easy to apply.The following summarizes the conclusions of the present work:(1)
It is found that a harbor with constant depth and variable width is more reflective than the case of harbors with constant width and variable depth.(2)
The results show that, for small values of the kinematic parameterκ2, the harbors are less reflective and that the maximum value of the transmission coefficient is CT~Oɛh-1/4ɛw-1/2, indicating that we can appreciate the motion of waves with large amplitudes for small values of ɛh and ɛw.(3)
The combination of exponentsm~O1 and n~O1 can take fractional or integer values with m≥1 and n≥1.(4)
As an example, for a harbor with variables width and depth,m=n=2 and ɛw=ɛh=0.25, the results show that the energy flux balance, (27), CR2+ɛwɛh1/2CT2=0.999 is fully satisfied.(5)
The geometrical transitions studied in this work can be considered as idealized representations of harbors. The applications of the present formulas for the coefficientsCR and CT are not limited to the studied cases in the present work, because the methodology used here can be extended to more complex harbor transitions.



---

*Source: 184147-2015-07-27.xml*