Next Article in Journal
Development and Validation of Overpressure Response Model in Steel Tunnels Subjected to External Explosion
Previous Article in Journal
Recommendations for Different Tasks Based on the Uniform Multimodal Joint Representation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long Wave Run-Up Resonance in a Multi-Reflection System

1
School of Mathematics and Statistics, Nanyang Normal University, Nanyang 473061, China
2
School of Mathematics and Statistics, University College Dublin, Dublin 4, Ireland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(18), 6172; https://doi.org/10.3390/app10186172
Submission received: 9 July 2020 / Revised: 10 August 2020 / Accepted: 2 September 2020 / Published: 4 September 2020
(This article belongs to the Section Applied Physics General)

Abstract

:
Wave reflection and wave trapping can lead to long wave run-up resonance. After reviewing the theory of run-up resonance in the framework of the linear shallow water equations, we perform numerical simulations of periodic waves incident on a linearly sloping beach in the framework of the nonlinear shallow water equations. Three different types of boundary conditions are tested: fully reflective boundary, relaxation zone, and influx transparent boundary. The effect of the boundary condition on wave run-up is investigated. For the fully reflective boundary condition, it is found that resonant regimes do exist for certain values of the frequency of the incoming wave, which is consistent with theoretical results. The influx transparent boundary condition does not lead to run-up resonance. Finally, by decomposing the left- and right-going waves into a multi-reflection system, we find that the relaxation zone can lead to run-up resonance depending on the length of the relaxation zone.

1. Introduction

Long wave run-up is difficult to observe in nature in real-time due to the catastrophic consequences it usually leads to. An illustration of a long wave is a tsunami, which has a wavelength far greater than the water depth and a small amplitude in the open ocean. When a tsunami reaches the nearshore region, its amplitude increases significantly, while its propagation speeds decreases. A large part of the theoretical results for run-up are based on shallow water theory. Analytical linear solutions for periodic long waves on a slope can be found, for example, in Green and Ferrers [1], Lamb [2], and Keller and Keller [3].
Tsunamis are known for their dramatic run-up height. From the data collection of post-tsunami surveys, extreme run-up values have been observed that exceed 20m and the most severe damages were not always caused by the first wave [4,5,6]. Extreme run-up values cannot be explained by classical analytical solutions. A possible explanation for these extreme run-ups is resonance. Miles [7] and Wilson [8] discussed the conditions for resonance in closed and open basins. The Helmholtz mode was also studied, which is the most important mode in natural basins and is normally observed in bays, inlets, and harbors with narrow entrance. Kajiura [9] described the possibility of resonances in bays. Rabinovich and Leviant [10] introduced the concept of shelf resonance, which is also called quarter-wavelength resonance. Neetu et al. [11] found that reflection and refraction from the nearshore topographic features can lead to the amplification of run-up values of non-leading tsunami waves. Resonance effects in the process of run-up have been studied numerically by Stefanakis et al. [12,13] and experimentally by Abcha et al. [14,15].
Obviously boundary conditions play an important role in numerical simulations: they can easily make the difference between successful and unsuccessful computations, or between fast and slow ones. There is still no clear understanding of how to impose both incoming and outgoing waves at a boundary. One generally uses an open boundary condition seaward, which is a numerical artifact that lies in a limited computational domain to solve a mathematical problem that has no boundary at all. Such open boundary conditions are usually called “radiation”, “absorbing”, “non-reflecting”, or “far-field” boundary conditions. They all allow waves to leave the truncated domain, thus avoiding spurious reflections that may pollute the solution in the interior of the computational domain of interest. They can be classified into two groups: non-reflecting boundary conditions and non-reflecting boundary layers. Non-reflecting boundary conditions absorb impinging waves on the artificial boundary. A classical example is the Sommerfeld boundary condition [16,17]. On the other hand, non-reflecting boundary layers have the property of absorbing waves that are traveling inside the layer. Open boundaries are considered in a widespread literature. Examples are Jensen [18], Carter and Merrifield [19], Ma and Madsen [20].
However, in a laboratory wave flume, not all reflections caused by a wave paddle can be eliminated. The resulting wave system is a multi-reflection system. In the framework of wave energy converters, Wei et al. [21] discussed the re-reflection effects on the wave field. In nature, waves on some topographies are easily trapped, such as in narrow bays. When the water depth changes abruptly in shelves with piece-wise geometry, a multi-reflection system can also occur. Such reflective boundary conditions are seldom discussed.
Run-up on a sloping beach with the use of an open boundary condition at the seaward boundary matches the linear analytical solution that shows no resonance at all. However, no matter if it is theory or experiment, there exist some resonant frequencies [2,13,22]. Wave trapping is one of the phenomena that lead to resonance [23,24]. Using reflective boundary conditions leads to the reflection of outgoing waves and leaves them in the computational domain, thus constructing a multi-reflection system. We discuss run-up in such a system.
In this paper, we provide the analytical solution of run-up on a bounded domain that consists of a sloping beach connected to a flat bottom (linear shallow water equations). Conditions for run-up resonance are given. Then, we investigate numerically the same problem (nonlinear shallow water equations) using three different boundary conditions: influx transparent boundary conditions, reflective boundary conditions, and relaxation zones. The outcome of the numerical simulations is shown. For the fully reflective boundary condition, it is found that resonant regimes do exist for certain values of the frequency of the incoming wave, which is consistent with theoretical results. The influx transparent boundary condition does not lead to run-up resonance. By decomposing the left- and right-going waves into a multi-reflection system, we find that the relaxation zone can lead to run-up resonance depending on the length of the relaxation zone. Finally, the applicability of the results to tsunami science is briefly discussed.

2. Materials and Methods

