Next Article in Journal
Numerical Study on Collisions of Solitons of Surface Waves in Finite Water Depth
Previous Article in Journal
Exact Solutions of Navier–Stokes Equations for Quasi-Two-Dimensional Flows with Rayleigh Friction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Flow around Circular and Rectangular Emergent Cylinders with Subcritical and Supercritical Conditions

by
Kannangara D. C. R. Dissanayaka
1 and
Norio Tanaka
2,*
1
Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura-Ku, Saitama-shi 338-8570, Japan
2
International Institute for Resilient Society, Saitama University, 255 Shimo-Okubo, Sakura-Ku, Saitama-shi 338-8570, Japan
*
Author to whom correspondence should be addressed.
Fluids 2023, 8(4), 124; https://doi.org/10.3390/fluids8040124
Submission received: 28 January 2023 / Revised: 27 March 2023 / Accepted: 29 March 2023 / Published: 3 April 2023

Abstract

:
There are multiple initiatives aimed at strengthening coastal communities against tsunami disaster risks, such as growing vegetation belts, construction of embankments, moats, and different hybrid alternatives. To find a solution for strengthening the coastal buildings themselves, we firstly reviewed the flow phenomena around a single emergent (circular and rectangular) cylinder (case C1), which was considered as a piloti-type column under different Froude conditions, and evaluated the formation of surface bow-waves, hydraulic jump detachment, and wall-jet-like bow-waves. Secondly, the flow characteristics were investigated under the same Froude conditions with side-by-side two-cylinder (case C2) and four-cylinder (case C4) arrays in an open channel. Surface bow-wave length ( L B w ) increased by 7–12% over the rectangular cylinders (RCs) compared to the circular cylinders (CCs) with a subcritical flow. For the supercritical flow with a 1/200 bed slope, hydraulic jump detachment was observed in relation to the Froude number. The observed length of the hydraulic jump detachment ( L j u m p ) varied between 3.1–8.5% and 4.2–12.9% for the CCs and RCs in the supercritical flow with a 1/200 bed slope. In addition, the wall-jet-like bow-wave height ( h j e t ) over the CCs was increased by 37% and 29% compared to the RCs with a supercritical flow and zero bed slope (orifice-type flow). For case C4, a hydraulic jump was observed for the supercritical flow over the horizontal channel bed. Finally, empirical equations were defined concerning the geometrical shape and arrangement based on the experiment data for the single and side-by-side configurations of the cylinders to validate the height of the wall-jet-like bow-wave as the most critical flow property.

1. Introduction

The Indian Ocean Tsunami (IOT) in 2004, Sulawesi Tsunami Indonesia (STI) in 2018, and Great East Japan Tsunami (GEJT) in 2011 destroyed numerous coastal structures. These were the most catastrophic failures ever recorded in human history [1,2,3]. After the GEJT, Japan’s Ministry of Land, Infrastructure, Transport, and Tourism (MLIT) divided the tsunami impact into two groups, referencing the recurrence interval: level 1 (recurrence interval of less than a hundred years) and level 2 (recurrence interval of several hundred to a thousand years) tsunamis [2]. Developing a mitigation strategy to safeguard coastal communities from level 2 tsunamis, which have the capability to overtop coastal embankments, has thus become essential. After the 2011 GEJT, the mitigation strategy and the proposed methodologies switched from a single defense system to a compound and/or hybrid mitigation system that integrates artificial and natural defenses against the hydrodynamic forces of the tsunami waves [3]. Numerous researchers have integrated different settings and numerical approaches into physical investigations. These include a combination of a coastal forest and moat [4]; double embankment systems [5]; a combination of an embankment, moat, and forest [6,7]; a mangrove forest [8]; a canal and dune [9]; a moat and/or canal behind the coastal embankment [10,11,12]; a multi-defensive line strategy [13]; vertical double-layer vegetation behind a coastal embankment [14]; and vegetation with different geometric conditions in the crest of a coastal embankment for energy reduction [15]. Moreover, coastal vegetation has been recognized as the most convenient natural approach that can be used to reduce the energy of tsunami waves [16,17,18]. The vegetation bio-shield can mitigate the hydrodynamic force of level 1 and level 2 tsunamis through a single mitigation approach. However, the only problem is the mitigation effect, as a single approach would be less beneficial during a disastrous situation [17]. Therefore, employing a hybrid defense system (HDS) or an alternative design for the coastal building structures must be investigated [1,2,3].
Cuomo et al. [19] examined the flow behavior around a building structure during a flood event. They highlighted the need for modern risk assessment techniques for research addressing critical infrastructures’ vulnerability. According to the computational and experimental assessments, the interaction of a tsunami traveling inland with the confined environment leads to a difference in flow depths and velocities, indicating the location affected from the Froude number [20,21,22,23]. Thus, an early study by Manawasekara et al. [23] highlighted the importance of an opening in the building frame for minimizing flood-induced loads, while Chock et al. [24] showed the need for a safe design against tsunami wave hydrodynamic loadings. Furthermore, the use of an internal break-away wall has been shown to be helpful in protecting against the collapse of buildings from increased loading during tsunamis and coastal flood catastrophes [25].
Consequently, during the hydrodynamic phase, when a break-away wall and internal partitions have eliminated the impact forces, the influences of the building design and its structural strength become critical [26]. Investigations in the aftermath of the 2018 STI also showed that elevated buildings, such as vertical parking garages and commercial buildings with multiple stories with porous construction, saved many lives [27]. Consequently, it would be advisable to elevate building structures by adding piloti-type columns up to a particular height, and the minimum height should be equal to the height of a single floor. Previous studies concentrated on the loading process over the building structures and the impact forces produced by tsunami-like surges. However, evaluating the flow behavior around these piloti-type columns in piloti-type building structures under different Froude conditions is compulsory to produce safer structures regarding hydrodynamic loading and the scouring process [23,24,28]. Yamamoto et al. [29] evaluated the flood damage throughout Japan and suggested that a new approach to develop a resilient environment in coastal and riverine regions by using piloti-type buildings would be the best solution for the adaptation measures. Dang-Vu et al. [30] also evaluated the influence of introducing additional shear walls, which can be constructed on the first-floor level, as a strengthening approach for piloti-type buildings subjected to earthquake loadings.
The hydrodynamic forces acting on an impermeable structural system in an open-channel flow are well-characterized under various flow conditions [27,31,32,33,34,35]. Mignot and Riviere [32] and Riviere et al. [35] focused on the bow waves’ key characteristics and the features of the flow behavior around a blunt rectangular obstacle in a supercritical flow. Mignot and Riviere [32] and Riviere et al. [35] also evaluated the formation of the detachment hydraulic jump and wall-jet-like bow-wave over the emergent obstacle’s front face while changing the flow depth, obstacle geometry, and channel slope. However, there is an unknown factor that is yet to be confirmed when the approach flow is a kind of tsunami inundation current with subcritical and supercritical characteristics through the coastal structures, especially with regard to piloti-type buildings, as Yamamoto et al. [29] and Samarasekara et al. [36] have explained.
Therefore, this study aimed to characterize the flow phenomena over single circular and rectangular emergent cylinders, assuming a piloti-type column, in a steady-state flow under different Froude number and channel conditions to address the subcritical and supercritical flow behaviors of an approaching tsunami flow current along the shoreline. Additionally, emergent circular cylindrical arrays were tested as piloti-type column groups, assuming they were connected to a coastal building under the same Froude conditions. The primary purpose of evaluating the flow fields over the piloti-type single cylinder or cylindrical arrays was to define the maximum possible spacing between the ground surface and the first ceiling level of the piloti-type building in order to propose more beneficial and safer designs for resilient structures as a resilience approach in the future.

2. Materials and Methods

2.1. Experimental Facility

