 Open Access
 Total Downloads : 54
 Authors : Karim Badr Hussein , M.I. Ibrahim
 Paper ID : IJERTV8IS070222
 Volume & Issue : Volume 08, Issue 07 (July 2019)
 Published (First Online): 26072019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Wave Interaction with Vertical Slotted Breakwaters
Karim Badr Hussein1 , M.I. Ibrahim1
1 Lecturer of Irrigation and Hydraulics, Faculty of Engineering, AlAzhar University, Egypt
Abstract: This study aims to introduce two types of an innovate breakwater with an economic feasibility. The first type consists of two vertical perforated walls, the first wall is permeable in lower (porosity =50%) and the upper part is impermeable. And second wall is permeable in upper (porosity =50%) and the lower part is impermeable. between them horizontal slotted wall. the second type is the same construction as on the first type but without horizontal slotted wall. The results indicate that the hydrodynamic performance of the first type is better than that of the second type in the percentage (1015%) because of the presence of the horizontal slotted wall. The effect of wave force on the first model is bigger than the second model in the range (10 15%). The wave force on the proposed models increases with increasing the relative length (h/L).The transmission coefficient (kt) decreases with increasing the relative length (h/L).The reflection coefficient (kr) increases with increasing the relative length (h/L).from FLOW3D results able to find values velocity front and behind barriers.
Key Words: Coastal structures – Permeable breakwater – perforated wall – numerical model – refraction – transmission – energy dissipation.
1. INTRODUCTION
Traditional breakwaters (i.e. rubblemound, vertical caissons and gravity wall) are widely used to provide a protected calm water area to accommodate vessels and allow loading and unloading processes. Such types possess a large width according to the water depth. Consequently, great amounts of construction material are required. Moreover; such breakwaters block the littoral drift leading to the occurrence of severe erosion or accretion. In addition, they dampen the water circulation leading to a deteriorated water quality and achieving unbalance to the ecosystem. Furthermore; traditional structures need skilled labor for their construction and certain foundation requirements. All the above leads to an uneconomic construction cost.
On the contrary, permeable breakwaters avoid the occurrence of the above sideeffects, at the same time they provide reasonable protection with economic construction cost. This research was thus initiated with the objective of proposing and investigating the hydrodynamic performance of an innovative economic breakwater, numerically.
Many journals, periodicals and researches in the field of breakwaters were assembled, reviewed and comprehended from which it was clear that many researchers were occupied with finding out innovative types of economic breakwaters. Among these researchers were the following:
Wiegel (1960) and Hayashi et al. (1966) investigated breakwaters in the form of a row of close piles. Herbich (1998) investigated double rows of close piles. Suh el at. (2006); K. Laju el at. (2007) stated that breakwaters in the form of thin, rigid, pilesupported vertical barriers or many rows of piles which is placed below the water surface would provide relatively greater protection. Rageh and Koraim (2010) examined the hydraulic performance of a vertical wall with horizontal slots. The upper part was impermeable but the lower part of model was horizontal slots. Ahmed et al. (2011) investigated the hydrodynamic characteristics of a vertical slotted wall breakwater. They further stated that for more protection and more dissipation of energy a pair of permeable barriers might be desired. Isaacson et al. (1999) examined a pair of thin vertical barriers placed below the water surface. Koraim et al. (2011) and (2014) investigated the hydrodynamic characteristics of double permeable breakwater under regular waves. The model consisted of double walls with horizontal slots. Suh et al. (1995) established an analytical model for predicting wave reflection from a perforatedwall caisson breakwater. They implemented the GalerkinEigen function method to predict the reflection coefficient of a perforated wall caisson mounted on a rubble mound foundation. Hsu and Wu (1999) developed a numerical solution based on the boundary element method and boundary value problem for linear and secondorder wave. Isaacson et al. (1998), Ahmed et al. (2011) and Suh el at. (2006) established the
Eigenfunction expansion method for linear waves, second order waves and nonlinear waves. Lin et al. (2007), Huang et al. (2003), Lara et al. (2006), (2008) and Karim et al. (2009) executed numerical solutions for wave interaction with structures. They provided interesting examples for wave interaction with porous and submerged structure are reported. Ahmed (2014) investigated regular wave interaction using a numerical model of (FLOW3D, VOF) with a single vertical perforated wall. From the reviewed literature, it was obvious that permeable breakwaters were meticulously investigated, even though; other types are required to be investigated. This research aims to propose and investigate the hydrodynamic performance of an innovative economic breakwater, numerically. The research also aims to compare the hydraulic performance for the first and the second model and the present study with the previous study.
2PROPOSING AN INNOVATIVE BREAKWATER
The two innovative economic breakwaters were proposed, as shown in figures (1) and (2). It is a permeable breakwater that
comprises a pair of identical vertical perforated walls. First wall is impermeable in the upper and second wall is impermeable in
The transmission coefficient kt was calculated directly from the wave transmitted profile by:
the lower. A porosity of a permeable part =50%. The first and second perforated walls are placed apart by a distance of 0.5 of the water depth. Model (1) is vertical perforated wall
kt Ht
Hi
where: Ht is the transmitted wave height.
(2)
breakwater with horizontal slotted wall and model (2) is
The energy dissipation coefficient kd is given
vertical perforated wall breakwater without horizontal slotted wall.
3INVESTIGATING THE PROPOSED BREAKWATER
kd 1 kr2 kt 2
3.4ANYLATICAL STUDY
(3)
NUMERICALLY
This section presents the implemented model and its theory. It
The velocity potential is assumed as periodic motion in time T and it can be expressed as follows:
presents the validation process of the model together with the executed numerical simulations to the proposed breakwater.
(x, z,t) Re

ighi
2
(x, z)
1
cosh kh
eit
(4)
3.1 THE IMPLEMENTED MODEL
The proposed breakwater was investigated. This was achieved via achieving numerical simulations using the commercial Computational Fluid Dynamics (CFD) code FLOW3D. This
Where Re is real part of a complex value, is wave angular frequency, g is gravitational acceleration, i 1 and
K is wave number (K=2/L).
is attributed to the fact that from the assembled literature, it was The dispersion relationship could be formulated as:
1
1
clear that CFD applications are common practice in all sectors of engineering and they are increasingly becoming important in maritime and coastal engineering. Therefore, the commercial CFD code (i.e. FLOW3D, Flow Science Inc.) was chosen to be implemented in this study.
gk tanh( Kh) The permeable boundary condition along
the structure has been developed on the basis of the formulation of Sollit and Cross (1972) and a adopted by Yu (1995) for a thin vertical barrier. This might be given by:
3.2 THEORY OF FLOW3D
Basically, FLOW3D applies the finite volume theory to solve
1 (x)
x
2 (x)
x
iG (
(x) 2
(x))
the threedimensional Reynolds Averaged Navier Stokes (RANS) equations. The model consists of a group of solid subcomponents, figures (3) and (4). The numerical model
at x a
for
h z D1
(5)
within FLOW3D represents the geometrical and hydraulic boundary conditions.
3.3NUMERICAL SIMULATIONS USING FLOW3D
(x)
2
2
x
(x)
3
3
x
iG
(2

(x))
3
3
Confident with the validation process, the model was implemented, varying the different parameters. Numerical
at x a for D2 z 0
G
(6)
replications were achieved to simulate the proposed breakwater.
The proportional constant G
, G is called the
b
In order to get a good compromise between precision/accuracy and computation time, two independent meshes with different cell sizes were used. Mesh cells are sized by 1 cm in each direction for waves of small frequencies and mesh cells are sized by 0.5 cm for waves of large frequencies. The time window for analyzing the wave height is carefully selected
permeability parameter and is expressed by: G
f is
where is the porosity of the structure given by the dimension
and spacing of the piles, f is the friction factor (empirical parameter) and s is the inertia coefficient and given
according to the wave length and adjusted to avoid any
by s 1 cm1
where cm is an added mass
reflection from the flume end or the wave paddle.
The reflection coefficient was calculated by the threeprobe method of Mansard and Funk (1980). The selected data are converted into frequency domain by Fast Fourier Transformation. Finally, the spectrum of the incident, transmitted and reflected wave height were calculated. Thereby, the reflection coefficient kr is calculated from extracted wave profiles by:
coefficient. Eigen function expansion solved the velocity potential in a series of infinite number of solutions as follow:
1 (x) i A1m cos[m (h z)]exp(m (x a))
m0
kr Hr (1)
Hi
Where: Hr is the reflected wave height, Hi is incident wave height.
at x a
(7)
2 (x)
m0
A2m
cos[m
(h z)]exp(m
(x a))
0
F * p(0, z)dz
d
(14)
A3m cos[m (h z)]exp( m (x a))
F * gH i
2k
(1 kr

kt
) tanh kh
(15)
m0
At and
a x a
(8)
F Re[F*eit ]
(16)
3 (x) A4m cos[m (h z)]exp(m (x a))
m0
At x a
(9)
Applying the matching conditions (i.e. combining (7), (8) and
(9) with (4) and (5) at the breakwater; the
coefficients A1m , A2m ,
A3m
and A4m could be determined
by the following matrix equation:
C( mn )
C( mn )
C( mn)
C( mn )
11
m0
12
m0
13
m0
14
m0
A
b
C( mn )
C( mn )
C( mn)
C( mn ) 1m 1n
21
22
23
24 A b
m0
m0
m0
m0 2 m
2 n
( mn )
( mn )
( mn)
( mn )
b
C31
C32
C33
C34
3m
3n
m0
m0
m0
m0 A b
A
A
C( mn )
C( mn )
C( mn)
C( mn ) 4 m 4 n
41
m0
42
m0
43
m0
44
m0
Figure (1) Definition sketch for a vertical perforated wall breakwater with horizontal slot
For n =1, 2, 3,.. (10)
Equation (10) could be solved by numerical tools. Consequently, kr and kt could be obtained;
kr =
kt =
A10
A40
The energy dissipation coefficient can be determined using equation (3).
3.5 HYDRODYNAMIC FORCE
The hydrodynamic pressure (p) exerted on the surfaces of a body can be expressed by linearizing Bernoullis equation as follows:
p
t
i (1
2 )
xa
x = a (11)
xa
p
t
i (2
3
) xa
x = a (12)
p gH i
2
xa
coshk(z h) (1 k
coshkh r

kt )
(13)
Figure (2) Definition sketch for a vertical perforated wall breakwater
Then, the hydrodynamic forces (F*) exerted on the breakwater can be evaluated by integrating the pressure around the bodys wetted surface as follows:
without horizontal slot
Figure (3) Breakwater model and Mesh blocks without horizontal slot
Figure (4) Breakwater model with horizontal slot
4. RESULTS AND DESICCATIONS
In this paper, the hydrodynamic performance of a new type breakwaters is studied using numerical models. The breakwater consists of a double vertical perforated walls (with and without horizontal slots).The first wall is impermeable in upper part and permeable in lower part. The second wall is permeable in upper part and impermeable in lower part. The transmission, the reflection, and the wave energy dissipation coefficients are presented for different waves and structure parameters. It Table
(1) represents the results obtained from the experiments that was carried out in the Irrigation and Drainage Engineering laboratory of the Faculty of Engineering, Zagazig University. The experimental work was carried out without proposed breakwater to determine incident wave height and wave periods.
Table (1) the wave periods, lengths and incident wave heights for different wave generator frequencies without models
Parameters
units
Ranges
Incident wave(Hi)
Cm
6.7
7.6
8.2
9
10.1
11.5
12.5
Wave period (T)
Sec.
1.4
1.35
1.3
1.2
1.1
1
0.9
Water depth(h)
M
0.4
0.4
0.4
0.4
0.4
0.4
0.4
Wave length(L)
M
2.42
2.30
2.17
1.93
1.68
1.44
1.21
h/L
D.L
0.17
0.17
0.18
0.21
0.24
0.28
0.33
Chamber width (2a)
0.5h
0.2
0.2
0.2
0.2
0.2
0.2
0.2
4.1 THE WAVE FLUME
The wave flume used in this work has a rectangular cross section with 2 m bed width and 1.2 m depth. The overall length of the flume is 12 m. All sides of the flume are made of reinforced concrete of 0.25 m thickess. This flume is divided into three parts (i.e the inlet, working section and the outlet parts). A gravel wave absorber with slope 3:1 is installed at outlet part of the flume to absorb the transmitted waves. A general view of the flume is photo (1).
Figure (5) General view of the wave flume.
4.2 ANALYSING THE RESULTS
Figure (6) shows the wave high resulted from the experimental test conducted on the laboratory without using breakwater at different frequencies. Figuer (7) clear shows that the transion coefficant (kt) . Decresses with the increase of the relative depth (h/L) and with the comparasion of the first model and the second model it is cleared that the energy dispation coefficant is better for the first model than the second in the range betwwen 1015%. Fig(8) Comparison of predicted(double perforated walls without horizontal slot) hydrodynamic coefficients (kt, kr and kl) with results from Ji and Suh (2010) as a function of d/L, when D/d = 0.5, B/d = 1,
= 0.5, f = 2. Fig. (9). Comparison of predicted hydrodynamic coefficients (kt, kr and kl) with results from Laju et al. (2011) as a function of d/L, when D/d = 0.35, B/d = 0.5, = 0.25 and f
= 1.2 Fig. (10) Comparison of predicted and (Flolw3D) hydrodynamic coefficients (kt and kr) with results from Laju et al. (2011) as a function of d/L, when D/d = 0.35, B/d = 0.5, =
0.25 and f = 1.2
From figure (11) the effect of the wave force on the first model is bigger than its effect on the second model. It also can be seen
that with the increase of the relative depth the variation between the two models is increased. Fig
(12) shows the model in CED (FLOW3D) and the location of measuring probes. Fig (13) shows the free surface elevation (cm)after 2.0 meters from breakwater by using (FLOW 3D) at wave period (T) = 0.9sec and 1.1 sec. Fig (14) shows surface elevation (cm) at wave period T =1.5 sec, wave translated at probe (1)and wave reflected at probe (2).
The numerical model results are used to detect the velocity field and the velocity vectors in the surrounding barrier.
Figure (15) shows the velocity vector and velocity field for one cycle of with time increment of 0.1 sec of wave frequency =1 Hz, calculated with FLOW3D
The higher velocities are observed at the crest of wave and around the slots.
The higher velocities are formed trough of slots as a result to effect of obstacle
The velocity magnitude is very high in front of the barrier and very low behind it
The flow in the area between the barriers is turbulence and the motion seems a vertically except the region that located near the slots.
The set of the experimental results used in order to examine the effect of the shape and chamber width in mean velocity fields at constant porosity = 50% and D/h=0.5 includes 2 tests. The first Test is a double perforated walls with horizontal slots for B/h=0.5 and B/h=1.0 the mean velocity field in the breakwater area is presented in figure (16) (a)t=10.25T, (b) t=10.50T and (c)t=10.75T over the 10th wave cycle. Maximum velocities 72.9 cm/s are observed in the region of the impermeable barriers (figure 16) and mainly in the seaward side of the structures where the main structureinteraction takes place and intense vortices are observed. A recirculating region is observed in the region beneath the body of each barrier with the first one being stronger and wider the mean velocity field for permeable barriers with porosity 0.50 is presented. observed in the impermeable barriers case under the same hydrodynamic conditions. The recirculating region this time is obvious only beneath the first barrier while the interaction of the water mass with the gaps on the body of the barriers is presented. In figure

Maximum velocities 56.7 cm/s observed velocity for B/h=1.0 bigger than velocity for B/h=0.5 at the same condition.
Figure (6) Waves profiles for wave periods: (a) T=1.40, (b) T=1.00 sec
Figure (7) Comparison between CFD (FLOW3D) and predicted results for a double perforated walls with horizontal wall slot and without horizontal wall slot as function of (h/L) at 2a = 0.5h, (a) kt, (b) kr and (c) kd
Figure (8) Comparison of predicted (double perforated walls without horizontal slot) hydrodynamic coefficients (kt, kr and kl) with results from Ji and Suh (2010) as a function of h/L, when D/d = 0.5, B/d = 1, = 0.5, f = 2
Fig.(9) Comparison of predicted (double perforated walls without horizontal slot) hydrodynamic coefficients (a)kt, (b)kr and (c)kl with results from Laju et al. (2011) as a function of h/L, when D/d = 0.35, B/d = 0.5, = 0.25 and f= 1.2
Fig. (12) wave direction, location wave reflection and wave translation
Fig.(10) Comparison of predicted and (Flolw3D) (double perforated walls without horizontal slot) hydrodynamic coefficients (a)kt and(b) kr with results from Laju et al. (2011) as a function of d/L, when D/d = 0.35, B/d = 0.5, = 0.25 and f = 1.2
Fig.(11) Comparison of dimensionless wave forces between a double perforated walls with horizontal slot and a double perforated walls without horizontal slot as function of (d/L = 0.5
Fig (13) Free surface elevation (cm) after 2.00 meters from breakwater by using (FLOW 3D).
Fig (14) surface elevation (cm) at wave period T =1.5 sec, wave translated at probe 1and wave reflected at probe 2
Figure (15) Velocity magnitude (cm/s) and velocity vector of FLOW3D
results for T=1.2 sec, Hi= 9 cm, dt = 0.1 sec 2a/h= 0.5 Figure (16) Detailed velocity field in the region of the barriers
(Porosity=0.50, B/h=0.5 and Hi=9 cmT=1.20 sec).
Figure (17) Detailed velocity field in the region of the barriers (Porosity=0.50, B/h=1.0 and Hi=9 cmT=1.20 sec).
5 CONCLUSIONS

The reviewed literature revealed that the breakwaters were physically modeled and investigated meticulously but breakwater numerical modeling has some discrepancies. It was also clear that extra investigations are needed. Among the reviewed available models, model Flow 3D was found to be capable of simulating the proposed breakwater.

Flow3D was validated against extensive laboratory investigations and theoretical model.

Flow 3D is capable of describing the wave interaction of a linear wave with double vertical perforated walls. Flow 3D is capable of reproducing most of the important features of the experimental data and semianalytical results. Flow 3D reproduced numerical results that are perfectly acceptable.

The wave force on the proposed models increases with increasing the relative length (h/L).

The effect of wave force on the first model is bigger than the second model in the range (1015%).

The transmission coefficient (kt) decreases with increasing the relative length (h/L).and the reflection coefficient (kr) increases with increasing the relative length (h/L).

The comparasion of the first and the second model it is cleared that the energy dispation coefficant is better for the first model than that for the seconde in the range betwwen 1015%.

The hydrodynamic performance of the seconed model is lower than that of the previous study in the rang of (37%).

FLOW3D results are able to find values velocity in front of and behind barriers.

Velocity magnitude for B/h=1.0 bigger than velocity for B/h=0.5 at the same condition.
NOTATIONS
The following symbols have been adopted for use in this paper:

A10 = complex reflection coefficient;

A40 = complex transmission coefficient;

A1n = complex unknown coefficients;

= half distance between the two walls;

b = thickness of the vertical wall;

C =added mass coefficient;

f = friction coefficient;
G = permeability parameter;

g = acceleration of gravity;

hi = incident wave height;

hr = reflected wave height;

Ht = transmitted wave height;

h = water depth;

k = incident wave number;

kl = energy dissipation coefficient;

kr = reflection coefficient;

kt = transmission coefficient;

L = wave length;

T = wave period;

t = time;


x , z = two dimensional axis;

1 = porosity of the lower part of the first wall;

p = total flow velocity potential;

1 = seaward velocity potential;

2 = velocity potential between the two walls;

3 = shoreward velocity potential;
= angular wave frequency.

F* = wave force

REFERENCES

Ahmed, H., 2011. Wave Interaction with Vertical Perforated wallsas a Permeable Breakwater, PhD. Thesis, Hydro Sciences (IGAW), Bergische University of Wuppertal, Germany, 2011.

Hayashi, T., & Kano, T., 1966. Hydraulic research on the closely space Pile breakwater. 10th Coastal Eng. Conf., ASCE, New York, Vol. 11, Chapter 50.

Herbich, J. B., 1989. Wave transmission through a doublerow Pile breakwater. Proc. 21st Int. Conf. on Coastal Eng., ASCE, Chapter 165, Torremolinos, Spain.

Hirt, C. W. and Nichols, B. D., 1981. Volume of Fluid (VOF) method for the dynamics of free boundaries, J. Computat. Phys., vol. 39, no. 1, pp. 201225.

Hsu, HH. & Wu, YC., 1999. Numerical solution for the second order wave interaction with porous structures. International Journal for Numerical Methods in Fluids, Vol. 29 Issue 3, pp. 265 288.

Huang, C. J.; Chang, H. H.; and Hwung, H. H., 2003. Structural permeability effects on the interaction of a solitary wave and a submerged breakwater, Coastal Engineering. Vol. 49, pp. 124.

Isaacson, M.; Baldwin, J.; Premasiri, S. and Yang, G., 1999. Wave interaction with double Perforatedbarrier,Applied ocean research, Vol. 21, pp. 8191.

Isaacson, M.; Premasiri, S. and Yang, G., 1998. Wave interaction with vertical Perforatedbarrier, J. of Waterway, Port, Coastal and Ocean Engineering, vol. 124, no. 3, pp. 118125.

Karim, M. F.; Tanimoto, K. and Hieu, P. D., 2009 Modelling and simulation of wave transformation in porous structures using VOF based twophase flow model, Applied Mathematical Modellin. vo 33, pp. 343360.

Koraim, A. S.; Heikal; E. M. and Rageh,O. S., 2011. Hydrodynamic characteristics of double permeable reakwater undregular waves, Coastal Eng., vol. 24, pp. pp.503527.

Koraim, A. S; Iskander, M. M. and Elsayed, W. R., 2014. Hydrodynamic performance odouble rows of piles suspending horizontal c shaped bars, Coastal Eng., vol. 84, pp. 8196.

Laju, K., Sundar, V. & Sundaravadivelu, R., 2007. Studies on Pile supported double skirt breakwater models. Journal of Ocean Technology, Vol. 2, No.1. pp. 3253.

Lara, J.; Garcia, LN. and Losada, I.J., 2006. RANS modelling applied to random wave interaction with submerged permeable structures, Coastal Engineering, vol. 53, pp. 395417.

Lara, J.L.; Losada, I.J. and .Guanche, R., 2008. Wave interaction with lowmound breakwaters using a RANS model, Ocean Engineering, vol. 35, pp. 1388 1400.

Lin, P.; and Karunarathna, .S.A., 2007. Numerical stud of solitary wave interaction with porous breakwaters, J. of waterway, port, coastal and ocean engineering. , pp. 352363.

Mansard, E .P. D. & Funke, E. R., 1980. The measurement of incident and reflected spectra using a least squares method. In Proc. 17th Coastal Eng. Conf., Sydney, Australia, pp 159174.

Rageh, O.S.; Koraim, A.S., 2010. Hydraulic performance of vertical walls with horizontal sloed as breakwater, Coastal Engineering , vol. 57, pp. 745756, 2010.

Sollit, C. K. & Cross, R.H., 1972. Wave transmission through permeable breakwaters. Proceedings of the 13th Coastal Eng. Conf., ASCE, Vancouver, pp. 18271846.

Suh, K. D., Shin, S. & Cox, D. T., 2006. Hydrodynamic characteristics of PileSupported vertical wall breakwaters. J. of Waterways, Port, Coastal and Ocean Engineering, Vol.132, No.2, pp.8396.

Suh, K.D., Park, W.S., 1995. Wave reflection from perforatedwall caisson breakwaters. Coastal Engineering 26, 177193.Suh, KD., Park, J.K., Park, W.S., 2006. Wave reflection from partially perforatedwall caisson breakwater. Coastal Engineering 33, 264 280.

Wiegel, R. L., 1960. Transmission of wave past a rigid vertical thin barrier J. Waterway, Port, Coastal and Ocean Eng., ASCE, Vol. 86, No.1.

Yu, X., 1995. Diffraction of water waves by porous breakwater.,
J. of Waterway, Port, Coastal and Ocean Engineering, Vol. 121, No. 6, pp. 275282.