We consider periodic waves of angular frequency ω incident on a linearly sloping beach with angle α (see Figure 1). The sloping beach starts at a distance L from the shore.
The one-dimensional (1D) nonlinear shallow water equations (NSWE) read
η t + x ( h u ) = 0 ,
u t + u u x + g η x = 0 ,
where x is the horizontal coordinate, t is time, η ( x , t ) is the free-surface elevation, h ( x , t ) = h 0 ( x ) + η ( x , t ) is the water depth with h 0 ( x ) the still water depth, u ( x , t ) is the depth-averaged horizontal velocity, and g is the acceleration due to gravity. The origin of the coordinate system is at the shoreline. The still water depth h 0 ( x ) is equal to L tan α for L t x L and x tan α for L x 0 .

2.1. Theoretical Model for Run-Up Resonance Based on the Linear Shallow Water Equations

Resonance phenomena play an important role in run-up amplification. They can explain extreme run-up heights and later arrival of waves with larger amplitude. For the geometry of Figure 1, Billingham and King [25] used linear shallow water theory to show that a resonance occurs when 2 ω L / ( g tan α ) is a zero of the Bessel function J 0 . This resonance can also be found both in the numerical results of Stefanakis et al. [13] and in the experimental results of Ezersky et al. [22], because of the presence of a reflective boundary condition seaward.
Even though the numerical simulations shown in the next section are based on the 1D NSWE (1) and (2), we use the 1D linear shallow water equations (LSWE) to explain the resonance phenomenon analytically. The 1D LSWE read
η t + x ( h 0 u ) = 0 ,
u t + g η x = 0 .
The phase velocity of the wave is denoted by c ( x ) = ω / k ( x ) = g h 0 ( x ) .
Equations (3) and (4) can be combined into a single equation:
2 η t 2 g x h 0 η x = 0 .
For a beach with a constant slope, we look for a harmonic solution of the form [26]
η ( x , t ) = A ( x ) e i ω t .
Inserting Equation (6) into Equation (5) leads to
d 2 A d σ 2 + 1 σ d A d σ + A = 0 ,
where σ = 0 x k ( x ) d x with k ( x ) the wave number. Equation (7) is the Bessel equation of the first kind, whose solution can be represented by
A = A 1 J 0 ( σ ) + A 2 Y 0 ( σ ) ,
where A 1 and A 2 are linear coefficients, and J 0 , resp. Y 0 , are the Bessel functions of the first and second kind. As Y 0 tends to infinity at the shore, the coefficient of Y 0 must be zero in the nearshore region if we assume that the amplitude of the free-surface elevation is finite. Therefore, A is expressed only in terms of the Bessel function of the first kind:
A = R m a x J 0 ( σ ) ,
where R m a x denotes the run-up. Billingham and King [25] found the resonance mentioned above by imposing a fixed amplitude at x = L . Indeed, σ ( x = L ) = 2 ω L / ( g tan α ) .
For a constant-depth region, the wave field can be described as a sum of incident and reflected waves:
η ( x , t ) = ( A i e i k x + A r e i k x ) e i ω t ,
where A i and A r are the amplitudes of the incident and reflected waves, respectively.
Next, we discuss run-up due to the more realistic bathymetry profile that consists of a constant depth region connected to a sloping beach. Hereafter, it will be referred to as the canonical case (see Figure 1). Following Kanoğlu and Synolakis [27] and Stefanakis et al. [13], we use the continuity of the free-surface elevation and of the horizontal fluxes to match surface elevations (Equation (11)) and particle velocities (according to Equation (4), it is equivalent to Equation (12)) at the toe of the sloping beach from two adjacent segments, and match the surface elevations at the left boundary Equation (13), respectively. We obtain
J 0 ( σ 1 ) R m a x J 0 ( σ 2 ) A 1 Y 0 ( σ 2 ) A 2 = 0 ,
J 1 ( σ 1 ) R m a x J 1 ( σ 2 ) A 1 Y 1 ( σ 2 ) A 2 = 0 ,
    J 0 ( σ 3 ) A 1 + Y 0 ( σ 3 ) A 2 = A i e i k x + A r e i k x ,
where
σ 1 = 2 ω L g tan α ,
σ 2 = ω L g L tan α ,
σ 3 = ω L t g L tan α .
The determinant of the system (11)–(13) is given by
d e t ( system ) = J 0 ( σ 1 ) J 0 ( σ 2 ) Y 0 ( σ 2 ) J 1 ( σ 1 ) J 1 ( σ 2 ) Y 1 ( σ 2 ) 0 J 0 ( σ 3 ) Y 0 ( σ 3 ) ·
When the determinant vanishes, run-up becomes resonant. Therefore, the resonance condition reads d e t ( system ) = 0 . Clearly, the distance L t from the shore to the seaward boundary that appears in σ 3 plays a crucial role in the resonance condition.

2.2. Numerical Integration of the Nonlinear Shallow Water Equations—The Boundary Conditions

We integrate the 1D NSWE (1) and (2) with a finite volume solver based on the Characteristic Flux scheme with UNO2 reconstruction for higher order terms. A third-order Runge–Kutta scheme is used for time discretization. The numerical model is described and validated in Dutykh et al. [28]. What is new here is the fact that we test three types of boundary conditions described below.

2.2.1. Influx Transparent Boundary Conditions

