Next Article in Journal
Cross Diffusion Effect on Linear and Nonlinear Double Diffusive Convection in a Viscoelastic Fluid Saturated Porous Layer with Internal Heat Source
Next Article in Special Issue
Generation of Gravity Waves by Pedal-Wavemakers
Previous Article in Journal
Near Stall Unsteady Flow Responses to Morphing Flap Deflections
Previous Article in Special Issue
A Computational Fluid Dynamics Model for the Small-Scale Dynamics of Wave, Ice Floe and Interstitial Grease Ice Interaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Resonant Periods of Seiches in Semi-Closed Basins with Complex Bottom Topography

Faculty of Mathematics and Natural Sciences, Bandung Institute of Technology, Bandung 40132, Indonesia
*
Author to whom correspondence should be addressed.
Submission received: 29 March 2021 / Revised: 25 April 2021 / Accepted: 29 April 2021 / Published: 9 May 2021
(This article belongs to the Special Issue Mathematical and Numerical Modeling of Water Waves)

Abstract

:
Seiches and resonances are two closely related phenomena that can cause damage to coastal areas. Seiches that occur in a basin at a distinct period named the resonant period may generate resonance when a wave induced by external forces enters the basin and has the same period as the seiches. Studying this period has become essential if we want to understand the resonance better. Thus, in this paper, we derive the resonant period in various shapes of semi-closed basin using the shallow water equations. The equations are then solved analytically using the separation of variables method and numerically using the finite volume method on staggered grid to discover the resonant period for each basin. To validate the numerical scheme, we compare its results against the analytical resonant periods, resulting in a very small error for each basin, suggesting that the numerical model is quite reliable in the estimation of the analytical resonant period. Further, resonant wave profiles are also observed. It is revealed that, in the coupled rectangular basin, the maximum wave elevation is disproportionate to the ratio of the length of the basin, while, in the trapezoidal basin, the ratio of the depth of the basin has no significant impact on the maximum wave elevation.

1. Introduction

Seiches are standing waves that arise in a closed basin (e.g., lake, pool, and reservoir) or semi-closed basin (e.g., bay, port, and harbor) with a unique period. Many theories surrounding seiches originate from Forel’s remarkable publications, which began in 1892 [1]. Seiches may be initiated by a number of exciting forces, such as seismic and atmospheric pressure disturbances, wind, or internal and surface gravity waves [2,3]. These oscillations occur at a particular period commonly called the natural (resonant) period, which strongly depends on the length and depth of the basin in which the seiches occur. The major issue is considered when an external wave enters the basin and has a period that matches the natural period, and therefore resonance happens. Such an occurrence is potentially destructive, as it can destroy moored ships and infrastructures in a port, disrupt coastal activities, and harm the structure located near lakes, reservoirs, or coastlines. The worst consequence that we might encounter is a tsunami wave accompanied by a resonance that ultimately leads to a hazard in the coastal region. Several researchers have documented tsunami cases that are associated with resonance, such as Munger and Cheung [4], Roeber et al. [5], and Yamazaki and Cheung [6], who confirmed the 2006 Kuril, 2009 Samoa, and 2010 Chilean tsunamis over the insular shelves and slopes of Hawaii, American Samoa, and Chile that triggered long oscillations in coastal waters.
In consideration of the negative impacts of tsunami waves on coastlines, we need to find an appropriate way to avoid resonance phenomena. One of the strategies to prevent resonance is to understand the characteristics of seiches over a certain type of basin, particularly its resonant period. Numerous studies have analytically examined the resonant period in the simpler geometric form of basin, which is a rectangular type [7,8,9,10,11]. In addition, many other researchers have discovered the resonant period of the triangular (isosceles) and parabolic shaped basin [9,11]. Magdalena, et. al [12] have extended these studies by involving bottom friction in the derivation of wave resonant period. More complicated bottom topographies, such as semi-closed quartic, triangular (right-angled), trapezoidal, and coupled rectangular type, are rarely studied analytically due to the complexities of the basin depth formula. One of the rare studies that have investigated the quartic basin using an analytical approach was published by Chrystal [1] and Magdalena et al. [13]. However, they only provided the resonant period for closed quartic basin. Meanwhile, for triangular (right-angled), trapezoidal, and coupled rectangular form, e.g., Wilson [14], Rabinovich [15], and Dean and Dalrymple [16] derived analytical solutions from potential theory to approximate the resonant period of seiches in closed basins of those types, while the subject analysis for semi-closed basins has not yet been published, experimentally, analytically, or computationally. Thus, the analytical and numerical solution presented in this paper are new approaches in calculating the resonant period for semi-closed quartic, triangular (right-angled), trapezoidal, and coupled rectangular basin. Moreover, these shapes are more realistic in approximating the structure of a waterbody (e.g., continental shelf, inlet, and gulf) than the simpler geometries mentioned above.
Therefore, in this paper, we propose a mathematical model to simulate seiches, along with corresponding resonance phenomena, in various types of basin, in particular semi-closed quartic, triangular (right-angled), trapezoidal, and coupled rectangular type, in which configurations are modification from Wilson [14]. To focus the study, we only examine the phenomena in semi-closed basins with constant width and no friction force is involved. The main objectives of this paper are to determine the fundamental resonant period of the seiches in each type of basin and analyze the wave profile as a result of the resonance phenomenon. Linear shallow water equations (LSWE) are used, which are to be solved analytically and numerically in order to obtain the resonant period. Separation of the variables method is applied to derive the periods analytically, while the finite volume method on the staggered grid is implemented to formulate a scheme to simulate the model and to acquire the periods numerically. The LSWE model have already been validated previously against the available experimental resonant periods for rectangular, triangular (isosceles), and parabolic cases [9,11]. However, in this work, because we use different basin configurations from the mentioned studies, the analytical approaches in those works are not applicable for our cases. Hence, we need to develop a new analytical procedure in finding the resonant period.
The numerical scheme used in this paper has been tested and validated [17,18,19,20] to emulate various fluid phenomena. In this paper, we apply the same numerical method to simulate the wave resonance phenomenon and obtain the natural period. We also confirm our numerical model, particularly on approximating the resonant period of quartic, triangular (right-angled), trapezoidal, and coupled rectangular semi-closed basin, by comparing the simulations results against the analytical solutions. The accuracy of the numerical method is tested along with the rate of convergence, which has not been provided yet by other works. Further, the profile of the resonant wave is also evaluated to understand the behavior of the wave produced by the resonance phenomenon, depending on the type of basin. From this paper, we expect that the developed model and computational scheme can be conveniently used to extend the topic of the resonant period in the future, so that we can examine seiches, resonance, and their periods in other types of basin that have not been discovered before without the need for a costly experimental model.
The paper is organized into six main sections. Section 1 provides a concise explanation of seiches, resonances, and why those events are dangerous and need to be investigated. Then, the LSWE model used to solve the problem stated in the previous section is introduced in Section 2. The model is then solved analytically in Section 3 and numerically in Section 4. Further, Section 5 presents the simulation results and discussion on the accuracy of our model. Finally, the conclusions of the findings of this study are addressed in Section 6.