Figure 1 shows the experimental setup arranged inside an open channel over a smoothly coated plywood bed, which was used to avoid frictional losses. The adjustable gate fixed upstream near the inlet regulated the flow with the different experimental characteristics to achieve supercritical flow conditions. The channel width ( B ) was 0.5 m, and the cylinders were selected with the parameters D 60 mm and W 30 mm to avoid lateral confinement effects from the side walls. Table 1 summarizes the preliminary characteristics for the experiments in this study.
The definitions of the notations shown in Figure 1a are as follows: L is the working length of the channel (=12.5 m); L u is the distance to the emergent cylindrical obstacle from the vertical gate (=2.8 m); L B is the length of the fixed bed installed inside the channel (=6.5 m); H B is the height of the bed fixed inside the channel (=0.1 m); Q i n  is the inlet flow discharge; the upstream water depth ( h o ) was measured at a distance equal to 10 D (or 10 W ) from the upstream face of the emergent cylinder; the mean velocity of the flow was obtained with U ( = Q i n / B h o ) ; H is the depth of the channel (=0.4 m); D is the diameter of the circular cylinder; and W is the width of the rectangular cylinder. Changing the inlet discharge, channel slope ( S o ), and cylinder dimensions ( D or W ) made it possible to independently modify the approach flow conditions h o / D (or h o / W ) and the Froude number ( F r ). The details for each experimental case considered in the experiments are presented in Table 2.
As illustrated in Figure 1b,c, the same procedure was followed with the side-by-side emergent cylinder array placed inside the experimental flume. Select side-by-side emergent cylinder arrays were produced using 32 mm diameter circular cylinders, which represented the piloti-type column pairs. The notations shown inside the cylinder arrays (case C2 and case C4) in Figure 1b,c indicate the points where water depths were observed at the defined locations along the center line C 0 , C 1 , and C 2 .

2.2. Instrumentation and Flow Observations

In the experiment facility, an electromagnetic flow meter (EFM) attached to the closed circuit was used to measure the channel’s inlet flow rate. Flow was controlled using a computer-driven software package called HYDRA (v2.02), which connected to the circulating pump of the experimental flume. Water depth upstream of the emergent cylinder was measured using a point gauge, as shown in Figure 1, at a distance of 10 D (or 10 W ) for both the circular and rectangular emergent cylinders. The uncertainty values regarding the discharge measurement ( Q i n ) of the flow meter in the channel and the water depth measurement upstream ( h o ) were 0.005 L/min and 0.01 mm, respectively. The pointwise x-directional velocity in each case was recorded using an electromagnetic flow current meter (EFCM), which had an uncertainty ( U ) of about 0.01 cm/s. Three high-resolution cameras were used to capture the flow variation inside the experimental flume, as shown in Figure 1a. Finally, an image analysis technique was used to analyze the images to calculate the average water surface elevation over the emergent cylinders in each case under both the subcritical and supercritical flow conditions. To calculate the average water depth value, images were captured over a period of 30 s by the cameras from the respective direction such that 30 pictures were used in each case (Figure 1a).

2.3. Dimensional Analysis

The experiment cases tabulated in Table 2 represent piloti-type columns with different geometries and arrangements under subcritical and supercritical flow conditions over horizontal and sloping beds with steady-state flow (Figure 1). The following dimensional variables were used to define the flow: upstream water depth ( h o ), upstream velocity ( U ), emergent cylinders’ dimensions (diameter D or width W ), properties of water (dynamic viscosity μ , density ρ , and the surface tension with air σ ), bed and obstacle roughness ( k s ), and gravitational acceleration ( g ). The obstruction length was ignored for the emergent rectangular cylinders due to its lower contribution to the flow structures, whether subcritical or supercritical [31,32,37]. This research focused on flow features, such as surface bow-waves, hydraulic jump detachment, and wall-jet-like bow-waves. Regarding the flow circumstances, these flow characteristics were dimensionally defined to develop an empirical relationship based on multiple linear regression analysis [32,35]. Equation (1) was used to characterize the primary parameters selected for the dimensional analysis of the cylindrical column.
( flow   form ,   wall - jet ) = f ( h o ,   h F ,   D ,   U ,   ρ ,   μ ,   g ,   σ ,   k s )
Vaschy–Buckingham’s pi theorem was used to develop the empirical relationship, taking into account the fundamentals of the flow behavior (Figure 1 and Table 1). The most important flow type observed in the entire experiment series was the wall-jet-like bow-wave along the front face of the cylinder [32,35] under the supercritical flow conditions. Yamamoto et al. [29] and Dang-Vu et al. [30] have highlighted the importance of determining the level of elevation and its influence on the seismic characteristics when evaluating the stability of a building structure. Hence, this paper describes an empirical relationship by considering only this wall-jet-like bow-wave flow condition.

2.3.1. Single Emergent Cylinder

For the supercritical flow, Equation (2) describes the identified non-dimensional relationship for case C1 concerning the wall-jet characteristics in the present experiment. The height of the wall-jet ( h j e t ) was made dimensionless relative to the kinetic height, and it finally became h j e t * = h j e t / U 2 / 2 g . Furthermore, the formation of the wall-jet-like bow-wave ( h j e t ) was identical to the supercritical flow with horizontal and sloping beds. The only difference was that the wall-jet height over the front surface in the sloping channel bed was lower than the supercritical flow with a horizontal channel bed [38,39,40].
h j e t * = f F r , R e , h o D

2.3.2. Side-by-Side Two-Sided and Four-Sided Emergent Cylinder Arrays

When a subcritical or supercritical flow passes through an emergent cylinder array (cases C2 and C4), the primary parameter is the water depth in front of the cylinder array. Referring to Equation (3) for the cylinder array, the non-dimensional wall-jet height ( h j e t * ) can be explained by the spacing ( S Y ) in the transverse direction (Figure 1). The Froude number ( F r ), upstream water depth ( h o ), and water depth at the center of the Y-Y axis at the first column pair ( h C ) were used to define the wall-jet height ( h j e t ) . The h C value was denoted as C o for case C2 and C 1 for case C4.
h j e t * = f F r , R e , h o S Y , h C S Y

3. Results and Discussion

3.1. Subcritical Flow through an Emergent Cylinder

In the subcritical flow, the Reynolds number ( R e ) varied in the range of 8900–23,600. A surface bow-wave was generated at the front surface of the cylinder when the subcritical flow passed through case C1. For the flows where F r < 1, referring to Figure 2a, only the surface bow-wave was apparent in the subcritical regime, and the flow deflection originated upstream from case C1, as shown in Figure 2b. When the h o / D or h o / W ratio became smaller, a larger surface bow-wave in terms of flow depth ( h o ) and geometry could be observed. Figure 2a shows both the CC and RC. The distance from the front face of the cylinder to the surface bow-wave toe ( L b w ) was visible with the streamlined separation. The L B w value increased with increasing water depth and the cylinder geometry ( D or W ). The non-dimensional surface bow-wave length L B w * (= L B w / D or L B w / W ) is tabulated in Table 3 for each geometry considered in the present experiment. When the subcritical flow passed through cases C2 and C4, similar behavior was observed over the front surface, such as the formation of a surface bow-wave. Figure 2c shows the surface bow-wave formation for case C2. As shown in Figure 2c,d, the parabolic pathway of this surface bow-wave merged at the center line along the flow direction when the flow passed through cases C2 and C4.
As shown in Figure 2b, when changing the geometry of the cylinder ( D or W ), the length of the surface bow-wave ( L B w ) remains proportional to the h o / D (or h o / W ) [32,35,41]. When the diameter or width of the emergent cylinder increases, the surface bow-wave’s length ( L B w ) increases due to the blockage ratio.
The importance of determining the interaction in the formation of the bow-wave is, as McDonald [42] has explained, due to the fact that the bow wave can be characterized as steady, unsteady, or quasi-steady. In the subcritical flow condition, due to the breaking dynamics of the incoming flow imposed by the surface tension, when the obstacle becomes larger, a larger bow-wave is observed [41,42,43]. Hence, reducing the interaction in the bow-wave would be beneficial for lowering the scouring and stability near the structure in terms of the Froude conditions of the incoming flow [29,30,32,36,41,42,43,44].

3.2. Supercritical Flow through Emergent Cylinders