Appropriate boundary conditions are needed to let the wave travel in and out of the domain in an undisturbed way. Influx transparent boundary conditions are boundary conditions that make waves on boundaries propagate into the computational domain and are transparent for wave propagation out of the computational domain. In the region of constant depth, the LSWE reduce to two wave Equations [2]:
η t t c 2 η x x = 0 ,
and
Φ t t c 2 Φ x x = 0 ,
where Φ is the time-integral of the displacement past the plane x, up to the time t: Φ ( x , t ) = t u d t .
Equations (18) and (19) have the general solutions
η ( x , t ) = η r g ( x c t ) + η l g ( x c + t ) ,
Φ ( x , t ) = Φ r g ( x c t ) + Φ l g ( x c + t ) ,
where the subscripts l g and r g denote the waves traveling to the left and to the right directions, respectively.
Differentiating Equation (21) with respect to x and t yields
Φ x = 1 c Φ r g + 1 c Φ l g ,
Φ t = Φ r g + Φ l g .
To build a left boundary at x = x l that is transparent for waves propagating to the left but allows waves to propagate to the right, we eliminate Φ l g :
Φ t c Φ x x = x l = 2 Φ r g x = x l .
This is called the left Influx Transparent Boundary Condition (ITBC).
As Φ t = u ,
u ( x l , t ) = c Φ x x = x l 2 Φ r g x = x l .
Moreover, from (3),
Φ x = t u x d t = 1 h 0 η = g c 2 η .
Therefore, the velocity at the left boundary can be written as
u ( x l , t ) = g η ( x l , t ) c ( x l ) + 2 g η r g ( x l , t ) c ( x l ) .
Similarly, the right ITBC ( x = x r ) that is transparent for waves propagating to the right but makes waves propagate to the left can be written as
Φ t + c Φ x x = x r = 2 Φ l g x = x r .
Therefore, the velocity at the right boundary is
u ( x r , t ) = g η ( x r , t ) c ( x r ) 2 g η l g ( x r , t ) c ( x r ) .
The head-on collision of two single waves on a flat bottom with the use of the left and right influx transparent boundary conditions is shown in Figure 2. Waves of the form η 0 ( t ) = 0.1 sech 2 ( 0.17 t 2 ) are generated at the boundaries. As the boundaries are transparent to the outgoing waves, no reflection occurs at the boundaries.

2.2.2. The Reflective Boundary Conditions

The reflective boundary conditions are boundary conditions that are prescribed on a bounded interval [ x l , x r ] and read
η x ( x l , t ) = η x ( x r , t ) = 0 ,
u ( x l , t ) = u ( x r , t ) = 0 .
For the classical NSWE system, we only need to impose the boundary conditions in one of the two dependent variables [29]. Therefore, it is sufficient to take u ( x l , t ) = u ( x r , t ) = 0 as reflective boundary condition.
In this paper, analogous reflective boundary conditions are considered. The incident wave η ( x l , t ) is coming from the left boundary, where the velocity can be expressed by
u ( x l , t ) = c ( x l , t ) + 1 h ( x l , t ) H ( x l , t ) c ( x l , t ) ,
where h ( x , t ) is the total depth of water and H ( x l , t ) = h 0 + η i ( t ) with η i ( t ) the free-surface elevation of the incident wave at the left boundary. This boundary condition was used by Stefanakis et al. [12]. A vertical wall is used at the right boundary, which means u ( x r , t ) = 0 . Figure 3 shows an incident single wave generated at the left boundary that propagates to the right and meets a vertical wall at the right boundary, where it is reflected. When the reflected wave propagates back towards the left boundary, it reflects again and the wave polarity is reversed. It makes a multi-reflection system.

2.2.3. The Relaxation Zone

The relaxation zone is a kind of non-reflecting boundary layer. It absorbs all the reflected waves in a region that lies at the inlet boundary. It is usually one wavelength long [26,30]. In practice, it can modify an incident wave field to any desired target function. This also extends the applicability of the sponge layer to inlet boundaries, as the target function could be that of the generated wave train. The relaxation zone works as a blend between a target solution and the computed solution. The modified velocity can be described by
u = c r e l a x u t a r g e t + ( 1 c r e l a x ) u m o d e l ,
where u t a r g e t is a target velocity, u m o d e l is the unaltered computed velocity, and c r e l a x is the relaxation coefficient. This technique has been discussed in many papers, such as Mayer et al. [31,32,33,34,35]. As the relaxation is applied explicitly as a postprocessing step, it is possible to apply it at the inlet. If there is no proper relaxation, the waves may propagate upstream toward the inlet boundary, be reflected back into the computational domain and pollute the flow field.
To obtain a smoothly modified internal flow field, it is required that “ c r e l a x ” and its spatial derivatives are equal to zero along the interface connecting the relaxation zone and the internal field of the computational domain. There are many options to build a proper function of the relaxation coefficient. In this paper, we use the form
c r e l a x ( β ) = b β 3 + ( 1 b ) β 6 ,
where b [ 0 , 2 ] is a free parameter, and β = 1 ( x x s t a r t ) / ( x e n d x s t a r t ) ( x s t a r t and x e n d are coordinates of the relaxation zone boundaries) is the nondimensional distance function.
The choice of the relaxation coefficient c r e l a x depends strongly on the problem and the length of the relaxation zone. What is important is that the choice is not unique for a given problem.
In numerical simulations, reflected waves from the beach will interact with the wave generation boundary, which will pollute the final results, especially for breaking waves. Therefore, we can test the applicability of the relaxation zone through the modeling of standing waves. If it can absorb reflected waves in front of the inlet, the incident wave profile is kept in the relaxation zone and the standing wave pattern is simulated in the non-relaxed part of the computational zone.
The standing wave is modeled as sketched in Figure 4. The relaxation zone at the left boundary is one wavelength long. The total horizontal length of the domain is three wavelengths. The incident wave amplitude is η 0 = 1 m, the wave period is T = 40 s, and the water depth is h 0 = 20 m. One can see that the relaxation zone can simulate a standing wave pattern well in the non-relaxed part of the computational domain.
The relaxation zone makes the boundary as “open”. Compared with the ITBC, which can only be applied to regular waves, it can be wildly used for any type of wave, regular or irregular, long or short.

3. Results for Waves Incident on a Linearly Sloping Beach

Having described the three types of boundary conditions, we now present their effect on the numerical simulations of waves incident on a linearly sloping beach. Periodic monochromatic incident waves η ( x l , t ) = η 0 sin ( ω t ) propagate towards a sloping beach, with sloping angle tan α = 0.02 (unless otherwise specified) and horizontal length L = 5000 m (unless otherwise specified). The value tan α = 0.02 is typical of nearshore conditions as seen in [36].

3.1. Results with the Influx Transparent Boundary Conditions