2. Mathematical Model

In this section, we introduce the governing equations that are used to model the propagation of seiches in a semi-closed basin with arbitrary bottom topography. It is presumed that the seiches oscillating within the basin follow the shallow water condition where the wavelength is much greater than the maximum depth of the basin (depth < λ / 20 ). Another assumption considered in this model is that the surface elevation, denoted by η ( x , t ) (see Figure 1), is relatively small compared to the total water depth h = η + d , where d ( x ) is the depth of the basin measured at z = 0 . Therefore, we can claim that h d . Thus, it is appropriate to apply the well-known linear shallow water equations (LSWE) as the governing equations of this model. The LSWE are illustrated in Figure 1 and consist of the two equations below.
η t + ( h u ) x = 0 ,
u t + g η x = 0 ,
As pictured in Figure 1, we assume that the seiches in the semi-closed basin are initiated by waves entering the open boundary at the right end of the basin. The waves travel to the left side of the basin at a horizontal velocity of u ( x , t ) , while the vertical velocity is ignored in this situation. Moreover, the g = 9.81 m/s 2 specified in Equation (2) represents the acceleration due to the force of gravity. In the next section, the model described in this part is used to investigate the period of seiches (resonant period) as well as resonance phenomena in several types of semi-closed basin, including quartic, triangular (right-angled), trapezoidal, and coupled rectangular basin.

3. Analytical Solution

Here, the LSWE model is implemented to determine the resonant period for quartic, triangular (right-angled), trapezoidal, and coupled rectangular semi-closed basin. Firstly, the model is solved using the separation of variables method to obtain the solution of η ( x , t ) and u ( x , t ) . We define the function η ( x , t ) and u ( x , t ) as monochromatic waves with constant angular frequency ω , written as
η ( x , t ) = F ( x ) e i ω t ,
u ( x , t ) = G ( x ) e i ω t ,
in which F ( x ) and G ( x ) are functions to be determined. Substituting Equations (3) and (4) into Equations (1) and (2) yields
ω 2 F + h g F x x + g F x h x = 0 ,
G ( x ) = i g ω F x .
Clearly, the function F ( x ) is the solution of the differential Equation (5), while G ( x ) can be obtained from Equation (6). Then, we can derive the fundamental resonant period of the given basin from the obtained solutions F ( x ) and G ( x ) as well as the boundary conditions of each basin. The fundamental resonant period itself is the seiche period in fundamental mode, which produces the biggest resonant period of the basin.

3.1. Quartic Semi-Closed Basin

