# An Extended Analytic Solution of Combined Refraction and

Authors: Yu-Shu Kuo; Tai-Wen Hsu; Chia-Cheng Tsai; Yu-Hsuan Huang
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

An analytic solution of long waves scattering by a cylindrical island mounted on a permeable circular shoal was obtained by solving the linear long wave equation (LWE). The solution is in terms of the Bessel function expressed by complex variables. The present solution is suitable for arbitrary bottom configurations described by a power function with two independent parameters. For the case of the paraboloidal shoal, there exists a singular point (α=2) which can be removed using Frobenius series, where α is a real constant. The present solution is reduced to Yu and Zhang’s (2003) solution for impermeable circular shoal. The numerical results show some special features of the combined effect of wave refraction and diffraction caused by a porous circular island. The effect of key parameters of the island dimension, the shoal slope, and permeability on wave scattering was discussed based on the analytic solution.


---

## Body

## 1. Introduction

The combined effect of wave refraction and diffraction is caused by the varying bottom topography or the presence of islands and structures. The mild-slope equation (MSE) was first derived by Berkhoff [1] and then was extended by Porter and Staziker [2], Hsu and Wen [3] to the rapidly varying sea bottom configuration. Up to now, the MSE still plays an important role in coastal and ocean engineering practice, as it can be widely implemented to solve the problems of combined refraction and diffraction based on linear wave theory.Although numerical models are commonly used in practical applications, it is still limited by the problems of accuracy and numerical stability. Numerical results are often needed to be verified by comparing with measured data or analytic solutions. In some cases, it takes time and expensive cost for laboratory experiments because of the complexity in completion of the model and facilities setup. An alternate is to look for analytic solution that could provide more clear expressions in which the physical properties of wave motion can be used to examine the predictability and validity of the numerical models.A typical example of the analytic study of combined refraction and diffraction is the case when the long-wave approximation is applied to the MSE in one dimension or axisymmetric bathymetry in two dimensions. Homma [4] is a pioneer to obtain the solution of long-wave equation (LWE) by considering long waves propagating over the top of a parabolic shoal. The solution was further extended and applied by Vastano and Reid [5] and Jonsson et al. [6]. Different geometries of circular or conical island were explored by Fujima and Goto [7], Fujima et al. [8], Zhang and Zhu [9], and Zhu and Zhang [10]. The research field is still under rapid developments. Details can be referred to recent articles of Jung et al. [11], Niu and Yu [12], Liu and Xie [13], and Jung and Lee [14]. Furthermore, there are also studies, such as Liu et al. [15] and Niu and Yu [16], using Hunt’s [17] approximate solution of the implicit dispersion equation to solve the MSE without the long-wave approximation.The present investigation is to extend Yu and Zhang’s [18] analytic solution to long wave propagation over a porous circular island with arbitrary topography. The combined effect of refraction and diffraction of long wave motion as well as key features of the intensity of wave ray focus behind the shoal is presented. A general solution was derived based on the mild slope equation (MSE). The wave motion with combined refraction and diffraction was described by the Bessel function (α≠2) and Frobenius series. The effect of key parameters on wave scattering is intensively discussed based on the analytic solution.


## 2. Theoretical Formulation