The left ITBC is used on the left boundary. Non-dimensional run-up with respect to nondimensional wavelength is shown in Figure 5. One can see that the numerical nondimensional run-up matches the results given by the analytical linear solution before wave breaking. The ITBC is only approximate when the wavelength is too large as can be seen in Figure 5 for λ 0 / L > 3 [23,37,38]. The left boundary is not entirely transparent to outgoing waves when the wavelength is large. The outgoing waves are partially reflected and trapped in the computational domain. However, the effect on run-up is slight.

3.2. Results with the Reflective Boundary Conditions

We construct a multi-reflection system for the canonical case of Figure 1, where the reflective boundary conditions are used. Waves generated at the left seaward boundary propagate forward and are reflected at the right boundary. The reflected waves propagate back to the left boundary and are re-reflected, then the re-reflected waves propagate forward again and the process is repeated. A reflection cycle is defined as a cycle during which the wave travels the total domain back and forth. The time it takes is called the reflection period T r e f l e c t . For the canonical case of Figure 1, it can be expressed by
T r e f l e c t = 2 L t L g h 0 + 2 L g h 0 .
In the simulations below, the total simulation time is always four reflection cycles, and the incident periodic wave at x = L t is expressed by
η i ( t ) = η 0 sin 2 π t T ,
where T is the wave period and η 0 = 0.1 m.
First, we considered three combinations of α , L, and L t , thus fixing T r e f l e c t for each combination. We tested the maximum run-up in a multi-reflection system for three cases by varying the nondimensional wave period within the range 0.25 < T / T r e f l e c t < 0.65 (see Figure 6). Resonance effects are clearly visible. The maximum nondimensional run-ups are sharply amplified for certain wave periods, which lead to three peaks in Figure 6a, i.e., the maximum nondimensional run-up is up to 27 when tan α = 0.02, L t = 8000 m, L = 5000 m, and T = 220.8 s ( T / T r e f l e c t = 0.2658). Second, we fixed α , L, and T instead of L t , the horizontal distance from the left boundary to the undisturbed shoreline. This time, we tested the maximum run-up as a function of L t , which in turn made T r e f l e c t and T / T r e f l e c t vary through Equation (35). The results are shown in Figure 7. Waves with the same period propagating on the same bathymetry ( tan α = 0.02, L = 2000 m) lead to different maximum run-ups with respect to various flat bottom lengths, which indicates that the wave behavior is affected by the location of the left boundary. This is different from the open boundary. However, resonance occurs when T / T r e f l e c t = 0.2658 , 0.3595 in the range of 0.25 < T / T r e f l e c t < 0.4 . From these two figures, monochromatic forcing is found to lead to resonant run-up by non-leading waves when the nondimensional wave periods T / T r e f l e c t are the roots of the determinant of system (17), which coincides with the resonance condition we gave above.

3.3. Results with the Relaxation Zone

Next, we investigate the wave behavior in the presence of a relaxation zone. Its coefficients are given in the equation for c r e l a x (34). However, we use the parameters b = 2 and β = ( x e n d x ) / λ 0 , where λ 0 is the wavelength of the incoming wave. When the full relaxation zone is applied, the value of “ c r e l a x ” decreases from 1 to 0 along the relaxation zone. When the relaxation zone is only partial, the value of “ c r e l a x ” decreases from a value that depends on the length of the zone (value less than 1) to 0.
Results are shown in Figure 8. Two different relaxation zones are considered. To emphasize the effect of the relaxation zone, we also show the free-surface elevation time series with the use of the other boundary conditions described above (reflective and left ITBC). One can observe that the partial relaxation zone only alters the wave amplitude, but maintains the other wave characteristics. In other words, it can be considered as a partially reflective boundary condition. When the length of the relaxation zone is one fifth of the incident wavelength, the results coincide with those of obtained with the left ITBC, which means that there is no reflection caused at the left boundary.

4. Conditions for Run-Up Resonance in a Multi-Reflection System

To investigate the influence of the boundary conditions on the resonance mechanism, we first consider the reflective boundary conditions. In Appendix A, we describe a technique to decompose a wave into left- and right-going waves in a multi-reflection system. Using that technique, we separated the recorded waves at a point x located 2000 m to the right of the left boundary, into left- and right-going waves with L t = 7000 , 8000 , 9000 m. The other parameters were fixed at L = 2000 m, tan α = 0.02 , T = 268.5 s, and η 0 = 0.1 m (see Figure 9). The amplitudes of left- and right-going waves are constant in each reflection cycle. When L t = 8000 m, the amplitude of right-going waves is the sum of amplitudes of the initial incident wave and of the left-going wave in the preceding reflection cycle, which are not observed in the other two cases (see Figure 9).
The wave motion could develop as a standing wave pattern with the use of open boundary condition [26]. For a multi-reflection system, it can be considered as a superposition of a standing wave and traveling waves. When the phase difference between the right-going and left-going waves is π , resonance can occur. When L t = 8000 m, one can observe that the phase difference between the right-going and left-going waves is close to π (see Figure 9b,c), which leads to resonance and amplifies the free surface elevation significantly (see Figure 9a). The period T = 268.5 s ( T / T r e f l e c t = 0.2658 ) used in Figure 9 is the resonant frequency that leads to the first peak in Figure 6a. Therefore, det(system) = 0 gives an analytical resonance condition, and Figure 9 reveals a wave motion resonance condition. We can conclude that the phase difference between the right-going and left-going waves must be π for the waves of the periods at the run-up peaks and det(system) = 0 in Figure 7.
Due to the new incident component coming from the reflected component in the previous reflection cycle, the reflection coefficient r l ( m ) in the mth reflection cycle at the left boundary can be defined by
r l ( m ) = A I ( 1 , m ) A R ( 1 , m 1 ) ,
where m > 1 . The average value is considered as the reflection coefficient of the left boundary.
Table 1 and Table 2 list reflection coefficients corresponding to cases of resonant and non-resonant regimes for the three types of boundary conditions. The amplitudes of the harmonic components being small, tiny errors may lead to deviations, and the reflection coefficients have small discrepancies even with the use of the same boundary condition. The reflection coefficient of the reflective boundary condition is up to 95 % , it can be regarded as “full reflection”. With the use of a partial relaxation zone ( L r e l a x = 0.2 λ 0 ) and the left ITBC, amplitudes of the harmonic components are equal to zero when m > 1 , which indicates that there is no reflection at the left. The left boundary is open.
In Figure 10, we tested the maximum run-up as a function of wave period and length of the partial relaxation zone, respectively. It can be observed that the maximum run-ups of the resonant regimes are determined by the length of the partial relaxation zone, which are sharply reduced by 40 % when L r e l a x = 0.1 λ 0 and hold constant when L r e l a x > 0.18 λ 0 . It indicates that the relaxation zone does not require a full wavelength to behave as an open boundary.