The first fundamental resonant period to be discovered is the one that belongs to the quartic semi-closed basin whose depth follows the formula h ( x ) = h 0 1 x 2 L 2 2 in the domain [ L , 0 ] . Notation h 0 is the maximum depth and L is the length of the basin. The illustration of this type of basin is shown in Figure 2. Now, to obtain the differential equation containing F ( x ) for this specific topography, we substitute the depth into Equation (5), and then we have
ω 2 F + h 0 1 x 2 L 2 2 g F x x g F x 4 x h 0 L 2 1 x 2 L 2 2 = 0
with hard-wall boundary condition at x = L which can be written mathematically as u ( L , t ) = 0 . According to Equation (6), this condition can be rewritten as
F x ( L ) = 0 .
In addition, there is another boundary condition that is adapted from the theory of the nodal lines in semi-closed basin for fundamental mode [15], which implies that, at the open end of the basin, the wave has minimum amplitude, or, in other words, η ( 0 , t ) = 0 , which can also be written as
F ( 0 ) = 0 .
Now, solving Equation (7) gives us the solution read as
F ( x ) = K 1 i x g h 0 + γ ( L 2 x 2 ) L x L + x i γ 4 L 2 g h 0 + K 2 i x g h 0 + γ ( L 2 x 2 ) L + x L x i γ 4 L 2 g h 0 ,
in which γ = L L 2 ω 2 g h 0 , while K 1 and K 2 are certain constants. The derivative of F with respect to x is
F x ( x ) = 1 / ( L 4 ω 2 ) ( γ 2 + g h 0 x 2 ) ( L 2 x 2 ) 3 K 1 i x g h 0 + γ L x L + x i γ 4 L 2 g h 0 ( φ ( x ) + L 8 ω 4 x ) + K 2 i x g h 0 + γ L + x L x i γ 4 L 2 g h 0 ( φ ( x ) L 8 ω 4 x ) ,
where φ ( x ) = i γ L 5 ω 2 L 2 g h 0 γ 2 2 g h 0 x 2 + L 2 g h 0 . Substituting Equation (11) into Equation (8) results in K 1 = 0 . Applying this result to Equation (10), and then substituting the new equation into Equation (9), yields:
L 2 ω 2 = g h 0 ,
where ω = 2 π T 1 . Hence, we obtain the fundamental resonant period of quartic semi-closed basin, which is
T 1 = 3.1416 2 L g h 0 .

3.2. Triangular (Right-Angled) Semi-Closed Basin

The second type of basin discussed is the triangular (right-angled) semi-closed basin with depth of h ( x ) = h 1 x L for domain [ 0 , L ] in which h 1 is the maximum depth. By substituting h ( x ) into Equation (5), we have a second-order ordinary differential equation written as
ω 2 F + h 1 x L g F x x + g F x h 1 L = 0 ,
with boundary condition defined according to the illustration in Figure 3, where we have hard-wall boundary condition at x = 0 . Similar to the previous case, here, we also consider the condition that is following the nodal lines theory, where the wave amplitude is minimum at the open end of the basin or at x = L . Therefore, we have the following conditions:
F x ( 0 ) = 0 ,
F ( L ) = 0 .
Now, solving the differential Equation (14), we obtain
F ( x ) = K 3 J 0 2 ω x L g h 1 + K 4 Y 0 2 ω x L g h 1 ,
where K 3 and K 4 denote certain constants, while J n ( x ) and Y n ( x ) are Bessel functions of the first and second kind of order n, respectively. These functions are two linearly independent solutions of the Bessel differential equation. Moreover, the derivative of F with respect to x is written as
F x ( x ) = K 3 J 1 2 ω x L g h 1 ω L g h 1 x K 4 Y 1 2 ω x L g h 1 ω L g h 1 x .
Next, to determine the fundamental resonant period, we apply Equation (15) to Equation (18) and obtain
K 3 J 1 2 ω x L g h 1 x = 0 + K 4 Y 1 2 ω x L g h 1 x = 0 = 0 .
Considering the behavior of Bessel functions of order one near x = 0 , where J 1 ( 0 ) = 0 and lim α 0 Y 1 ( α ) = , then K 4 = 0 must be satisfied so that Equation (19) is also satisfied. Now, using the second condition, we have
F ( L ) = K 3 J 0 2 ω L g h 1 = 0 .
For the nontrivial solution to be achieved, K 3 must not equals to zero. Thus, to satisfy Equation (20), we need to find the root of the first-kind Bessel function of order zero. In addition, since we aim to obtain the fundamental resonant period, we are looking for the first root of the function, which, from the Wolfram Language website, is found to be 2.4048 . Therefore, 2 ω L g h 1 = 2.4048 must be satisfied. Hence, the sought period is
T 1 = 2.6128 2 L g h 1 .

3.3. Trapezoidal Semi-Closed Basin

In the third case, we determine the fundamental resonant period for trapezoidal semi-closed basin with depth h ( x ) = h 0 + ( h 1 h 0 ) x L in domain [ 0 , L ] , as depicted in Figure 4. Here, the maximum depth is denoted by notation h 1 while h 0 is the minimum depth. Three configurations of h 0 and h 1 are used in this paper to evaluate the effect of the difference in the ratio between the minimum and maximum depth on the value of the resonant period and wave profile: h 0 h 1 = 1 2 , h 0 h 1 = 1 3 , and h 0 h 1 = 1 4 . The detailed derivation is only provided for the first configuration, i.e., h 0 h 1 = 1 2 .
The substitution of h and h 1 = 2 h 0 into Equation (5) gives us
ω 2 F + g h 0 L F x + g h 0 + h 0 x L F x x = 0 ,
with boundary conditions
F x ( 0 ) = 0 ,
F ( L ) = 0 .
Solving Equation (22) gives us
F ( x ) = K 5 J 0 2 ω L ( L + x ) g h 0 + K 6 Y 0 2 ω L ( L + x ) g h 0 ,
in which K 5 and K 6 are constants. The first derivative of F is
F x ( x ) = ω L g h 0 ( L + x ) K 5 J 1 2 ω L ( L + x ) g h 0 + K 6 Y 1 2 ω L ( L + x ) g h 0 .
By substituting the boundary conditions (23) and (24) into Equations (26) and (25), respectively, we obtain two equations read as
K 5 J 1 2 ω L g h 0 + K 6 Y 1 2 ω L g h 0 = 0 ,
K 5 J 0 2 2 ω L g h 0 + K 6 Y 0 2 2 ω L g h 0 = 0 .
Eliminating K 5 from the equations and ignoring the trivial solution, we have
J 0 2 2 ω L g h 0 Y 1 2 ω L g h 0 J 1 2 ω L g h 0 Y 0 2 2 ω L g h 0 = 0 .
The above result is also obtained if we eliminate K 6 instead of K 5 . Now, Equation (29) is only satisfied if
ω L g h 0 = 2.0308 .
Because ω = 2 π / T 1 , the fundamental resonant period for trapezoidal semi-closed basin with h 0 h 1 = 1 2 is
T 1 a = 1.5470 2 L g h 0 .
Furthermore, using the exact same steps and method, we determine the fundamental resonant period for trapezoidal basin with the other two ratios, h 0 h 1 = 1 3 and h 0 h 1 = 1 4 , which are
T 1 b = 1.3157 2 L g h 0 .
and
T 1 c = 1.1674 2 L g h 0 ,
respectively.