The long wave equation for combined refraction and diffraction of waves propagating over porous media is express as(2.1)∇h((hw+δhp)∇hϕ)+k2(hw+δhp)ϕ=0,
where ϕ is the velocity potential, ∇h=(∂/∂x,∂/∂y) the horizontal gradient operator, hw the thickness of the water, hp the thickness of porous layer, k is the wavenumber and δ is defined as n0/(S+ifp) with n0 the porosity, S=n0+(1-n0)CM, the inertial coefficient, CM the virtual mass coefficient, i=-1 the unit complex number, and fp the friction coefficient. In the derivation of (2.1), the time harmonic eiωt is assumed with ω being the angular frequency of the incident wave. Here, the angular frequency ω is related to the wavenumber k by the following dispersion relation:
(2.2)ω2g=k(δkhp+khw).
Equations (2.1) and (2.2) can be derived from the mild-slope equation of Rojanakamthorn et al. [19] by neglecting terms of orders O(k2hw2,k2hp2,k2hwhp). Note that (2.1) and (2.2) are only valid for long waves, that is, the relative depth is restricted in the shallow water region, that is, hw/L≤1/20 and hp/L≤1/20 where L is the wave length. More details can be found in the Appendix.In order to study the wave scattering over a circular porous shoal, we assume that the shoal is symmetric and of the following form:(2.3)hw=hw(r)=h-wr0αrα∝rα,r≤r0,hp=hp(r)=h-pr0αrα∝rα,r≤r0,
where r is the distance of any point on the shoal from the center, r1 is the radius of the circular island, r0 is the radius of the impermeable circular shoal, h-w and h-p are the depths of water and porous layers, respectively, for r>r0 and α is the power parameter of bottom shape. The definition sketch of the permeable circular island is presented in Figure 1.Figure 1
Definition sketch of a circular island.Using the method of separation variables, the velocity potential is written as(2.4)ϕ=∑n=0∞Rn(r)Θn(θ)
in which Θn is the eigenfunction given by
(2.5)Θn(θ)=C1ncosnθ+C2nsinnθ(n=0,1,2,…),
where C1n and C2n are coefficients to be determined. C2n is set to be zero due to the symmetric property of the present problem with respect to the coordinate x. For convenience of computational procedure, the coefficient C1n is specified as C1n=1. Substituting (2.4) and (2.5) into (2.1) gives
(2.6)r2d2Rndr2+r(1+rh~dh~dr)+(ω2r2gh~-n2)Rn=0,
where h~ is an equivalent water depth defined by
(2.7)h~=hw+δhp.
The following dimensionless variables are introduced for convenience of mathematical formulation:
(2.8)ρ=ω2rg,ξw2=ω2hwg,ξp2=ω2hpg.
Substitution of ξw2 and ξp2 from (2.8) into (2.3) leads to the resultant equations.
(2.9)ξw2=ξ-w2ρ0αρα,ρ≤ρ0,(2.10)ξp2=ξ-p2ρ0αρα,ρ≤ρ0.
Equation (2.10) is multiplied by δ and added with (2.9). The following equation is thus obtained, that is
(2.11)ξ2=ξw2+δξp2=ξ-w2ρ0αρα+δξ-p2ρ0αρα=(ξ-w2ρ0α+δξ-p2ρ0α)ρα=1β2ρα.And, (2.6) is rearranged with the help of the dimensionless variable ρ=ω2r/g (2.8) and the result reads
(2.12)ρ2d2Rndρ2+(α+1)ρdRndρ+(β2ρ2-α-n2)Rn=0.
For the special case of α=2, (2.12) is solved by making use of the Frobenius series, the solution is given by
(2.13)Rn={C3nρ-1+μn+C4nρ-1-μnif1+n2-β2≠0C3nρ-1+C4nρ-1lnρif1+n2-β2=0,
where C3n and C4n are constants, and μn is represented by
(2.14)μn=|1+n2-β2|.
We introduce two dimensionless parameters Qn=ρα/2Rn and χ=(2β/|2-α|)ρ1-α/2 for α≠2, (2.6) is rewritten in the form
(2.15)χ2d2Qndχ2+χdQndχ+(χ2-νn2)Qn=0,
where νn=(α2+4n2)/(|2-α|).The general solution of (2.15) is the Bessel function expressed as
(2.16)Qn=C3nJνn(χ)+C4nYνn(χ)(n=0,1,2…),
where Jνn and Yνn are the Bessel functions of the first and second kind of vn order. Substitution the corresponding boundary conditions in (2.16) results in the particular solution for α≠2 given by
(2.17)Rn(ρ)=ρ-α/2[C3nJνn(2β|2-α|ρ1-(α/2))+C4nYνn(2β|2-α|ρ1-(α/2))](n=0,1,2,...).
For an arbitrary bathymetry of a circular island, the solution is thus written as follows
(2.18)Rn(ρ)=C3nR1n(ρ)+C4nR2n(ρ),
where
(2.19)R1n(ρ)={ρ-1+μnifα=2&1+n2-β2≠0,ρ-1ifα=2&1+n2-β2=0,ρ-α/2Jνn(2β|2-α|ρ1-(α/2))ifα≠2,R2n(ρ)={ρ-1-μnifα=2&1+n2-β2≠0,ρ-1lnρifα=2&1+n2-β2=0,ρ-α/2Yνn(2β|2-α|ρ1-(α/2))ifα≠2.
It is noted that the LWE of (2.1) has been solved for long waves propagating over a permeable circular island. The velocity potential is finally expressed as
(2.20)ϕ=∑n=0∞[C3nR1n(ρ)+C4nR2n(ρ)]cosnθ.
Considering an incident wave expressed in a nondimensional form:
(2.21)ϕ0=AIeik0x=AI∑n=0∞inεnJn(k0r)cosnθ=AI∑n=0∞inεnJn(κ0ρ)cosnθ,
where AI is the reference wave height defined as AI=ϕ0(ρ=0). In addition, the wavenumber k0 for r>r0 can be obtained from the dispersion relation (2.2), and the dimensionless wave number κ0 is defined as k0g/ω2. In addition, εn is the Jacobian parameter given by
(2.22)εn={1n=0,2n≥1.
The physical problem is presented in Figure 1. Following Yu and Zhang [18], the method of eigenexpansion is adopted to solve the problem. The study domain is divided into two subregions on the horizontal plane: the finite region with varying depth (r1<r<r0) and semi-infinite far region with a porous constant depth (r>r0). The general solution of the complex velocity potential subject to the Somerfield condition as well as the symmetric condition is written in the form
(2.23)ϕ1=ϕ0+∑n=0∞DnΨn(ρ),(2.24)Ψn(ρ)≡Hn(κ0ρ),
where Hn is the Hankel function of the first kind of the nth order, Dn are complex coefficients to be determined.In the finite region, the permeable bottom configuration is represented by (2.3). The complex velocity potential is given by
(2.25)ϕ2=∑n=0∞[C3nR1n(ρ)+C4nR2n(ρ)]cosnθ.
Because the present investigation assumes the symmetric condition, the corresponding terms of the odd function sinnθ are negligible.The following boundary conditions are used for solving the problem:(2.26)∂ϕ2∂ρ=0forρ=ρ1(orr=r1).
Using (2.25), the relation between two coefficients C3n and C4n is given by
(2.27)Cn=C3nR2n′(ρ1)=C4nR1n′(ρ1).
The prime denotes the first derivative with respect to ρ. The rearrangement of (2.25) yields
(2.28)ϕ2=∑n=0∞Cn[R2n′(ρ1)R1n(ρ)+R1n′(ρ1)R2n(ρ)]cosnθ.
Two matching boundary conditions of continuity of flow mass and pressure are again adopted which are given by
(2.29)ϕ1=ϕ2forρ=ρ0,∂ϕ1∂ρ=∂ϕ2∂ρforρ=ρ0.
Substituting (2.29) into (2.23) and (2.25), we have
(2.30)CnΓn=AIinεnJn(κ0ρ0)+DnΨn(ρ0),CnΛn=AIinεnJn′(κ0ρ0)+DnΨn′(ρ0),
where
(2.31)Γn=R2n′(ρ1)R1n(ρ0)-R1n′(ρ1)R2n(ρ0),Λn=R2n′(ρ1)R1n′(ρ0)+R1n′(ρ1)R2n′(ρ0).The coefficientsCn and Dn are solved by (2.30) and are expressed by
(2.32)Cn=-AIinεnΨn(ρ0)Jn′(κ0ρ0)-Ψn′(ρ0)Jn(κ0ρ0)ΓnΨn′(ρ0)-ΛnΨn(ρ0),Dn=-AIinεnΓnJn′(κ0ρ0)-ΛnJn(κ0ρ0)ΓnΨn′(ρ0)-ΛnΨn(ρ0).


## 3. Results and Discussions

In the following numerical results, we have typically setω2/g=1, the inertial coefficient S=1 and h-p=h-w=ξ-p2=ξ-w2=0.02. The infinite series in the solution of the velocity potential should be truncated, indicating that the number of N must be large enough to represent the infinity so that the infinite series in (2.23) and (2.28) can effectively accounts for the wave motion by a prescribed accuracy. For the case of α=1 and α=2, respectively. Figure 2 presents the variation of the relative wave runup height (A/AI) around the cross-section along the y=0. The result demonstrates the wave height distribution at different values of α for different permeable materials with variant porosities. In Figure 2(a), the wave height distribution before the island increases monotonically with the increase of n0, but the variation of wave height distribution behind the island decreases monotonically with the increasing n0 for α=1, r0=0.2π. Note that the variation of wave height distribution at α=2 shows the same feature shown in Figure 2(b).Effects of the porosity of the permeable material (ρ0=0.2π, ρ1=0.1π, fp=1).
(a)
α
=
1
(b)
α
=
2Figure3 demonstrates the variation of wave height distribution at different values of α as a linearized friction coefficient changes. We notice that the wave height increases with the increase of the friction factor fp in the front of the island but opposite property behind the island.Effects of linearized friction coefficient (ρ0=0.2π, ρ1=0.1π, n0=0.5).
(a)
α
=
1
(b)
α
=
2Figures4 and 5 show different bottom slopes for different values of α, ρ0, ξ0 and ρ1. It is interesting to note that a higher value could produce a steeper bottom geometry, the intensity of wave ray focusing increases with the decrease of the bottom slope. Consequently, a significant increase of wave runup height is evident.Effects of island on wave scattering forα=1 (ρ0=0.2π, fp=1, n0=0.5).
(a)
ρ
1
=
0.02
π
(b)
ρ
1
=
0.15
πEffects of island on wave scattering forα=2 (ρ0=0.2π, fp=1, n0=0.5).
(a)
ρ1=0.02π
(b)
ρ
1
=
0.15
π


## 4. Conclusions

An analytic solution of the long-wave equation for waves propagating over a porous circular island was derived. The solution is in terms of the Bessel functions. The singularity was removed by the Frobenius series for the case of a parabolic island. The solution is able to account for the combined effect of refraction and diffraction around a circular island on the top of a shoal for arbitrary geometry described by two independent parameters. The effect of key parameters on waves scattering was investigated based on numerical calculations from the analytic solution.



---

*Source: 742705-2012-11-13.xml*