5. Conclusions

Long wave run-up resonance in multi-reflection systems is studied in the paper. A simplified geometry involving long periodic waves incident on a linearly sloping beach is considered. By decomposing the surface wave into right- and left-going waves, we are able to obtain a better understanding of the conditions leading to resonance. It is found that wave trapping and wave–wave interactions can lead to run-up resonance. Resonance occurs when a condition obtained by linearizing the equations is satisfied. This condition is written as a determinant that must vanish. This determinant involves the angle of the slope, the length of the slope and the location of the artificial seaward boundary where boundary conditions are applied. The maximum run-up for the resonant regime is largely determined by the intensity of wave trapping.
Meanwhile, the fully and partially reflective boundary conditions are discussed. The reflection coefficient is determined by the length of the partial relaxation zone. The partial relaxation zone can also behave as an open boundary, even if it is less than one wavelength long.
The present study is mathematical. What it says is that the wrong choice of boundary conditions can lead to a misinterpretation of what happens in the real-world. However, there are several physical phenomena that can lead to wave trapping, and therefore to run-up resonance if the bathymetry has the right properties. In recent tsunamis, areas that were close to each other have experienced different run-ups, sometimes with an order of magnitude difference. The results in the present paper could provide a partial explanation.

Author Contributions

S.X. coordinated the scientific effort together with F.D. S.X. performed all numerical simulations. The overall supervision was provided by F.D. Both authors participated in the interpretation of results and the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Nanyang Normal University grant number 2018ZX016.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript.
ITBCInflux Transparent Boundary Condition
SWEShallow Water Equations

Appendix A. Decomposition of Left- and Right-Going Waves in a Multi-Reflection System

A multi-reflection system can be regarded as a superposition of waves propagating in the positive and negative directions of the coordinate x. We call them the right- and left-going waves, respectively. To understand the characteristics of the wave field in the computational domain, decomposing the composite water surface elevation into the right- and left-going waves is often required. The reconstruction of surface elevation by superposition of harmonic waves can be found, for example, in Wei et al. [21,39,40,41]. The wave profile in a multi-reflection system can be expressed by
η ( x , t ) = m 1 A I ( 1 , m ) cos k x ω t + κ I ( 1 , m ) + m 1 A R ( 1 , m ) cos k x + ω t + κ R ( 1 , m ) + m 1 n 2 A I , B ( n , m ) cos n k x n ω t + κ I , B ( n , m ) + m 1 n 2 A R , B ( n ) cos n k x + n ω t + κ R , B ( n , m ) + m 1 n 2 A I , F ( n , m ) cos k ( n ) x n ω t + κ I , F ( n , m ) + m 1 n 2 A R , F ( n , m ) cos k ( n ) x + n ω t + κ R , F ( n , m )
where A is the amplitude and κ is the phase difference. The subscripts I and R denote the incident and reflected waves, respectively. The second subscripts B and F denote the bounded and free waves, respectively. The superscript n denotes the nth harmonic, and the superscript m denotes the mth reflection cycle.
The waves are decomposed into the incident and reflected waves in terms of the harmonic bound and free waves. The bound waves are in phase-locked modes that travel at the same phase velocity as the first harmonic waves and are bound to them. The free waves propagate at their respective phase velocities, which can be presented by the dispersion relation
( n ω ) 2 = g k ( n ) tanh k ( n ) h 0 .
Let us verify the reconstruction of the wave motion provided by Equation (A1) on an example. Considering the canonical problem shown in Figure 1, together with the use of reflective boundary conditions, we reconstruct the recorded waves on the flat segment at x = 4000 m for a resonant and a non-resonant regimes (see Figure A1). A good agreement is found between the reconstructed and recorded waves. A new component is added in each reflection cycle. A jump appears in the time series of the free-surface elevation when it arrives (see Figure A1b). The jumps can vanish when we smooth the amplitude numerically. As the result shown in Figure A1 proves that Equation (A1) can describe the multi-reflection system, we consider that the smoothing process can be ignored.
Figure A1. A comparison between the recorded and reconstructed wave at x = 4000 m for the canonical case with the use of reflective boundary conditions, where L t = 8000 m, L = 2000 m, tan α = 0.02 , and η 0 = 0.1 m. The reconstructed wave is given by Equation (A1), where n = 1 and m = 4 . (a) T = 268.5 s ( T / T r e f l e c t = 0.2658 ); (b) T = 299.2 s ( T / T r e f l e c t = 0.2962 ).
Figure A1. A comparison between the recorded and reconstructed wave at x = 4000 m for the canonical case with the use of reflective boundary conditions, where L t = 8000 m, L = 2000 m, tan α = 0.02 , and η 0 = 0.1 m. The reconstructed wave is given by Equation (A1), where n = 1 and m = 4 . (a) T = 268.5 s ( T / T r e f l e c t = 0.2658 ); (b) T = 299.2 s ( T / T r e f l e c t = 0.2962 ).
Applsci 10 06172 g0a1aApplsci 10 06172 g0a1b