For the supercritical flow, the primary concern was the formation of the wall-jet-like bow-wave. When evaluating the wall jets’ behavior, there was a scale effect with respect to the field-measured data as the entire experimental setup considered the field-measured data after the post-tsunami survey to finalize the experiment scale [1,2,27,45,46,47]. Furthermore, based on the overall objective of the experiment, deciding on a minimum elevation for the building structure was a secondary solution following the defensive measures on the coastline in the case of a tsunami flow approaching [1,2,3,16,29,30,48]. Hence, the laboratory-produced wall-jet-like bow-waves seemed to be extremely similar to those seen in nature and represented in this study. However, when the upstream water depth ( h o ) or the geometry of the cylinder ( D or W ) decreased significantly, scale effects and the capillary effects became apparent. Therefore, the scale effect in the experimental setup was checked with the supercritical flow, both with the channel bed kept horizontal and with a slope, to confirm that pressure affected the wall-jet height and its characteristics, as discussed by Riviere et al. [35]. Furthermore, in order to prevent any lateral-confinement effects from the side walls of the channel, flows with high h o / D or h o / W were created in the channel by employing small-size CCs and RCs, as tabulated in Table 1 [35,44,49]. Moreover, extended droplets from the wall jet’s tip were seen in the vicinity of the front face of the emergent cylinder in the present experiment due to the capillary effects. These minor capillary and wall-jet changes did not affect the recorded time-averaged water depth in front of the emergent cylinders [32,35]. Furthermore, the capillary effect does not affect the flow transition from a surface bow-wave or detached hydraulic jump to a wall jet, even in laboratory-scale experiments [32,35,37,41,42,43].

3.2.1. Supercritical Flow over the Sloping Channel Bed

The bed slope of the channel was set to 1/200 to achieve the supercritical flow condition, replicating the tsunami flow over sloping ground [27,32,35,36,48]. Under the sloping bed condition, the Reynold’s number ranged between 10,700 and 35,700. Each physical model (case C1, case C2, and case C4) was kept inside the experimental flume to identify the flow structures and characteristics. As shown in Figure 3, when F r > 1 and the channel had a 1/200 bed slope, wall jets were demonstrated over the front face of the cylindrical columns in case C1 with increasing upstream water depth ( h o ). In case C1, for the high h o / D or h o / W ratios corresponding to smaller cylinder size, the wall jet could be evacuated as a parallel jet where the flow was blocked by the CC or RC columns with a relatively small detached hydraulic jump [31,32,35,41,43]. However, for the low h o / D or h o / W ratios, which corresponded to larger cylinder dimensions in relation to the upstream water depth ( h o ), wall jets occurred with a more significant detached hydraulic jump due to the slight increase in water depth near the upstream face resulting from the bluntness of the emergent cylinder.
As shown in Figure 4a,b for case C1 over a 1/200 sloping bed, with high h o / D and h o / W ratios, wall-jet-like bow-waves ( h j e t ) with a minor detachment length occurred (Figure 3). Furthermore, with low h o / D ratios (with 60 mm dia. CCs), small wall jets with a more significant detached hydraulic jump were observed [41,43]. The formation of the hydraulic jump detachment and the variations in its length ( L j u m p ) for case C1 relative to the geometry ( D or W ) for the 1/200 sloping bed are tabulated in Table 4. As shown in Table 4, the non-dimensional hydraulic jump detachment length L j u m p / D (or L j u m p / W ) relative to the geometry demonstrated higher values for the RCs than the CCs. From a general perspective, when F r > 1 for supercritical flows, the wall-jet ( h j e t ) height is the most critical parameter for the channel horizontal and 1/200 sloping bed cases as uplift forces can be created over the ceiling level of a piloti-type building [29,30,35].
Table 4 summarizes the non-dimensional wall-jet height h j e t / D (or h j e t / W ) relative to the geometry for the RCs and CCs. For cases C2 and C4, the hydraulic jump detachment resembled a surface bow-wave due to the increased water depth at the front. Figure 4c,d show the side views of the flow profile for cases C2 and C4, respectively. In cases C2 and C4, a relatively low wall-jet-like bow-wave height was observed, and the surface bow-wave was not similar to that in case C1 over the 1/200 sloping channel bed, as shown in Figure 4. Furthermore, the height of the wall jet ( h j e t ) slightly increased with an increase in the S Y spacing, as shown in Table 5.
Table 5 shows the average changes in the non-dimensional surface bow-wave length L B w * (= L B w / S Y ) and the non-dimensional wall-jet bow-wave height h j e t * (= h j e t / S Y ) calculated relative to the S Y spacing. Arnaud et al. [50] explained that, when a tsunami flows through a cylinder array (porous structure), due to the refraction and diffraction effects resulting from the adjacent cylinders in the array, the wall-jet-like bow-waves demonstrate minor changes, as seen in the results obtained for case C1.

3.2.2. Supercritical Flow over the Horizontal Channel Bed

The supercritical flow over the horizontal bed was achieved using a gate fixed at the channel’s inlet, as shown in Figure 1. The flow reached a supercritical condition after exiting from the gate, following a vena contracta at the gate exit, without any transition to the front face of the cylinder in each case. The water depth ( h 0 ) and Froude number ( F r ) used under this flow condition were selected in order to maintain the flow transition represented by the formation of a hydraulic jump downstream in all cases in the present experiment and to achieve the described fluid conditions [31,32,34,35]. This flow replicated the flow immediately after a tsunami crosses a shoreline and flows through the frontline buildings following the mitigation measures, attacking the frontal structures on its flow path. In this flow condition, the Reynolds number ( R e ) varied in the range of 20,400–44,800.
Figure 5 shows the results related to flow forms when F r > 1 and when the channel slope was kept horizontal with flow controlled by a sluice gate, as shown in Figure 1. Due to the supercritical condition and the kinetic energy of the flow, for both the high and low h o / D and h o / W ratios, wall jets were observed over the front face for case C1. When the flow reached the obstacle front, it was converted to a potential head from a kinetic head when Q o u t > Q i n , as it was blocked by the CC or RC, as shown in Figure 6a [35,37,38]. Figure 6 shows that, when the supercritical flow reached the cylinder front in case C1, case C2, and case C4, its kinetic energy was converted to a wall-jet-like bow-wave and dissipated as it skirted the emergent cylinder through the formation of a lateral jet [32,35]. The height of the wall-jet bow-wave was several times that of the upstream water depth ( h o ). Therefore, the non-dimensional wall-jet height h j e t * ( = h j e t / h o ) was evaluated relative to the upstream water depth ( h o ). Table 6 shows the non-dimensional wall-jet height ( h j e t / h o ) relative to the upstream water depth ( h o ) for the supercritical flow over the horizontal channel bed.
Furthermore, we observed a slight, abrupt water depth increase over the RC compared to the CC at the toe of the wall jet, as shown in Figure 6a. For the CCs, when increasing the cylinder diameter, the contact area for the flow also increased, but, due to its shape, no abrupt increase in water depth occurred at the front except for the wall jet. As shown in Figure 6a, the supercritical flow in case C1 was skirted with a smaller aerated region over the CCs compared to the RCs due to its geometry [35]. A similar flow phenomenon occurred at the front face for the supercritical flow in case C2, as shown in Figure 6b. With the supercritical flow in case C4, a small bow-wave was formed over the cylinders’ front face in the second row, as shown in Figure 6c. Figure 6d shows the skirting of the wall jet around the circular cylindrical column surface and the dissipation from both sides as lateral jets in case C4.

3.3. Variation in the Depth of the Flow through the Side-by-Side Cylinder Array

Concerning the results of the present experiment, significant differences were observed in the flow characteristics inside the cylinder array, including in case C2 and case C4, for both the subcritical and supercritical flows. Among the Froude characteristics selected in the present experiment, supercritical flow with a horizontal channel bed was given extra attention due to the behavior of the flow in case C4. The flow depth in case C4 changed along the center line for the supercritical flow over the horizontal channel bed. In case C4, water depth drastically changed at the points C o , C 1 , and C 2 in relation to the spacing in the y-direction ( S Y ) and x-direction ( S X ), as shown in Figure 1c. The supercritical flow formed a hydraulic jump inside the cylinder array. This hydraulic jump formation would create an additional upthrust over the bottom floor of piloti-type structures during a supercritical overland tsunami flow [28,29]. Figure 7 shows the variation in the non-dimensional water depth relative to the upstream water depth measured at each point considered in case C4 (see Figure 1c).
The percentage increment in the water depth for the supercritical flow over the horizontal channel bed in case C4-1 at point C 1 ranged between 28.8 and 34.5%, while at point C o , it ranged between 36.4 and 45.7%, and at point C 2 , it ranged between 30.5 and 37.5%. Further, for case C4-2 with supercritical flow over the horizontal channel bed, the percentage increments in the water depth at points C 1 , C o , and C 2 were 27.7–30.5%, 32.4–42.0%, and 26.7–33.0%, respectively. This scenario showed the opposite results as for the supercritical flow over a 1/200 sloping channel bed. Due to this phenomenon over the sloping bed, the flow became hypercritical downstream following the cylinder arrays in cases C2 and C4 [39,40,51,52].

