Numerical Model of Wave Effects on Permeable Vertical Barriers above the Seabed By Gang Yang B.E . , Dalian University of Technology, P.R.China, 1992 M . A . S c , Dalian University of Technology, P.R.China, 1995 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENT FOR THE DEGREE OF Master of Applied science in The Faculty of Graduate Studies Department of Civil Engineering We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA December, 1996 © Gang Yang, 1996 In presenting this thesis in partial fulfillment of the requirements for an advance degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of the department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Civil Engineering, The University of British Columbia, 2324 Main Mall, Vancouver, B.C., Canada, V6T 1Z4 Date: 3>e&. 4 , WS 11 Abstract This thesis describes a numerical model used to predict the interaction of a regular small amplitude wave train with a thin permeable vertical barrier extending from the water surface to some distance above the seabed. The case of wave interactions with a pair of such barriers is also treated. The approach used is based on an eigenfunction expansion method and utilizes a boundary condition at the barrier surface which accounts for energy dissipation within the barrier. Comparisons of results based on the method have been carried out with those of previous numerical studies for related situations, and close agreement has been obtained in all such cases. A selection of results based on the method are presented for the transmission, reflection, and energy dissipation coefficients, the wave runup, and the maximum horizontal force on the barrier. These exhibit various features of interest which are discussed. The numerical model is written in the FORTRAN language. With appropriate input, it can be used to calculate wave interactions with one or two barriers under a wide range of conditions. It is confirmed that the results of the numerical model are consistent for various limiting conditions including those of a solid barrier, a fully transparent barrier, and a barrier extending to the seabed. Overall, the numerical model is found to provide a reasonably flexible and reliable means of predicting wave effects on thin permeable vertical barriers. iii Table of Contents page ABSTRACT ii TABLE OF CONTENTS iii LIST OF FIGURES v LIST OF SYMBOLS viii ACKNOWLEDGMENTS xii CHAPTER 1 INTRODUCTION 1 1.1 General 1 1.2 Literature review 1 1.2.1 Wave effects on impermeable barriers 2 1.2.2 Wave effects on permeable structures 3 1.3 Scope of present work 6 CHAPTER 2 THEORETICAL DEVELOPMENT 7 2.1 Assumptions 7 2.2 Governing equations '. 8 2.3 Flow within barrier 11 2.4 Solution by eigenfunction expansion method 14 2.4.1 Eigenfunction expansion method 14 2.4.2 Transmission and reflection coefficients 17 2.4.3 Energy balance 18 2.4.4 Wave runup 20 2.4.5 Wave loads 21 IV 2.4.6 Limiting conditions 22 2.5 Extension to the two-barrier problem 23 2.5.1 Governing equations and boundary conditions 23 2.5.2 Solutions by eigenfunction expansion method 25 CHAPTER 3 DEVELOPMENT OF NUMERICAL MODEL 33 3.1 General 33 3.2 Parameter selection 33 CHAPTER 4 RESULTS AND DISCUSSION 35 4.1 One-barrier problem 35 4.1.1 Comparisons with previous predictions 35 4.1.2 Transmission and reflection coefficients 36 4.1.3 Energy dissipation coefficient 37 4.1.4 Wave loads 37 4.1.5 Wave runup 38 4.2 Two-barrier problem 39 4.2.1 Transmission and reflection coefficients 39 4.2.2 Energy dissipation coefficient 40 4.2.3 Wave forces 41 CHAPTER 5 CONCLUSIONS AND RECOMMENDATIONS 42 5.1 Conclusions 42 5.2 Recommendations for further study 43 REFERENCES 45 V List of Figures Figure 1 Definition sketch of one-barrier problem. Figure 2 Definition sketch of two-barrier problem. Figure 3 Transmission and reflection coefficients shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). Figure 4 Transmission coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 5 Reflection coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 6 Transmission coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 7 Reflection coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 8 Energy dissipation coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). Figure 9 Energy dissipation coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 10 Energy dissipation coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 11 Force coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). Figure 12 Force coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 13 Force coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 14 Overturning moment coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). vi Figure 15 Overturning moment coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 16 Overturning moment coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 17 Wave runup coefficients shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). Figure 18 Wave runup (upwave) coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 19 Wave runup (downwave) coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). Figure 20 Wave runup (upwave) coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 21 Wave runup (downwave) coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 22 Transmission and reflection coefficients shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5, X/d- 1.0 (two-barrier problem). Figure 23 Transmission coefficient as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S (two-barrier problem). Figure 24 Reflection coefficient as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S (two-barrier problem). Figure 25 Transmission coefficient as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). Figure 26 Reflection coefficient as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). Figure 27 Energy dissipation coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5, X/d = 1.0 (two-barrier problem). Figure 28 Energy dissipation coefficient as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S (two-barrier problem). Figure 29 Energy dissipation coefficient as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). Figure 30 Force coefficients of the upwave barrier shown as contours in the complex S plane for the h/d = 0.5, kh = 0.5, X/d = 1.0 (two-barrier problem). Figure 31 Force coefficients of the downwave barrier shown as contours in the complex S plane for the case h/d = 0.5, kh = 05, X/d = 1.0 (two-barrier problem). Figure 32 Force coefficient of the upwave barrier as a function of kh for h/d = 0.5, X/d - 1.0 and for various values of S (two-barrier problem). Figure 33 Force coefficient of the downwave barrier as a function of kh for h/d = 0.5, X/d — 1.0 and for various values of S (two-barrier problem). Figure 34 Force coefficient of the upwave barrier as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). Figure 35 Force coefficients of the downwave barrier as function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). List of Symbols The following symbols are used in this thesis: a = barrier height, see Fig. 1 A + = see Eqs. 2.63 and 2.64 A- = see Eqs. 2.63 and 2.64 Aim = unknown coefficients, see Eqs. 2.92-2.94 A m = unknown coefficients, see Eq. 2.39 b barrier thickness, see Fig. 1 b n = see Eqs. 2.45 and 2.47 Bm = unknown coefficients, see Eq. 2.40 C = wave speed C G = group speed C = constant, see Eq. 2.2 Cm = added mass coefficient Cmn = matrix coefficients, see Eqs. 2.45 and 2.46 c F = dimensionless maximum horizontal force coefficient, see Eq. 3.1 CM = dimensionless maximum overturning moment coefficient, see Eq. 3.2 C R = complex reflection coefficient C T = complex transmission coefficient d = still water depth E = energy flux f linearized resistance coefficient f see Eq. 2.48 F maximum horizontal wave force F s = maximum horizontal force of standing waves on a solid barrier, see Eq. 2.79 g gravitational constant ix G = complex porous-wall-effect parameter, see Eq. 2.35 G r = real part of G Gj = imaginary part of G h = draft, see Fig. 1 H = wave height i = V ^ I K r = reflection coefficient K t = transmission coefficient Kd = energy dissipation coefficient k = wave number m = evanescent mode M = maximum overturning moment M s = maximum overturning moment of standing waves on a solid barrier, see Eq. 2.80 R = wave runup T = wave period u = horizontal fluid velocity u = spatial velocity vector U = general velocity vector p = spatial pressure P = general pressure S = S'd, dimensionless complex barrier permeability parameter S r = real part of S S; = imaginary part of S S' = G/b, complex barrier permeability parameter t = time x = horizontal coordinate, see Figs. 1 and 2 y = horizontal coordinate along the barrier axis X z = vertical coordinate measured upwards from the seabed Ap = pressure difference across barrier 8 l m = see Eq. 2.139 5 2 m = see Eq. 2.142 e = porosity of barrier r) = wave elevation X = half distance between two barriers, see Fig. 2 (I rn = evanescent mode wave numbers for m > 1 p = water density 0 = incident wave angle O = velocity potential, see Eq. 2.1 <|) = potential function, see Eq. 2.1 <|>i = upwave potential function, see Eqs. 2.39 and 2.92 (J>2 = potential function, downwave in case of one barrier, intermediate in case of two barriers, see Eqs. 2.40 and 2.93 § 3 = downwave potential function in case of two barriers, see Eq. 2.94 <|)w = incident wave potential function, see Eq. 2.7 co = wave angular frequency Q = dimensionless porous impedance, see Eq. 2.24 The following subscripts are used in this thesis + = downwave side (one-barrier case) = upwave side (one-barrier case) 0+ = within the barrier at x = 0 0- = outside the barrier at x = 0 1+ = downwave side of the upwave barrier (two-barrier case) 1- = upwave side of the upwave barrier (two-barrier case) 2+ = downwave side of the downwave barrier (two-barrier case) xi 2- = upwave side of the downwave barrier (two-barrier case) 1 = fluid region number or the upwave barrier, see Figs. 1 and 2 2 = fluid region number or the downwave barrier, see Figs. 1 and 2 3 = fluid region number in case of two barriers, see Fig. 2 b+ = outside the barrier at x = b b- = within the barrier at x = b r = reflected wave s = standing wave condition for solid barrier t = transmitted wave w = incident wave Xll Acknowledgments The author would like to thank his supervisor, Dr. Michael Isaacson, for his guidance throughout the preparation of this thesis. The author would like to express his gratitude to Prof. Dahong Qiu and Prof. Yucheng Li of the Dalian University of Technology for their brilliant teaching and help. Finally, the author is grateful to Mr. Shankar S. Bhat and Mr. Sundaralingam Premasiri for their kind help and encouragement. Chapter 1 Introduction Chapter 1 INTRODUCTION 1 1.1 General Breakwaters have been used extensively to provide economical protection from waves at many coastal locations. Generally, floating breakwaters are effective only for more modest wave conditions, whereas rubble mound structures are widely used for relatively severe waves. For intermediate wave conditions, floating structures may not be sufficiently effective or robust, whereas rubble mound structures may be expensive and may create environmental problems associated with restricted water circulation. As an alternative to these two categories of breakwater, breakwaters in the form of thin, rigid vertical barriers have been considered or adopted, especially for marinas and small fishing harbours. A vertical barrier extending from the water surface to some distance above the seabed has the advantages of allowing water circulation, fish passage and sediment transport beneath the breakwater. Such barriers may be impermeable or porous. Parameters such as reflection and transmission coefficients, wave runup, wave force and overturning moment, are very important in the design of such barriers. However, very little is known about wave effects on permeable barriers above the seabed. Thus, wave interactions with thin vertical porous barriers is the subject of this thesis. The author hopes to extend knowledge on this topic by solving the corresponding linearized boundary value problem in order to obtain various parameters of engineering interest. 1.2 Literature review Commonly, two approaches are used to predict wave interactions with coastal and offshore structures. One is based on the use of Morison equation, originally proposed by Morison, Johnson, O'Brien and Schaaf (1950). The application of this method carries the implicit assumption that the structure size is small relative to the wave length so that the incident flow is virtually uniform in the vicinity of the structure. When the structure is relatively large, this is no longer true and an alternative formulation is required. However, since flow separation is now usually unimportant, a theoretical approach may be used in which the corresponding potential flow is treated taking account of boundary conditions introduced by the presence of the structure. The wave forces predicted by solving this diffraction problem to a first approximation were first obtained by MacCamy and Fuchs (1954) for the special case of a large vertical circular cylinder. Within the last few decades, structures of more general shape have been treated by numerical schemes, and this approach generally forms the basis of present prediction methods for large offshore structures. For certain restricted geometric conditions, the eigenfunction expansion method is widely used as an alternative numerical approach to diffraction problems. Although, problems relating to wave effects on solid barriers and permeable structures have been the subject of numerous theoretical and experimental studies for many years, numerical predictions of wave effects on permeable barriers have been given attention only more recently. In the following sections, a brief review of available theoretical approaches and experimental investigations relating to wave effects on solid and permeable barriers is presented. 1.2.1 Wave effects on impermeable barriers The prediction of wave effects on fixed, impermeable vertical barriers on the basis of linear wave diffraction theory has been considered by a number of authors. Ursell (1947) presented an early analytical solution for the transmission coefficient of a barrier extending down from the water surface in deep water waves. For the case of finite depth, Liu and Abbaspour (1982) developed a numerical solution to the linear diffraction problem using a boundary element method. Losada et al. (1992) and Abul-Azm (1993) developed a numerical solution based on the eigenfunction expansion method, and Losada et al. (1993) extended this to the case of random waves. Chapter 1 Introduction 3 Wiegel (1960) had earlier proposed a simple prediction of the transmission coefficient based on the hypothesis that the wave energy transmitted past the barrier corresponds to the proportion of incident energy flux associated with the gap beneath the barrier. However, Losada et al. (1992) and Abul-Azm (1993) found notable differences between this approximation and their own accurate solutions. Recently, Isaacson (1996) developed a variant of an eigenfunction method to provide accurate results over a comprehensive range of conditions, and used this solution to develop simplified expressions for various parameters of engineering interest. These enable such parameters to be estimated with relatively little effort. 1.2.2 Wave effects on permeable structures Wave interactions with permeable structures have been studied theoretically by a number of authors. These structures include rubble mound structures, wave absorbers, as well as vertical permeable barriers. Le Mehaute (1971) developed a theory for progressive wave absorbers. The basic principle of the theory is that the energy absorbed by the absorbers is proportional to the area of the plates in contact with water—the plates cause partial reflection of the incident wave due to the sudden narrowing the channel. Experimental results demonstrated qualitatively the validity of this theory. Sollitt and Cross (1972) derived a theory to predict wave reflection and transmission at a permeable breakwater of rectangular cross section extending from the seabed to the water surface. The theory solves for a damped wave component within the breakwater and matches boundary conditions at the upwave and downwave breakwater faces to predict the reflected and transmitted wave components. Madsen (1974) derived a simple solution for the reflection and transmission coefficients for a rectangular homogeneous porous structure based on the assumption of Chapter 1 Introduction 4 relative long normally incident waves. The solution was based on the linearized form of the governing equations and a linearized form of the flow resistance formula governing the flow within the porous structure. Macaskill (1978) considered the linearized problem of wave reflection by a thin barrier of arbitrary permeability with the restriction that the flow be two-dimensional. The formulation included the special case of transmission through one or more gaps in an otherwise impermeable barrier. He reduced the general problem to a set of integral equations using standard techniques. A representative range of solutions was obtained numerically for both the finite and infinite depth problems. Chwang (1982) proposed a porous-wavemaker theory to analyze small-amplitude surface waves on water of finite depth produced by horizontal oscillations of a porous vertical plate. Closed-form solutions were obtained for the surface-wave profile, the hydrodynamic-pressure distribution and the total force on the wavemaker. The influence of a wave-effect parameter and a porous-effect parameter, both dimensionless, on the surface waves and on the hydrodynamic pressure was discussed in detail. Madsen (1983) presented a theoretical solution for the reflection of linear shallow-water waves from a vertical porous wave absorber on a horizontal bottom. Periodic solutions were matched at the front face of the absorber by assuming continuity of pressure and mass flux. The friction term describing the energy loss inside the absorber was linearized by using Lorentz principle of equivalent work, and, the reflection coefficient was determined as function of parameters describing the incoming waves and the absorber characteristics. Dalrymple et al. (1990) studied wave reflection and transmission from a porous breakwater using eigenfunction expansions for the water regions in front of the structure, within the porous medium, and behind the structure. Chapter 1 Introduction 5 Twu and Lin (1990) investigated a highly effective wave absorber containing a finite number of porous plates with various porous-effect parameters. The relationships of the spacing between the adjacent porous plates as well as the alignment of the plates were analyzed and a formula for evaluating the reflection coefficient was derived. Also, experiments were carried out and the results agree well with the theoretical solutions. Mallayachari and Sundar (1992) set up a numerical model based on the methodology proposed by Ijima et al. (1976) to determine the reflection characteristics of permeable vertical seawalls and sloping permeable walls. Comparisons of the numerical results with previous analytical results were given. Losada et al. (1993) used linear theory for water waves impinging obliquely on dissipative multi-layered media to evaluate reflection and transmission coefficients. They examined the case of a periodic medium consisting of alternating layers of upright porous sections and water regions of equal or different thickness. Yu (1995) studied the diffraction of water waves by thin porous breakwaters extending from the water surface down to the bottom of the seabed based on the linear potential wave theory. A relation for the fluid flow through thin porous structures in addition to the conventional governing equation and boundary conditions for small-amplitude waves in an ideal fluid was derived. A two-dimensional problem with long-crested waves propagating in the normal direction of such an infinite porous wall was first solved and the solution was compared with available experimental data. Generally, as summarized above, previous theoretical research on permeable barriers are related to barriers or breakwaters extending from the seabed to above the wave surface, whereas little research has been conducted for permeable wave barriers extending from the water surface to above the seabed. Furthermore, past research on permeable wave barriers or breakwaters can be grouped into two categories, depending on the type of flow in the medium. One is based on the assumption of Darcy flow within the medium. In this case, the barriers function as wave absorbers and are often very thin and have a very fine homogeneous porosity. The other is based on non-Darcy flow within the medium. In this case, the barriers function as rubble mound breakwaters. For the Darcy flow problem, the porous flow velocity is linearly proportional to the pressure difference between the two sides of the porous plate. For the non-Darcy flow problem, the Lorentz principle of equivalent work is often used to represent the velocity field within the barriers. 1.3 Scope of present work The primary aim of the present investigation is to develop a numerical model to predict the interaction of a regular small-amplitude wave train with a thin permeable vertical barrier extending from the wave surface to some distance above the seabed. The case of wave interactions with a pair of such barriers is also studied. The approach used is based on an eigenfunction expansion method and linear diffraction theory. Fluid mechanics within a permeable medium is studied and the boundary conditions at the barrier surfaces are thereby developed. Numerical solutions to the problems are presented and the corresponding numerical models are setup. Comparisons of the numerical results are carried out with those of previous numerical studies for related situations. A selection of numerical results based on the method are presented for the transmission, reflection and energy dissipation coefficients, the wave runup, and the maximum horizontal force on the barrier(s). It is hoped that this numerical model will help provide designers of coastal barriers with a simple, effective and comprehensive method for wave analysis and the prediction of wave interactions with such barriers. Chapter 2 Theoretical development Chapter 2 THEORETICAL DEVELOPMENT 7 A definition sketch of a single thin permeable barrier is shown in Fig. 1. It is noted that in practice the barrier is fixed to the seabed by piles, and these are not shown in the figure. As mentioned earlier, the purpose of the present study is to develop a method for describing regular wave interactions with a thin vertical permeable barrier, extending down from the wave surface to above the seabed (shown as Fig. 1), and then use this to develop expressions for various quantities of engineering interest, including the reflection and transmission coefficients, the wave runup, the wave force and overturning moment. 2.1 Assumptions In order to simplify the problem, the following assumptions have been made: i . The fluid is inviscid and incompressible. Thus the flow is irrotational and may be described by a velocity potential which is governed by the Laplace equation. ii . Linear wave theory is considered suitable for the problem at hand. Thus, the seabed is horizontal and impermeable and wave amplitudes are small. iii. The barrier is considered to be thin and fixed. Any effect caused by support piles is neglected. iv. The problem is assumed to be two-dimensional, with the incident waves approaching the barrier normally. v. There is no overtopping of the waves. Chapter 2 Theoretical development 2.2 Governing equations A regular normally incident small-amplitude wave train of height H and angular frequency co propagates in water of constant depth d past a thin vertical barrier as shown in Fig. 1. As indicated in the figure, the fluid field is divided into regions 1 and 2, corresponding to the upwave and downwave sides respectively. A Cartesian coordinate system (x, z) is defined with x measured in the direction of wave propagation from the location of the barrier, and z measured upwards from the seabed. Based on the assumptions mentioned above, the fluid motion can be described by a velocity potential O which satisfies the Laplace equation within the fluid region, i.e. V 2 0 = 0. In addition, since the wave height is sufficiently small for linear wave theory to apply, <P is subject to the usual boundary conditions, linearized where appropriate, at the seabed, free surface and far field. Finally, O is subject to appropriate boundary conditions at the structure surface. The velocity potential may thus be expressed, using complex notation, in the form: 0(x,z,t) = Re[C(|>(x,z)exp(-icot)] (2.1) where C = - ' ^ - — L - (2.2) 2co cosh(kd) Also, Re[] denotes the real part of the argument, i = ^ 4 , t is time, k is the wave number, and g is the gravitational constant. The potential <j) in Eq. 2.1 is denoted <\>{ and (|)2 in regions 1 and 2, corresponding to the upwave and downwave sides of the barrier respectively (see Fig. 1), and different expressions for <(), and (|)2 are developed. Chapter 2 Theoretical development Along the seabed, the vertical velocity of water particle is zero, conditions at the seabed are 9 Thus the boundary d^_ dz ^ 2 dz 0 = 0 at z = 0 (2.3) Along the wave surface, the free-surface boundary conditions include the kinematic and dynamic boundary conditions, which can be combined into: d<t>, CO dz dz ^,=0 CO 0 2 = 0 at z = d (2.4) where co is angular frequency of the wave motion. At large distances from the barrier, 0 must satisfy a radiation condition, namely 3(4>i-4>w) dx = -ik(<l>i -<t>w) atx = - o 302 dx = ik())2 at x = + o o (2.5) (2.6) where (j)w is incident wave potential which is given as: <j)w = cosh(kz)exp(ikx) (2.7) Along the matching boundary within the fluid, x = 0 at 0 < z < a, the velocity potential and horizontal velocity are equated: 01 =4>2 dx dx forx = 0± 0<z<a (2.8) Chapter 2 Theoretical development 10 Finally, in order to develop suitable boundary conditions along the barrier surface, the fluid flow within porous media needs to be studied and this is considered in the following section. This is shown to give rise to boundary conditions on the barrier surface as: ^ = ^ = -iS'(<|>2-<!>,) forx = 0± a<z<d (2.9) where S' is a complex constant (with unit m"1) which depends on the porosity and the thickness of the barrier. The preceding equations may be collected as follows: V\=0 V2(|)2 =0 dz dz = 0 = 0 d(|), CO 4>,=o dz g d4>2 co2 §2 =° in ( - oo < x < 0; 0 < z < d) in (o < x < +oo; 0 < z < d) at z = 0 at z = d dz g ^ = 3 ^ =-iS'(<|>2 -(!>,) forx = 0± a<z<d <t>l =<|>2 dx dx forx = 0± 0<z<a (2.10) (2.11) (2.12) (2.13) (2.14) (2.15) Chapter 2 Theoretical development 11 9(4>.-4»w) = - i k (0 ! -0 w ) forx = - (2.16) • o o d(j)2 = ik())2 for x = +< (2.17) • o o dx A solution to the above well-posed boundary-value problem is required. 2.3 Flow within barrier The development of Eq. 2.9 is now considered. The small-amplitude wave motion in an undeformable porous medium may be described by the continuity equation and the convection-neglected Euler equation proposed by Sollitt and Cross (1972) as follows: where U is velocity vector; P is dynamic pressure; p is fluid density; e, f and C m are the porosity, the linearized resistance coefficient, and the added-mass coefficient of the porous medium respectively. The last two terms in Eq. 2.19 arise from modeling the porous effect on the fluid motion by two body forces: a linear resistance, which is proportional to the velocity, and an inertia force, which is proportional to the local acceleration. Since the equation arises from a linearization of the friction term which is approximately proportional to velocity squared, it follows that the linearized resistance coefficient f is itself expected to depend on the overall velocity amplitude, and therefore on wave steepness. This aspect is not considered, and here, f is treated as a constant for any given barrier. It is understood that f and C m , which can only be determined experimentally, depend on the viscosity of the fluid and the porosity of the medium, as well as the geometric shape and the surface properties of the medium pores. Thus f and C m are treated as constants for any given barrier. V U = 0 (2.18) m 1-eoTJ e dt (2.19) Chapter 2 Theoretical development 12 Considering only harmonic solutions, U and P can be expressed in the following form: U = u exp(-icot) (2.20) P = p • exp(-icot) (2.21) where u and p are complex amplitudes which are spatial functions only. Substituting Eqs. 2.20 and 2.21 into Eqs. 2.18 and 2.19 yields: V u = 0 Vp + pcoQu = 0 (2.22) (2.23) where Q = f - i 1 + C r l - e (2.24) is a dimensionless porous impedance, whose real and imaginary parts are related respectively to the resistance and the inertial effects of the medium. For the two-dimensional problem at hand, Eqs. 2.22 and 2.23 can be written as: 9p 3x + pcoQu = 0 (2.25) (2.26) where u is velocity component in the x-direction. Integrating Eq. 2.26 with respect to x over the width b of the barrier, one obtains: u = -(Pb--Po+)/(p«t>Q) (2.27) where p^ and pb_ denote the pressure within the porous medium at x = 0 and x = b, respectively. Chapter 2 Theoretical development 13 Equation 2.27 states that the x-component of the velocity within the porous medium is proportional to the pressure decrease in the same direction, but with a phase lag. Since Q has a positive real part and a negative imaginary part, the phase lag is valued between 0 and 7i/2. For a medium in which the resistance dominates the inertial effect, the last term of Eq. 2.24 vanishes and Q. - f . The phase lag then tends to zero and, consequently, Eq. 2.27 corresponds to the Darcy's equation: ^ + pcofu = 0 (2.28) dx In order to relate the velocity and the pressure within the porous medium with those outside, we consider matching conditions at x = 0 and x = b. These include the continuity of the mass flux and pressure expressed as: u^ = ub_ = eu (2.29) Po- = Po+ (2-30) P b-=Pb + (2-31) where the subscripts 0+ and 0- denote the values inside and outside the porous medium at x = 0, and b- and b+ denote the corresponding values at x = b. Substituting Eqs. 2.29, 2.30 and 2.31 into Eq. 2.28, we obtain (p b _-p 0 + ) atx = 0 - a<z<d (2.32) pcobQ ^ ^ ( p b _ - p 0 + ) atx = b+ a<z<d (2.33) Without contradicting the small thickness assumption of the porous barrier, the barrier can be treated approximately as a boundary along x = 0 and, accordingly, Eqs. 2.32 and 2.33 are applied at x = 0 - and x = 0 + , respectively. Therefore, the boundary condition reduces to: Chapter 2 Theoretical development 14 u 0 + = u 0 _ = - ^ ( p b + - p 0 _ ) atx = 0+ a<z<d (2.34) where G = e/Q = G r -r-iG; (2.35) is a complex porous-wall-effect parameter, with G r being its real part and Gj its imaginary part. It is obvious that G = 0 when either e = 0 (i.e. solid wall limit) or Q. —» oo (both the real and the imaginary parts of Q tend to infinity). An increase of the value of the real and/or the imaginary part of G corresponds to an increase of the transparency of the wall. For a medium in which the resistance dominates the inertial effect, the porous-wall-effect parameter G becomes real and Eq. 2.34 reduces to that obeying Darcy's law. Substituting into Eq. 2.34 the following relationships u = - ^ (2.36) dx p = - - ^ = ipaxj> (2.37) the boundary conditions along the surface of the barrier may eventually be expressed as: dp_ = dp_ = , v forx = 0 ± a < z < d ( 2 3 8 ) dx ox where S' = G/b . This corresponds to Eq. 2.9 which has been indicated earlier. 2.4 Solution by eigenfunction expansion method 2.4.1 Eigenfunction expansion method For the barrier above the seabed, appropriate solutions for the potentials (J), and <|>2 that satisfy Eqs. 2.10-2.13, 2.16, and 2.17 may be written as Chapter 2 Theoretical development 15 <t>i = <l>w" S A m cos(ixmz)exp(|imx) (x < 0) (2.39) m=0 OO <t>2 = K + S B m cos(nmz)exp(-p,mx) (x > 0) (2.40) m=0 where the summation terms represent the scattered wave modes. The variables A m and B m are unknown complex coefficients to be determined. In the summation terms, H 0 = - i k (2.41) and the remaining values of | i m are the positive roots of the following equation, taken in ascending order: 2 u.mdtan(p:md) = -—— form>l (2.42) g Thus, each summation in Eqs. 2.39 and 2.40 represents a superposition of a propagating mode (m = 0) and a series of evanescent modes (m > 1) which decay with distance away from the barrier. Equations 2.39 and 2.40 satisfy all the relevant boundary conditions, except the conditions of pressure continuity along the matching boundary. The boundary condition at the barrier surface are needed to determine the coefficients A m and B m . Through the application of the boundary conditions along the surface of the barrier given by Eq. 2.14, and the matching boundary conditions given by Eq. 2.15, and noting that u,0 = - ik and thus cos(|X0z) = cosh(kz), it can be shown that A m = B m for all m, and oo £ A m c o s ( n m z ) = 0 for0<z<a (2.43) m=0 X A m cos(|imz)(|im - 2iS') = - | i 0 cos(n0z) for a < z < d (2.44) m=0 Chapter 2 Theoretical development 16 In previous work based on the eigenfunction expansion method (Losada et al., 1992 and Abul-Azm, 1993), a method involving a least squares fit to the required conditions was used as a basis for estimating A m . Isaacson (1996) used a simplified integration method to solve for A m . Here, using Isaacson's method, Eqs. 2.43 and 2.44 are integrated in an appropriate way and then added to develop a suitable matrix equation for A m . Thus, each equation is first multiplied by cos(pnz), then integrated with respect to z over the appropriate domain of z (i.e. from z = 0 to a, or from z = a to d, as appropriate), and the resulting two equations are then added. Note that in order to make the units consistent in this addition, Eq. 2.44 is first multiplied by the water depth d. This gives rise to the following matrix equation for A m : £ C m n A m - b n m=0 for n = 0,1,...,oo (2.45) where Cmn=f^{0,a) + (iimd-2is)fmn{a,d) b„ = - M f 0 „ ( a , d ) (2.46) (2.47) P f ™ (a. P) = j cos(p. mz) cos(|i n z)dz sin(p:J) ^ sin(|i+(3) sin(p:_a) sin((i+a) 4n — {[2fxmp + sinljij)] - [2[ima + sin(u,ma)]} (m* n) (m=n) (2.48) For convenience, the dimensionless parameter S = S'd is used throughout. Also =^m-lin.and [i+ =p.m+p. n. In a numerical solution to the problem, Eq. 2.45 is truncated to a finite number of terms N, and thus becomes a complex matrix equation of rank N which can be solved for the first Chapter 2 Theoretical development 17 N unknown coefficients A m . Once these have been obtained, the various quantities of engineering interest may readily be obtained. 2.4.2 Transmission and reflection coefficients According to linear wave theory, the wave surface elevation is given as: 1 r| = — g V dt J z = d (2.49) Substituting Eqs. 2.1, 2.2, 2.39 and 2.40 into the above equation, the complex wave elevation can be obtained as: T i l TI2 H H 2 e x P ( ^ ) - c o s h ( k d ) Z A m cos(timz)exp(p,mx) exp(- icot) 1 0 0 exp(ikx) + c o s h ( k d ) X A m cos(ixmz)exp(- n m x) exp(- icot) (2.50) (2.51) where the subscripts 1 and 2 denote the upwave and downwave sides respectively. Since the evanescent wave modes (m> 1) in Eqs. 2.50 and 2.51 decay with distance away from the barrier, only propagating wave modes (m = 0 ) remain at locations far from the barrier. Noting that | l 0 = - ik and thus cos(|l0z) = cosh(kz), the complex transmission and reflection coefficients, denoted here by CT and C R respectively, are obtained as: C T = 1 + A r (2.52) (2.53) It is obvious that: (2.54) Chapter 2 Theoretical development 18 For engineering purposes, one usually defines instead real transmission and reflection coefficients, denoted here as K t and K r respectively, as the appropriate ratios of wave heights: K t = H t/H and K r = H r/H, where H t and H r are the transmitted and reflected wave heights respectively. The relationship between the complex and engineering definitions of transmission and reflection coefficients is simply: K t = C T = 1 + A 0 K r — C^ - A 0 (2.55) (2.56) 2.4.3 Energy balance An energy loss due to the presence of the barrier should also be considered. The energy flux, which is the average rate of transfer of energy per unit width across a plane of constant x, is defined as E= J p +—p(u2 + w 2 ) + pgz udz 0 L. (2.57) where an overbar denotes a time average. On the basis of the unsteady Bernoulli equation, and replacing the limit z - d + T| by z = d, since linear theory is used, Eq. 2.57 becomes: u E = -pJ dt dx dz (2.58) Recalling the form of Eq. 2.1, E can be expressed as: E = pco J Re[iC()) exp(- icot)] Re d§ , . C—exp(-icot) dx dx (2.59) Initially, the energy flux of the incident wave train alone is considered. Thus one takes initially, ())(x,z) = ())w(x,z) = cosh(kz)exp(ikx) (2.60) 19 Chapter 2 Theoretical development Substituting Eq. 2.60 into 2.59, the energy flux of the incident wave is eventually given by the well-known result (e.g. Sarpkaya and Isaacson, 1981): E w = ^ p g H 2 c G (2.61) where 2kd c r = — 2 1 + -sinh(2kd) and c is the wave speed. (2.62) We consider now the more general case of the wave field at large distances upwave or downwave of the barrier. Far from the barrier, the evanescent wave components are absent and only propagating terms are present. Thus, referring to Eqs. 2.39 and 2.40, in general the velocity potential may be written in the form: (|)(x, z) = cosh(kz)[A+ exp(ikx) + A_ exp(- ikx)] (2.63) where A+ and A. are the complex dimensionless amplitudes of the waves propagating in the + X and -x directions, respectively. Substituting Eq. 2.63 into Eq. 2.59, and carrying out the necessary integration, one eventually obtains: E = ^ p g H 2 c G ( | A + | 2 - | A _ f ) (2.64) On the upwave side, A + = 1 and A_ = - A 0 , and on the downwave side, A + = 1 +A 0 and A = 0 . Thus one obtains: E ^ E ^ I - I A J 2 ) (2.65) E 2 = E w ( | l + A 0 | 2 ) (2.66) Chapter 2 Theoretical development Consequently the energy flux disipated at the barrier, E d, is given by 20 Ed - Ej E 2 (2.67) Substitute Eqs. 2.65 and 2.66 into Eq. 2.67, finally, E d = E w ( l - | A 0 | 2 - | l - A 0 | 2 ) An energy loss coefficient can be defined as (2.68) K - = ] f - = H A o l - | 1 - A 0 | From Eqs. 2.55 and 2.56, K , 2 + K R 2 = | A 0 | 2 + | I - A 0 | 2 It is obvious that (2.69) (2.70) K t 2 + K / + K d = l (2.71) 2.4.4 Wave runup When (()[ and <|)2 are used in the representation of <E> (see Eq. 2.1), and x = 0, z = d are taken, then Eq. 2.49 provides the runup on the upwave (subscript -) and downwave (subscript +) sides of the barrier, respectively. The runup can thereby be expressed as H R = H 1 -1 + s5 )S A - c o s ( ^ d ) j j i 5 j J > . c « < M ) (2.72) (2.73) Chapter 2 Theoretical development 2.4.5 Wave loads 21 On the basis of the linearized, unsteady Bernoulli equation, the hydrodynamic pressure along the upwave and downwave sides of the barrier surface can be expressed respectively as: ao, Po- =-p - p -at ao x = - 0 at x=+0 Hence"the maximum pressure difference across the barrier, Ap(z), is given as: Ap(z) - max(p0_ - p 0 + ) = max 1 d ( * 2 - * . y at = pgH cos h(kd) XAmCos(|imz) m=0 (2.74) (2.75) (2.76) The maximum horizontal force per unit width on the barrier, F, can be obtained by integrating Ap with respect to z. This leads to: F = p g H ^ i S d T Am[sin(^md)-sin(itma)] m=0 (2.77) Likewise, the overturning moment per unit width on the barrier about the seabed, M , can be obtained by an appropriately weighted integration of Ap with respect to z. This leads to: Chapter 2 Theoretical development 22 d M = } Ap(z)zdz a ; H — - — cosh(kd) (2.78) ~ Am[|imdsin(|imd)-\imasin(p:ma) + cos(p:md)- cos(p.ma)] It is noted that the terms m = 0 in the above equations may readily be expressed directly in terms of real quantities, if required. 2.4.6 Limiting conditions Various limiting conditions to the preceding solution can be established and are discussed in the following. Solid barrier (e = 0) For a solid barrier (Sr = S; =0j, various limits to the above solution are known. As kh —> oo (short wave limit) or as h/d —> 1 (full barrier extending to the seabed), the solution should correspond to full reflection on the upwave side, and no wave activity on the downwave side, so that K t ^ 0 , K r - > 1 , K d - > 0 , R_/H->1, R + / H - > 0 , Ap(z) —» Ap s(z), F—>FS, and M—>M S ; where the subscript s corresponds to the quantities associated with a standing wave train on the upwave side and no wave activity on the downwave side. Expressions for F s and M s are easily derived and given respectively as: pgH sinh(kd) - sinh(ka) (2.79) k [ cosh(kd) 2 k 2 | cosh(kd) 1 pgH J [kd sinh(kd) - cosh(kd)] - [ka sinh(ka) - cosh(ka)] j (2.80) Chapter 2 Theoretical development 23 On the other hand, as kh —» 0 (long wave limit), or as h/d —> 0 (no barrier), the solution should correspond to the undisturbed incident wave field on either side of the barrier, so that K t - > 1 , K r - > 0 , K d - > 0 , R_/H->l /2 , R+/H->l/2, Ap(z)-»0, F->0,and M - » 0 . Fully transparent barrier (e = 1) A fully transparent barrier (Sr = +°oJ should correspond to the undisturbed incident wave field on either side of the barrier, and so will lead to the same results as above for kh —» 0 (the long wave limit), or as h/d —> 0 (no barrier). Thus K t —» 1, K r —»0, K d —»0, R_/H->l /2 , R + / H - » l / 2 , Ap(z)->0, F —> 0, and M -» 0. 2.5 Extension to the two-barrier problem The preceding method can be extended to the two-barrier problem in which an incident wave train interacts with two identical barriers a distance 2\ apart. A definition sketch of the problem is given in Fig. 2. As indicated in the figure, the flow field is now divided into three regions, denoted 1, 2 and 3, corresponding to the regions upwave of the barriers, between the barriers, and downwave of the barriers, respectively. 2.5.1 Governing equations and boundary conditions The velocity potential <P within the three fluid regions may be taken as being expressed by the same form as Eq. 2.1. Here, the potential 0 in Eq. 2.1 is denoted 0t , 02 and 03 in regions 1, 2 and 3. The potentials satisfy the Laplace equation within the fluid region: V 2 0 t =O (2.81) V 2 0 2 =0 (2.82) V > 3 =0 (2.83) Chapter 2 Theoretical development 24 The boundary conditions along the seabed, free surface, barrier surface, together with the matching conditions and radiation conditions are given respectively by: 3z dz dz dz dz 0 = 0 = 0 ^ 0 ) 2 A n 30 2 CO " 0 2 = 0 co' • — (t)3=0 0, =02 90, 302 " 3x 02 302 303 dx 3x 30, 302 3x 3x 302 303 3x 3x 3(0, -0w) at z = 0 at z = d at (x = -X±, 0 < z < a) at (x = +X±, 0<z<a) 3x 303 3x = ik03 -iS'(02-0,) at(x = -X±, a<z<d) -iS'(03 -0 2) at(x = +X±, a<z<d) -ik(0,-0w) for x = - ° o for x = + o o (2.84) (2.85) (2.86) (2.87) (2.88) (2.89) (2.90) (2.91) Chapter 2 Theoretical development 25 2.5.2 Solutions by eigenfunction expansion method The solutions satisfying the above governing equations and boundary conditions except Eqs. 2.86 - 2.89 can be expressed as: <i>i =<l>w +X A imcos(p.mz)exp[p.m(x + ?0] (2.92) m=0 oo oo <t>2 =X A 2 m cos(p: m z)exp[-H r a (x + ?i)]+XA3mCos(p:mz)exp[ixm(x-X)] (2.93) m=0 m=0 <t>3 = XA 4 m cos ( | i m z)exp[ - | i m (x -X) ] (2.94) m=0 where ^ is incident wave potential which is given by Eq. 2.7. The summation terms in the above equations represent the scattered wave modes. Of these, each summation term represents a superposition of a propagating mode (m = 0) and a series of evanescent modes (m > 1) which decay with distance away from the barriers. The variables A i m , A2m, A3m, and A4m are unknown complex coefficients to be determined. The eigenvalues | j , m are determined by Eqs. 2.41 and 2.42. By applying the boundary conditions expressed in Eqs. 2.86 and 2.87, the following equations are obtained for 0 < z < a: CO CO oo X A i m cos(nmz) - X A 2 m cos(|imz)- X A 3 m cos(|imz)exp(- 2|imX) m=0 m=0 m=0 (2.95) = -cos(p:0z)exp(|a,0X.) X A i m l t m cos(|imz)+ X A 2 J r a cos(|imz) m=0 m=0 " X A 3 m ^ m cos(fimz)exp(- 2\iml) = \i0 cos(p:0z)exp(|i0^) m=0 oo oo oo X A 2 m cos(p:mz)exp(-2ixm^)+ X A 3 m c o s ^ z ) - X A 4 m cos(pimz) = 0 (2.97) m=0 m=0 m=0 Chapter 2 Theoretical development 26 X A 2 m H m cos(nmz)exp(- 2\Lj) - X A 3m[im cos(pmz) - X A 4 m ^ m C O s ( ( I m z ) = 0 m = ° . m = 0 (2.98) m=0 By applying the boundary conditions along the barrier's surface (expressed as Eqs. 2.88 and 2.89), the following equations are obtained for a < z < d : X A l m cos(nmz)(pm - iS') + i S ' X A 2 m cos(|imz) m=0 m=0 + i S ' X A 3 m cos(pmz)exp(-2p:mX) = ((i0-i-iS/)cos(N.0z)exp(p.0^) (2.99) m=0 i S ' E A l m cos(p.mz)+ X A 2 m cos(pmz)(pm - iS') m=0 m=0 " Z A 3 m cos(p.mz)exp(- 2|imX)(pm + iS') = -iS'cos(p0z)exp(|i0?i) (2.100) m=0 XA 2 mcos(u.mz)exp(-2u.m? L)( |^m+iS')-XA 3 mcos(|Llmz)(n.m-iS') m=0 m=0 - i S ' Z A 4 m c o s ( | l m z ) = 0 (2.101) m=0 i S ' S A 2 m cos(p.mz)exp(-2pmX) + i S ' £ A 3 m cos(pmz) m=0 m=0 oo + X A 4 m c o s ( u m z ) ( | i m - i S / ) = 0 (2.102) m=0 As an extension to the solution for the single barrier problem, Eqs. 2.95 and 2.99, 2.96 and 2.100, 2.97 and 2.101, 2.98 and 2.102 are integrated with respect to z over 0 to a or over a to d as appropriate, and each pair of resulting equations is then added to obtain four sets of equations for Aj m . (Some equations are first multiplied by d to make the units consistent.) These are: Chapter 2 Theoretical development 27 X A l m [ f m n fa, a) + (u.md - iS)f „„, (a, d)] m=0 + EA 2 m [ - f m n (0 , a ) + iSfmn(a,d)] m=0 + X A 3 m exp(- 2iimX)[- fm{0,a) + iSfmn(a,d)] m=0 = - exp(mX)f0n (0, a) + (n0d + is)exp(p.0X)f0n (a, d) XA l m [ ix m df m n (0 ,a) + iSfmn(a,d)] m=0 + X A 2 m [ixmdfmn (0, a) + (u.md - is)f m (a, d)] m=0 + X A 3 m exp(- 2|im^)[- |J.mdfmn (0, a) - (jimd + isjf^ (a,d)] m=0 = M exp(|i0X)fOn (0, a) - iS exp(p:0?i)f0n (a, d) X A 2 m exp(- 2M-m^)[fmn (0, a) + (|i Bd + is)f ^ (a, d)] m=0 + X A 3 m [f „„, (0, a) - ( j L L m d - is)f m (a, d)] m=0 + XA 4 m [-f m n (0 ,a)-iSf m n (a ,d)] = 0 m=0 (2.103) (2.104) (2.105) (2.106) . X A 2 m exp(- 2\Lj.i\Lm& ^ (0, a) + iSf „„, (a, d)] m=0 + X A 3 m [- | i m df „„ (0, a) + iSf^ (a, d)] m=0 + X A 4 m [- a mdf „ (0, a) + (u n d - is)f m n (a, d)] = 0 m=0 where the function f^Q is defined by Eq. 2.48. Equations 2.101 - 2.104 can be expressed as the following matrix equation for A i m , A 2 m , A 3 m and A ^ : Chapter 2 Theoretical development 28 I c r m=0 Z c g - } m=0 m=0 ^ C ( 7 ) m=0 m=0 m=0 I c ( 2 r } m=0 £c<™> m=0 m=0 m=0 w m=0 m=0 m=0 S c s r > m=0 m=0 S c r m=0 A l m " V A 2 m b2„ A 3 m A4m. b 4 n . (2.107) where p(mn) M 2 C |r ) = fmn(0,a) + (nJ md-iS)fm n(a,d) -fmn(0,a) + iSfmn(a,d) exp(- 2u.m (0, a) + iSf „„ (a, d)] C<7) = u.mdfmn(0,a) + iSfmn(a,d) C ^ = pmdf m n (0, a) + (|imd - is)f (a,d) exp(- 2pmk)[- u.mdf „„, (0,a) - (pmd + i s ) ^ (a,d)] exp(- 2pmX)[f m (0, a) + (|imd + iS)f „„, (a, d)] fM1(0,a)-(jimd-iS)f11111(a,d) Cl7 ,=-fn-(0,a)-iSfB B(a,d) C42™ > = exp(- 2u m A,Xlindf „ (0, a) + iSf M (a, d)] c!r)=-M.B1dfII11(0,a) + iSfI11I1(a,d) -p:mdfmil(0,a) + (nind-is)f11111(a,d) p(mn) ^23 p(mn) *-32 i(mn) '33 i(mn) '44 C(mn) p(mn) p(mn) p(mn) /-> '24 " ^31 41 (2.108) (2.109) (2.110) (2.111) (2.112) (2.113) (2.114) (2.115) (2.116) (2.117) (2.118) (2.119) (2.120) Chapter 2 Theoretical development 29 b l n =exp(fi0A.)[-f0n(0,a) + (|i 0d + is)fon(a,d)] (2.121) b 2 n =exp(|i0?i)[|i0df0n(0,a)-iSf0n(a,d)] (2.122) b 3 „ = b 4 n = 0 (2.123) In a numerical solution to the problem, Eq. 2.107 is truncated to a finite number of terms N, and thus becomes a complex matrix equation of rank 4N which can be solved for the first N unknown values of each set of coefficients A i m , A2 m , A 3 m and A 4 m . Once these have been calculated, the various quantities of engineering interest may readily be obtained. The transmission and reflection coefficients are given in terms of A i m and A4m by: K t = | A j (2.124) K r =|A 1 0 | (2.125) Considering energy conservation, at large distances from the two barriers, x - X and x + X are equivalent to x. Therefore, Eq. 2.63 can still be used to express the upwave and downwave wave potentials ()), and (j)3. Here, on the upwave side, A + = 1, A_ = A 1 0 , and on the downwave side, A_ = A 4 0 , A_ =0. Substituting these values into Eq. 2.64, the energy flux on the upwave and downwave sides, Ei and E3 respectively, can be expressed as: = E w ( l - | A 1 0 | 2 ) (2.126) E 3 = E W | A 4 0 | 2 (2.127) Consequently the energy flux dissipated at the barriers, E d , is given by E d = E, - E 3 = E w ( l - | A 1 0 | 2 - | A 4 0 | 2 ) (2.128) Chapter 2 Theoretical development Therefore the energy dissipation coefficient can be obtained as: 30 Kd - g — 1 | A 1 0 | A 4 0 (2.129) From considerations of energy conservation, substituting Eqs. 2.124 and 2.125 into Eq. 2.129, these are related to the energy dissipation coefficient Kj: K 2 + K 2 + K d = l The wave runup on both sides of each barrier can be obtained as: (2.130) H R, = — 1 _ 2 expi (H<A)- cosh(kd) ~ 0 SAlmCOs(p.md) (2.131) H R l + - — 1 + 2 H R, = — 2" 2 H R l t — 2 + 2 cos h(kd) X A 2 m cos(p.md)+ X A 3 m cos(p,md)exp(- 2u.mA,) ,m=0 m=0 cos h(kd) X A 2 m C 0 S ( l J 'm a 1 )eXp(- 2(1 mA,) + X A 3 m C O S ( M ) m=0 m=0 ^ ^ l A ^ e o s ^ d ) (2.132) (2.133) (2.134) where Ri_ and Ri+ denote the wave runup on the upwave side and downwave side of the barrier 1, and R 2. and R 2 + denote the wave runup on the upwave side and downwave side of the barrier 2, respectively. The maximum pressure difference across each barrier can be obtained by: Chapter 2 Theoretical development 31 Ap1(z) = max(p_JL_-p_x+) = max 3t pgH X A i m cos(p:mz)- X A 2 m cos(p: m-0 m-0 2cosh(kd) X A 3 m cos((imz)exp(- 2\imX) + cos(|i0z)exp(|i0?i)| (2.135) m - 0 Ap 2 (z) = max(px_ - p u ) = max d ( * 3 - * 2 y at PgH 2cosh(kd) X A 2 m cos(|imz)exp(- 2\Lm\) + m-0 X A 3 m cos((tmz) - X A 4 m cos(^mz)| (2.136) m-0 m-0 Integrating Eqs. 2.135 and 2.136 with respect to z, result in the maximum horizontal force per unit length on each barrier: PgH F l 2cosh(kd) X A i m ^ i m X A 2 m § i m m=0 m=0 e x p ( - 2ixm^) + 510exp(^0?i)| m=0 (2.137) PgH F ' 2cosh(kd) 2^ A 2 m 8 l m exp(-2\imX)+ X A 3 m 5 l m - X A 4 m 5 lm m=0 m=0 m=0 (2.138) where sin(^md)-sin(p:ma) (2.139) The overturning moment per unit width on each barrier about the seabed can be obtained as: Chapter 2 Theoretical development 32 PgH M = 1 2cosh(kd) XA l m ^ 2 m X A 2 m ^ 2 m m=0 m=0 - X A 3 m 5 2 m exp(- 2\imX) + 5 2 0 exp(n,0A,) m=0 (2.140) M , -PgH 2cosh(kd) _2LrfA2m^2m exp(- 2 p m ? i ) + X A 3 m 5 2 m - X A 4 m S 2m m=0 m=0 m=0 (2.141) where S 2m = M-md sin(u,md) - p m a sin(p.ma) + cos(pmd) - cos(|ima) (2.142) Chapter 3 Development of numerical model Chapter 3 DEVELOPMENT OF N U M E R I C A L M O D E L 33 3.1 General Based on the theoretical analysis outlined in the preceding chapter, a numerical model written in the FORTRAN language has been developed. The objective of setting up the numerical model is to obtain various quantities of engineering interest, such as the transmission and reflection coefficients, the dimensionless wave runup, and the wave forces acting on the barriers. As mentioned in Chapter 2, all these quantities depend on the coefficients A m (one-barrier problem) or the coefficients A i m , A 2 m , A3 m and Aim (two-barrier problem). Thus the task becomes one of calculating the various coefficients C m n and solving the matrix equations given by Eq. 2.45 (one-barrier problem) and Eq. 2.107 (two-barrier problem). Once these coefficients are obtained, the required quantities can readily be calculated in the manner indicated. 3.2 Parameter selection Before setting up the numerical model, it is appropriate to consider initially the selection of suitable input and output parameters. The input parameters include the incident wave conditions (wave height, wave period, water depth), barrier conditions (the draft h, the permeability parameter S'). The output parameters include the transmission and reflection coefficients, the wave runup, the wave forces and the overturning moments. In presenting results of the numerical model, it is useful to carry out a dimensional analysis. Bearing in mind that the independent parameters may be considered as H, d, h, k and S' (also X for the two-barrier problem), a dimensional analysis indicates that the independent parameters can be selected as, for example, kh (or kd), h/d, X/d, and S (S = S'd). It is unnecessary to specify the wave height H since the dimensionless output parameters are obtained in dimensionless form which are themselves independent of H. As Chapter 3 Development of numerical model 34 indicated in Chapter 2, the permeability parameter S' can be only obtained by experiments, thus S is treated as an input value in the numerical model. The output parameters are all expressed in dimensionless form. These parameters include transmission and reflection coefficients, wave runup coefficients, wave force coefficients and overturning moment coefficients. The wave runup coefficients are defined as R/H, and the dimensionless force coefficients and the overturning moment coefficients are defined as: F C F = pgHd M C m ~ PgHd2 (3.1) (3.2) Chapter 4 Results and discussion Chapter 4 RESULTS AND DISCUSSION 35 A range of numerical results have been obtained in order to examine the validity of the numerical model and to provide predictions for practical situations. Also, comparisons with previous predictions have been carried out. 4.1 One-barrier problem Previous numerical results for wave interactions with a permeable vertical barrier fixed above the seabed appear not to be available. Thus, in order to examine the effectiveness of the numerical model, comparisons are made with previous predictions for the following two limiting cases: a solid barrier above the seabed (Abul-Azm, 1993, Losada et al., 1992, and Isaacson, 1996), and a permeable barrier extending down to the seabed (Yu, 1995). Subsequently, results for a single permeable barrier above the seabed are presented and discussed. These provide various quantities of engineering interest as functions of the dimensionless wave frequency kh, the relative penetration of the barrier h/d, and the dimensionless permeability parameter S. 4.1.1 Comparisons with previous predictions For a solid barrier above the seabed (S = 0), the problem becomes exactly that treated previously by several authors, including Losada et al. (1992), Abul-Azm (1993) and Isaacson (1996). For this limiting case, the present numerical model coincides exactly with that of Isaacson (1996), and is slightly different from that used by Abul-Azm and Losada et al. in that the coefficients A m are obtained here by a simplified integration method rather than a least squares method. The present numerical results have been compared with those of previous studies for a wide range of conditions, and in all cases show very close agreement. Differences among the results are not distinguishable on ordinary plots and hence are not shown here. Chapter 4 Results and discussion 36 For a permeable barrier extending from the free surface down to the bottom of the seabed (h/d = l), the problem becomes exactly that solved by Yu (1995). The transmission and reflection coefficients as functions of the permeability parameter S have been calculated by the present method for a range of relative water depths kd and compared with the results given by Yu. Once more, the agreement has been found to be very close. 4.1.2 Transmission and reflection coefficients Figure 3 shows the transmission coefficient K t and the reflection coefficient K r as functions of S for the case kh = 0.5 and h/d = 0.5. The results are shown as contours of each parameter plotted in the complex S plane. S r and Sj denote the real and imaginary parts of S respectively. The waviness in some of the contours is due to the contour fitting program that was used. The results for S = 0, for which K t = 0.93 and K r = 0.37, correspond to the case of an impermeable barrier. The results along the abscissa (Sr) correspond to the Darcy flow assumption, with increasing S r corresponding to increasing porosity, decreasing friction coefficient, and/or decreasing barrier width (see Eqs. 2.24 and 2.35). This figure shows that K r decreases with increasing values of both the real and imaginary components of S; that K t varies quite strongly with Sr, initially decreasing and then increasing as the value of S r increases; and that K t varies quite weakly with Sj. Figures 4 and 5 show K t and K r as functions of kh for the case h / d = 0.5 and for various values of S. The selected values of S correspond to an impermeable barrier (S = 0), the Darcy flow assumption with a smaller porosity (S = 0.5) and a larger porosity (S = 2.0), and two cases when both frictional and inertial effects are equally important (S = 0.5 + 0.5i, S = 2.0 + 2.0i). Figure 4 shows how K t decreases as kh increases. The case of an impermeable barrier (S = 0) is seen to give a lower transmission coefficient than a porous barrier, as expected, while the case of the more porous barrier (S = 2.0 + 2.0i) generally gives the highest transmission coefficient. The two cases S = 2.0 and S = 2.0 + 2.0i indicate that Si has a relatively small effect on K t and K r in comparison to Sr. The opposite trends for the reflection coefficient are indicated in Fig. 5. Chapter 4 Results and discussion 37 Figures 6 and 7 show K t and K r as functions of h/d for kd = 1.0 and for various values of S. The figures indicate that K t decreases as h/d increases, and K r increases as h/d increases, as expected. It is noted that the differences of K t and K r among the barriers with different S values become more significant as h/d increases. It is also noted that for a solid barrier (S = 0), as kh —> °° (short wave limit) or as h/d —> 1 (full barrier extending to the seabed), the numerical results correspond to the expected limit of full reflection on the upwave side, and no wave activity on the downwave side: K t —> 0, K r —> 1. 4.1.3 Energy dissipation coefficient Figure 8 shows the energy dissipation coefficient Kd as contours in the complex S plane for the case kh = 0.5, h/d = 0.5 corresponding to Fig. 3. Figure 9 shows Kd as a function of kh for h/d = 0.5 and for various values of S. As indicated in these figures, Kd increases initially and then decreases as S r increases, and decreases quite weakly with Sj. It is seen that the energy dissipation coefficient reaches a maximum for a certain value of S under a specified wave condition, or under a certain wave condition for a specified barrier. For example, in Fig. 8, when h/d = 0.5 and kh = 0.5, a barrier with approximately S = 0.7 gives rise to a maximum energy dissipation coefficient K d ~ 0.26, and as shown in Fig. 9, when S = 0.5, a maximum K d = 0.45 is reached when kh = 1.0. Figure 10 shows the energy dissipation coefficient as a function of h/d for kd = 1.0 and for various values of S. This figure indicates how Kd increases as h/d increases. A solid barrier does not give rise to any energy dissipation so that a solid line cannot be observed. 4.1.4 Wave loads Figure 11 shows the force coefficient C F as contours in the complex S plane for the case h/d = 0.5, kh = 0.5; Fig. 12 shows CF as a function of kh for h/d = 0.5 and for various values of S; and Fig. 13 shows C F as a function of h/d for kd = 1.0 and for various values of S. These figures indicate that the wave force coefficient first increases with a decrease Chapter 4 Results and discussion 38 of Si or Sr, and reaches the maximum as S = 0; that the wave force coefficient increases sharply and then decreases more gradually as kh increases giving a peak value under certain wave conditions; and that the wave force coefficient increases as h/d increases, as expected. Figures 14, 15 and 16 show the overturning moment coefficient C M as a function of S, kh, and h/d respectively for various conditions. These figures indicate exactly the same trends as those of the force coefficient CF-4.1.5 Wave runup Figure 17 shows the wave runup coefficients for both sides of the barrier, R_/H and R + / H , as contours in the complex S plane for the case h/d = 0.5, kh = 0.5. The wave runup on the upwave side of the barrier (shown as solid lines) initially increases, then decreases as Sr increases. The opposite trend is seen for the wave runup on the downwave side of the barrier (shown as broken lines). The wave runup on both sides vary slowly as Sj increases. These phenomena coincide with that of the wave energy dissipation coefficient. The wave runup coefficients on the upwave side and the downwave side as functions of kh for h/d = 0.5 and for various values of S are shown in Figs. 18 and 19. These figures show that, for a porous barrier, the wave runup on the upwave side increases as kh increases, and the magnitude ranges from H/2 to H. On the other hand, the wave runup on the downwave side decreases from H/2 to a small value (0 for a solid barrier) as kh increases. The wave runup coefficients on the upwave side and the downwave side as functions of h/d for kd = 1.0 and for various values of S are shown in Figs. 20 and 21. It is seen that the wave runup on the upwave side increases as h/d increases, and reaches a constant value for a more porous barrier. The wave runup on the downwave side, on the other hand, decreases as h/d increases, and also reaches a constant value for more porous barrier. Chapter 4 Results and discussion 39 It is noticed that for a solid barrier (S = 0), as kh —> <*> (short wave limit) or as h/d —> 1 (full barrier extending to the seabed), the numerical results correspond to full reflection on the upwave side, and no wave activity on the downwave side, so that, R_/H-> 1, R + / H - > 0 . 4.2 Two-barrier problem Previous results for wave interactions with two permeable vertical barriers fixed above the seabed appear not to be available. Thus no comparison can be made to examine the effectiveness of the numerical model for two-barrier problem. However, for some limiting conditions, some results are evident and can be examined by the numerical model. Bearing in mind the difficulty of checking the numerical model, it is nevertheless useful to examine the numerical prediction of wave effects on two permeable vertical barriers. 4.2.1 Transmission and reflection coefficients Figure 22 shows the transmission and reflection coefficients as contours in the complex S plane for the case h/d = 0.5, kh = 0.5, X/d = 1.0. It is seen that, similar to the one-barrier case, K r decreases with increasing values of both the real and imaginary components of S; K t varies quite strongly with Sr, initially decreasing and then increasing as the value of S r increases; K t varies quite weakly with Sj. It is noted that both the values of K t and K r are smaller than the corresponding values for the one-barrier case. Figures 23 and 24 show the transmission and reflection coefficients as functions of kh for h/d = 0.5, ?i/d = 1.0 and for various values of S. Generally, K t decreases and K r increases as kh increases. However, under certain wave conditions, phenomena of resonance occur, where K t increases sharply, and K r decreases sharply. The figures show that for solid barriers there are two points of significant resonance occurring at kh = 0.907 and kh = 1.585. For permeable barriers, additional phenomena of resonance can be seen. Thus K t and K r vary periodically with kh. Chapter 4 Results and discussion 40 Figures 25 and 26 show the transmission and reflection coefficients as functions of X/d for h/d = 0.5, kd = 1.0 and for various values of S. It is seen that K t and K r increase and decrease periodically as X/d increases, and the magnitudes of changing values of K t and K r decrease as S increases. It is found that the peak (or trough) values of K t and/or K r repeatedly occur when the distance between the two barriers is increase by nL/4 (n=l,2,...), where L is the wave length. Since the phase angles are dependent on the wave conditions (wave period, water depth) and barrier conditions (barrier length under the water, the width, the porosity parameter), the points of X/d at which the peak (or trough) values of K t (or K r) occur are different for different values of S as h/d = 0.5 and kd = 1.0, but still can be calculated. 4.2.2 Energy dissipation coefficient The energy dissipation coefficient IQ as contours in the complex S plane for the case h/d = 0.5, kh = 0.5, X/d - 1.0 is shown in Fig. 27. As indicated in the figure, Ka increases initially and then decreases as S r increases, and decreases quite weakly with Sj. The maximum energy dissipation coefficient occurs for a certain value of S under a specified wave condition, or under a certain wave condition for a specified barrier. As shown in Fig. 27, when h/d = 0.5, kh = 0.5 and X/d = 1.0, the barriers with approximately S = 0.7 give rise to a maximum wave energy dissipation coefficient K d ~ 0.52. Compared with the one-barrier case shown in Fig. 8 (a maximum wave energy dissipation coefficient K d = 0.26), it's seen that more energy is dissipated. Figure 28 shows the energy dissipation coefficient as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S. As shown in the figure, generally Kd increases initially and then decreases as kh increases, but varies periodically as kh changes in small ranges. The energy dissipation coefficient for solid barriers appears to be zero, as expected. For permeable barriers, the Kd value reaches its maximum under certain wave conditions. Chapter 4 Results and discussion 41 Figure 29 shows the energy dissipation coefficient as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S. This figure indicates that, Kd varies periodically as the distance between the two barriers increases. The magnitudes of changing values of Kd decrease as S increases. It is also found that the maximum (or minimum) values of Kd repeatedly occur as the distance between the two barriers is increased by nL/4 (n = 1,2,...), where L is the wave length. 4.2.3 Wave forces Figures 30 and 31 show the wave force coefficients CF of the upwave barrier and the downwave barrier as contours in the complex S plane for h/d = 0.5, kh = 0.5 and X/d = 1.0. The force coefficients Cp of the upwave barrier and the downwave barrier as functions of kh are shown in Figs. 32 and 33 for h/d = 0.5, X/d = 1.0 and for various values of S. Figures 34 and 35 show the CF values for the upwave and the downwave barriers as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S. The figures indicate that the wave forces acting on the barriers decrease as S r and/or Si increases. The wave force acting on the downwave side barrier is smaller than that acting on the upwave barrier. Some resonance phenomena occur under certain wave conditions, which give rise to relatively high or low values of CF- The CF values of both barriers vary periodically as the distance between the barriers increases by nL/4 (n=l,2,...), where L is the wave length. As the barriers are more porous (S increases), the forces acting on the barriers become more stable and vary slightly as X/d increases. Chapter 5 Conclusion and recommendations Chapter 5 CONCLUSIONS AND RECOMMENDATIONS 42 5.1 Conclusions A numerical model is set up based on the linear diffraction theory and the eigenfunction expansion method to predict the interaction of a regular small amplitude wave train with a thin permeable vertical barrier. The fluid mechanics within a permeable barrier is studied and the boundary conditions along the barrier's surface are thereby developed. Expressions are obtained for parameters of engineering interest, including the transmission and reflection coefficients, the wave runup on the upwave and downwave sides of the barrier, and the maximum horizontal force and overturning moment on the barrier. Comparisons have been carried out with previous numerical studies for the limiting cases of an impermeable barrier and a permeable barrier extending down to the seabed, and very close agreement has been obtained in all such cases. A selection of numerical results based on the present method have been presented. These show the transmission, reflection and energy dissipation coefficients, the maximum horizontal force and overturning moment coefficients, as well as the wave runup on the upwave and downwave sides of the barrier as functions of the dimensionless wave frequency kh, the relative penetration of the barrier h/d, and the permeability parameter S. These results exhibit various features of interest, some of which are expected. It is concluded that, the various parameters of interest depend on kh and h/d for a permeable barrier in a similar manner as for a solid barrier, while they also depend significantly on the permeability parameter. Under a specified wave condition, a barrier with a given permeability will give rise to a maximum energy dissipation coefficient. It is also seen that a certain wave condition may produce a maximum horizontal force and overturning moment on the barrier. As an extension to the present method, another numerical model is set up to predict the wave interactions with a pair of thin permeable vertical barriers fixed above the seabed. The two barriers are identical in permeability and penetration. A selection of numerical Chapter 5 Conclusion and recommendations 43 results are presented. The results show the transmission, reflection and energy dissipation coefficients, the maximum horizontal force acting on both the upwave and downwave barriers as functions of the complex barrier permeability parameter S, the wave frequency kh and the dimensionless distance between the barriers X/ d. Under the same wave conditions, the transmission and reflection coefficients in the two-barrier case are smaller than those in the one-barrier case, and the wave energy dissipation coefficient is greater. Thus compared to one-barrier problem, two barriers generate relatively smaller reflection and transmission and dissipate more wave energy. Some resonance occurs under certain wave conditions, which K t , K r , Kd and CF are extraordinarily high or low. Under a certain wave condition, all the values of K t , K r , Kd and CF vary periodically as the distance between the barriers is increased by a quarter of the wave length, and the amplitudes decrease as the permeability of the barriers increases. The features of the parameters of engineering interest discovered for one or two thin permeable vertical barriers are considered to be meaningful and are hoped to be used to predict wave effects on barrier(s) in practical applications. 5.2 Recommendations for further study Only the numerical model for the one-barrier problem has been validated against previous numerical results for some limiting conditions. At present no results for the two-barrier problem are available for comparision. To further examine the effectiveness of the numerical models, a wide range of experiments need to be carried out. The present approach assumes that the barrier is thin and that the fluid velocities along both the barrier surfaces are the same. However, the thickness of the barrier affecting the fluid mechanics within the barrier needs to be investigated more properly. Finally, further studies in this area are needed to develop numerical models of a single wave barrier under random wave conditions, oblique waves and combined wave-current Chapter 5 Conclusion and recommendations conditions. For two or more barriers, different drafts and barrier widths, as well permeability parameters could be taken into account. 45 References Abul-Azm, R.G. (1992). "Water diffraction through submerged breakwaters." Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, 93(6), 587-605. Chwang, A.T. (1983). "A porous-wavemaker theory." Journal of Fluid Mechanics, 132, 395-406. Dalrymple, R.A., Losada, M.A., and Martin, P.A. (1991). "Reflection and transmission from porous structures under oblique wave attack." Journal of Fluid Mechanics, 224, 625-644. Dalrymple, Robert A and Martin, P.A. (1990). "Wave diffraction through offshore breakwaters." Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, 116, 727-741. Gilman, J.F., and Nottingham, D. (1992). "Wave barriers: an environmentally benign alternative." Proceeding of the Conference of Coastal Engineering Practice '92 , ASCE, 479-486 Isaacson, M . (1996). "Simplified estimation of wave effects on vertical barriers." Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, submitted. Kondo, H. (1979). "Analysis of breakwaters having two porous walls." Coastal Structures, ASCE, 79 ,962-977. Le Mehaute, B. (1972). "Progressive wave absorber." Journal of Hydraulic Research, IAHR, 10(2), 153-169. Liu, P. L-F., and Abbaspour, M . (1982). "Wave scattering by a rigid thin barrier." Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, 108(4), 479-491. 46 Losada, I.J., Losada, M.A., and Baquerizo, A. (1993). "An analytical method to evaluate the efficiency of porous screen as wave dampers." Applied Ocean Research, 15, 207-215. Losada, I.J., Losada, M.A., and Roldan, A.J. (1992). "Propagation of oblique incident waves past rigid vertical thin barriers." Applied Ocean Research, 14, 191-199. Macaskill, C. (1979). "Reflection of water waves by a permeable barrier." Journal of Fluid Mechanics, 95, 141-157. Madsen, O.S. (1974). "Wave transmission through porous structures." Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, 100, 169-188. Madsen, P.A. (1983). "Wave reflection from a vertical permeable wave absorber." Coastal Engineering, 7, 381-396. Mallayachari, V., and Sundar, V. (1994). "Reflection characteristics of permeable seawalls." Coastal Engineering, 23, 135-150. Mei, C.C.(1985). "Resonant reflection of surface water waves by periodic sand bars." Journal of Fluid Mechanics, 152, 315-335. Naciri, M. , and Mei, C C . (1988). "Bragg scattering of water waves by a doubly periodic seabed." Journal of Fluid Mechanics, 192, 51-74. Sollitt. C.K., and Cross, R.H. (1972). "Wave transmission through permeable breakwaters." Proceedings of the 13th Coastal Engineering Conference, ASCE, Vancouver, B.C, Canada, 1827-1846. Twu, S.W., and Lin, D.T. (1990). "On highly effective wave absorber." Coastal Engineering, 15, 389-405. 47 Twu, S.W., and Lin, D.T. (1990). "Wave reflection by a number of thin porous plates fixed in a semi-infinitely long flume." Proceedings of the 22nd Coastal Engineering Conference, ASCE, Delft, The Netherlands, 1046-1059. Yu, Xiping (1995). "Diffraction of water waves by porous breakwaters." Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, 121, 275-282. 48 wave direction region 1 „ t x j region 2 Figure 1 Definition sketch of one-barrier problem. |«i X ^ wave direction . x z . barrier 1 «—X H barrier 2 region 1 region 2 region 3 Figure 2 Definition sketch of two-barrier problem. 49 Figure 3 Transmission and reflection coefficients shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). , K r ; , IQ. 50 1.0 0.8 0.6 0.4 \ / s=( S = ( 0.0,0.0) 0.5,0.0) 2.0,0.0) 0.5,0.5) 2.0,2.0) J -\ \ \\ \ \\\ \ V S-( S = ( — - S = ( \ \ \ \ \ 1 \ 1 \ \ ' \ \ \ \ \ \ X \ \ • \ 1 1 I 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 4 Transmission coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). 51 1.0 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 5 Reflection coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). 52 1.0 0.8 0.6 K t 0.4 0.0 - ~ ~ -\ --f l -\ S=(0.0,0.0) S=(0.5,0.0) S=(2.0,0.0) S = (0.5,0.5) S=(2.0,2.0) J I \ -I i i 0.0 0.2 0.4 0.6 0.8 1.0 h/d Figure 6 Transmission coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 7 Reflection coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). Figure 8 Energy dissipation coefficient shown as contours in the complex S plane for the case h/d = 05, kh = 0.5 (one-barrier problem). 55 1.0 0.6 0.4 0.2 -f S = (0.0,0.0) S = (0.5,0.0) S = (2.0,0.0) -S = (0.5,0.5) S = (2.0,2.0) J -V / / , / - ¥r Vi P.i l.'i i i i l 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 9 Energy dissipation coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). 56 0.5 0.4 0.3 0.2 0.1 -f V S = (0.0,0.0) S = (0.5,0.0) S = (2.0,0.0) S = (0.5,0.5) S = (2.0,2.0) J -/ / • / - ' ' -- S i / 1 / ' / / / , / ' / ' / ,' , -7 ' S ' 1 / ' / / ' / ' ^ / / — / / / / / / / / ' / / / j/^ 1 1 1 0.0 0.2 0.4 0.6 0.8 1.0 h/d Figure 10 Energy dissipation coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). 57 Figure 11 Force coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). 58 59 60 Figure 14 Overturning moment coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). 61 62 Figure 16 Overturning moment coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). 63 Figure 17 Wave runup coefficients shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5 (one-barrier problem). , R7H (upwave); ,R+/H (downwave). 64 1.2 1.0 0.6 0.4 / / ' — ^ ~~ ; / / ,<''' _ / / / s s -=(0.0,0.0) =(0.5,0.0) =(2.0,0.0) = (0.5,0.5) =(2.0,2.0) ) — r 1 / V f// i i i s s — — s V 1 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 18 Wave runup (upwave) coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). 65 0.6 0.4 R+/H 0.2 0.0 0.0 2.0 r S = (0.0,0.0) S = (0.5,0.0) S = (2.0,0.0) S = (0.5,0.5) S = (2.0,2.0) 4.0 6.0 kh 8.0 10.0 Figure 19 Wave runup (downwave) coefficient as a function of kh for h/d = 0.5 and for various values of S (one-barrier problem). 66 1.2 R . / H 0.8 0.6 -r V S = (0.0,0.0) S=(0.5,0.0) S = (2.0,0.0) S=(0.5,0.5) S=(2.0,2.0) J -/ -- -I I I 0.0 0.2 0.4 0.6 0.8 1.0 h/d Figure 20 Wave runup (upwave) coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). 67 0.6 0.4 R + / H 0.2 0.0 ^ \ — — " — — ' — -f S=(0.0,0.0) """" \ — "' S=(0.5,0.0) S=(2.0,0.0) S = (0.5,0.5) V 1 S = (2.0,2.0) ) 1 1 1 1 0.0 0.2 0.4 0.6 0.8 1.0 h/d Figure 21 Wave runup (downwave) coefficient as a function of h/d for kd = 1.0 and for various values of S (one-barrier problem). 68 Figure 22 Transmission and reflection coefficients shown as contours in the complex S plane for the case h/d = 05, kh = 0.5, X/d = 1.0 (two-barrier problem). , K r ; , Kt'. 69 1.0 0.6 K, 0.4 0.2 0.0 - 4 l l \ f _ _ _ _ _ _ V S=(0.0,0.0) S=(0.5,0.0) S=(2.0,0.0) S=(0.5,0.5) S=(2.0,2.0) J \ \ l\ i-\ \ -f, \ -\ I I / . i -lj \ \ I 1 / 1 \ ' M n / \ v < M w r i 11 ~ ; \ / . •'. - ;<\ /i 'i ii ; 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 23 Transmission coefficient as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S (two-barrier problem). 70 1.0 0.8 0.6 0.4 0.0 1 1 * 4 J. ? \ i\ | <l /: -I -, / ! . ' ;; ' 1 i 1/ , f> • ft /: ' i n ' i • i ' h l\ /' • //I1 1 V / \ ' \ i: / l ( V 1, i ; 7 -ji ,i ^ -, li ; ki-ll i II '; | \ ! ' i / f i 1 -I 1 i ; i" r \ir I i V 1 S = (0.0,0.0) S=(0.5,0.0) S=(2.0,0.0) S=(0.5,0.5) S = (2.0,2.0) J 1 If ') 1 i i 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 24 Reflection coefficient as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S (two-barrier problem). 71 1.0 0.8 K, 0.6 •-• m / \\ // \ / / \ r m \ f 1 N \ s f ^ 1 s V / \ 7 N N V ^ I ^ \' 1 ^ \ J 1 V V / s V / S = (0.0,0.0) S = (0.5,0.0) S = (2.0,0.0) S = (0.5,0.5) — - S = (2.0,2.0) L , — ) 1 I 0.0 2.0 4.0 6.0 8.0 10.0 X/d Figure 25 Transmission coefficient as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). 72 Figure 26 Reflection coefficient as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). 73 Figure 27 Energy dissipation coefficient shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5, X/d = 1.0 (two-barrier problem). 1.0 0.6 0.4 0.2 S=(0.5,0.0) S=(2.0,0.0) S=(0.5,0.5) S=(2.0,2.0) ' ) -n J 1 ' ' i r 1 i Mi ' \ i 1, 1 ri I. ii '/1\ \ ]\ } , 1 i.'ii >i J Ii ! 1! it t! i 1 i' * <\ r i i •,'\ \ i 1 \ '. M M !i 1 i 1 1 i •\ V. V, ' V, i \[ v! t ii i! //' • ii ' 1 i .1 f i i { t i 1 0.0 2.0 4.0 6.0 8.0 10.0 kh Figure 28 Energy dissipation coefficient as a function of kfi for h/d = 05, X/d and for various values of S (two-barrier problem). 0.6 0.4 0.0 » i -: • . ' ' x ', • ./ \ \ -/ > . * \ / ' ; X J ' ' ~ \ ' V r / 1 1 w /a \ < ,/ ' \ ' ' ' ~ \ '. ' \ ' ' \ ', \ ' \ ', \ i \ 1 i 1 / / >/ / \ , : ' \ • < i \ ', ;/ \ . ,/ "\ / \ ; / \ i V i \ / W \ ! \ 1 !, -Y \ ^ ^ \ y S = (0.5,0.0) S = (2.0,0.0) S = (0.5,0.5) S = (2.0,2.0) J — 1 1 1 0.0 2.0 4.0 6.0 8.0 10.0 X/d Figure 29 Energy dissipation coefficient as a function of X/d for h/d = 05, kd = 1.0 and for various values of S (two-barrier problem). 76 77 Figure 31 Force coefficients of the downwave barrier shown as contours in the complex S plane for the case h/d = 0.5, kh = 0.5, X/d = 1.0 (two-barrier problem). • 78 1.0 10.0 Figure 32 Force coefficient of the upwave barrier as a function of kh for h/d = 0.5, X/d = 1.0 and for various values of S (two-barrier problem). 79 0.8 0.6 0.2 0.0 _ . _ 5=(0.0,0.0) a=(0.5,0.0) s=(2.0,0.0) 3=(0.5,0.5) 3 = (2.0,2.0) J --r V 1 i 0.0 2.0 4.0 6.0 8.0 10.0 Figure 33 Force coefficient of the downwave barrier as a function of kh for h/d = 0.5, A/d = 1.0 and for various values of S (two-barrier problem). 80 0.6 0.4 0.2 S = (0.0,0.0) S = (0.5,0.0) S = (2.0,0.0) S = (0.5,0.5) S = (2.0,2.0) ) \ r / \ , y v y ^ - — - — 1 0.0 2.0 4.0 6.0 8.0 10.0 X/6 Figure 34 Force coefficient of the upwave barrier as a function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem). 81 0.4 0.3 -r > S = (0.0,0.0) S = (0.5,0.0) S = (2.0,0.0) S = (0.5,0.5) — S = (2.0,2.0) L J - A J \ A \ J 1 1 I • - - - - - — _ ^ i 0.0 2.0 4.0 6.0 8.0 10.0 . A,/d Figure 35 Force coefficients of the downwave barrier as function of X/d for h/d = 0.5, kd = 1.0 and for various values of S (two-barrier problem).
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Numerical model of wave effects on permeable vertical...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Numerical model of wave effects on permeable vertical barriers above the seabed Yang, Gang 1996
pdf
Page Metadata
Item Metadata
Title | Numerical model of wave effects on permeable vertical barriers above the seabed |
Creator |
Yang, Gang |
Date Issued | 1996 |
Description | This thesis describes a numerical model used to predict the interaction of a regular small amplitude wave train with a thin permeable vertical barrier extending from the water surface to some distance above the seabed. The case of wave interactions with a pair of such barriers is also treated. The approach used is based on an eigenfunction expansion method and utilizes a boundary condition at the barrier surface which accounts for energy dissipation within the barrier. Comparisons of results based on the method have been carried out with those of previous numerical studies for related situations, and close agreement has been obtained in all such cases. A selection of results based on the method are presented for the transmission, reflection, and energy dissipation coefficients, the wave runup, and the maximum horizontal force on the barrier. These exhibit various features of interest which are discussed. The numerical model is written in the FORTRAN language. With appropriate input, it can be used to calculate wave interactions with one or two barriers under a wide range of conditions. It is confirmed that the results of the numerical model are consistent for various limiting conditions including those of a solid barrier, a fully transparent barrier, and a barrier extending to the seabed. Overall, the numerical model is found to provide a reasonably flexible and reliable means of predicting wave effects on thin permeable vertical barriers. |
Extent | 2834303 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
FileFormat | application/pdf |
Language | eng |
Date Available | 2009-03-09 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0050366 |
URI | http://hdl.handle.net/2429/5763 |
Degree |
Master of Applied Science - MASc |
Program |
Civil Engineering |
Affiliation |
Applied Science, Faculty of Civil Engineering, Department of |
Degree Grantor | University of British Columbia |
GraduationDate | 1997-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
AggregatedSourceRepository | DSpace |
Download
- Media
- 831-ubc_1997-0086.pdf [ 2.7MB ]
- Metadata
- JSON: 831-1.0050366.json
- JSON-LD: 831-1.0050366-ld.json
- RDF/XML (Pretty): 831-1.0050366-rdf.xml
- RDF/JSON: 831-1.0050366-rdf.json
- Turtle: 831-1.0050366-turtle.txt
- N-Triples: 831-1.0050366-rdf-ntriples.txt
- Original Record: 831-1.0050366-source.json
- Full Text
- 831-1.0050366-fulltext.txt
- Citation
- 831-1.0050366.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
https://iiif.library.ubc.ca/presentation/dsp.831.1-0050366/manifest