References

  1. Green, G.; Ferrers, N. On the motion of waves in a variable canal of small depth and width. Trans. Camb. Philos. Soc. 1838, 6, 457. [Google Scholar]
  2. Lamb, H. Hydrodynamics; Cambridge University Press: Cambridge, UK, 1932. [Google Scholar]
  3. Keller, J.; Keller, H. Water Wave Run-Up on a Beach; Technical Report; Department of the Navy: Washington, DC, USA, 1964.
  4. Liu, P.; Lynett, P.; Fernando, H.; Jaffe, B.; Fritz, H.; Higman, B.; Morton, R.; Goff, J.; Synolakis, C. Observations by the international tsunami survey team in Sri Lanka. Science 2005, 308, 1595. [Google Scholar] [CrossRef]
  5. Fritz, H.M.; Kongko, W.; Moore, A.; McAdoo, B.; Goff, J.; Harbitz, C.; Uslu, B.; Kalligeris, N.; Suteja, D.; Kalsum, K.; et al. Extreme runup from the 17 July 2006 Java tsunami. Geophys. Res. Lett. 2007, 34, L12602. [Google Scholar] [CrossRef]
  6. Mori, N.; Takahashi, T.; Yasuda, T.; Yanagisawa, H. Survey of 2011 Tohoku earthquake tsunami inundation and run-up. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef]
  7. Miles, J. Resonant response of harbours: An equivalent-circuit analysis. J. Fluid Mech. 1971, 46, 241–265. [Google Scholar] [CrossRef]
  8. Wilson, B. Seiches. Adv. Hydrosci. 1972, 8, 1–94. [Google Scholar]
  9. Kajiura, K. Local behaviour of tsunamis. In Waves on Water of Variable Depth; Lecture Notes in Physics; Provis, D., Radok, R., Eds.; Springer: Berlin/Heidelberg, Germany, 1977; Volume 64, pp. 72–79. [Google Scholar]
  10. Rabinovich, A.; Leviant, A. Influence of seiche oscillations on the formation of the long-wave spectrum near the coast of the Southern Kuriles. Oceanology 1992, 32, 17–23. [Google Scholar]
  11. Neetu, S.; Suresh, I.; Shankar, R.; Nagarajan, B.; Sharma, R.; Shenoi, S.S.C.; Unnikrishnan, A.S.; Sundar, D. Trapped waves of the 27 November 1945 Makran tsunami: Observations and numerical modeling. Nat. Hazards 2011, 59, 1609–1618. [Google Scholar] [CrossRef]
  12. Stefanakis, T.S.; Dias, F.; Dutykh, D. Local run-up amplification by resonant wave interactions. Phys. Rev. Lett. 2011, 107, 124502. [Google Scholar] [CrossRef] [Green Version]
  13. Stefanakis, T.; Xu, S.; Dutykh, D.; Dias, F. Run-up amplification of transient long waves. Q. Appl. Math. 2015, 73, 177–199. [Google Scholar] [CrossRef] [Green Version]
  14. Abcha, N.; Ezersky, A.; Pelinovsky, E. Physical modeling of resonance phenomena in the long wave dynamics. La Houille Blanche 2016, 2016, 57–64. [Google Scholar] [CrossRef] [Green Version]
  15. Abcha, N.; Pelinovsky, E.; Didenkulova, I. Laboratory modeling of resonance phenomena in the long wave dynamics. La Houille Blanche 2018, 2016, 58–65. [Google Scholar]
  16. Sommerfeld, A. Partial Differential Equations in Physics; Academic Press: Cambridge, MA, USA, 1949; Volume 1. [Google Scholar]
  17. Romate, J. Absorbing boundary conditions for free surface waves. J. Comput. Phys. 1992, 99, 135–145. [Google Scholar] [CrossRef] [Green Version]
  18. Jensen, T. Open boundary conditions in stratified ocean models. J. Mar. Syst. 1998, 16, 297–322. [Google Scholar] [CrossRef]
  19. Carter, G.; Merrifield, M. Open boundary conditions for regional tidal simulations. Ocean Model. 2007, 18, 194–209. [Google Scholar] [CrossRef]
  20. Ma, P.; Madsen, O.S. An open boundary condition for numerical coastal circulation models. J. Phys. Oceanogr. 2011, 41, 2363–2380. [Google Scholar] [CrossRef]
  21. Wei, Y.; Abadie, T.; Henry, A.; Dias, F. Wave interaction with an Oscillating Wave Surge Converter. Part II: Slamming. Ocean Eng. 2016, 113, 319–334. [Google Scholar] [CrossRef]
  22. Ezersky, A.; Abcha, N.; Pelinovsky, E. Physical simulation of resonant wave run-up on a beach. Nonlinear Proc. Geoph. 2013, 20, 35–40. [Google Scholar] [CrossRef] [Green Version]
  23. Agnon, Y.; Mei, C. Trapping and resonance of long shelf waves due to groups of short waves. J. Fluid Mech. 1988, 195, 201–221. [Google Scholar] [CrossRef]
  24. Grataloup, G.; Mei, C. Localization of harmonics generated in nonlinear shallow water waves. Phys. Rev. E 2003, 68, 026314. [Google Scholar] [CrossRef]
  25. Billingham, J.; King, A. Wave Motion; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  26. Madsen, P.; Fuhrman, D. Run-up of tsunamis and long waves in terms of surf-similarity. Coast. Eng. 2008, 55, 209–223. [Google Scholar] [CrossRef]
  27. Kânoğlu, U.; Synolakis, C. Long wave runup on piecewise linear topographies. J. Fluid Mech. 1998, 374, 1–28. [Google Scholar] [CrossRef]
  28. Dutykh, D.; Katsaounis, T.; Mitsotakis, D. Finite volume schemes for dispersive wave propagation and runup. J. Comput. Phys. 2011, 230, 3035–3061. [Google Scholar] [CrossRef] [Green Version]
  29. Fokas, A.; Pelloni, B. Boundary value problems for Boussinesq type systems. Math. Phys. Anal. Geom. 2005, 8, 59–96. [Google Scholar] [CrossRef]
  30. Bingham, H.; Agnon, Y. A Fourier—Boussinesq method for nonlinear water waves. Eur. J. Mech. B/Fluids 2005, 24, 255–274. [Google Scholar] [CrossRef]
  31. Mayer, S.; Garapon, A.; Soerensen, L. A fractional step method for unsteady free-surface flow with applications to non-linear wave dynamics. Int. J. Numer. Meth. Fluids 1998, 28, 293–315. [Google Scholar] [CrossRef]
  32. Fuhrman, D. Numerical Solutions of Boussinesq Equations for Fully Nonlinear and Extremely Dispersive Water Waves. Ph.D. Thesis, Technical University of Denmark, Department of Mechanical Engineering, Maritime Enineering, Lyngby, Denmark, 2004. [Google Scholar]
  33. Engsig-Karup, A. Unstructured Nodal DG-FEM Solution of High-Order Boussinesq-Type Equations. Ph.D. Thesis, Technical University of Denmark, Lyngby, Denmark, 2006. [Google Scholar]
  34. Jacobsen, N. A Full Hydro-and Morphodynamic Description of Breaker Bar Development. Ph.D. Thesis, Technical University of Denmark, Lyngby, Denmark, 2011. [Google Scholar]
  35. Seng, S. Slamming and Whipping Analysis of Ships. Ph.D. Thesis, Department of Naval Architecture and Offshore EngineeringInstitut for Skibs-og Havteknik, Technical University of Denmark Tekniske Universitet, Lyngby, Denmark, 2012. [Google Scholar]
  36. Athanasiou, P.; van Dongeren, A.; Giardino, A.; Vousdoukas, M.; Gaytan-Aguilar, S.; Ranasinghe, R. Global distribution of nearshore slopes with implications for coastal retreat. Earth Syst. Sci. Data 2019, 11, 1515–1529. [Google Scholar] [CrossRef] [Green Version]
  37. Sohankar, A.; Norberg, C.; Davidson, L. Low-Reynolds-number flow around a square cylinder at incidence: Study of blockage, onset of vortex shedding and outlet boundary condition. Int. J. Numer. Methods Fluids 1998, 26, 39–56. [Google Scholar] [CrossRef]
  38. Espinoza, H.; Codina, R.; Badia, S. A Sommerfeld non-reflecting boundary condition for the wave equation in mixed form. Comput. Methods Appl. Mech. Eng. 2014, 276, 122–148. [Google Scholar] [CrossRef]
  39. Goda, Y.; Suzuki, T. Estimation of incident and reflected waves in random wave experiments. Coast. Eng. Proc. 1976, 1, 828–845. [Google Scholar] [CrossRef] [Green Version]
  40. Mansard, E.; Funke, E. The measurement of incident and reflected spectra using a least squares method. Coast. Eng. Proc. 1980, 1, 154–172. [Google Scholar]
  41. Lin, C.; Huang, C. Decomposition of incident and reflected higher harmonic waves using four wave gauges. Coast. Eng. 2004, 51, 395–406. [Google Scholar] [CrossRef]