3.4. Coupled Rectangular Semi-Closed Basin

The last type of basin is the coupled rectangular semi-closed basin, the topography of which is shown in Figure 5 and defined in a piece-wise function
h ( x ) = h 1 , L 1 x 0 , h 2 , 0 < x L 2 .
h 1 is the water depth in the shallower area with length of L 1 , while h 2 denotes the water depth in the deeper region with length of L 2 . Unlike the previous case, here, the depth ratio between h 1 and h 2 is fixed to h 1 h 2 = 1 4 [14]. However, the ratio between the length of the shallower and deeper domain, L 1 and L 2 , respectively, is varied to observe its effect on the resonant period and wave profile. There are four configurations of the lengths: L 1 L 2 = 1 2 , L 1 L 2 = 1 3 , L 1 L 2 = 1 4 , and L 1 L 2 = 1 8 . Only the case with length ratio of L 1 L 2 = 1 2 is explained in detail below.
To start the derivation, we assume that the function F ( x ) from Equation (3) can be rewritten as
F ( x ) = F 1 ( x ) , L 1 x 0 , F 2 ( x ) , 0 < x L 2 .
Then, we substitute h ( x ) with the mentioned depth ratio and length criteria to Equation (5) with F ( x ) following the piece-wise function (35), which gives us
0 = ω 2 F 1 + g h 2 4 F 1 x x , L 1 x 0 , ω 2 F 2 + g h 2 F 2 x x , 0 < x L 2 .
The solution of Equation (36) is obtained using the characteristic method, read as
F ( x ) = F 1 ( x ) = A cos 2 ω g h 2 x + B sin 2 ω g h 2 x , L 1 x 0 , F 2 ( x ) = C cos ω g h 2 x + D sin ω g h 2 x , 0 < x L 2 .
where A, B, C, and D are certain constants.
Next, we need to make sure that the solution is continue at the discontinue point, which is x = 0 , by applying the matching conditions. According to Dean and Dalrymple [16], there are two matching conditions based on the continuity of surface elevation and mass flow rate at x = 0 , written as
F 1 ( 0 ) = F 2 ( 0 ) ,
h 1 F 1 x ( 0 ) = h 2 F 2 x ( 0 ) ,
respectively. From the first condition (Equation (38)), we find that A = C , while Equation (39) gives us D = B / 2 . Hence, we have a new solution for both F 1 and F 2 read as
F 1 ( x ) = A cos 2 ω g h 2 x + B sin 2 ω g h 2 x ,
F 2 ( x ) = A cos ω g h 2 x + B 2 sin ω g h 2 x .
By taking the derivative of F 1 and F 2 , we obtain
F 1 x ( x ) = A sin 2 ω g h 2 x 2 ω g h 2 + B cos 2 ω g h 2 x 2 ω g h 2 ,
F 2 x ( x ) = A sin ω g h 2 x ω g h 2 + B 2 cos ω g h 2 x ω g h 2 .
Similar to the previous cases, for this type of basin, we have two boundary conditions:
F 1 x ( L 1 ) = F 1 x ( L 2 / 2 ) = 0 ,
F 2 ( L 2 ) = 0 .
Substituting the boundary conditions to Equations (42) and (41), we obtain
A sin ω L 2 g h 2 + B cos ω L 2 g h 2 = 0 ,
A cos ω L 2 g h 2 + B 2 sin ω L 2 g h 2 = 0 .
By eliminating constant A and ignoring the trivial solution, we have
cos 2 ω L 2 g h 2 1 2 sin 2 ω L 2 g h 2 = 0 ,
which is also obtained if we eliminate the other constant. Solving Equation (48) yields
T 1 d = 1.4370 2 L 2 g h 2 ,
which is the fundamental resonant period we are looking for. Moreover, a similar procedure is applied to determine the fundamental resonant period for the coupled rectangular basin with length ratio of L 1 L 2 = 1 3 , L 1 L 2 = 1 4 , and L 1 L 2 = 1 8 . Those resonant periods, respectively, are written as
T 1 e = 1.1626 2 L 2 g h 2 ,
T 1 f = 1.0000 2 L 2 g h 2 ,
and
T 1 g = 0.7776 2 L 2 g h 2 .
All the results found in this section are listed in Table 1, along with the fundamental resonant periods obtained using numerical approach that are explained in the next section.

4. Numerical Scheme

