Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2014, Article ID 341385, 15 pages http://dx.doi.org/10.1155/2014/341385
Research Article Second-Order Time-Dependent Mild-Slope Equation for Wave Transformation
Ching-Piao Tsai,1 Hong-Bin Chen,2 and John R. C. Hsu3,4
1 Department of Civil Engineering, National Chung Hsing University, Taichung 402, Taiwan 2 Department of Leisure and Recreation Management, Chihlee Institute of Technology, New Taipei 220, Taiwan 3 Department of Marine Environment & Engineering, National Sun Yat-Sen University, Kaohsiung 804, Taiwan 4 School of Civil & Resources Engineering, University of Western Australia, Crawley, WA 6009, Australia
Correspondence should be addressed to Ching-Piao Tsai; [email protected]
Received 16 October 2013; Accepted 1 May 2014; Published 4 May 2014
Academic Editor: Efstratios Tzirtzilakis
Copyright Β© 2014 Ching-Piao Tsai 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.
This study is to propose a wave model with both wave dispersivity and nonlinearity for the wave field without water depth restriction. A narrow-banded sea state centred around a certain dominant wave frequency is considered for applications in coastal engineering. A system of fully nonlinear governing equations is first derived by depth integration of the incompressible Navier-Stokes equation in conservative form. A set of second-order nonlinear time-dependent mild-slope equations is then developed by a perturbation scheme. The present nonlinear equations can be simplified to the linear time-dependent mild-slope equation, nonlinear longwave equation, and traditional Boussinesq wave equation, respectively. A finite volume method with the fourth-order Adams-Moulton predictor-corrector numerical scheme is adopted to directly compute the wave transformation. The validity of the present model is demonstrated by the simulation of the Stokes wave, cnoidal wave, and solitary wave on uniform depth, nonlinear wave shoaling on a sloping beach, and wave propagation over an elliptic shoal. The nearshore wave transformation across the surf zone is simulated for 1D wave on a uniform slope and on a composite bar profile and 2D wave field around a jetty. These computed wave height distributions show very good agreement with the experimental results available.
1. Introduction mild-slope equation was in an elliptic type, it was not efficient in performing direct numerical calculation and also Accurate estimation of design waves, especially wave heights in treating the open boundary conditions. The conventional and their spatial distributions, is essential for the plan- schemes also required large computer storage and computing ning and installation of coastal structures, calculation of time. Later, numerical models, based on an elliptic form of longshore sediment transport, and prediction of nearshore the mild-slope equation, were proposed to compute the wave topographical changes. In practical engineering applications, field for a large coastal area2 [ β5]. These models obtained a narrow-banded sea state centred on a certain dominant the solution for the governing equations by iterations using wave frequency has remained popular for several decades, powerful and sophisticated algorithm. despite the use of random wave theories in more recent time. The mild-slope equation was also cast in a parabolic Because wave transformation from deep water into shallow form to promote computational efficiency [6] and to handle region involves nonlinear process due to shoaling, combined the problem with weak wave reflection7 [ ]. However, the refraction and diffraction, and reflection and breaking, var- parabolicapproximationmethodwasfoundinadequatefor ious forms for mild-slope equations have been developed. determining the wave field in the presence of a highly reflec- Previously, wave transformations involving refraction and tive structure. Later, Copeland [8] expressed the mild-slope diffraction had been calculated separately for simple beach equation as a pair of first-order hyperbolic equations, without topography, until Berkhoff [1] who presented a linear mild- thelossofthereflectedwave.WatanabeandMaruyama[9] slope equation for the combined effects. Because the original extended the set of time-dependent mild-slope equations 2 Mathematical Problems in Engineering including a wave energy dissipation factor for modelling 2. Nonlinear Wave Equations the nearshore wave field. Tsai et al. [10]incorporatedboth the nonlinear shoaling factor and the breaking-wave energy 2.1. Governing Equations. Consider a three-dimensional dissipation factor into the time-dependent mild-slope equa- wave motion in an incompressible, inviscid, and homoge- tionstoimprovetheaccuracyofthecalculationofwave neous fluid, the principle of conservation of mass gives the transformation across the surf zone. Taking large bottom following continuity equation: gradient into the mild-slope equation, the so-called extended ππ€ mild-slope equation or the modified mild-slope equation ββ π’+β =0 (1) were also developed for waves over rapid variation of bottom ππ§ [11, 12]. On the other hand, the complementary mild-slope equation was derived in terms of a stream function vector and the incompressible Navier-Stokes equations in conserva- rather than in terms of a velocity potential [13, 14]. All of these tive form lead to Eulerβs equations as are in the family of linear mild-slope equations. ππ’β ππ’β 1 Based on the second-order Stokes wave theory, the + π’β ββ π’+π€β =β βπ, (2) nonlinear mild-slope equation with parabolic type was devel- ππ‘ ππ§ π oped [15, 16]. These theories are adaptable to the relatively ππ€ ππ€ 1 ππ deeper water. Belibassakis and Athanassoulis [17]proposed + π’β βπ€+π€β =β βπ ππ‘ ππ§ π ππ§ (3) a coupled-mode theory, by using a local mode expansion of the wave potential in conjunction with Luke variational principle, to extend the second-order Stokes theory for the in the horizontal and vertical directions of a Cartesian π’β = (π’, V) π€ case of a generally shaped bottom profile connecting two half- coordinates system, respectively. In (1)β(3), and strips of constant depths; their model was then extended with are the Eulerian velocity components in the horizontal plane π π application to nonlinear water waves propagating in variable and vertical direction, respectively; is the pressure, is the π π‘ bathymetry regions [18]. On the other hand, Boussinesq- density of the fluid, is the acceleration due to gravity, is the β = (π/ππ₯, π/ππ¦) type equations developed by Peregrine [19] were extended to time, and is the horizontal vector operator, π₯ π¦ π§ cover varying depth adaptable to the shallow water [20]. The where , ,and are the Cartesian coordinates. improved Boussinesq-type equations can provide the wave The fluid motion described by(1)β(3)hastosatisfy dispersion and shoaling characteristics from shallow water to boundary conditions at the free surface and on the seabed. water depth with πβ =40,whereπ isthewavenumberand First, the dynamic and kinematic boundary conditions at the β is the water depth [21]. Recently, Li [22] derived nonlinear free surface are wave equations based on combining the second-order Stokes π=0 π§=π(π₯, π¦,) π‘ , wave theory with the Boussinesq theory. (4) Many investigators have tried to develop nonlinear wave ππ equations adaptable to arbitrary water depth. For instance, a π€= + π’β βπβ π§=π(π₯,π¦,π‘) (5) ππ‘ weakly nonlinear and fully dispersive model was introduced by Nadaoka et al. [23], who generated a set of equations and the bottom boundary condition on an impermeable that rendered a single time-dependent nonlinear mild-slope seabed is equation for water waves [24]. The numerically simulated cnoidal waves and the second-order Stokes waves were found π€=βπ’β βββ π§=ββ(π₯,π¦) (6) in excellent agreement with theoretical solutions. However, their nonlinear model was only simulated to unidirectional π(π₯, π¦, π‘) β(π₯, π¦) waves, rather than a two-dimensional wave field involving in which is the water surface elevation and waves combining refraction-diffraction and breaking. the local water depth. In this paper, a new nonlinear wave equation is developed for modelling wave field covering the entire range from 2.2. Derivation of the Nonlinear Wave Equation. To develop deep to shallow water and involving wave breaking. It starts a nonlinear wave equation, we perform depth-integration on from the depth-integration of the continuity equation and the governing equations (1)β(3) separately and then substitute nonlinear momentum equation. A set of fully dispersive the boundary conditions (4)β(6) into the appropriate resul- and weakly nonlinear wave equations is then developed by tant forms. First, upon integrating (1)verticallyfromseabed introducing a perturbation technique. The theoretical deriva- tothefreesurfaceandsubstitutingtheboundaryconditions tion is outlined in Section 2.InSection 3,weintroducethe given by (5)and(6), it yields finite volume method with a fourth-order Adams-Moulton π predictor-corrector technique, rather than the usual finite ππ +ββ β« π’ππ§=0.β (7) difference and finite element schemes, for directly computing ππ‘ ββ the nonlinear equations. Finally, the present model is verified against Stokes wave and cnoidal wave propagation on uni- Second, after integrating the momentum equation (3)in form depth and nonlinear wave shoaling on a sloping beach. the vertical direction, from an arbitrary depth π§ to the Examples for wave propagation on an elliptic shoal and wave free surface, and inserting the boundary condition for the transformation across the surf zone are also presented. pressure at the free surface, (4), it produces Mathematical Problems in Engineering 3
π π ππ€ π ππ€ =π(πβπ§) + β« ( ) ππ§ + β« (π’β βπ€+π€β ) ππ§. The flow velocity π€ in the vertical direction π§ can be π π§ ππ‘ π§ ππ§ obtained upon depth integration of (1)fromββ to π§ and the (8) bottom boundary condition equation (6), rendering
π (πβ,) π§ Substituting (8)into(2) renders the momentum equation π€(π₯,π¦,π§,π‘)=βββ [ π(π₯,π¦,π‘)],β π for the π₯-andπ¦-axis in the horizontal directions, respectively, (12) as follows: π (β+π§) π (πβ,) π§ = sinh . cosh πβ πβσ³¨π’ ππ’β +πβπ+(π’β ββ π’+π€β ) ππ‘ ππ§ Inserting (11)into(7) gives the following depth-integrated (9) continuity equation: π ππ€ π ππ€ =βββ« ( )ππ§βββ« (π’β βπ€+π€β )ππ§. π§ ππ‘ π§ ππ§ ππ π +ββ [(β« πΉππ§)π]=0.β (13) ππ‘ ββ It is worth noting that the convective terms, (π’β ββ π’+β π€ππ’/ππ§)β and (π’β βπ€+π€ππ€/ππ§)β in (8)and(9), are transformed Substituting (11)and(12)into(8), and after laborious manip- 2 2 to β(π’β β π’+π€β ) and π(π’β β π’+π€β )/ππ§ in Nadaoka et al. [23], ulations, we obtain the pressure expression as based on the irrotational assumption before the integration process described above. However, such transformations π ππβ ππβ =π(πβπ§)βπΜ (β β )βπΜ (ββ β ) do not completely satisfy the irrotationality for variable π π ππ‘ π ππ‘ depth, that is, π’βπ§ ββπ€=0ΜΈ ,whichwasmentionedin[23]. Subsequently, the present equation (9)isretainedinthe β πΜ [(πβ β)(ββ β π)β β (β β π)β (β β π)]β (14) original forms of the nonlinear convective terms; that is, we π preserve, completely, the irrotational condition. It is noted Μ β β 2 that Belibassakis and Athanassoulis [18] preserved this flow β ππ [(β β π) (ββ β π)] + π [(ββ) ], irrotationality by working with a potential formulation; they derived a nonlinear coupled-mode system of equations on the where horizontal plane, with respect to the free surface elevation and the unknown modal amplitudes, in conjunction with Luke π π πΜπ = β« ( )ππ§, variational principle. π§ π The horizontal velocity vector can be expressed in terms of a vertical distribution function and a corresponding π π π πΜπ = β« ( )ππ§, horizontal velocity vector at the still water level, in the form π§ πβ π of (15) π π πΜπ = β« πΉ( )ππ§, π π§ π βσ³¨ cosh ππ (β+π§) β π’(π₯,π¦,π§,π‘)= β ππ (π₯,π¦,π‘), (10) π β π π=1 cosh π π π ππΉ π πΜπ = β« [πΉ ( )β ( )] ππ§. π§ πβ π πβ π β where ππ(π₯, π¦, π‘) is the πth component of the horizontal The last term on the right-hand-side (RHS) of(14)isthe velocity vector at π§=0,andππ(π₯, π¦) is the local wave number. The full dispersivity for a broad-banded wave field higher order term of the bottom slope which may be omitted can be attained by taking only a few terms [23]. However, under the assumption of mild-slope. Multiplying (9) by the vertical distribution function πΉ considering the fact that a narrow-banded sea state centred π§=ββ π around a certain dominant wave frequency may be described in (11) and then integrating it from to ,aswell with sufficient accuracy by the single-component model as after complex mathematical manipulations, produces a equations [24], that is, π=1,itbecomes depth-integrated expression for the momentum equation in the form of
βσ³¨π’(π₯,π¦,π§,π‘)=πΉ(πβ,) π§ π(π₯,π¦,π‘),β ππβ πΌΜ + πΌΜ πβπ + πΌΜ (πβ β)β πβ 1 ππ‘ 2 3 π (β+π§) (11) πΉ (πβ,) π§ = cosh . β β β β cosh πβ + πΌΜ4 (β β π) π+πΌΜ5 (ββ β π) π
ππβ ππβ This single-component model has remained popular and = πΌΜ β(ββ )+πΌΜ (β β )ββ useful in practical coastal engineering applications. 6 ππ‘ 7 ππ‘ 4 Mathematical Problems in Engineering
ππβ + πΌΜ (β β )βπ Substituting (18)into(13)and(16)andtakingTaylorseries 8 ππ‘ expansion on the still water level π§=0,weobtaintheconti- nuity equation and the momentum equation, respectively, as ππβ follows: + πΌΜ9 (ββ β βπ) ππ‘ 0 ππ 3 +ββ [(β« πΉππ§+π)π]=π(πβ ), (19) β β β β ππ‘ ββ + πΌΜ10β[(πβ β)(ββ π) β (β β π) (β β π)] 0 β 0 2 ππ + πΌΜ [(β β π)β (ββ β β) π]β (β« πΉ ππ§) +(β« πΉππ§)πβπ 11 ββ ππ‘ ββ
0 + πΌΜ [(πβ β)(ββ β π)β β (β β π)β (β β π)]β βπ 1 2 2 12 + [ β« πΉ(πΉ βπ )ππ§ 2 ββ + πΌΜ [(β β π)β (ββ β π)β βπ]( +π[ ββ)2] 13 0 π2π4 β2 β« πΉππ§ ]β(πβ β π)β 2 (16) ββ π
0 πΉβπΉ2 ππβ in which =(β« ππ§) β (ββ ) 2 ββ π ππ‘ π π π πΌΜ = β« πΉ2ππ§, πΌΜ = β« πΉππ§, πΌΜ = β« πΉ3ππ§, 0 π 1βπΉ ππβ 1 2 3 +[β« πΉ ( )ππ§](ββ )ββ+π(π3). ββ ββ ββ 2 ββ πβ π ππ‘ π 2 (20) πΌΜ4 =ββ« πΉπ ππ§, ββ These two equations are accurate to the second-order of π. π π πΉ π 2 ππΉ π π Substituting the expressions of and givenin(11)and πΌΜ5 = β« [πΉ βππΉπ ( )] ππ§, πΌΜ6 = β« πΉπΜπππ§, ββ πβ πβ π ββ (12)into(19)and(20), we obtain the final expression of the continuity equation and momentum equation as follows: π 2 πΌΜ7 =2β« πΉπΜπππ§, ππ π ββ +ββ [( +π)π]β = 0, ππ‘ π (21) π ππΜ π ππΜ π π 2 4 πΌΜ8 = β« πΉ ππ§, πΌΜ9 = β« πΉ ππ§, ππβ 1 π π ββ ππ ββ ππ π +πβπ+ (1 β 3 )β(πβ β π)β ππ‘ 2 π2 π (22) 2 πΌΜ10 = β« πΉπΜπππ§, 1 π (1βπ) ππβ ββ = β[ (β β )] π2 π2 ππ‘ π π ππΜ Μ π πΌΜ11 =2β« πΉππππ§, πΌΜ12 = β« πΉ ππ§, in which ββ ββ ππ 2 π π π = tanh πβ, (23) ππΜπ π πΌΜ13 = β« πΉ ππ§. ββ ππ 1 2πβ (17) π= [1+ ]. (24) 2 sinh (2πβ) β β The momentum equation, (16), retains a full nonlinear The coefficient to the third term β(πβ π) on the left- form under the mild-slope condition. As it is difficult to solve hand-side of (22) consists of two parts. The first part, 1/2,is derived from the nonlinear convective terms of the horizontal (16) directly by a numerical scheme; a perturbation technique 2 4 2 is then proposed. The surface displacement and velocity field velocity;andthesecondpart,β3/2(π π /π ), is the nonlinear 2 4 2 are expressed by effect of the vertical velocity. The value of (π π /π ) tends to 1 in deep water, implying that the nonlinear variation of vertical π=ππ +π2π +β β β , velocityisinthesameorderasthehorizontalvelocityandit 1 2 cannot be neglected, whereas it approximates reasonably to (18) π π=πβ πβ +π2πβ +β β β , zero in the very shallow water. It is worth noting that the and 1 2 π can be computed for a prescribed wave period and a local depth based on the dispersion relationship shown in (23). The where π is the dimensionless nonlinearity parameter of wave dispersion characteristics are essentially the same as in [23] amplitude, and the subscriptions π (π = 1, 2, ..) in (18) denote because both models are approximated to the second-order the πth order approximation. of π. Mathematical Problems in Engineering 5
2 Equations (21)and(22) together become the so-called The term of π[(ββ) ] canbeomittedbytheassumption second-order time-dependent mild-slope equation, because of mild-slope bottom. Upon inserting (29)into(28), the theycanbeusedtorepresentthewavefieldsthatareweakly nonlinear equation developed in this study is reduced to nonlinear and fully dispersive and on the basis of the mild- βσ³¨ slopeassumption.Differingfrom[23], the present equations ππ π2 + β(ππ)=0. (30) (21)and(22) can be directly used to calculate the wave trans- ππ‘ π formationofwavefieldinplanebecausetheverticalvelocity has been embedded in the nonlinear convective terms into Equations (25)and(30) are identical to the so-called time- a simpler wave model. The consistency of the equation can dependent mild-slope equation in Copeland [8]. Eliminating β be simplified into several well-established theories, such as flow rate π from (25)and(30) renders the elliptic mild-slope linear mild-slope equation, nonlinear long wave equation, equation in Berkhoff1 [ ]. and Boussinesq equation. The process is demonstrated below. 2.4. Nonlinear Long Wave Equation. For long wave theory in 2.3. Linear Mild-Slope Equation. Upon dropping the nonlin- shallow water, wave celerity (π),groupvelocity(ππ), and flow ear term in (21) for the sake of linear (Airy) wave theory, this rate (πβ)areusuallygivenby equationcanberewrittenas π=π = βπβ, πβ = (β+π) πβ, (31) ππ βσ³¨ π +ββ π=0, (25) ππ‘ where πβ is the horizontal velocity identical to mean flow πβ velocity in shallow water condition. By dropping the nonlin- where is the linear flow rate of the waves defined by ear terms in the vertical direction in (22) and introducing (31), 0 2 a system of nonlinear long wave equations can be written in π β πβ = β« π’ππ§=β πβ. (26) terms of the mean flow velocity π,givenby ββ π ππβ 1 +πβπ+ β(πβ β π)β = 0. (32) On the other hand, (22)canberewrittenforthelinearwaves ππ‘ 2 as For an irrotational flow, the term β(πβ β π)/2β canbereplacedby ππβ 1 π2 (1βπ) ππβ π +πβπ= β[ (β β )] . (πβ β)β πβ, thus rendering the usual expression for the nonlinear ππ‘ π2 π2 ππ‘ (27) long wave equation.
Substituting (26)into(27), then we obtain 2.5. Boussinesq Equations. Boussinesq equations account for the effect of wave nonlinearity but weak dispersion, that the ππβ π2 1 (1βπ) ππβ (1βπ) ππβ + βπ = β [ ββ β βπ2 β ] . wave celerity and group velocity may be approximated by the 2 2 2 ππ‘ π π π ππ‘ π π ππ‘ first two terms in the Taylor series, (28) 1 π2 =πβ(1β π2β2), In the right-hand side of (28), 3 π (33) π 1 2 2 1 (1βπ) ππβ 1 (1βπ) π2π π= =(1β π β ). β[ (β β )] = β β[ ] π 3 π π2 ππ‘ π π2 ππ‘2 2 Substituting (33)into(21)and(22) and replacing π πβ by 1 = β[π2 (1βπ) π] ββ(ββ π)β for the irrotational flow in shallow water condition, π the present nonlinear mild-slope equation gives (1βπ) = πβπ2 ππ β3 π +ββ [(β+π)π]β + β2 (β β π)β = 0, ππ‘ 3 π2 π2 (34) + βπ β β(ππ), ππβ 1 π π +πβπ+ β(πβ β π)β = 0. ππ‘ 2 1 (1βπ) ππβ (1βπ) ππβ β β[ βπ2 β ]=β βπ2 (β β ) π π2π2 ππ‘ ππ2π2 ππ‘ Equations (34) are the standard form of traditional Boussi- nesq equations. +π[(ββ)2] 3. Numerical Methods (1βπ) =β πβπ2 +π[(ββ)2]. π 3.1. Finite Volume Method. A finite volume method is used to (29) calculate the wave heights based on the governing equations 6 Mathematical Problems in Engineering given by (21)and(22). These two equations are first proceeded And the equations of the corrector-step are expressed by using a volume integration over a specific control volume β, as follows: π+1 π Ξπ‘ π+1 π πβ1 πβ2 ππ π2 π =π + (9Ξ¦ + 19Ξ¦ β5Ξ¦ +Ξ¦ ) , (39) β« ( )πβ+β« [β β ( +π)π]β πβ =0, 24 β ππ‘ β π π+1 π {[π΄] π}β ={[π΄] π}β ππβ 1 3π2π4 β« (π )πβ+β« [πβπ + (1 β )β(πβ β π)]β πβ 2 β ππ‘ β 2 π Ξπ‘ + (9Ξ¨π+1 + 19Ξ¨π β5Ξ¨πβ1 +Ξ¨πβ2), 24 1 π2 (1βπ) ππβ = β« [β ( ββ )] πβ. (40) 2 2 β π π ππ‘ (35) where By the divergence theorem, (35) renders ππ 1 π2 1 2 π2 + π2 β ( ) =β β« [πβ (β +π)π]β ππ +β β ββ β Ξ¦=β β {ππΎ β [( +π)π] β ππΎ β [( +π)π] }, ππ‘ β π π β π π π πΎ=1 πΎ πΎ π 1 1 π2 (1βπ) {ππββ ( ) β« πβ [β (β β π)β πΌ]β ππ} π 2 ππ‘ β π2 π2 Ξ¨=β β [π+β β (π)+ πΌββ πββ β (π)β πΌ]β π π π πΎ πΎ πΎ πΎ β πΎ=1 1 =β π β« πβ (ππΌ)β ππ β 1 1 3π2π4 π β [ (1 β )] β 2 π2 1 1 π2π4 β [ (1 β 3 )] β« πβ [(β πβ β π)β πΌ]β ππ 2 β 2 π π 2 π + β Γ β [π+β β (πβ β π)β πΌββ πββ β (πβ β π)β πΌ]β , (36) πΎ πΎ πΎ πΎ πΎ=1 in which subscription π denotes the mean quantities refer- β π πΌβ 2 2 + ring to the centre of the volume , is the control surface, is 1 1 π (1 β π) + πβ π [π΄] π=πβ πββ ( ) β [( ) β (ββ π)β πΌβ the unit vector, and is the normal vector of .Byknowing β π2 π2 πΎ πΎ=1 πΎ the quantities at a starting time step (denotes as πΊπ), the integral form of (36) can be used to calculate the change in 2 β π (1βπ) β πΊπ within a time increment, and thus update the πΊπ values β( ) β (ββ π)β πΌ]β , π2 πΎ of the next time step. πΎ
Regarding the computational mesh, it is envisaged that 2 1 + β thegridsizeshouldbefineenoughtoresolveimportant ββ π=β β [π+β β (π)β β πββ β (π)β ] πΎ πΎ πΎ πΎ characteristic details yet coarse enough to keep storage and β πΎ=1 computational time requirements reasonable. For instance, a (41) fine mesh is used in the area where the wave transformation is significant, such as the region around the structure or near wavebreaking.Fortheefficientandaccuratecomputationfor in which superscript π represents the computed values of π, an irregular geometric boundary, the boundary-fitted mesh πβ, Ξ¦,andΞ¨ at the πth time iteration step, and superscript system is obtained in the present computation according to a π+1denotes the predicted values at the π+1th time method based on Poisson equations [25]. step. The subscript πΎ in (41)representsthetwodirectionsin π₯-π¦ plane, in which superscript β+β denotes the quantities 3.2. Adams-Moulton Predictor-Corrector Schemes. The intheupstreamofcontrol-volumeandβββforthosein +β ββ fourth-order Adams-Moulton predictor-corrector technique the downstream. The section vectors denoted as ππΎ and ππΎ is introduced to compute (36). The discrete equations of the represent opposing faces of the middle volume. The definition predictor-step are given by of variables in the finite volume method is illustrated in Ξπ‘ Figure 1. In the discrete continuity equations, that is, (37) ππ+1 =ππ + (55Ξ¦π β 59Ξ¦πβ1 + 37Ξ¦πβ2 β9Ξ¦πβ3), 24 and (39), the new values can be computed explicitly. But the (37) implicit process is required to obtain the new values for the momentum equations, (38)and(40), because the dispersive π+1 π {[π΄] π}β ={[π΄] π}β term involves the time and spatial derivatives simultaneously. That is, the surface elevation π canbesolvedexplicitly,but Ξπ‘ πβ + (55Ξ¨π β 59Ξ¨πβ1 + 37Ξ¨πβ2 β9Ξ¨πβ3). the velocity is solved by an implicit numerical scheme. The 24 successive over-relaxation iteration method (SOR) is adopted (38) for the implicit scheme. Mathematical Problems in Engineering 7
β β where + + + S2 , U2 , π2
β β β π΄1 ={[(π1) β (π π+1,π+1 + π π+1,π) β β π+1,π+1 + + + S1 , U1 , π1 +(πβ ) β (π β + π β )] 1 π+1,π π+1,πβ1 π+1,π β β β Sβ, Uβ, πβ U, π 1 1 1 β[(πβ ) β (π β + π β ) 1 πβ1,π+1 πβ1,π+1 πβ1,π
+(πβ ) β (π β + π β )] 1 πβ1,π πβ1,πβ1 πβ1,π β β β β β S2 , U2 , π2 +[(πβ ) β (π β + π β ) 2 π+1,π+1 π+1,π+1 π,π+1 Figure 1: Definition sketch of the variables in finite volume method. +(πβ ) β (π β + π β )] 2 π,π+1 πβ1,π+1 π,π+1
β[(πβ ) β (π β + π β ) β β 2 π+1,πβ1 π+1,πβ1 π,πβ1 iβ1, j+1 ( ) ( ) S2 i, j+1 S2 i+1, j+1 i,j+1 i+1,j+1 +(πβ ) β (π β + π β )]} , 2 π,πβ1 πβ1,πβ1 π,πβ1 β ( ) β S1 ( ) π΄ ={[(πβ ) β (π +π ) πΌβ iβ1,j+1 S1 2 1 π+1,π+1 π+1,π+1 π+1,π i+1,j+1 β U, π i, j +(πβ ) β (π +π ) πΌ]β iβ1, j i+1, j 1 π+1,π π+1,πβ1 π+1,π β β β β[(π1) β (ππβ1,π+1 +ππβ1,π) πΌ β (S ) πβ1,π+1 ( ) 1 S1 i+1,j iβ1,j +(πβ ) β (π +π ) πΌ]β 1 πβ1,π πβ1,πβ1 πβ1,π iβ1, jβ1 i+1, jβ1 +[(πβ ) β (π +π ) πΌβ 2 π+1,π+1 π+1,π+1 π,π+1 β i, jβ1 ( ) β S2 ( ) i,jβ1 S2 i+1,jβ1 +(πβ ) β (π +π ) πΌ]β 2 π,π+1 πβ1,π+1 π,π+1
β[(πβ ) β (π +π ) πΌβ Figure 2: Mesh arrangement and corresponding variables around 2 π+1,πβ1 π+1,πβ1 π,πβ1 the centre of a control volume. +(πβ ) β (π +π ) πΌ]}β , 2 π,πβ1 πβ1,πβ1 π,πβ1 β Μ Μ β π΄3 ={[(π1) β (ππ+1,π+1 + ππ+1,π) πΌ As shown in Figure 2, the quantities of variables π and π+1,π+1 πβ are set in the centre of the control volume. Then the flux +(πβ ) β (πΜ + πΜ ) πΌ]β terms, Ξ¦ and Ξ¨,in(37)to(40),canbeexpressedas 1 π+1,π π+1,πβ1 π+1,π β Μ Μ β β[(π1) β (ππβ1,π+1 + ππβ1,π) πΌ 1 πβ1,π+1 Ξ¦π,π =β β π΄1, 2(β +β +β +β ) +(πβ ) β (πΜ + πΜ ) πΌ]β π,π π+1,π π,π+1 π+1,π+1 1 πβ1,π πβ1,πβ1 πβ1,π π Ξ¨π,π =β β π΄2 +[(πβ ) β (πΜ + πΜ ) πΌβ 2 π+1,π+1 π+1,π+1 π,π+1 2(βπ,π +βπ+1,π +βπ,π+1 +βπ+1,π+1) 1 (42) +(πβ ) β (πΜ + πΜ ) πΌ]β β 2 π,π+1 πβ1,π+1 π,π+1 4(βπ,π +βπ+1,π +βπ,π+1 +βπ+1,π+1) β Μ Μ β β[(π2) β (ππ+1,πβ1 + ππ,πβ1) πΌ 3π2π4 π+1,πβ1 Γ(1β ) β π΄ , π2 3 π,π +(πβ ) β (πΜ + πΜ ) πΌ]}β , 2 π,πβ1 πβ1,πβ1 π,πβ1 8 Mathematical Problems in Engineering
π2 Incident waves π β =( π,π +π ) πβ , π,π π π,π π,π
σ³° πΌinc Ξr πΜ = πβ β πβ =π2 +π2 . B B π,π π,π π,π π,π π,π πΌout (43) Outgoing waves 3.3. Initial and Converge Conditions. The condition of zero flow or cold start has been set as the initial condition for most numerical simulations. In this case, an unsteady fluctuation Figure 3: Sketch of the offshore open boundary. will commonly be produced in the initial computing stage, but it will become negligible when the computation time is long enough, and finally a steady state can be obtained. In the explicit iteration for solving the surface elevation, the where the subscriptions inc and out represent the quantities of incoming and outgoing waves, respectively, Ξπ is the length Courant stability condition has to be satisfied at each time σΈ σΈ of π΅π΅ , π΅ is the phantom point of π΅,andπΌ is the angle step, that is, σΈ out between outgoing wave and π΅π΅ ,asshowninFigure 3. β The second type is for the side or lateral boundary of Ξπ‘ β€π π β min (σ΅¨ σ΅¨) (44) σ΅¨πβ β πσ΅¨β the computation domain, which is generally a virtual open σ΅¨ σ΅¨ boundary. It is treated as having no wave reflection at this boundary, hence, allowing waves to pass freely through them. in which ππ is the safe coefficient (commonly taken as 0.6). In This can be represented by the computations, several iteration cycles are usually required toachieveasteadyresult.Thecomputedresultsareregarded π‘ π‘βπ (π) σΈ =(π) , as reaching a steady state when the iterative errors of the π΅ π΅ π‘ π‘βπ (48) computed wave height are less than 2%. As for the velocity (π)β =(π)β . in the horizontal plane, πβ, each time interval in the SOR π΅σΈ π΅ iteration method needs to satisfy the following converge And the third is the internal boundary. If an impermeable condition: structure exists, where full reflection from the structure σ΅¨ β π+1 β π σ΅¨ surface is assumed for the internal boundary associated with σ΅¨π β π σ΅¨ β4 Ξπβσ³¨ = (σ΅¨ σ΅¨)β€10 π max σ΅¨ σ΅¨ (45) wave reflection, it can be set as σ΅¨ πβ π+1 σ΅¨ π‘ π‘ (π) σΈ =(π) , Ξπβσ³¨ π΅ π΅ in which π denotes the absolute maximum error. π‘ π‘ (πβ ) =(πβ ) , ππ‘ π΅σΈ ππ‘ π΅ (49) 3.4. Boundary Specifications. Three types of boundary con- ditions are specified for computing the wave field. These are β π‘ β π‘ (πππ) σΈ =β(πππ) , appliedusingphantompointexteriorfromthecomputational π΅ π΅ mesh. The first one is an offshore boundary which isthe where subscriptions St and Sn represent the tangent and incident wave conditions. In general, there exist not only normal components at the boundary, respectively. incident waves but also outgoing wave components across the offshore boundary due to the reflection from structures or the shore. Referring to [9, 26], such outgoing waves should be 4. Numerical Results transmitted freely through the boundary. Thus the offshore Validation of the present model is made by comparing the boundary is set as numerical results with analytical solutions and/or experi- (π)π‘ =(π )π‘ +(π )π‘βπ mental data for Stokes wave and cnoidal wave propagation π΅σΈ inc π΅σΈ out π΅ on uniform depth, one-dimensional wave shoaling, and two- π‘ π‘ π‘βπ (46) (π)β =(πβ ) +(πβ ) , dimensional wave transformation across an elliptic shoal. π΅σΈ inc π΅σΈ out π΅ Wave transformation involving energy dissipation due to wave breaking is finally simulated for wave propagation on where auniformslopebeachandabar-typebeachandwavearound (πΌ ) ajetty. π=Ξπβ cos out , π 4.1. Stokes Waves. The first case of validation is the numerical (π )π‘ =(π)π‘ β(π )π‘ , (47) out π΅ π΅ inc π΅ simulation of Stokes waves to the second-order. The com- ππ» π‘ π‘ π‘ puted results for two different values in intermediate (πβ ) =(π)β β(πβ ) , β/πΏ = 0.2 π» out π΅ π΅ inc π΅ depth ( ) are depicted in Figure 4,inwhich is Mathematical Problems in Engineering 9
πΏ π =π»/2 kH = 0.40, h/L = 0.2 the wave height, is the wavelength and π .The 1.5 temporal and spatial resolutions are Ξπ‘ = π/60 and Ξπ₯ = πΏ/30, respectively, where π is the wave period. It can be seen 1.0 that the numerical results are in good agreement with the 0.5 o a theoretical solution of the Stokes waves to a second-order. The / 0.0 π minor discrepancy in phase for large π₯/πΏ is possibly caused β0.5 by the truncation error in the discrete numerical process β1.0 which affects the solution of the phase velocity. Nevertheless, β1.5 the verification indicates the present model is capable of 012345 producing identical results to the Stokes waves in water of x/L kH = 0.50, h/L = 0.2 finite depth. 1.5 The vertical distribution of the horizontal velocity can be 1.0 calculated from (10) once the solution of πβ is obtained. The 0.5
maximum nondimensional horizontal velocity calculated o a π§/β / 0.0 from (10) is plotted against dimensionless depth in π Figure 5. Good agreement is found upon comparing the β0.5 results with the experimental data of Kim et al. [27]. β1.0 β1.5 0 12345 4.2. Cnoidal Waves. It is well known that cnoidal wave x/L theory is preferred than Stokes wave theory in shallow water. Figure 6 demonstrates the simulated wave profiles of cnoidal Figure 4: Simulated Stokes waves in intermediate water depth waves in relative depth of β/πΏ = 0.05 for two different π»/β compared with analytical solution (circles: present solution, line: values. The temporal and spatial resolutions are Ξπ‘ = π/40 2nd-order Stokes wave theory). and Ξπ₯ = πΏ/20 in this computation. The present numerical solution is also found in good agreement with the theoretical 0.4 solution of cnoidal waves. 2 In very shallow water region, cnoidal waves approach the h/(gT ) = 0.050 limit of a solitary wave. Since it is not possible to define a H/h = 0.277 finite wavelength for a solitary wave, as it is infinite in theory, 0.2 the simulation of the solitary wave is often made by taking a very small value of β/πΏ.Hereβ/πΏ is taken as 1/50. Figure 7 β/πΏ illustrates the simulated wave profiles in relative depth of 0 =0.02forπ»/β =0.20and0.30,respectively.
From the comparisons above for simulations of nonlinear z/h waves on uniform depth, it reveals that the present model is β0.2 adequate for the entire range of water depth from deep to very shallow. The following simulations will cover nonlinear wave shoaling and two-dimensional wave field involving combined wave refraction-diffraction, even wave breaking. β0.4
4.3. Wave Shoaling on Uniform Slope. Shuto [28]proposedan β0.6 equation for nonlinear wave shoaling prior to breaking along 0 0.1 0.2 0.3 0.4 a uniform slope based on a K-dV formulation. His results U/(gh)0.5 were found to agree reasonably well with the experimental Present solution dataforawiderangeofwavesteepness[29]. The validity Experiment, Kim et al. (1992) of the present nonlinear model is further tested against the results of Shuto [28] by comparing the shoaling coefficients Figure 5: Vertical distribution of the computed maximum horizon- π»/π»π versus relative water depth β/πΏπ for wave steepness talvelocitycomparedwithexperimentaldata. π»π/πΏπ ranging from 0.001 to 0.08, on uniform beach slope of 1/30, in which π»π and πΏπ are the wave height and wavelength in deep water. As seen in Figure 8,theresults computed from the present model (in solid circles) are in the shoaling coefficient increases as wave steepness increases. good agreement with the nonlinear shoaling theory of Shuto However, the solution computed by the linear wave shoaling [28], despite small fluctuations about the simulated curves is independent of the wave steepness π»π/πΏπ which apparently within some range of β/πΏπ.Tsaietal.[10]incorporated underestimates any shoaling wave in shallow water compared a nonlinear shoaling factor into the linear time-dependent with the present nonlinear shoaling model, but has the same mild-slope equations; their results (circles in Figure 8)are value between the nonlinear and linear wave shoaling as β/πΏπ also in good agreement with Shuto [28]. At a particular β/πΏπ, greater than 0.15. 10 Mathematical Problems in Engineering
H/h = 0.10, h/L = 0.05 3.5 2.0 1.5 3.0 Ho/Lo 1.0 0.001
o 0.5 0.002 a
/ 2.5 π 0.0 0.005
β0.5 o 2.0 β1.0 0.01 H/H β1.5 0.02 012345 1.5 x/L 0.04 0.08 H/h = 0.20, h/L = 0.05 1.0 2.0 1.5 0.5 1.0 0.005 0.01 0.02 0.05 0.1 0.2
o 0.5 a h/Lo /
π 0.0 β0.5 Shuto (1988) Tsai et al. (2001) β1.0 Linear solution Present solution β1.5 012345 Figure 8: Comparison on the wave height from the present nonlinear wave shoaling solution with other solutions. x/L
Figure 6: Simulated wave profiles for cnoidal waves compared favourably with analytical solution (circles: present solution, line: better resolve the wave transformation in the computation. Cnoidal wave theory). The simulation results for wave height distribution behind the shoal (seen in Figure 10) exhibit significant contribution of
ao/ho = 0.20, h/L = 0.02 the combined effect of refraction-diffraction. 0.5 The simulation results for wave height distribution along 0.4 the eight transects indicated in Figure 11 are compared with 0.3 the experimental data and the numerical model of Li and o h
/ 0.2 Anastasiou [3], who derived an improved solution using an π 0.1 elliptic form of the linear mild-slope equation. The com- 0.0 parison in Figure 11 demonstrates that the present nonlinear model can yield the wave height oscillations spatially dis- β0.1 012345tributed due to combined refraction-diffraction effect behind π΅ πΈ x/L an elliptical shoal, particularly on transects to . a /h = 0.30, h/L = 0.02 0.5 o o 4.5. 1D and 2D Wave Transformation across Surf Zone. More- 0.4 over, the momentum equation used in the present model can 0.3 be applied to simulate wave transformation across the surf o h / 0.2 zone, upon incorporating a mixing term, ππ·.Inthisway, π 0.1 (22)becomes 0.0 ππβ π 1 π2π4 β0.1 + βπ + (1 β 3 )β(πβ β π)β 012345 ππ‘ π 2π π2 x/L (50) 1 π2 (1βπ) ππβ Figure7:Simulatedwaveprofileforsolitarywavesinveryshallow = β[ (β β )] + π . ππ2 π2 ππ‘ π· water compared with analytical solution (circles: present solution, line: solitary wave theory). Following Kabiling and Sato [31], the mixing term due to energy dissipation can be expressed by
4.4. Wave Transformation across an Elliptic Shoal. The 2 π = ] β π,β (51) hydraulic model test [30] for wave propagation over a 3D π· π elliptical shoal on a plane beach with uniform slope 1 : 50 has where ]π is the eddy viscosity related to an energy dissipation usuallybeenemployedforwavemodelverification.Thewave factor, ππ·,duetobreaking, height distribution was measured along eight cross sections indicated in Figure 9.Anincidentwavewithaperiodof πβ π π βπ 1.0 sec and wave height of 2.32 cm was used in the simulation. ] = π ;π=πΌ π½β ( π π ). π 2 π· π· π· tan (52) The finer rectangular grids are set in the shoal region to π β ππ βππ Mathematical Problems in Engineering 11
Outside shoal Shoal boundary Shoal thickness h = 0.45 β 0.02(5.84 β yσ³°) (xσ³°/4)2 +(yσ³°/3)2 =1 β0.3 + 0.5[1 β (xσ³°/5)2 β(yσ³°/3.75)2]0.5 10
5 Wave Shoal xσ³° T = 1.0 s Hi = 2.32 cm yσ³° H
(m) 0 G y
F
ABC DE β5 m 5
0.4 h = 0.4 0.3 0.2 0.15 0.25 0.35 0.1 0.05
β10 β10 β50 51015 x (m)
Figure 9: Geometry showing wave propagation over an elliptic shoal.
T=1s, H = 2.32 cm In (52), ππ· is the energy dissipation factor, π (= 2π/π) 10 i is the wave frequency, tanπ½ is an average bottom slope 2 πΌ 2 from shoreline to wave breaking point, D is a modification coefficient (usually taken as 2.5 as suggested in[32]), ππ is the 3 2 π π 5 2 amplitude of flow rate, and π and π are determined from 2 2 the following equations: 1 2 1 2 3 π = 0.4 (0.57 + 5.3 π½) βπβ3, 3 π tan (m) 0 4 2
y 5 1 2 3 (53) 3 1 2 π = 0.135βπβ . 3 π 2 β5 2 1 In (52), the dissipation factor ππ· is set equal to zero outside 2 3 the surf zone and in any region where ππ <ππ. On the basis of shallow water assumption for wave 2 β 2 β β10 breaking, we can simplify β π=βππ and render ππ· = β10 β50 51015 β βππ·π,thus(50)becomes x (m) ππβ π 1 π2π4 H/H + βπ + (1 β 3 )β(πβ β π)β + π πβ i ππ‘ π 2π π2 π· (7) 2.4 (3) 1.2 (6) 2.1 (2) 0.9 (54) 1 π2 (1βπ) ππβ (5) 1.8 (1) 0.6 = β[ (β β )] . (4) 1.5 ππ2 π2 ππ‘
Figure 10: Simulation results on wave height distribution behind an This modified momentum equation can now be used for the elliptic shoal. wave transformation after wave breaking. To estimate the 12 Mathematical Problems in Engineering
3.0 3.0 ) ) i i 2.5 2.5 H/H H/H 2.0 2.0 1.5 1.5 1.0 1.0 0.5 0.5 Relative wave height ( height wave Relative Relative wave height ( height wave Relative 0.0 0.0 β5 β4 β3 β2 β1 0 1 2 3 4 5 β5β4β3β2β1012345 Distance y (m) Distance y (m) (a) (b) 3.0 3.0 ) ) i i 2.5 2.5 H/H H/H 2.0 2.0 1.5 1.5 1.0 1.0 0.5 0.5 Relative wave height ( height wave Relative Relative wave height ( height wave Relative 0.0 0.0 β5β4β3β2β1012345 β5β4β3β2β1012345 Distance y (m) Distance y (m) (c) (d) 3.0 3.0 ) ) i i 2.5 2.5 H/H H/H 2.0 2.0 1.5 1.5 1.0 1.0 0.5 0.5
Relative wave height ( height wave Relative 0.0 ( height wave Relative 0.0 β5 β4 β3 β2 β1 0 1 2 3 4 5 012345678910 11 Distance y (m) Distance x (m) (e) (f) 3.0
3.0 ) ) i i 2.5 2.5 H/H H/H 2.0 2.0 1.5 1.5 1.0 1.0 0.5 0.5 Relative wave height ( height wave Relative Relative wave height ( height wave Relative 0.0 0.0 012345678910 11 012345678910 11 Distance x (m) Distance x (m)
Experiment, Berkhoff 1982( ) Experiment, Berkhoff 1982( ) Li and Anastasiou (1992) Li and Anastasiou (1992) Present solution Present solution (g) (h)
Figure 11: Simulation results of wave height distribution along eight transects shown in Figure 9 andcomparisonwithexperimentaldataand other numerical models. Mathematical Problems in Engineering 13
T = 1.19 s Ho = 2.0 cm H = 6.00 cm 5.0 T = 1.2 s o πΌ =60β Incident waves 10 o 4.0
(m) 3.0
x 2.0 5
(cm) 2.0 3.0 H 2.5 1.5 2.0 2.5 2.5 2.5 2.5 1.0 2.0 1.5 1.5 1.0 1.0 0 0.5 0.5 0.0 1 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 20 y (m) β5 Figure 14: The computed wave field in the vicinity of a jetty. 0123456 x (m)
Experiment, Nagayama (1983) Tsai et al. (2001) instead of (22), is then used for calculating the wave height Watanabe and Maruyama (1986) Present solution transformation after wave breaking. Figure 12: Comparisons among numerical solutions and laboratory Numerical verification on one-dimensional wave trans- experiment for wave transformation on uniform slope. formation across the surf zone is now considered for two different bottom configurations, a uniform slope in 1/20 and a composite profile with a bar in 1/20 slope, as shown T = 0.94 s in Figures 12 and 13, respectively. The numerical results of Ho = 7.00 cm the present nonlinear model are compared with the linear 10 solution of wave transformation [9] and the experimental data of Nagayama [34].TheresultsofTsaietal.[10]arealso included in these figures. The local wave heights calculated using the present nonlinear model are in better agreement 5 with the experimental data than the linear model. In Figure 12, the computed wave heights increases as it (cm) shoals, which ultimately approaches a maximum at breaking, H then decreases rapidly until reaching the shoreline. For the 0 resultantwaveheighttransformationonacompositeslope, 1 1 20 as in Figure 13,waveheightreachesamaximumpriorto 1 20 20 the bar crest where first breaking occurs on the leading β5 uniform slope of 1/20. Wave height recovers soon after it 01234567 passesthebartroughconnectedtothetrailinguniformslope x (m) on which the second breaking appears. Finally, the surviving Experiment, Nagayama (1983) Tsai et al. (2001) wave diminishes at the shoreline. For the region before the Watanabe and Maruyama (1986) Present solution breaking on a uniform slope (Figure 12)andpriortothe first breaking on a composite slope (Figure 13), the present Figure 13: Comparisons among numerical solutions and laboratory nonlinear solution and linear model agree with each other, experiment for wave transformation on a composite profile of bar otherwise the linear model is found to deviate from the type. experimental results. This confirms that the present nonlinear model is adequate for wave decay, recovery, and secondary breaking on a composite beach profile. location of wave breaking, a breaking criterion proposed by The example of two-dimensional wave transformation Goda [33]isusedasfollows: was conducted for wave field in the vicinity of a jetty, for which experimental and numerical results by Watanabe and π»π βπ 4/3 Maruyama [9] were available. The jetty in the experimental = 0.17 {1βexp [β1.5π (1+15tan π½)]} (55) πΏπ πΏπ tank was placed perpendicular to the shoreline on a beach with a slope of 1/50. The model jetty was extended to an in which π»π and βπ are the wave height and water depth offshore distance of 4 m, where the water depth was 8 cm. The β of the breaking wave, respectively, and πΏπ is the deep water waves were incident with angle 60 to the contour lines, and wavelength. The computed wave height at any location based the wave height and period in deep water were 2.0 cm and on (21)and(22)isalwayscomparedwiththevaluecalculated 1.2 s, respectively. from the right-hand side of (55). If the computed wave height The computed wave height distribution around the jetty is reaches the breaking criterion, the momentum equation (54), shown in Figure 14.Aclapotis gaufreΜ wave system is visible in 14 Mathematical Problems in Engineering
5 The resultant nonlinear mild-slope equation is computed y=3m by a finite volume method with a fourth-order Adams- 4 Moulton predictor-corrector technique, rather than the usual 3 finite difference or finite element scheme. The validity ofthe
(cm) present nonlinear model is also satisfactorily verified with H 2 examples that cover a wide range of practical problems, such 1 as Stokes wave to a second-order, cnoidal wave propagation on uniform depth, nonlinear wave shoaling on a sloping 0 012345 6beach,andwavetransformationonanellipticshoal.By x (m) introducing an extra term representing energy dissipation 5 duetowavebreakingtothemomentumequation,thepresent y=4.8m nonlinear model is suitable for the wave transformation 4 withinthesurfzoneonauniformslopeandonacomposite 3 profile with a bar on which second wave breaking occurs after
(cm) wave recovering. The comparisons of computed wave height H 2 distributions around a jetty indicate the present nonlinear 1 solutions are in very good agreement with the experimental results. 0 012345 6 x (m) Conflict of Interests 5 y=5.2m The authors declare that there is no conflict of interests 4 regarding the publication of this paper. 3 (cm)
H 2 Acknowledgment 1 Thismaterialwasbasedupontheworksupportedbythe 0 Ministry of Science and Technology, Taiwan, under Grant no. 012345 6NSC 99-2221-E-005-116-MY3. x (m)
Measured (Watanabe and Maruyama 1986) References Linear solution (Watanabe and Maruyama 1986) Present solution [1] J. C. W. Berkhoff, βComputation of combined refraction- diffraction,βin Proceedings of 13th International Confererence on Figure 15: Comparisons among numerical solutions and experi- Coastal Engineering (ASCE β72),vol.1,pp.10β1972,1972. mental results for distributions of wave height at sections π¦ = 3.0 m, [2] V. G. Panchang, B. R. Pearce, G. Wei, and B. Cushman-Roisin, 4.8 m, and 5.2 m shown in Figure 14. βSolution of the mild-slope wave problem by iteration,β Applied Ocean Research,vol.13,no.4,pp.187β199,1991. [3] B. Li and K. Anastasiou, βEfficient elliptic solvers for the mild-slope equation using the multigrid technique,β Coastal frontofthejettybythereflectionfromit.Ontheotherhand, Engineering,vol.16,no.3,pp.245β266,1992. the wave diffraction was observed in the lee of the jetty. The [4] Y. Zhao and K. Anastasiou, βModelling of wave propaga- comparisons of the wave height distributions at cross-shore tion in the nearshore region using the mild slope equation sections π¦ = 3.0 m, 4.8 m, and 5.2 m are shown in Figure 15, with GMRES-based iterative solvers,β International Journal for which indicates the present nonlinear solutions are in better Numerical Methods in Fluids,vol.23,no.4,pp.397β411,1996. agreement with the experimental results. [5]F.S.B.F.OliveiraandK.Anastasiou,βAnefficientcomputa- tional model for water wave propagation in coastal regions,β Applied Ocean Research, vol. 20, no. 5, pp. 263β271, 1998. 5. Conclusions [6] A. C. Radder, βOn the parabolic equation method for water- wave propagation,β Journal of Fluid Mechanics,vol.95,no.1,pp. In this paper, the development of a second-order time- 159β176, 1979. dependent mild-slope equation is reported, based on the [7] P. L.-F. Liu and T.-K. Tsay, βOn weak reflection of water waves,β fundamental principles of conservation of mass and momen- Journal of Fluid Mechanics,vol.131,pp.59β72,1983. tum, together with a depth-integration procedure. After some [8]G.J.M.Copeland,βApracticalalternativetotheβmild-slopeβ appropriate mathematical simplifications, the present nonlin- wave equation,β Coastal Engineering,vol.9,no.2,pp.125β149, ear equation can yield the usual form of linear mild-slope 1985. equation, nonlinear long wave equation, and Boussinesq [9]A.WatanabeandK.Maruyama,βNumericalmodelingof equation, respectively. We have also demonstrated that the nearshore wave field under combined refraction, diffraction present nonlinear model is applicable to the entire range of and breaking,β Coastal Engineering Journal,vol.29,pp.19β39, water depth from deep to very shallow, even wave breaking. 1986. Mathematical Problems in Engineering 15
[10] C.-P. Tsai, H.-B. Chen, and J. R.-C. Hsu, βCalculations of wave [30] J. C. W. Berkhoff, βRefraction and diffraction of water waves; transformation across the surf zone,β Ocean Engineering,vol.28, wave deformation by a shoal, comparison between computa- no. 8, pp. 941β955, 2001. tions and measurements,βReport on Mathematics Investigation [11] P. G. Chamberlain and D. Porter, βThe modified mild-slope W 154-III, Delft Hydraulics Laboratory, Delft, The Netherland, equation,β Journal of Fluid Mechanics,vol.291,pp.393β407, 1982. 1995. [31] M. B. Kabiling and S. Sato, βTwo-dimensional nonlinear disper- [12] T.-W. Hsu, T.-Y. Lin, C.-C. Wen, and S.-H. Ou, βA comple- sive wave-current and three-dimensional beach deformation mentary mild-slope equation derived using higher-order depth model,β Coastal Engineering in Japan,vol.36,pp.195β212,1993. function for waves obliquely propagating on sloping bottom,β [32] A. Watanabe and M. Dibajnia, βNumerical model of wave Physics of Fluids,vol.18,no.8,ArticleID087106,2006. deformation in surf zone,β in Proceedings of 21st International [13]J.W.KimandK.J.Bai,βAnewcomplementarymild-slope Conference on Coastal Engineering (ASCE β88),pp.578β587,June equation,β Journal of Fluid Mechanics,vol.511,pp.25β40,2004. 1988. [14] Y. Toledo and Y. Agnon, βThree dimensional application of the [33] Y. Goda, βIrregular wave deformation in the surf zone,β Coastal complementary mild-slope equation,β Coastal Engineering,vol. Engineering in Japan,vol.18,pp.13β26,1975. 58,no.1,pp.1β8,2011. [34] S. Nagayama, Study on the change of wave height and energy in [15] D. K. P.Yue and C. C. Mei, βForward diffraction of Stokes waves thesurfzone[Bachelorthesis], Yokohama National University, by a thin wedge,β Journal of Fluid Mechanics,vol.99,no.1,pp. Yokohama, Japan, 1983, (Japanese). 33β52, 1980. [16] P. L.-F. Liu and T.-K. Tsay, βRefraction-diffraction model for weakly nonlinear water waves,β Journal of Fluid Mechanics,vol. 141, pp. 265β274, 1984. [17]K.A.BelibassakisandG.A.Athanassoulis,βExtensionof second-order Stokes theory to variable bathymetry,β Journal of Fluid Mechanics,vol.464,pp.35β80,2002. [18] K. A. Belibassakis and G. A. Athanassoulis, βA coupled-mode system with application to nonlinear water waves propagating in finite water depth and in variable bathymetry regions,β Coastal Engineering,vol.58,no.4,pp.337β350,2011. [19] D. H. Peregrine, βLong waves on a beach,β Journal of Fluid Mechanics,vol.27,pp.815β827,1967. [20] P. A. Madseon and H. A. Schaffer, βA review of Boussinesq- type equations for surface gravity waves,β in Advanced Coastal and Ocean Engineering, P. L.-F. Liu, Ed., vol. 5, Word Scientific, Singapore, 1999. [21]P.A.Madsen,H.B.Bingham,andH.Liu,βAnewBoussinesq method for fully nonlinear waves from shallow to deep water,β Journal of Fluid Mechanics,vol.462,pp.1β30,2002. [22] B. Li, βA mathematical model for weakly nonlinear water wave propagation,β Wave Motion,vol.47,no.5,pp.265β278,2010. [23] K. Nadaoka, S. Beji, and Y. Nakagawa, βAfully dispersive weakly nonlinearmodelforwaterwaves,βProceedings of the Royal Society A,vol.453,no.1957,pp.303β318,1997. [24] S. Beji and K. Nadaoka, βA time-dependent nonlinear mild- slope equation for water waves,β Proceedings of the Royal Society A,vol.453,no.1957,pp.319β332,1997. [25] J. F. Thompson, F. C. Thames, and C. Wayne Mastin, βTOMCATβa code for numerical generation of boundary- fitted curvilinear coordinate systems on fields containing any number of arbitrary two-dimensional bodies,β JournalofCom- putational Physics,vol.24,no.3,pp.274β302,1977. [26] H.-B. Chen and C.-P. Tsai, βComputations of nearshore wave transformation using finite-volume method,β Applied Ocean Research,vol.38,pp.32β39,2012. [27] C. H. Kim, R. E. Randall, S. Y.Boo, and M. J. Krafft, βKinematics of 2-D transient water waves using laser Doppler anemometry,β JournalofWaterway,Port,Coastal,&OceanEngineering,vol. 118, no. 2, pp. 147β165, 1992. [28] N. Shuto, βNonlinear long waves in a channel of variable section,β Coastal Engineering in Japan,vol.17,pp.1β12,1974. [29] K. Horikawa, Nearshore Dynamics and Coastal Processes, Uni- versity of Tokyo Press, 1988. Copyright of Mathematical Problems in Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.