3.4. Development of the Conceptual Models and Explanation of the Flow Transition

Out of the three flow characteristics considered in the present experiment, the most critical flow type in all the cases was the formation of the wall-jet-like bow-wave over the front face of the emergent cylinders. Furthermore, in an elevated building, the wall-jet-like bow-wave produces an additional uplifting pressure force over the superstructure due to the water depth increase inside the cylinder array, as explained in Section 3.3 for the supercritical flow [28]. This increase in water depth inside the cylinder array (case C4) further increases the risk affecting the stability of the building structure [23,24,27,28,35,45,46,53,54]. Hence, developing an empirical equation to predict this wall-jet height would be necessary for sustainable structural design in the future.

3.4.1. Single Emergent Cylinders in the Supercritical Flow

A similar phenomenon as with the supercritical flow over the horizontal bed was observed for the supercritical flow with a 1/200 bed slope, but the height of the wall jet was lower than in the horizontal bed case due to the flow transition in front of the cylinder. In both the supercritical flows considered, the wall jet dissipated by skirting the cylinder surface and was converted to a lateral jet after reaching its maximum possible height. Figure 8a represents the time-average relative wall-jet height h j e t * = h j e t / U 2 / 2 g as a function of the upstream Froude number observed over the emergent cylinders in case C1 for the supercritical flow with both horizontal and 1/200 sloping beds. As shown in Figure 8a, most of the values of the h j e t *   w e r e   >   1 . This was because the surface velocity over the front face of the cylinder was more significant than the mean upstream velocity ( U ). Utilizing Graf and Altinakar’s [52] explanation, the upper bound of the experimental values for the non-dimensional wall-jet height ( h j e t * ) was equal to 2.071 ≈ (1.15/0.8)2 in the present experiment, as shown in Figure 8a by the dashed line. In Figure 8a, open symbols refer to the supercritical flow over a 1/200 bed slope, and closed symbols refer to the supercritical flow with a horizontal bed in case C1. Equation (4) was developed to predict the dimensionless wall-jet height ( h j e t * ) based on the present experimental results for the supercritical flow with a horizontal bed in case C1. Also, for the rectangular cylinders, term D in Equation (4) is replaced as term W . Equation (4) was subsequently validated for the supercritical flow data observed in the 1/200 bed-slope case.
h j e t * = 1.72 F r 0.38 h o D 0.33
According to the t-test results, the p-values of the selected dimensionless parameters F r and ( h o / D ) in Equation (4) were lower than 0.05, which meant they satisfied the null hypothesis. The values for the experimental dimensionless wall-jet height versus those for the predicted wall-jet height obtained with Equation (4) are presented in Figure 8b. They indicated that the selected parameters in Equation (4) explained all the response data around the mean. For the circular cylinders, the average percentage errors in the predictions of the non-dimensional wall-jet height using Equation (4) were −0.25% and −6.79% for the supercritical flow over a horizontal bed and over a 1/200 sloping bed, respectively. Each case’s correlation coefficient is also shown in Figure 8b.
Moreover, Equation (4) was tested to predict the wall-jet height over the emergent rectangular cylinders, and the correlation plots for the predicted versus experimental values for each case, along with the correlation coefficients, are shown in Figure 8c for both the horizontal and 1/200 sloping bed supercritical flow conditions in case C1. The applicability of Equation (4) for the prediction of the dimensionless wall-jet height for supercritical flow over horizontal and sloping beds was evaluated in terms of the percentage of error. The average error percentages were 0.27% and −7.51% for the supercritical flow over a horizontal bed and the supercritical flow over a 1/200 sloping bed, respectively.
Additionally, the empirical equation developed by Riviere et al. [35], adopted and as Equation (5) was tested with modifications in the present experiment for a circular and rectangular emergent cylinders in case C1. In Equation (5), h is upstream water depth, F is Froude number, and R is obstacle width. To apply the present experiment, the term R in Equation (5) is replaced by D for the circular cylinders and replaced by W for the rectangular cylinders in case C1, respcetively.
h j e t * = 1.7 h R 0.24 F [ 0.24 h / R 0.17 ]
Referencing Equation (5), the correlation plot for the experimental and predicted dimensionless wall-jet height is shown in Figure 8d. The goodness of fit ( R 2 ) value for the calculated dimensionless wall-jet height was 0.977 for the supercritical flow when the channel bed was kept horizontal and 0.954 for the supercritical flow when the channel bed sloped by 1/200. Furthermore, the average percentage errors in the prediction of the dimensionless wall-jet height using Equation (5) for the circular and rectangular emergent cylinders in case C1 were 3.24% and −3.76%, respectively.

3.4.2. Side-by-Side Emergent Cylinder Array in the Supercritical Flow

Similar to the explanation in Section 3.2 and Section 3.3, when water flowed through the cylinder arrays in cases C2 and C4, the formation of the wall jet depended on the parameters given in Equation (6). The parameters are also shown in Equation (4) for the formation of the wall-jet-like bow-wave for the supercritical flow in cases C2 and C4. The water depth between the cylinders ( h C ) along the centerline was considered in relation to the spacing between the cylinders ( S Y ), as shown in Figure 1b,c. Equation (6) was formulated using the dataset for the wall jet with a supercritical flow passing through the cylinder array on the horizontal channel bed. Then, Equation (6) was validated with the 1/200 sloping channel bed dataset to predict the dimensionless wall-jet height.
h j e t * = 3.89 F r 2.15 h o D 0.85 h C S Y 0.05
Figure 9a,b show the correlation plots for the predicted wall-jet height compared to the experimental results for the supercritical flow with horizontal and 1/200 sloping channel beds for both cases C2 and C4. For each case, the goodness of the decision coefficient ( R 2 ) is shown on the graph. According to the plotted graphs shown in Figure 9 for cases C2 and C4, the non-dimensional wall-jet height over the sloping bed was slightly underestimated. The average percentage errors for case C2 in the results predicted by Equation (6) were 2.1% and −7.8% for the horizontal and 1/200 sloping beds, respectively. Furthermore, the average percentage errors for case C4 in the results predicted by Equation (6) were 3.4% and −4.5% for the horizontal and 1/200 sloping beds, respectively.

4. Summary