In this section, we aim to observe the propagation of waves in a semi-closed basin when a resonance phenomenon arises using a numerical approach. Hence, we construct a numerical scheme to simulate the occurrence accurately. In addition, this computational scheme is also used to investigate the numerical fundamental resonant period for each type of basin, which is tested and verified against analytical solutions from the previous section.
Firstly, we revisit our mathematical model, i.e., the LSWE model, represented by Equations (1) and (2). Consider [ 0 , L ] as the spatial domain of a semi-closed basin where the waves propagate. Then, the domain is partitioned in a staggered way into x 1 2 = 0 , x 1 , , x j 1 2 , x j , x j + 1 2 , , x N x + 1 2 = L . The numerical scheme is formulated based on the staggered grid using the finite volume method, as illustrated in Figure 6. Using this method, the cells centered at full grid points ( x j , j = 1 , 2 , , N x ) are used to approximate the mass conservation Equation (1), while their adjacent cells are used to approximate the momentum balance Equation (2). In addition, the value of the horizontal velocity u is stored only in half grid points ( x j + 1 / 2 , j = 1 , 2 , , N x ), whereas η and h are calculated only at the full grid points.
According to the explanation above, Equations (1) and (2) are, respectively, approximated by
η j n + 1 = η j n Δ t Δ x ( ( h * u ) j + 1 2 n ( h * u ) j 1 2 n ) ,
u j + 1 2 n + 1 = u j + 1 2 n g Δ t Δ x ( η j + 1 n + 1 η j n + 1 ) ,
where
h j + 1 2 * = h j + 1 + h j 2 .
Using the von Neumann stability analysis, it is found that the Courant–Friedrichs–Lewy stability condition of our numerical scheme is g h max Δ t Δ x 1 . For the following section, we simulate the wave propagation and resonance phenomenon in each type of basin using the scheme (Equations (53) and (54)) in order to obtain the numerical fundamental resonant period. Concerning these simulations, the computational scheme is validated by comparing the numerical results with the analytical solutions obtained previously.

5. Results and Discussion

Now, the numerical scheme established in the previous section is applied to simulate the resonance phenomenon in each type of semi-closed basin discussed. The results are presented to determine the fundamental resonant period numerically as well as to observe the resonant wave profile. The numerical resonant period for each basin is then compared to its analytical counterpart in order to evaluate the reliability of our numerical scheme on simulating resonance phenomena.
For the computation, we set the total length of the basin to be L = 25 m and the maximum depth to be H = 1.5 m in order to fulfill the shallow water requirements. The spatial domain is then divided into cells with length of Δ x = 0.1 m. Moreover, the observation time is set to be within the range of [ 150 , 500 ] s, so that we can capture the occurrence of resonance clearer. As for the initial condition, we assume that the water is calm and undisturbed, which means that η ( x , 0 ) = 0 and u ( x , 0 ) = 0 . Then, we presume that an incoming wave with an amplitude of A i = 0.1 m enters the basin and disturbs the calm water, generating seiches inside the basin. In addition, we use the same boundary conditions explained and illustrated in Section 3.
To calculate the numerical resonant period, we use the trial and error procedure, where we input a certain value of wave period into the model and simulate the wave propagation, then observe whether or not the wave undergoes resonance. These steps are carried out numerous times until the resonance-generating period is found. The resonance itself is indicated by a gradual rise of wave amplitude throughout the observation. Thus, we obtain a numerical fundamental resonant period for each basin with the specific parameters specified above. In addition, we compare the results against the analytical solutions, as listed in Table 1.
In Table 1, the numerical fundamental resonant period for each basin is compared to the corresponding analytical solution to examine the accuracy of our numerical scheme. The error presented in the table is calculated using the formula e r r o r = N u m e r i c a l T 1 A n a l y t i c a l T 1 A n a l y t i c a l T 1 × 100 % . Notice that relative errors in nearly all cases are less than 0.65 % , with the exception of a quartic basin case which resulted in an error of around 2.07 % . The larger error shown in quartic basin is most likely caused by the dry area at the end of the domain. In simulations, this dry area is illustrated by the small part at the end of the domain that is significantly close to z = 0 and is considered as z = 0 by the numerical system. The existence of this dry area eventually leads to the emergence of wave non-linearity effect that is not captured well by our linear numerical scheme. Therefore, the error in this particular basin is larger than the others.
On the other hand, our analytical model seems to estimate the real fundamental resonant period quite well. This is supported by the observation done by Darbyshire [21] in Slangkop Point, South Africa. In this observation, the fundamental resonant period was measured for a quartic-like coastal area with length of L = 37.04 km and depth of H 0 = 274.32 m. It was found that the observed resonant period is T 1 = 73 s. Using our analytical model and numerical scheme, we obtain that the respective T 1 are 74.77 and 76.32 s. From those findings, we can calculate the relative errors between analytical and observed T 1 as well as between numerical and observed T 1 , which are 2.425 % and 4.548 % , respectively. Both errors are considerably small for model comparison against real measurement, since there are many natural factors that are ignored in the linear model, such as friction, non-linearity, and wind. Therefore, it is safe to say that both models successfully estimate the real T 1 measured. However, the numerical error is still almost twice as big as the analytical error, which confirms that, in this specific type of basin, our numerical scheme needs to be improved by involving the wave non-linearity effect on the dry area.
Additionally, the accuracy of our computational scheme is determined by calculating the rate of convergence ( τ ). Since our model involves two independent variables, x (spatial) and t (time), we evaluate its accuracy from both perspectives. We perform numerous simulations in trapezoidal basin ( H 0 / H 1 = 1 / 2 ) with the same parameters but different Δ x and Δ t values in each simulation. In terms of time, we use a fixed Δ x = 0.05 m, while Δ t varies between 0.003093 and 0.013034 s. The error between analytical and numerical T 1 is then calculated for each time step using the formula E T j = | A n a l y t i c a l T 1 N u m e r i c a l T 1 | . Each step’s rate of convergence is computed using τ t j = l o g ( E T ( j + 1 ) / E T j ) / l o g ( Δ t j + 1 / Δ t j ) . The results are listed in Table 2.
In Table 2, we can see that all values of | τ t j | are within the range of 0 | τ t j | 1 ; thus, the order of accuracy of our numerical scheme in term of time is 1. Then, using the same calculations, we simulate the numerical scheme with fixed Δ t = 0.013034 s and various Δ x of range 0.05 Δ x 0.50625 m. Table 3 presents the calculation results for each spatial step. In this case, the values of | τ x j | are in the range of 0 to 2, indicating that the order of accuracy in term of spatial is 2. As a result, the order of accuracy of our numerical scheme is O ( Δ t , Δ x 2 ) .
Moreover, the same scheme is used to observe the wave profile when resonance phenomenon occurs and how different the behavior of the resonant wave depends on the type of the basin where it happens. Using the same parameters, we simulate wave resonance in each basin with the incoming wave period follows the obtained fundamental resonant period. The results of the simulations for quartic, triangular (right-angled), trapezoidal, and coupled rectangular semi-closed basin are shown in Figure 7, Figure 8 and Figure 9, respectively.
As shown in Figure 7, Figure 8 and Figure 9, the wave elevation within each basin increases over time, which implies that a resonance occurs when an external wave with a period of T 1 enters the basin. Note that resonance occurrences in all cases produce waves with amplitudes that continue to increase and eventually exceed 30 times the initial amplitude (3 m) in the first 200 s. In real life, this could mean that, in just a few minutes after the wave entered the basin, the amplitude would rise dramatically and, as a result, there would be a risk of major destruction of the facility or environment in the coastal region. The worse case would happen if the initial wave has bigger amplitude. In addition, from the results in all types of basin, we are able to conclude that the worse resonance happens in the quartic basin, as the wave amplitude in this basin rises the most, almost seven times the initial amplitude, compared to the other three at only 3–4 times. This outcome might happen because of the same factor explained above, which is the existence of the dry area that eventually leads to the arise of wave non-linearity effect. Because our linear scheme cannot handle the non-linearity very well, the wave amplitude rises more quickly than in other basins. This means, for coastal areas or semi-closed basins in general with quartic type of bottom topography, we need to be more attentive on dealing with resonances and might be need major modification in the numerical scheme itself.
Furthermore, we observe the differences of the resonant wave profile in trapezoidal and coupled rectangular semi-closed basin with different configurations. Figure 10 shows that, considering the variations in configuration of the minimum and maximum depth of the topography in the trapezoidal basin, there are no significant differences in the amplitude increases. This means that the change in the depth ratio will not significantly affect the resonance phenomenon in the trapezoidal basin. On the other hand, the different configurations of the length ratio of the coupled rectangular basin have a significant impact on the resonant wave amplitude. The lowest wave amplitude is produced by the greatest ratio, that is L 1 L 2 = 1 2 , or when the length of the deeper region is twice as long as the shallower domain. With the same manner, the smallest ratio L 1 L 2 = 1 8 , which is when the basin has much longer deep domain than the shallow one, gives the highest amplitude. In real life, this discovery might be helpful in deciding the optimum design of a human-made coupled rectangular harbor’s topography, that is keeping the length ratio to be as large as possible. This finding may also be used as a basic information on the construction of a protection structure in this specific type of harbor or other coastal areas.

6. Conclusions

The seiches and resonance phenomena over quartic, triangular (right-angled), trapezoidal, and coupled rectangular semi-closed basin were successfully simulated using the linear shallow water equations, which are solved analytically and numerically to obtain the fundamental resonant period of each basin. The separation of variables method is used to derive the resonant period from the equations analytically, whereas the numerical resonant periods are determined by applying the finite volume method on a staggered grid to the equations. Comparisons were conducted between the analytical and numerical fundamental resonant period in each basin to verify the formulated numerical scheme, resulting in errors of less than 2.1 % , which means that our scheme is highly accurate on estimating the resonant period. Moreover, the numerical simulations successfully simulated the resonance phenomena in each basin as a result of an external wave, whose period is the same as the derived fundamental resonant period, enters the basin mouth. In the simulation results, this phenomenon is indicated by the gradual increase of the amplitude of the wave. Further, from the numerical simulations, we can draw several noticeable outcomes. In general, since the amplitude of the resonant wave in quartic basin increased the most in a certain period of time, we need to be more attentive on handling the resonance phenomenon in this particular basin. Another perceptible finding is the negative correlation between the ratio of the basin’s length and the wave elevation for coupled rectangular basin. The smaller is the ratio of the basin’s length, the higher is the wave amplitude it produces, and vice versa. On the other hand, we found no significant effect of the various depth ratios of the trapezoidal basin on the increasing of wave amplitude. Furthermore, since this paper only focuses on basins with constant width and no friction force is involved, for further research, we suggest to do some modifications on the basin width and consider including the friction term into the model. The effect of the varying width and friction force to resonant wave amplitude might be an interesting case to study.

Author Contributions

Conceptualization, I.M.; model and methodology, I.M.; analytical solution, N.K.; numerical scheme, I.M.; simulation program, I.M.; numerical simulation, N.K.; investigation, I.M., H.Q.R.; writing—original draft preparation, N.K.; writing—review and editing, I.M. and H.Q.R.; supervision, I.M.; project administration, I.M.; and funding acquisition, I.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the ITB and The Indonesian Ministry of Research, Technology and Higher Education Research Grant.

Acknowledgments

This work was supported by the ITB and The Indonesian Ministry of Research, Technology and Higher Education Research Grant.

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:
LSWELinear Shallow Water Equations

References

  1. Chrystal. Hydrodynamical theory of seiches. Trans. R. Soc. Edinb. 1906, 41, 620–630. [Google Scholar]
  2. Chapman, D.; Giese, G. Seiches. In Encyclopedia of Ocean Sciences, 2nd ed.; Elsevier: New York, NY, USA, 2001; Volume 5, pp. 344–350. [Google Scholar]
  3. Rozas, C.; de la Fuente, A.; Ulloa, H.; Davies, P.; Nino, Y. Quantifying the effect of wind on internal wave resonance in Lake Villarrica, Chile. Environ. Fluid Mech. 2014, 14, 849–871. [Google Scholar] [CrossRef]
  4. Munger, S.; Cheung, K.F. Resonance in Hawaii waters from the 2006 Kuril Islands Tsunami. Geophys. Res. Lett. 2008, 35, L07605. [Google Scholar] [CrossRef]
  5. Roeber, V.; Yamazaki, Y.; Cheung, K.F. Resonance and impact of the 2009 Samoa tsunami around Tutulia, American Samoa. Geophys. Res. Lett. 2010, 37, L21604. [Google Scholar]
  6. Yamazaki, Y.; Cheung, K.F. Shelf resonance and impact of near-field tsunami generated by the 2010 Chile earthquake. Geophys. Res. Lett. 2011, 38, L12605. [Google Scholar] [CrossRef]
  7. Hwang, L.S.; Tuck, E.O. On the oscillations of harbours of arbitrary shape. J. Fluid Mech. 1970, 42, 447–464. [Google Scholar] [CrossRef]
  8. Kabiri-Samani, A.; Ataie-Ashtiani, B. Free Water Surface Oscillations in a Closed Rectangular Basin with Internal Barriers. Sci. Iran. 2008, 15, 315–322. [Google Scholar]
  9. Magdalena, I.; Rif’atin, H.Q.; Reeve, D.E. Seiches and Harbour Oscillations in A Porous Semi-Closed Basin. Appl. Math. Comput. 2020, 369, 124835. [Google Scholar] [CrossRef]
  10. Martinez, F.M.; Naverac, V.S. An experimental study of harbor resonance phenomena. In Proceedings of the Twenty-First International Conference on Coastal Engineering, Costa del Sol-Malaga, Spain, 20–25 June 1988; pp. 1280–1292. [Google Scholar]
  11. Rif’atin, H.Q.; Magdalena, I. Seiches in A Closed Basin of Various Geometric Shapes. J. Phys. Conf. Ser. 2019, 1245, 012061. [Google Scholar] [CrossRef]
  12. Magdalena, I.; Rif’atin, H.Q.; Matin, A.M.A. Analytical and numerical studies for Seiches in a closed basin with bottom friction. Theor. Appl. Mech. Lett. 2020, 10, 429–437. [Google Scholar] [CrossRef]
  13. Magdalena, I.; Karima, N.; Rif’atin, H.Q. A Mathematical Model for Investigating The Resonance Phenomenon in Lakes. Wave Motion 2021, 100, 102669. [Google Scholar] [CrossRef]
  14. Wilson, B. Chapter: Seiches in Advances in Hydroscience; Academic Press: Cambridge, MA, USA, 1972; Volume 8, pp. 1–93. [Google Scholar]
  15. Rabinovich, A. Chapter 9: Seiches and Harbor Oscillations in Handbook of Coastal and Ocean Engineering; Kim, Y.C., Ed.; World Scientific Publ.: Singapore, 2009; pp. 193–236. [Google Scholar]
  16. Dean, R.; Dalrymple, R. Water Wave Mechanics for Engineers and Scientist; World Scientific: Singapore, 1991. [Google Scholar]
  17. Magdalena, I.; Adityawan, M.; Jonathan, C. Numerical model for dam break over a movable bed using finite volume method. Int. J. Geomate 2020, 19, 98–105. [Google Scholar] [CrossRef]
  18. Magdalena, I.; Erwina, N.; Pudjaprasetya, S. Staggered Momentum Conservative Scheme for Radial Dam break Simulation. J. Sci. Comput. 2015, 65, 867–874. [Google Scholar] [CrossRef]
  19. Mungkasi, S.; Magdalena, I.; Pudjaprasetya, S.R.; Wiryanto, L.H.; Roberts, S.G. A Staggered Method for The Shallow Water Equations Involving Varying Channel Width and Topography. Int. J. Multiscale Comput. Eng. 2018, 16, 231–244. [Google Scholar] [CrossRef]
  20. Pudjaprasetya, S.; Magdalena, I. Momentum conservative schemes for shallow water flows. E. Asian J. Appl. Math. 2014, 4, 152–165. [Google Scholar] [CrossRef]
  21. Darbyshire, M. Long Waves on the Coast of the Cape Peninsula. Dtsch. Hydrogr. Z. 1963, 16, 167–185. [Google Scholar] [CrossRef]
Figure 1. Description of the shallow water equations model in a semi-closed basin with arbitrary bottom topography.
Figure 1. Description of the shallow water equations model in a semi-closed basin with arbitrary bottom topography.
Fluids 06 00181 g001
Figure 2. Illustration of the quartic semi-closed basin.
Figure 2. Illustration of the quartic semi-closed basin.
Fluids 06 00181 g002
Figure 3. Triangular (right-angled) semi-closed basin.
Figure 3. Triangular (right-angled) semi-closed basin.
Fluids 06 00181 g003
Figure 4. Trapezoidal semi-closed basin.
Figure 4. Trapezoidal semi-closed basin.
Fluids 06 00181 g004
Figure 5. Coupled rectangular semi-closed basin.
Figure 5. Coupled rectangular semi-closed basin.
Fluids 06 00181 g005
Figure 6. Illustration of the finite volume method on a staggered grid.
Figure 6. Illustration of the finite volume method on a staggered grid.
Fluids 06 00181 g006
Figure 7. Simulation result for quartic semi-closed basin.
Figure 7. Simulation result for quartic semi-closed basin.
Fluids 06 00181 g007
Figure 8. Simulation result for triangular (right-angled) semi-closed basin.
Figure 8. Simulation result for triangular (right-angled) semi-closed basin.
Fluids 06 00181 g008
Figure 9. Simulation results for trapezoidal semi-closed basin.
Figure 9. Simulation results for trapezoidal semi-closed basin.
Fluids 06 00181 g009
Figure 10. Simulation results for coupled rectangular semi-closed basin.
Figure 10. Simulation results for coupled rectangular semi-closed basin.
Fluids 06 00181 g010
Table 1. The fundamental resonant period obtained for each type of semi-closed basin using analytical and numerical approach.
Table 1. The fundamental resonant period obtained for each type of semi-closed basin using analytical and numerical approach.
Basin Type h ( x ) Analytical T 1 Numerical T 1 Error (%)
Quartic h 0 1 x 2 L 2 2 3.1416 2 L g h 0 3.2066 2 L g h 0 2.07
Triangular (right-angled) h 1 x L 2.6128 2 L g h 1 2.6212 2 L g h 1 0.32
Trapezoidal h 0 + ( h 1 h 0 ) x L h 0 / h 1 = 1 / 2 1.5470 2 L g h 0 1.5494 2 L g h 0 0.16
h 0 / h 1 = 1 / 3 1.3157 2 L g h 0 1.3181 2 L g h 0 0.18
h 0 / h 1 = 1 / 4 1.1674 2 L g h 0 1.1694 2 L g h 0 0.17
Coupled rectangular 1 4 h 2 , ( L 1 x 0 ) ; h 2 , ( 0 < x L 2 ) L 1 / L 2 = 1 / 2 1.4370 2 L 2 g h 2 1.4404 2 L 2 g h 2 0.24
L 1 / L 2 = 1 / 3 1.1626 2 L 2 g h 2 1.1652 2 L 2 g h 2 0.22
L 1 / L 2 = 1 / 4 1.0000 2 L 2 g h 2 1.0064 2 L 2 g h 2 0.64
L 1 / L 2 = 1 / 8 0.7776 2 L 2 g h 2 0.7788 2 L 2 g h 2 0.15
Table 2. Rates of convergence of our numerical scheme in terms of time (t).
Table 2. Rates of convergence of our numerical scheme in terms of time (t).
j Δ t j ( s ) E Tj | τ tj | j Δ t j ( s ) E Tj | τ tj |
10.0130340.0011710.11474540.0054990.0012320.611683
20.0097760.0012100.24495150.0041240.0010340.450214
30.0073320.0012990.18216160.0030930.001176
Table 3. Rates of convergence of our numerical scheme in terms of spatial (x).
Table 3. Rates of convergence of our numerical scheme in terms of spatial (x).
j Δ x j ( s ) E Xj | τ xj | j Δ x j ( s ) E Xj | τ xj |
10.050000.0011710.77054740.168750.0025071.152100
20.075000.0008571.93562950.337500.0055710.630723
30.112500.0018780.71199660.506250.004318
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Magdalena, I.; Karima, N.; Rif’atin, H.Q. Resonant Periods of Seiches in Semi-Closed Basins with Complex Bottom Topography. Fluids 2021, 6, 181. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6050181

AMA Style

Magdalena I, Karima N, Rif’atin HQ. Resonant Periods of Seiches in Semi-Closed Basins with Complex Bottom Topography. Fluids. 2021; 6(5):181. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6050181

Chicago/Turabian Style

Magdalena, Ikha, Nadhira Karima, and Hany Qoshirotur Rif’atin. 2021. "Resonant Periods of Seiches in Semi-Closed Basins with Complex Bottom Topography" Fluids 6, no. 5: 181. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6050181

Article Metrics

Back to TopTop