Figure 1. The geometry of a sloping beach connected to a region of constant depth. The still water depth is given by h 0 ( x ) , with h 0 ( x ) = h 0 = L tan α for L t x L and h 0 ( x ) = x tan α for L x 0 .
Figure 1. The geometry of a sloping beach connected to a region of constant depth. The still water depth is given by h 0 ( x ) , with h 0 ( x ) = h 0 = L tan α for L t x L and h 0 ( x ) = x tan α for L x 0 .
Applsci 10 06172 g001
Figure 2. Head-on collision of two single waves on a flat bottom with the use of the left and right influx transparent boundary conditions, where profiles of incident waves both have the form η 0 ( t ) = 0.1 sech 2 ( 0.17 t 2 ) . The red arrow shows the direction of wave propagation.
Figure 2. Head-on collision of two single waves on a flat bottom with the use of the left and right influx transparent boundary conditions, where profiles of incident waves both have the form η 0 ( t ) = 0.1 sech 2 ( 0.17 t 2 ) . The red arrow shows the direction of wave propagation.
Applsci 10 06172 g002
Figure 3. Free surface elevation of a single wave on a flat bottom with the use of the reflective boundary conditions, when the space interval is [ 1000 , 400 ] and a wall is located at x = 400 . The incident wave comes from the left, whose profile has the form η 0 ( t ) = 0.1 sech 2 ( 0.17 t 2 ) . The red arrow shows the direction of wave propagation.
Figure 3. Free surface elevation of a single wave on a flat bottom with the use of the reflective boundary conditions, when the space interval is [ 1000 , 400 ] and a wall is located at x = 400 . The incident wave comes from the left, whose profile has the form η 0 ( t ) = 0.1 sech 2 ( 0.17 t 2 ) . The red arrow shows the direction of wave propagation.
Applsci 10 06172 g003
Figure 4. (a) The geometry of a flat bottom connected to a vertical wall. (b) The modeled numerical solution to standing wave problem on the geometry of panel (a).
Figure 4. (a) The geometry of a flat bottom connected to a vertical wall. (b) The modeled numerical solution to standing wave problem on the geometry of panel (a).
Applsci 10 06172 g004
Figure 5. Maximum nondimensional run-up ratio as a function of nondimensional wavelength on a sloping beach. Line with stars: numerical solution of the nonlinear shallow water equations with the use of the left ITBC; solid line: analytical linear solution; dashed line: wave breaking criterion. The incident wave is a periodic monochromatic wave.
Figure 5. Maximum nondimensional run-up ratio as a function of nondimensional wavelength on a sloping beach. Line with stars: numerical solution of the nonlinear shallow water equations with the use of the left ITBC; solid line: analytical linear solution; dashed line: wave breaking criterion. The incident wave is a periodic monochromatic wave.
Applsci 10 06172 g005
Figure 6. (a) Maximum nondimensional run-up of monochromatic waves as a function of nondimensional wave periods for the canonical case with the use of reflective boundary conditions. (b) Plot of the determinant of the system (17). The circles show nondimensional resonant periods, which correspond to the zeros of the determinant.
Figure 6. (a) Maximum nondimensional run-up of monochromatic waves as a function of nondimensional wave periods for the canonical case with the use of reflective boundary conditions. (b) Plot of the determinant of the system (17). The circles show nondimensional resonant periods, which correspond to the zeros of the determinant.
Applsci 10 06172 g006
Figure 7. (a) Maximum nondimensional run-up as a function of L t for the canonical case with the use of the reflective boundary condition when tan α = 0.02 , L = 2000 m and T = 268.5 s. Note that we still use T / T r e f l e c t along the horizontal axis as L t makes T r e f l e c t and T / T r e f l e c t vary through Equation (35). (b) Plot of the corresponding determinant of the system (17). The circles show the zeros of the determinant, which indicate nondimensional resonant periods.
Figure 7. (a) Maximum nondimensional run-up as a function of L t for the canonical case with the use of the reflective boundary condition when tan α = 0.02 , L = 2000 m and T = 268.5 s. Note that we still use T / T r e f l e c t along the horizontal axis as L t makes T r e f l e c t and T / T r e f l e c t vary through Equation (35). (b) Plot of the corresponding determinant of the system (17). The circles show the zeros of the determinant, which indicate nondimensional resonant periods.
Applsci 10 06172 g007
Figure 8. Free-surface elevation time series of the recorded wave at x for the canonical case, when L = 2000 m, L t = 8000 m, tan α = 0.02, η 0 = 0.1 m, and the period T = 268.5 s ( T / T r e f l e c t = 0.2658 ). The distance from x to the toe of the slope is 2000 m. The results with the relaxation zone are compared with those obtained with the other types of boundary conditions.
Figure 8. Free-surface elevation time series of the recorded wave at x for the canonical case, when L = 2000 m, L t = 8000 m, tan α = 0.02, η 0 = 0.1 m, and the period T = 268.5 s ( T / T r e f l e c t = 0.2658 ). The distance from x to the toe of the slope is 2000 m. The results with the relaxation zone are compared with those obtained with the other types of boundary conditions.
Applsci 10 06172 g008
Figure 9. Free-surface elevation time series. (a) Recorded waves; (b) Right-going waves; (c) Left-going waves at x (the distant from x to the left boundary is 2000 m) for the canonical case with the use of reflective boundary condition, when L = 2000 m, tan α = 0.02 , the period T = 268.5 s, and η 0 = 0.1 m.
Figure 9. Free-surface elevation time series. (a) Recorded waves; (b) Right-going waves; (c) Left-going waves at x (the distant from x to the left boundary is 2000 m) for the canonical case with the use of reflective boundary condition, when L = 2000 m, tan α = 0.02 , the period T = 268.5 s, and η 0 = 0.1 m.
Applsci 10 06172 g009
Figure 10. Maximum nondimensional run-up for the canonical case with the use of partial relaxation zone as a function of (a) wave period; (b) relaxation zone length, when L t = 8000 m, L = 2000 m, and tan α = 0.02 .
Figure 10. Maximum nondimensional run-up for the canonical case with the use of partial relaxation zone as a function of (a) wave period; (b) relaxation zone length, when L t = 8000 m, L = 2000 m, and tan α = 0.02 .
Applsci 10 06172 g010
Table 1. Reflection coefficients at the left boundary with four boundary conditions for the canonical case, when L t = 8000 m, L = 2000 m, tan α = 0.02 , and T = 268.5 s ( T / T r e f l e c t = 0.2658 ).
Table 1. Reflection coefficients at the left boundary with four boundary conditions for the canonical case, when L t = 8000 m, L = 2000 m, tan α = 0.02 , and T = 268.5 s ( T / T r e f l e c t = 0.2658 ).
Reflective Boundary ConditionPartial Relaxation Zone ( L relax = 0.1 λ 0 )Partial Relaxation Zone ( L relax = 0.2 λ 0 )Left ITBC
m = 1
m = 2 0.97180.7419
m = 3 0.97250.7380
m = 4 0.95000.7219
reflection coefficients0.96480.7339no reflectionno reflection
Table 2. Reflection coefficients at the left boundary with four boundary conditions for the canonical case, when L t = 7000 m, L = 200 0 m, tan α = 0.02 , and T = 268.5 s ( T / T r e f l e c t = 0.2953 ).
Table 2. Reflection coefficients at the left boundary with four boundary conditions for the canonical case, when L t = 7000 m, L = 200 0 m, tan α = 0.02 , and T = 268.5 s ( T / T r e f l e c t = 0.2953 ).
Reflective Boundary ConditionPartial Relaxation Zone ( L relax = 0.1 λ 0 )Partial Relaxation Zone ( L relax = 0.2 λ 0 )Left ITBC
m = 1
m = 2 0.92410.7246
m = 3 0.94060.7485
m = 4 0.96340.7606
reflection coefficients0.94270.7446no reflectionno reflection

Share and Cite

MDPI and ACS Style

Xu, S.; Dias, F. Long Wave Run-Up Resonance in a Multi-Reflection System. Appl. Sci. 2020, 10, 6172. https://doi.org/10.3390/app10186172

AMA Style

Xu S, Dias F. Long Wave Run-Up Resonance in a Multi-Reflection System. Applied Sciences. 2020; 10(18):6172. https://doi.org/10.3390/app10186172

Chicago/Turabian Style

Xu, Shanshan, and Frédéric Dias. 2020. "Long Wave Run-Up Resonance in a Multi-Reflection System" Applied Sciences 10, no. 18: 6172. https://doi.org/10.3390/app10186172

APA Style

Xu, S., & Dias, F. (2020). Long Wave Run-Up Resonance in a Multi-Reflection System. Applied Sciences, 10(18), 6172. https://doi.org/10.3390/app10186172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop