1. Introduction
Tsunamis, oceanic waves with small amplitudes (usually less than 1 m) in deep water and long wavelengths (usually hundreds of kilometers), are generated by the rapid displacement of water mainly due to earthquakes, underwater landslides, and volcanic eruptions. In shallow areas or dry lands, tsunami waves break into a number of solitons and usually form a turbulent bore; see [
1,
2,
3,
4]. A train of breaking waves and bore are formed, impacting coastal structures. Significant damage has been caused by tsunamis’ impact on coastal structures, especially coastal bridges. In particular, the failure of coastal bridges significantly impacts the transportation system, given that the connections to the disaster-affected areas are cut off. This may result in a delayed rescue and recovery of areas affected by the destructive events.
In the failure of a bridge, a body of fluid, i.e., wave or bore, impacts the bridge, resulting in loads larger than the bridge’s capacity. In almost every major oceanic event in the past few decades, coastal bridge failures have been observed. The tsunami caused by the 2004 Sumatra earthquake impacted hundreds of bridges on the west coast of Sumatra, Indonesia; see [
5]. More than 250 coastal bridges were damaged by the 2011 Tōhoku tsunamis; see [
6]. Studies on the bore impact on bridges are traditionally divided into two categories, namely (i) the generation of a tsunami-like bore by using solitary waves, breaking solitary waves, and breaking the dam of a reservoir (namely dam-break), and (ii) bore impact on the bridge. Both of them were studied in this work.
In most tsunami events, tsunamis are caused by earthquakes. Many bridges can survive the earthquake but fail to resist the impact of tsunamis; see [
7]. In these failures, the bridge pipes remain whereas the bridge deck collapses. The impact of fluid on the bridge deck can be subdivided into horizontal and vertical force. If the peak value of vertical force equals or exceeds the vertical resistance of the connection part between the bridge deck and pipes, even a small horizontal force could push the deck off, or failure happens directly due to the large peak value of horizontal force. Hence, strengthening the resistance of the connections is suggested to reduce this kind of failure, where the wave–structure interaction is considered. On the other hand, failure still happens when the tsunami-induced impact is too large, even if the connection was designed to be strong enough according to the industrial standard; see [
7]. To this extent, the generation mechanism of a tsunami-like bore directly affects the profile, propagation speed, and flow field. The peak value of the wave–structure interaction of a tsunami-induced bore is suggested to be studied.
When a tsunami approaches the shore, its wave speed decreases according to
in shallow water, where
c is the wave speed,
g is the gravitational acceleration, and
h is the water depth, and its skewness, asymmetry, and wave height increase; see [
3,
8]. On dry lands and shallow areas, tsunami waves usually form a turbulent bore, and the effects of frequency dispersion balance non-linearity, while some certain huge waves break directly; see [
1,
2,
3,
4]. The un-breaking waves transform into a bore. As the formed bore shoaling, a train of breaking waves is formed, impacting the coastal and adjacent structures. Due to the instantaneity of the occurrence of a tsunami, there is little record about the Grimshaw-derived [
9] equations used to describe the amplitudes of solitary waves with slowly varied floors by energy conservation. Benilov et al. [
10] related the amplitude and velocity of an envelope soliton (packet of surface gravity waves or a series of solitary) to the depth in the topography-modified Schrödinger equation. They reported that the amplitude of a shoaling soliton decreases as it propagates. Grimshaw and Annenkov [
11] used a higher-order nonlinear Schrödinger equation to describe the deformation of a soliton. In their equation, the soliton can penetrate a shallow depth smaller than
. Rajan et al. [
12] used the Djordjević and Redekopp equation—see [
13]—to study a soliton. The maximum amplitude of the soliton decreases or increases as the water depth decreases or increases. Nik et al. [
14] used the Korteweg–de Vries (KdV) equation—see [
15]—to study the internal solitary wave propagating on variable topography, such as a constant-slope bottom and specific bottom profile following the previous studies. A secondary trailing wave packet was found after the wave passes the topographies. Wang et al. [
16] carried out laboratorial experiments to study O-solitons (solitary waves whose crest is shaped like the letter “o”). The phase of this kind of solitary wave remains after the collapse for a short period. Quezada and Ketcheson [
17] derived a KdV-type equation to study the propagations of solitary waves over a periodically variated bottom. The effective dispersion of solitary waves on the bottom was reported in detail. According to these studies, the breaking method of a solitary wave could usually be the bottom variation, such as a gentle slope.
Many experiments have been conducted to study the impact on deck. Bradner et al. [
18] conducted experiments to study the wave forces on causeway-type coastal highway bridges using a wide range of regular and random wave conditions. In order to gain an insight into the tsunami forces that various types of bridges have to withstand, ref. [
19] conducted experiments to test the tsunami bore impact on 1:20-scale bridge specimens with different shapes. A series of laboratory experiments and CFD calculations were conducted by [
20] to study the horizontal and vertical forces of solitary waves on a horizontal plate. Istrati et al. [
21] conducted the largest experiments, on a 1:5 scale, of a tsunami impact on a straight single-span bridge. The total forces on the bridge and the distribution in individual bearing connections were reported. Xiang et al. [
22] employed a solitary wave to represent a tsunami. The horizontal and vertical fluid impact forces on a bridge superstructure, on a 1:5 scale, were studied through large-scale laboratory experiments and numerical simulations. Istrati and Buckle [
23]) conducted the largest-scale experiments of tsunami impact on an open-girded skewed bridge specimen with a 45° angle.
For the analytical simulations of the impact on deck, ref. [
24] investigated the interaction between a train of surface gravity wave and a horizontal plate by assuming that the wavelength and plate length are more significant than the water depth. An analytical solution was provided under that assumption. Hayatdavoodi et al. [
25] studied the horizontal and vertical wave forces of cnoidal waves on a horizontal flat plate in shallow water. Based on Level I Green–Naghdi equations, a nonlinear shallow-water model was provided to study and compare with laboratorial experimental data.
Two main types of mesh are used in the computational simulations based on the Navier–Stokes equations. One type is the mesh-based method. Azadbakht [
26] conducted a 2D numerical simulation of tsunami impact on simplified models of bridges located on the West Coast of the US by using finite element models. Istratiand and Buckle [
27] investigated, via 2D and more realistic 3D CFD modeling, the tsunami impact on prototype box-girder bridges located in Oregon, US, and developed appropriate design procedures that were included in the first tsunami design guidelines developed in the US. Hayatdavoodi and Ertekin [
28] and Hayatdavoodi et al. [
29]) employed the Navier–Stokes equations and Green–Naghdi equations to study the solitary wave impact on a horizontal deck. The second type is the particle-based method. Mohammad and Ali [
30] used the smoothed particle hydrodynamics (SPHs) method to investigate the tsunami wave loads on bridge decks. Hasanpour et al. [
31] used the SPHs method to study the impact of tsunami debris on bridge decks.
For the generation mechanism of a tsunami-like bore, a solitary wave has been frequently used to study tsunamis (see [
32,
33,
34,
35]) at their beginning. However, these studies do not explain why solitary waves or cnoidal waves are used to study tsunamis, as pointed out by [
3]. Many observations and much research suggest that the length and profile of tsunamis are different from those of solitary waves (see [
3,
36]). Moreover, most coastal structures are near or after the shore, where the waves become bores. The phenomenon of the disintegration of long waves into shorter waves has been observed in many tsunami events, such as the Indian Ocean tsunami in 2004, when the tsunami approached the shore; see [
37]. When the wavefront becomes sufficiently steep due to the change in water depth, it breaks into small waves and an undular bore. Hence, a dam-break-generated bore is employed as a tsunami-like bore in many tsunami studies (see [
38,
39]).
Madsen et al. [
3] investigated the propagation of a tsunami by considering a symmetric smooth wave propagating over a gentle ramp, which could induce the breaking of the wave. The conclusion is that the breaking has an additional effect, where the breaking happens, on bore impact on coastal structures. Some researchers have studied the solitary wave break during the propagation over a ramp and studied the impact on the structure in their studies, e.g., [
40]. Manoochehr et al. [
41] carried out laboratorial experiments to study the solitary breaking wave impact, where a solitary wave passes over a ramp with different coastal slopes, breaks, and propagates on a flat platform. The selection of the coastal slopes is due to the slope range tested by [
42]. Istrati [
43] broke the solitary wave with a ramp, whose slope is 1:12, in their large wave flume to study the impact on the box-girder. The onset of the breaking point and the location of the bore formation were identified. Yang et al. [
44] broke the solitary wave in the same way by a fixed slope to study the impact on a box-girder bridge. Neither of them discussed the relationship between the slopes and the breaking of the solitary wave but only focused on the impact, while the solitary experiments of [
45,
46] show that the ramp slope plays an important role in the breaking of a solitary wave.
In this study, we are interested in investigating which mechanism, whether a solitary wave, dam-break wave, or breaking solitary wave, should be more suitable as a tsunami-like bore. Only with the proper tsunami-like bore can the interaction between the wave and structure be confirmed with confidence.
This study concerns the impact of breaking solitary waves on a horizontal deck. As shown by [
3], the breaking has an additional effect on the bore impact on coastal structures. Hence, we studied two breaking waves, namely the solitary breaking wave and dam-break wave. A solitary wave propagates over a submerged ramp and breaks over the shelf, where a horizontal deck is located. In order to compare the breaking solitary wave with the dam-break bore, various slopes from 1:20 to 1:0 (step) were used to find a proper way to break the solitary wave.
The Navier–Stokes equations used to study this problem and the setup of the numerical experiments are introduced in
Section 2. The details of the method of breaking are presented in
Section 3. Two comparisons with existing laboratory experiments are provided, followed by the results of the simulations of the breaking solitary wave cases and the dam-break cases. The paper is closed with concluding remarks.
3. Bore Generated by Breaking of a Solitary Wave
In this study, we considered the impact of a bore on a deck. The bores are generated by (i) the breaking of a solitary wave and (ii) the breaking of the dam of a reservoir (dam-break).
As the initial shape and form of the solitary wave and the dam-break are different, a reasonable way to compare these two breaking mechanisms is introduced in this section.
Figure 1a shows a schematic of the numerical tank where the solitary-wave-breaking bore is studied. The solitary wave breaks when the downwave water depth decreases sufficiently. The schematic of the numerical tank of the dam-break case is shown in
Figure 1b.
In this study, is the upwave water depth, is the distance from the top ramp to the leading edge of the deck, is the slope of the ramp, where is the vertical height from the bottom of the upwave tank to the bottom of the downwave tank, and is the horizontal length from the trail of the upwave tank to the head of the downwave tank. In all solitary wave cases, the ramp top position is set at .
At the time
, the solitary wave propagates from upwave to downwave (left to right), and the solitary wave crest is at
, where
is the effective wavelength of the solitary wave.
is used as the criteria to determine the upwave and downwave boundary of
, which means the wave starts from the upwave point where
and ends at the downwave point where
. After propagating over the ramp, the solitary wave breaks into a bore with long tail waves. The shape of the initial solitary wave varies with
A and
; see [
52]. The generated bore varies according to the slope of the ramp. Appropriate values of
A,
, and
R are discussed in the following sections.
3.1. Equivalent Solitary Wave
According to Equation (
3), where
is the area between the surface elevation of the solitary wave and SWL (still-water level), the shape of the solitary wave is related to two variables, namely
A and
. Hence, two conditions are required to find the most appropriate solitary wave to break, and the shape of the generated bore should be as close as possible to the dam-break bore.
The first condition is that the mass of the initial water reservoir, , and the solitary wave, , should be the same, i.e., . As the density of water, , and the width of the tank, W (equal to the unit length in two dimensions), is constant in this study, and , , where is the cross-section area of the reservoir.
The area of the initial dam is
, where
is the initial height of the water reservoir measured from the SWL and
is the initial length of the reservoir, as shown in
Figure 1b.
The second condition is to set the overlapping area between the water reservoir and solitary wave,
, to be at maximum, which means that
shown in
Figure 3 should be minimum.
This results in obtaining the closest form between the solitary wave and dam-break.
is given as
With these two conditions, the two variables of the initial solitary wave are determined, and the closest agreement between the initial solitary wave and the water reservoir is obtained. This will allow for a fair comparison between bores generated by these two methods.
By considering these two conditions, i.e., using
and Equation (
4), three cases with proper
A and
, where
, were chosen for each solitary-wave-breaking case to compare with three different dam-break cases. These are listed in
Table 3.
3.2. Determination of the Ramp Slope
This section studies the solitary wave bore generation and propagation over the ramp. Twelve cases with different slopes were considered to find a proper ramp slope, listed in
Table 4.
The surface elevation is measured at G1 at
(top of the ramp), G2 at
, and G3 at
. Drastic changes between cases 10 and 8 are found. The surface elevations of cases 10 to 12 have many fluctuations after the bore front, whereas those of cases 1 to 8 are smoother. The results of four (sample) ramp-slope solitary wave cases, Cases 1, 2, 7, 10, with slopes
1:0, 1:1, 1:3, 1:5, and the corresponding dam-break case are shown in
Figure 4.
Figure 4 shows that the peak elevations of all solitary wave cases are higher than those of the dam-break case. Surface elevations of the case with
1:5 are the farthest from that of the other three cases. There are two peaks of the surface elevation in this case, whereas only one peak is found in the other solitary wave cases. Surface elevations of cases with
1:3, 1:1, 1:0 are close to each other. Based on these results, the case with
1:0 (i.e., step) agrees best with the dam-break case. Hence,
1:0 was chosen for the rest of the study.
3.3. Determination of the Deck Location
Bores generated by the breaking of a solitary wave or an initial reservoir are unsteady and change shape as they propagate downwave of the ramp. In this section, we discuss an appropriate location for the deck, where the profile of the two bores are close, such that a fair comparison between the solitary wave and reservoir case can be made. To find the position, comparisons were made between surface elevations recorded by gauges at in the dam-break cases and at in the solitary-wave-breaking cases. Three dam-break cases with were considered in this section, with corresponding three solitary-wave-breaking cases with , , respectively. There is no deck in this section.
In total, six cases were considered to investigate the appropriate deck location. The solitary wave surface elevations with various
,
A, and
were studied in these six cases. The results were compared with the dam-break surface elevation at
. For all three dam-break cases with
, and 8, the surface elevation of all three solitary wave cases with
is far from that of the respective dam-break cases, and that of the solitary wave cases with
is the closest one. The results of the solitary wave case with
and
and the dam-break case with
are shown in
Figure 5 for comparison.
Figure 5 shows that the solitary-wave surface elevation at
gives the closest agreement with the dam-break surface elevation before
, whereas the solitary-wave surface elevation at
is higher than dam-break before
. It is concluded that the solitary wave at
is more suitable for comparison with the dam-break bore.
Snapshots of the solitary wave propagating over a submerged step, (a) without a deck, or with a downwave deck whose leading edge is at (b)
, or (c)
, are shown in
Figure 6. The two dashed lines show the position of the leading edge of the deck at
in
Figure 6b and at
in
Figure 6c, respectively.
Figure 6a shows the solitary wave profile changes as it propagates over the submerged step at
. Beyond
, the solitary wave disperses into several small waves. The solitary wave profile changes as it propagates.
In
Figure 6b,c, the evolution of the solitary wave profile shows little change due to the presence of the submerged deck before meeting with the deck, which can be proved by the snapshots of the solitary wave taken between
and
in these sub-figures. The wave profiles of
Figure 6a,c are closer to each other than that of
Figure 6b. The wave profiles of
Figure 6a,c show some differences from each other after
when the wave passes over the deck, which means that the deck plays a critical role in the evolution of the solitary wave profile.
Consider the solitary wave shape, the ramp slope, and the deck location. Three groups of solitary-wave-breaking cases are compared against dam-break cases. The information on the groups is listed in
Table 5.
4. Comparisons with Experiments
Comparisons of the model’s results with available laboratory experiments are presented here. First, the surface elevations of a solitary wave propagating over a submerged step were considered, followed by the surface elevations and impacts of a solitary wave on a submerged deck.
Seabra-Santos et al. [
53] carried out experiments to study the propagation of solitary waves over a downwave step in the long channel of the Institut de Mécanique de Grenoble. The channel is 36 m long,
m wide, and
m deep. In the case considered here,
,
, and
1:0.
The computed and measured wave surface elevation at wave gauges
at
and
at
of [
53] were compared. The comparisons of the time series of surface elevation measured at different gauges are shown in
Figure 7. Excellent agreement is observed between the model’s results and the laboratory measurements, indicating that the solitary wave propagation over a step is captured correctly in our model.
Seiffert et al. [
20] carried out experiments to study the propagation and impact of a solitary wave on a flat deck. The flume is
m long,
m wide, and
m high. A horizontal deck of width
m (where the width is defined in the direction into the page), length
m (where the length is defined in the direction of wave propagation), thickness
m, and deck height
m was placed downwave. Three wave gauges,
(
m upwave of the leading edge of the deck),
(
m downwave the leading edge of the deck), and
(
m downwave the leading edge of the deck), were used to measure the wave elevations. The horizontal force on the deck was measured by three load cells, whereas the vertical force was measured by one cell. The tank floor was flat in this case.
The computed and measured surface elevations and forces of [
20] were compared with the result of this study. The results are shown in
Figure 8.
Figure 8a–c show the surface elevations of the solitary wave propagating over the submerged deck, and the results of NS agree well with the experimental data of [
20]. Overall, an excellent agreement is observed between the results.
The same analysis was performed on the vertical force and horizontal force results shown in
Figure 8d,e, and good agreement is observed.
5. Comparisons between Dam-Break and Breaking Solitary Wave
This section present the results of bore impact on a horizontal deck for a combination of two bore generation methods and five deck elevations from the seafloor. Solitary waves with three initial amplitudes were considered, and the results were compared with the equivalent dam-break bore. These include I: , II: , III: .
As discussed in the previous section, two locations of the deck were considered, namely and , for breaking solitary wave cases. Five vertical positions of the deck were considered, including the submerged positions (, and ), when the deck is located at the free surface, , and when the deck is above the free surface, . The solitary wave-generated bore impact results were compared with those of the dam-break bore.
Hence, fifteen (
) cases were considered in this study, as shown in
Table 5. They are solitary wave cases with the leading edge of the deck located at (i)
and (ii)
, and (iii) dam-break cases with the leading edge of the deck at
. Five different deck elevations were considered for each of the three cases.
5.1. I:
The time series of solitary wave horizontal and vertical forces on the deck is shown in
Figure 9 and
Figure 10, along with the dam-break horizontal and vertical forces for
, respectively. The peak values of the forces are listed in
Table 6.
As expected for both horizontal and vertical forces, the maximum value of the solitary-wave bore is larger when the deck is closer to the step. The maximum value due to the solitary wave (regardless of the deck’s location) is larger than that of the dam-break for all submerged cases, where . The difference between the solitary waves becomes smaller when the deck is closer to the free surface. In the elevated case, where , the value force of the solitary-wave bore is significantly smaller than that of the dam-break bore.
The previous work of [
54], which investigated the dam break bore impact on a bridge superstructure, found significant impulsive uplift forces followed by downward forces. The impulsive uplift component is non-existent or has a minor role in
Figure 9. However, it had a governing role in [
54].
The bore’s uplift force on the deck is tightly related to the bore shape. We can compare the bore shape of these two studies between
Figure 5 in this study and Figure 2 in [
54].
In
Figure 5 in this study, the bore’s front foot arrives at the deck at
, and the peak arrives at
, where the duration is 6. In Figure 2 in [
54], the bore’s front foot arrives at the deck at
, and the peak at
, where the duration is 1. The duration in this study is longer than that in [
54], which means that the bore’s front in [
54] is much steeper than that in this study. Hence, it is easy to find the impulsive uplift component.
Moreover, in this study, the rate between the reservoir and downstream water depth (
) and the vertical position of the deck (
) are variables. The vertical force of cases with
, and 1 show that no or small impulsive uplift component is reasonable, as the deck’s top is under or on the free surface. The rate between the reservoir and downstream water depth (
) of cases in this study is
, and
, which is smaller than the
of the case shown in Figure 2 in the study in [
54].
5.2. II:
The time series of solitary wave horizontal and vertical forces on the deck is shown in
Figure 11 and
Figure 12, along with the dam-break horizontal and vertical forces for
, respectively. The peak values of the forces are listed in
Table 7.
The differences between solitary waves and dam-break waves are smaller than that for
, especially for vertical force in
Figure 12. In the elevated case, the peak value of the horizontal force of solitary waves is larger than the dam-break wave when the deck is further away from the step.
5.3. III:
The time series of solitary wave horizontal and vertical forces on the deck is shown in
Figure 13 and
Figure 14, along with the dam-break horizontal and vertical forces for
, respectively. The peak values of the forces are listed in
Table 8.
The horizontal force of the solitary-wave bore is in closer agreement with the dam-break force in the submerged case, . In the unsubmerged case, the peak value of the horizontal force of solitary waves is larger than the dam-break wave when the deck is further away from the step, also found in II: .
The vertical force of the solitary-wave bore is smaller than the dam-break force, which is different from that in I: and II: .
6. Concluding Remarks
The impact of bores generated by breaking solitary waves over a submerged step was studied using computational fluid dynamics. The proper method for breaking a solitary wave is discussed. The tsunami-like bores generated by different methods were compared. In order to study the effect of the downwave slope on breaking solitary waves, a wide range of seafloor slopes, from 1:20 to 1:0 (step), was used. In the step case, the shape of the generated bore is closer to that of the dam-break bore. The bore is unsteady. It changes its form as it propagates. The propagation distance is important in the bore shape and impact on the downwave structures. An appropriate location for the deck is defined by considering a range of locations downwave of the step. The results show that a breaking solitary wave takes a longer distance to form the same bore shape as the dam-break bore. A comparable solitary wave to dam-break was studied by determining appropriate values for the wave amplitude and water depth. Given that the bore generated by the solitary wave is unsteady, the discussion is provided on the appropriate deck location to allow for a fair comparison.
For submerged structures, the force of the solitary-wave bore is larger than the dam-break bore. For structures above the free surface, the horizontal force of the solitary-wave case is smaller, whereas the vertical force of the solitary-wave bore is larger. By increasing the initial amplitude, an increase in impacts on the deck for solitary-wave bores and dam-break bores is expected. A more significant increase is observed for the impacts of solitary-wave cases when the deck is further away from the step and dam-break cases. The maximum impact of the solitary-wave bore is larger when the deck is closer to the step. The impact of the solitary wave is smoother when the deck is closer to the step, which shows that the breaking solitary wave is less turbulent at that position. In all submerged cases, where , the maximum value due to the solitary wave (regardless of the deck’s location) is larger than that of the dam-break. The difference becomes smaller when the deck is closer to the free surface. The force of the solitary-wave bore is in closer agreement with the dam-break force when the deck is further from the step. The results are complicated in the unsubmerged cases, where , as the free surface varies for different initial amplitudes. For a small initial amplitude, , the value force of the solitary-wave bore is significantly smaller than that of the dam-break bore. For a large initial amplitude, , the maximum horizontal force due to the solitary wave with a farther deck is larger than that of the dam-break case. In comparison, the maximum vertical force due to the solitary wave with a farther deck is smaller than that of the dam-break case.
Hence, a breaking solitary wave is suggested to be used as a tsunami-like bore in studying the tsunami impact on structures, as the peak value is larger than the dam-break wave when the amplitude is small. For large-amplitude tsunami conditions, dam-break appears to provide more conservative loading on the deck.
There are many ways to find equivalent solitary waves and dam-break bores. We considered the initial velocity, height, and even equivalent length. There is much work that can be carried out and discussed on this topic. Velocity is a good option, but the velocity of the dam-break bore is not constant during propagation. The velocity of a solitary wave after breaking is complex to determine. Velocity in this study was one of the most complex variables that we dealt with. It would be more complex if velocity were picked. This study tries to distinguish between the breaking solitary wave and the dam-breaking bore. Hence, we fixed the volume of the water, which was easier to be fixed in this study, and discuss the complicated variables such as velocity and force. This study will be carried on. Our further study will discuss all parameters, such as velocity and height.
Many works could be carried out in the future. More details about the difference between the dam-break bore and breaking solitary wave could be discussed by comparing the particle velocity around the horizontal deck. The overturning moment on the deck also plays an important role in the bridge’s destruction; see [
22,
55]. Simplified force equations could be derived for the impact of a bore generated by dam-break and a breaking solitary wave, such as [
56], which related the length of the deck with the solitary-wave-induced uplift forces. The dynamic fluid–structure interaction of bridges and the effect of structural flexibility could be considered, such as [
57,
58], which proved that this kind of structure can reduce the tsunami forces significantly. The various length scales and time scales involved in the bore’s evolution also need to be studied in the future.