This study primarily examined the characteristics of subcritical flow, supercritical flow with a horizontal bed, and supercritical flow with a sloping bed around single circular and rectangular emergent cylinders in an open channel. Furthermore, side-by-side circular cylinder arrays (Figure 1) with varying spacings, as shown in Table 1, were tested to evaluate the flow characteristics with the same Froude number conditions. Three types of flow were observed during the experiment: surface bow-waves, detached hydraulic jumps, and wall-jet-like bow-waves [32,35,41,43,44]. The surface bow-wave appeared over the front face of the cylinders due to the subcritical flow condition with the single (case C1) and side-by-side emergent cylinder arrays (cases C2 and C4). This was the reason for the streamlined curvature of the flow, which flowed around the cylinders (Figure 2), and the length of the surface bow-wave ( L B w ) changed in relation to the geometry of the emergent cylinder. The development of the surface bow-wave length ( L B w ) could be altered by changing the h o / D (or h o / W ) ratio in relation to the Froude condition in the subcritical flow with a single emergent cylinder or the h C / S X ratio for side-by-side cylinder arrays.
Among these three flow types, wall-jet-like bow-waves were given extra attention in terms of the Froude number and the channel conditions. The wall-jet-like bow-wave was observed in the supercritical flow when the channel bed was kept horizontal and when it was sloped. The supercritical flow with the horizontal bed reached the emergent cylinder without any transition, and on the sloping bed, due to the bluntness of the cylinder, there was a water depth increase upstream resulting from the detached hydraulic jump in addition to the wall jet observed [28,32,41,43]. A wall-jet-like bow-wave was observed for the supercritical flow when it passed through the side-by-side cylinder arrays [35,38,41,51,55]. However, the shape of the emergent cylinder may be changed to correct the existing findings quantitatively. Placing more streamlined obstacles in contact with a supercritical flow, such as bridge piers or piloti-type columns with circular or rectangular cross-sections, may enable flow evacuation with the presence of a wall-jet-like bow-wave, as previously explained.
Following the post-tsunami surveys, the experimental setup employed field-measured data to determine the minimum elevation for building structures above ground level that would maintain a free surface beneath them as a secondary solution following defensive measures on the coastline in the case of a tsunami approach [3,56]. Hence, the laboratory-produced wall-jet-like bow-waves seemed to be extremely similar to those seen in nature and represented in Figure 4 and Figure 6. However, when the upstream water depth ( h o ) or cylinder geometry parameters ( D or W ) decreased significantly, scale and capillary effects became apparent. The scale effect in the experimental setup was checked to confirm that the pressure affected the wall-jet height and its characteristics, as discussed by Riviere et al. [35]. Due to the capillary effects, extended droplets from the wall jet’s tip were seen in the vicinity of such obstacles. These minor capillary and wall-jet changes did not affect the recorded average water depth in front of the cylinders. Such capillary effects are insignificant at the actual scale relevant to hydraulic engineering applications. They do not affect the form of the transition from a wall jet to a surface bow-wave or a detached hydraulic jump, even in laboratory-scale experiments [32,35,44].
As seen in the sloping roadways during urban floods [26,53] or the overland flow following a tsunami catastrophe as a tsunami run-up, different obstacles may respond in different ways to a supercritical flow [32,35]. Inland flow eventually reaches a supercritical regime after a tsunami and interacts with concurrently emergent man-made (buildings, automobiles, etc.) and natural (hills, etc.) obstructions [45,46,54,56]. The supercritical inflow must, therefore, skirt these obstacles on all sides since they are impermeable obstructions mounted on the base/bed and emerging from the free surface. This demonstrates that the primary contradiction may be explained by one aspect of the open-channel flow: its vertical confinement between the free surface and the bottom. In subcritical and supercritical flows, the flow may evacuate the constrained region between the bed and the free surface through a wall jet or a surface bow-wave/detached hydraulic jump [28,29,32,35]. Furthermore, without the need for a horizontal streamline curve, a supercritical flow workaround may occur beyond the free stream. However, sometimes, particularly when the cylinder/obstacle diameter ( D ) or width ( W ) are much greater than the upstream water depth ( h o ), a wall jet may not always be able to clear the whole volume of the water blocked by the obstacle [28]. This phenomenon replaces the wall-jet bow-wave with a hydraulic jump detachment. However, to form a hydraulic jump detachment, the cross-stream width of the obstruction should be larger than the upstream water depth ( h o ) in the supercritical flow [35]. For the supercritical flows, by altering the h o / D (or h o / W ) ratio, the Froude characteristics of the formation of the detached hydraulic jump and the wall-jet-like bow-wave can be adjusted.
However, the current findings may be quantitatively adjusted by modifying the geometry of the emerging cylindrical obstacles, including the single and side-by-side arrangements. The transition’s threshold ratio h o / D or h o / W is reduced with the Froude number, confirming the theoretical hypothesis that the change is triggered by mass conservation [32,34,35]. Concentrating on the wall-jet-like bow-wave, the empirical equations presented in this article show that mass conservation also impacts the height of the wall jet above the free surface flow, which was determined using F r and h o / D or h o / W . Lastly, the wall jet showed height oscillations, and these oscillations were likely due to the reversed spillage over the upstream face of the obstacle and related to the rising and lowering water cycle. From a general perspective, the determination of the wall-jet height and its prediction using the empirical equations defined in this article will be beneficial for the identification of relationships between wall jets and the upstream water depth ( h o ) in the construction of piloti-type buildings as a resilience approach against tsunami destructiveness in the future.

5. Conclusions

In the present study, we conducted a model-scale experiment to investigate the flow characteristics over emergent single and side-by-side cylinder arrays in an open channel with different Froude numbers and channel slopes. The Froude numbers lay within the subcritical to supercritical zones and different cylinder sizes and water depth ratios ( h o / D or h o / W ) were used, including a side-by-side cylinder arrangement with different spacings in the x and y directions ( S X and S Y ). The following conclusions were derived after completing a series of experiments in the present study:
  • The height of the wall-jet ( h j e t ) was altered with the changes in the h o / D (or h o / W ) ratio in case C1 and with the changes in the h C / S X ratio in the side-by-side cylinder arrays (cases C2 and C4) for both the supercritical conditions considered;
  • The scale effects did not influence the wall-jet height ( h j e t ) in the laboratory-scale experiment, and the extended droplets observed at the wall jets’ tips due to the capillary effect also did not affect the time-averaged wall-jet height;
  • For the supercritical flows, when the h o / D (or h o / W ) ratio was > 1, a wall-jet-like bow-wave could develop in front of the obstacle. Furthermore, when the h o / D (or h o / W ) ratio was <1, a distinct detached hydraulic jump could occur in front of the obstacle with different geometrical characteristics;
  • For the supercritical flow, the flow depth could increase within the side-by-side square cylinder array (case C4) in relation to the x and y directional spacing in the cylinder array and the sloping characteristics of the experimental flume;
  • The derived equations were used to predict the non-dimensional wall-jet height with the supercritical flow over the single and side-by-side cylinder arrays, which can be helpful in designing piloti-type building structures in urban coastal environments as a resilience measure against future tsunami risks.
In the present experiment, due to the lack of facilities, the S X and S Y spacings were limited to two cases for the side-by-side cylinder arrays to avoid sidewall effects from the experimental flume on the flow structure. Hence, additional research is needed with variations in the S X and S Y spacings to evaluate the flow phenomena inside the cylinder arrays further. Moreover, flow following a coastal embankment due to tsunami overtopping creates downstream subcritical and supercritical flow conditions. In addition, emergent single and side-by-side cylinder arrays can be tested with a flow following the overtopped coastal embankment.

Author Contributions

Conceptualization, K.D.C.R.D. and N.T.; Methodology, K.D.C.R.D.; Model preparation, K.D.C.R.D.; Investigation, K.D.C.R.D.; Data processing and analysis, K.D.C.R.D.; Writing—original draft, K.D.C.R.D.; Visualization, K.D.C.R.D.; Resources, N.T.; Writing—review and editing, N.T.; Supervision, N.T.; Project administration, N.T.; Funding acquisition, N.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors appreciate the support from the Japanese Ministry of Education, Culture, Sports, Science, and Technology (Monbukagakusho–MEXT Scholarship). The authors also acknowledge the anonymous reviewers for their valuable comments to improve this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suppasri, A.; Mas, E.; Koshimura, S.; Imai, K.; Harada, K.; Imamura, F. Developing Tsunami Fragility Curves from the Surveyed Data of the 2011 Great East Japan Tsunami in Sendai and Ishinomaki Plains. Coast. Eng. J. 2012, 54, 1250008-1–1250008-16. [Google Scholar] [CrossRef]
  2. Okumura, N.; Jonkman, S.N.; Esteban, M.; Hofland, B.; Shibayama, T. A Method for Tsunami Risk Assessment: A Case Study for Kamakura, Japan. Nat. Hazards 2017, 88, 1451–1472. [Google Scholar] [CrossRef] [Green Version]
  3. Dissanayaka, K.D.C.R.; Tanaka, N.; Vinodh, T.L.C. Integration of Eco-DRR and Hybrid Defense System on Mitigation of Natural Disasters (Tsunami and Coastal Flooding): A Review. Nat. Hazards 2021, 110, 1–28. [Google Scholar] [CrossRef]
  4. Fadly, U.; Murakami, K. Study on Reducing Tsunami Inundation Energy by the Modification of Topography Based on Local Wisdom. J. Jpn. Soc. Civ. Eng. Ser. B3 Ocean Eng. 2012, 68, I_66–I_71. [Google Scholar] [CrossRef] [Green Version]
  5. Yoshiya, I.; Norio, T. Evaluation of the Effect of Level 2 Tsunami on Reducing Energy Based on Changes in Flow Conditions When Overflowing a Double Levee Structure. J. Jpn. Soc. Civ. Eng. B1 Hydraul. Eng. 2017, 73, I_1009–I_1014. [Google Scholar] [CrossRef]
  6. Kimiwada, Y.; Tanaka, N.; Zaha, T. Differences in Effectiveness of a Hybrid Tsunami Defense System Comprising an Embankment, Moat, and Forest in Submerged, Emergent, or Combined Conditions. Ocean Eng. 2020, 208, 107457. [Google Scholar] [CrossRef]
  7. Zaha, T.; Tanaka, N.; Kimiwada, Y. Flume Experiments on Optimal Arrangement of Hybrid Defense System Comprising an Embankment, Moat, and Emergent Vegetation to Mitigate Inundating Tsunami Current. Ocean Eng. 2019, 173, 45–57. [Google Scholar] [CrossRef]
  8. Teh, S.Y.; Koh, H.L.; Liu, P.L.; Ismail, A.I.M.; Lee, H.L. Analytical and Numerical Simulation of Tsunami Mitigation by Mangroves in Penang, Malaysia. J. Asian Earth Sci. 2009, 36, 38–46. [Google Scholar] [CrossRef]
  9. Rahman, M.M.; Schaab, C.; Nakaza, E. Experimental and Numerical Modeling of Tsunami Mitigation by Canals. J. Waterw. Port Coast. Ocean Eng. 2017, 143, 04016012. [Google Scholar] [CrossRef]
  10. Ryo, M.; Gozo, T.; Fumihiko, Y. A Study on Scouring by Overflow from Seawall and Application of an Artificial Trench. J. Jpn. Soc. Civ. Eng. Ser B2 Coast. Eng. 2013, 69, I_791–I_795. [Google Scholar] [CrossRef] [Green Version]
  11. Tokida, K.; Tanimoto, R. Lessons for Countermeasures Using Earth Structures against Tsunami Obtained in the 2011 Off the Pacific Coast of Tohoku Earthquake. Soils Found. 2014, 54, 523–543. [Google Scholar] [CrossRef]
  12. Tsujimoto, G.; Mineura, R.; Yamada, F.; Kakinoki, T.; Uno, K. Scouring Mechanism Behind Seawall from Tsunami Overflow and Optimum Conditions to Reduce Tsunami Energy with an Artificial Trench. Coast. Eng. Proc. 2014, 1, 38. [Google Scholar] [CrossRef] [Green Version]
  13. Strusińska-Correia, A. Tsunami Mitigation in Japan after the 2011 Tōhoku Tsunami. Int. J. Disaster Risk Reduct. 2017, 22, 397–411. [Google Scholar] [CrossRef]
  14. Rashedunnabi, A.H.M.; Tanaka, N. Energy Reduction of a Tsunami Current through a Hybrid Defense System Comprising a Sea Embankment Followed by a Coastal Forest. Geosciences 2019, 9, 247. [Google Scholar] [CrossRef] [Green Version]
  15. Dissanayaka, K.D.C.R.; Tanaka, N.; Hasan, M.K. Effect of Orientation and Vegetation over the Embankment Crest for Energy Reduction at Downstream. Geosciences 2022, 12, 354. [Google Scholar] [CrossRef]
  16. Tanaka, N. Effectiveness and Limitations of Coastal Forest in Large Tsunami: Conditions of Japanese Pine Trees on Coastal Sand Dunes in Tsunami Caused by Great East Japan Earthquake. J. Jpn. Soc. Civ. Eng. Ser. B1 Hydraul. Eng. 2012, 68, II_7–II_15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Tanaka, N. Vegetation Bioshields for Tsunami Mitigation: Review of Effectiveness, Limitations, Construction, and Sustainable Management. Landsc. Ecol. Eng. 2009, 5, 71–79. [Google Scholar] [CrossRef] [Green Version]
  18. Tanaka, D.S.; Istiyanto, D.D.C.; Kuribayashi, M.D. Planning and Design of Tsunami-Mitigative Coastal Vegetation Belts; ICHARM Publ. No18, Technical note of PWRI no. 4177; International Centre for Water Hazard and Risk Management, United Nations Educational, Scientific and Cultural Organization: Paris, France, 2010.
  19. Cuomo, G.; Shams, G.; Jonkman, S.; van Gelder, P. Hydrodynamic Loadings of Buildings in Floods. In Coastal Engineering 2008; World Scientific Publishing Company: Singapore, 2009; pp. 3744–3756. ISBN 978-981-4277-36-5. [Google Scholar]
  20. Cox, D.; Tomita, T.; Lynett, P.; Holman, R. Tsunami Inundation with Macro-Roughness in the Constructed Environment. In Coastal Engineering 2008; World Scientific Publishing Company: Singapore, 2009; pp. 1421–1432. ISBN 978-981-4277-36-5. [Google Scholar]
  21. Goseberg, N. Reduction of Maximum Tsunami Run-up due to the Interaction with Beachfront Development—Application of Single Sinusoidal Waves. Nat. Hazards Earth Syst. Sci. 2013, 13, 2991–3010. [Google Scholar] [CrossRef] [Green Version]
  22. Park, H.; Cox, D.T.; Lynett, P.J.; Wiebe, D.M.; Shin, S. Tsunami Inundation Modeling in Constructed Environments: A Physical and Numerical Comparison of Free-Surface Elevation, Velocity, and Momentum Flux. Coast. Eng. 2013, 79, 9–21. [Google Scholar] [CrossRef]
  23. Manawasekara, C.; Mizutani, N.; Aoki, S. Influence of Openings and Orientation on Tsunami Generated Forces on Buildings. J. Disaster Res. 2016, 11, 670–679. [Google Scholar] [CrossRef]
  24. Chock, G.; Robertson, I.; Kriebel, D.; Francis, M.; Nistor, I. Tohoku, Japan, Earthquake and Tsunami of 2011; American Society of Civil Engineers: Reston, VA, USA, 2013; ISBN 978-0-7844-1249-7. [Google Scholar]
  25. American Society of Civil Engineers. Minimum Design Loads for Buildings and Other Structures; American Society of Civil Engineers: Reston, VA, USA, 2013; ISBN 978-0-7844-1291-6. [Google Scholar]
  26. Lukkunaprasit, P.; Ruangrassamee, A. Tsunami Loading on Buildings with Openings. Sci. Tsunami Hazards 2009, 28, 303–310. [Google Scholar]
  27. Mikami, T.; Shibayama, T.; Esteban, M.; Takabatake, T.; Nakamura, R.; Nishida, Y.; Achiari, H.; Rusli; Marzuki, A.G.; Marzuki, M.F.H.; et al. Field Survey of the 2018 Sulawesi Tsunami: Inundation and Run-up Heights and Damage to Coastal Communities. Pure Appl. Geophys. 2019, 176, 3291–3304. [Google Scholar] [CrossRef]
  28. Honda, T.; Oda, Y.; Ito, K.; Watanabe, M.; Takabatake, T. An Experimental Study on the Tsunami Pressure Acting on the Piloti-Type Buildings. Coast. Eng. Proc. 2014, 34, 41. [Google Scholar] [CrossRef]
  29. Yamamoto, T.; Kazama, S.; Touge, Y.; Yanagihara, H.; Tada, T.; Yamashita, T.; Takizawa, H. Evaluation of Flood Damage Reduction throughout Japan from Adaptation Measures Taken under a Range of Emissions Mitigation Scenarios. Clim. Chang. 2021, 165, 60. [Google Scholar] [CrossRef]
  30. Dang-Vu, H.; Shin, J.; Lee, K. Seismic Fragility Assessment of Columns in a Piloti-Type Building Retrofitted with Additional Shear Walls. Sustainability 2020, 12, 6530. [Google Scholar] [CrossRef]
  31. Benjamin, T.B. On the Flow in Channels When Rigid Obstacles Are Placed in the Stream. J. Fluid Mech. 1956, 1, 227–248. [Google Scholar] [CrossRef]
  32. Mignot, E.; Riviere, N. Bow-Wave-like Hydraulic Jump and Horseshoe Vortex around an Obstacle in a Supercritical Open Channel Flow. Phys. Fluids 2010, 22, 117105. [Google Scholar] [CrossRef] [Green Version]
  33. Mignot, E.; Moyne, T.; Doppler, D.; Riviere, N. Clear-Water Scouring Process in a Flow in Supercritical Regime. J. Hydraul. Eng. 2016, 142, 04015063. [Google Scholar] [CrossRef] [Green Version]
  34. Raju, K.G.R.; Rana, O.P.S.; Asawa, G.L.; Pillai, A.S.N. Rational Assessment of Blockage Effect in Channel Flow Past Smooth Circular Cylinders. J. Hydraul. Res. 1983, 21, 289–302. [Google Scholar] [CrossRef]
  35. Riviere, N.; Vouaillat, G.; Launay, G.; Mignot, E. Emerging Obstacles in Supercritical Open-Channel Flows: Detached Hydraulic Jump versus Wall-Jet-like Bow Wave. J. Hydraul. Eng. 2017, 143, 04017011. [Google Scholar] [CrossRef] [Green Version]
  36. Samarasekara, R.S.M.; Sasaki, J.; Esteban, M.; Matsuda, H. Assessment of the Co-Benefits of Structures in Coastal Areas for Tsunami Mitigation and Improving Community Resilience in Sri Lanka. Int. J. Disaster Risk Reduct. 2017, 23, 80–92. [Google Scholar] [CrossRef]
  37. Potts, D.A. Hydrodynamics of Vertical Surface-Piercing Cylinders. Master’s Thesis, University of Tasmania, Tasmania, Australia, 2019. [Google Scholar]
  38. Lam, K.; Zou, L. Experimental Study and Large Eddy Simulation for the Turbulent Flow around Four Cylinders in an In-Line Square Configuration. Int. J. Heat Fluid Flow 2009, 30, 276–285. [Google Scholar] [CrossRef]
  39. Wang, X.K.; Gong, K.; Liu, H.; Zhang, J.-X.; Tan, S.K. Flow around Four Cylinders Arranged in a Square Configuration. J. Fluids Struct. 2013, 43, 179–199. [Google Scholar] [CrossRef]
  40. Sumner, D.; Wong, S.S.T.; Price, S.J.; Païdoussis, M.P. Fluid Behaviour of Side-by-Side Circular Cylinders in Steady Cross-Flow. J. Fluids Struct. 1999, 13, 309–338. [Google Scholar] [CrossRef]
  41. Waniewski, T.A.; Brennen, C.E.; Raichlen, F. Bow Wave Dynamics. J. Ship Res. 2002, 46, 1–15. [Google Scholar] [CrossRef]
  42. McDonald, A.K. Experimental Investigation of Small-Scale Breaking Waves: Flow Visualization across the Air-Water Interface; Massachusetts Institute of Technology: Cambridge, MA, USA, 2005. [Google Scholar]
  43. Elhanafi, A.; Duffy, J.; Macfarlane, G.; Binns, J.; Keough, S.J. Inline Forces and Bow Wave Height on a Vertical Cylinder Moving in Waves—Experimental Study and CFD Validation. J. Fluids Struct. 2021, 107, 103387. [Google Scholar] [CrossRef]
  44. Yarnell, D.L. Bridge Piers as Channel Obstructions; USDA: Washington, DC, USA, 1934; Volume 52, p. 8.
  45. Synolakis, C.E.; Kong, L. Runup Measurements of the December 2004 Indian Ocean Tsunami. Earthq. Spectra 2006, 22, 67–91. [Google Scholar] [CrossRef]
  46. Goff, J.; Liu, P.L.-F.; Higman, B.; Morton, R.; Jaffe, B.E.; Fernando, H.; Lynett, P.; Fritz, H.; Synolakis, C.; Fernando, S. Sri Lanka Field Survey after the December 2004 Indian Ocean Tsunami. Earthq. Spectra 2006, 22, 155–172. [Google Scholar] [CrossRef]
  47. Yeh, H.; Sato, S. Tsunami Effects on Buildings and Coastal Structures. J. Disaster Res. 2016, 11, 662–669. [Google Scholar] [CrossRef]
  48. Tappin, D.R.; Evans, H.M.; Jordan, C.J.; Richmond, B.; Sugawara, D.; Goto, K. Coastal Changes in the Sendai Area from the Impact of the 2011 Tōhoku-Oki Tsunami: Interpretations of Time Series Satellite Images, Helicopter-Borne Video Footage and Field Observations. Sediment. Geol. 2012, 282, 151–174. [Google Scholar] [CrossRef]
  49. Peakall, J.; Warburton, J. Surface Tension in Small Hydraulic River Models—The Significance of the Weber Number. J. Hydrol. N. Z. 1996, 35, 199–212. [Google Scholar]
  50. Arnaud, G.; Rey, V.; Touboul, J.; Sous, D. Wave Propagation Through Dense Vertical Cylinder Arrays: 3D Experimental Study. Int. J. Ocean Coast. Eng. 2019, 2, 1950001. [Google Scholar] [CrossRef] [Green Version]
  51. Zdravkovich, M.M. Conceptual Overview of Laminar and Turbulent Flows Past Smooth and Rough Circular Cylinders. J. Wind Eng. Ind. Aerodyn. 1990, 33, 53–62. [Google Scholar] [CrossRef]
  52. Graf, W.H.; Altinakar, M.S. A Review of “Hydraulique Fluviale Ecoulement Permanent Uniforme et Non Uniforme”. Eur. J. Eng. Educ. 1994, 19, 378–379. [Google Scholar] [CrossRef]
  53. Wüthrich, D.; Pfister, M.; Nistor, I.; Schleiss, A.J. Experimental Study on the Hydrodynamic Impact of Tsunami-like Waves against Impervious Free-Standing Buildings. Coast. Eng. J. 2018, 60, 180–199. [Google Scholar] [CrossRef]
  54. Nandasena, N.A.K.; Sasaki, Y.; Tanaka, N. Modeling Field Observations of the 2011 Great East Japan Tsunami: Efficacy of Artificial and Natural Structures on Tsunami Mitigation. Coast. Eng. 2012, 67, 1–13. [Google Scholar] [CrossRef]
  55. Qi, Z.X.; Eames, I.; Johnson, E.R. Force Acting on a Square Cylinder Fixed in a Free-Surface Channel Flow. J. Fluid Mech. 2014, 756, 716–727. [Google Scholar] [CrossRef] [Green Version]
  56. Dias, P.; Dissanayake, R.; Chandratilake, R. Lessons Learned from Tsunami Damage in Sri Lanka. Proc. Inst. Civ. Eng. Civ. Eng. 2006, 159, 74–81. [Google Scholar] [CrossRef]
Figure 1. Schematic of the experimental setup and side-by-side cylinder arrangement: (a) plan view and side view of the experimental setup inside the channel, (b) two-cylinder array (case C2 (C2-1 or C2-2)) set up inside the channel, (c) four-cylinder array (case C4 (C4-1 or C4-2)) set up inside the channel.
Figure 1. Schematic of the experimental setup and side-by-side cylinder arrangement: (a) plan view and side view of the experimental setup inside the channel, (b) two-cylinder array (case C2 (C2-1 or C2-2)) set up inside the channel, (c) four-cylinder array (case C4 (C4-1 or C4-2)) set up inside the channel.
Fluids 08 00124 g001aFluids 08 00124 g001b
Figure 2. Flow forms with non-dimensional water depth relative to the cylinder geometry and the formation of the surface bow-wave around an impervious emergent cylinder in the subcritical open-channel flow: (a) flow forms for case C1 and surface bow-wave formation (b) for case C1-RC, (c) for case C2-1, and (d) for case C4-1. The numbers after D  and  W indicate the diameter and width, respectively. The subscript “sub” means “subcritical”.
Figure 2. Flow forms with non-dimensional water depth relative to the cylinder geometry and the formation of the surface bow-wave around an impervious emergent cylinder in the subcritical open-channel flow: (a) flow forms for case C1 and surface bow-wave formation (b) for case C1-RC, (c) for case C2-1, and (d) for case C4-1. The numbers after D  and  W indicate the diameter and width, respectively. The subscript “sub” means “subcritical”.
Fluids 08 00124 g002aFluids 08 00124 g002b
Figure 3. Flow forms of the non-dimensional water depth relative to the geometry for case C1 with the supercritical flow with a sloping bed. The subscript “sup-1” indicates supercritical flow over the 1/200 sloping bed.
Figure 3. Flow forms of the non-dimensional water depth relative to the geometry for case C1 with the supercritical flow with a sloping bed. The subscript “sup-1” indicates supercritical flow over the 1/200 sloping bed.
Fluids 08 00124 g003
Figure 4. Side view of the supercritical flow passing through emergent cylinders kept inside the experimental flume with a 1/200 bed slope: (a) case C1-CC, (b) case C1-RC, (c) case C2–1, (d) case C4-1.
Figure 4. Side view of the supercritical flow passing through emergent cylinders kept inside the experimental flume with a 1/200 bed slope: (a) case C1-CC, (b) case C1-RC, (c) case C2–1, (d) case C4-1.
Fluids 08 00124 g004
Figure 5. Flow forms of the non-dimensional water depth relative to the geometry for case C1 for the supercritical flow. The subscript “sup-2” indicates the supercritical flow over the horizontal channel bed.
Figure 5. Flow forms of the non-dimensional water depth relative to the geometry for case C1 for the supercritical flow. The subscript “sup-2” indicates the supercritical flow over the horizontal channel bed.
Fluids 08 00124 g005
Figure 6. Formation of the wall-jet-like bow-wave for the supercritical flow on the horizontal channel bed: (a) wall-jet in case C1, (b) front view of case C2-1, (c) front view of case C4-1, (d) back view of case C4-1.
Figure 6. Formation of the wall-jet-like bow-wave for the supercritical flow on the horizontal channel bed: (a) wall-jet in case C1, (b) front view of case C2-1, (c) front view of case C4-1, (d) back view of case C4-1.
Fluids 08 00124 g006aFluids 08 00124 g006b
Figure 7. Non-dimensional water depth relative to the upstream water depth ( h o ) along the centerline in case C4 for the supercritical flow over the horizontal channel bed: (a) at point C 1 , (b) point C o , and (c) point C 2 .
Figure 7. Non-dimensional water depth relative to the upstream water depth ( h o ) along the centerline in case C4 for the supercritical flow over the horizontal channel bed: (a) at point C 1 , (b) point C o , and (c) point C 2 .
Fluids 08 00124 g007
Figure 8. (a) Relative jet height in relation to time h j e t *   ( = h j e t / U 2 / 2 g recorded as a function of the upstream Froude number for case C1, (b) correlation plot for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (4) for the circular cylinders, (c) correlation plot for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (4) for the rectangular cylinders, (d) correlation plot for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (5) “Adapted with permission from Ref. [35]. 2017, Riviere et al.”.
Figure 8. (a) Relative jet height in relation to time h j e t *   ( = h j e t / U 2 / 2 g recorded as a function of the upstream Froude number for case C1, (b) correlation plot for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (4) for the circular cylinders, (c) correlation plot for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (4) for the rectangular cylinders, (d) correlation plot for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (5) “Adapted with permission from Ref. [35]. 2017, Riviere et al.”.
Fluids 08 00124 g008
Figure 9. Correlation plots for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (6) (a) for the case C2 array and (b) for the case C4 array.
Figure 9. Correlation plots for the experimental dimensionless wall-jet height ( h j e t * ) and that predicted according to Equation (6) (a) for the case C2 array and (b) for the case C4 array.
Fluids 08 00124 g009
Table 1. Water depths, geometric characteristics, and Froude numbers considered in the experiment.
Table 1. Water depths, geometric characteristics, and Froude numbers considered in the experiment.
DescriptionValue Froude   Number   ( F r )
Diameter of the circular cylinder  D (mm)32, 47, 60
Width of the rectangular cylinder  W (mm)30, 20
Water depth  h o (mm) (subcritical)30.7, 37.0, 45.5, 49.1, 55.70.59, 0.63, 0.65, 0.67, 0.70
Water depth  h o (mm) (supercritical with a 1/200 bed slope)20.5, 26.7, 30.7, 38.1, 44.21.27, 1.37, 1.40, 1.41, 1.45
Water depth  h o (mm) (supercritical with zero bed slope)53.2, 45.5, 38.0, 30.7, 23.51.45, 1.52, 1.70, 1.74, 1.99
Table 2. Experimental cases considered.
Table 2. Experimental cases considered.
CaseCase Name
Case C1-CCSingle emergent circular cylinder
Case C1-RCSingle emergent rectangular cylinder
Case C2-1Side-by-side two-cylinder array (center-to-center spacing  S Y was 200 mm)
Case C2-2Side-by-side two-cylinder array (center-to-center spacing  S Y was 100 mm)
Case C4-1Side-by-side four-cylinder array (spacings of  S X   and   S Y were 100 mm and 100 mm)
Case C4-2Side-by-side four-cylinder array (spacings of  S X   and   S Y were 200 mm and 200 mm)
Table 3. Non-dimensional surface bow-wave length.
Table 3. Non-dimensional surface bow-wave length.
F r L B w *   ( = L B w / D   or   L B w / W )
32 mm dia.47 mm dia. 60 mm dia.30 × 40 mm20 × 20 mm
0.590.560.420.370.911.03
0.630.570.490.421.071.18
0.650.700.550.501.251.30
0.670.790.580.521.371.43
0.700.870.690.611.541.53
Table 4. Non-dimensional hydraulic jump detachment length ( L j u m p * ) and non-dimensional wall-jet bow-wave height ( h j e t * ) relative to the geometry for case C1 for the supercritical flow when the channel bed slope equaled 1/200.
Table 4. Non-dimensional hydraulic jump detachment length ( L j u m p * ) and non-dimensional wall-jet bow-wave height ( h j e t * ) relative to the geometry for case C1 for the supercritical flow when the channel bed slope equaled 1/200.
Non-Dimensional PropertyUpstream Water Depth  ( h o ) —1/200 Sloping Channel Bed (mm)Cylinder
Dimension
20.526.730.738.144.2
L j u m p / D 0.450.530.680.650.7432 mm dia.
0.310.370.450.520.5847 mm dia.
0.320.330.390.480.5960 mm dia.
L j u m p / W 0.690.810.931.081.1130 × 40 mm
0.830.921.071.201.3920 × 20 mm
h j e t / D 0.640.931.331.291.4732 mm dia.
0.390.620.730.910.9647 mm dia.
0.330.490.600.730.8160 mm dia.
h j e t / W 0.630.971.151.361.5730 × 40 mm
1.011.521.622.122.4520 × 20 mm
Table 5. Non-dimensional values for surface bow-wave length ( L B w * ) and wall-jet height ( h j e t * ) relative to the spacing of S Y .
Table 5. Non-dimensional values for surface bow-wave length ( L B w * ) and wall-jet height ( h j e t * ) relative to the spacing of S Y .
Non-Dimensional PropertyUpstream Water Depth  ( h o ) —1/200 Sloping Channel Bed (mm)Case S Y   ( mm )
20.526.730.738.144.2
L B w * 0.170.190.220.240.26C2-1100
h j e t * 0.210.330.370.480.58
L B w * 0.070.090.100.120.13C2-2200
h j e t * 0.130.170.190.260.32
L B w * 0.170.180.200.230.26C4-1100
h j e t * 0.190.330.370.500.61
L B w * 0.080.100.120.130.14C4-2200
h j e t * 0.110.180.210.270.33
Table 6. Non-dimensional wall-jet height h j e t *  ( = h j e t / h o ) relative to the upstream water depth ( h o ) for the supercritical flow over the horizontal channel bed for case C1.
Table 6. Non-dimensional wall-jet height h j e t *  ( = h j e t / h o ) relative to the upstream water depth ( h o ) for the supercritical flow over the horizontal channel bed for case C1.
Non-Dimensional PropertyUpstream Water Depth  ( h o ) —1/200 Sloping Channel Bed (mm)Cylinder
Dimension
23.530.738.045.553.2
h j e t *
( = h j e t / h o )
2.922.191.881.481.2932 mm dia.
3.232.392.161.721.4647 mm dia.
3.482.662.281.791.5060 mm dia.
2.942.161.911.471.2630 × 40 mm
2.431.821.611.241.0620 × 20 mm
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dissanayaka, K.D.C.R.; Tanaka, N. Comparison of the Flow around Circular and Rectangular Emergent Cylinders with Subcritical and Supercritical Conditions. Fluids 2023, 8, 124. https://doi.org/10.3390/fluids8040124

AMA Style

Dissanayaka KDCR, Tanaka N. Comparison of the Flow around Circular and Rectangular Emergent Cylinders with Subcritical and Supercritical Conditions. Fluids. 2023; 8(4):124. https://doi.org/10.3390/fluids8040124

Chicago/Turabian Style

Dissanayaka, Kannangara D. C. R., and Norio Tanaka. 2023. "Comparison of the Flow around Circular and Rectangular Emergent Cylinders with Subcritical and Supercritical Conditions" Fluids 8, no. 4: 124. https://doi.org/10.3390/fluids8040124

APA Style

Dissanayaka, K. D. C. R., & Tanaka, N. (2023). Comparison of the Flow around Circular and Rectangular Emergent Cylinders with Subcritical and Supercritical Conditions. Fluids, 8(4), 124. https://doi.org/10.3390/fluids8040124

Article Metrics

Back to TopTop