Experimental Investigation on Changes of Water Surface Profile with Gaussian Shaped Bottom and Side Roughness

Similar documents
OPEN CHANNEL FLOW WORKSHEET 3 WATER SURFACE PROFILES

Exercise (3): Open Channel Flow Rapidly Varied Flow

ANSWERS TO QUESTIONS IN THE NOTES AUTUMN 2018

Experiment (13): Flow channel

Effect of Fluid Density and Temperature on Discharge Coefficient of Ogee Spillways Using Physical Models

HYDRAULIC JUMP AND WEIR FLOW

EXAMPLES (OPEN-CHANNEL FLOW) AUTUMN 2018

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 3, No 1, 2012

3. GRADUALLY-VARIED FLOW (GVF) AUTUMN 2018

LABORATORY EXPERIMENTS ON WAVE OVERTOPPING OVER SMOOTH AND STEPPED GENTLE SLOPE SEAWALLS

Exercise (4): Open Channel Flow - Gradually Varied Flow

International Journal of Scientific & Engineering Research, Volume 5, Issue 1, January ISSN

Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati

Effect of channel slope on flow characteristics of undular hydraulic jumps

Wind Directional Effect on a Single Storey House Using Educational Wind Tunnel

Scour Investigations Behind a Vertical Sluice Gate without Apron.

Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati

HY-8 Version 7.2 Build Date January 17, Federal Highway Administration.

FREE OVERFALL IN A HORIZONTAL SMOOTH RECTANGULAR CHANNEL

Evaluation of step s slope on energy dissipation in stepped spillway

Numerical investigation of transition between free surface flow and pressurized flow for a circular pipe flowing full upstream

Determination of the wind pressure distribution on the facade of the triangularly shaped high-rise building structure

Transactions on Ecology and the Environment vol 12, 1996 WIT Press, ISSN

Discharge Coefficient in Oblique Side Weirs

Experimental Investigation of Clear-Water Local Scour at Pile Groups

Wind Pressure Distribution on Rectangular Plan Buildings with Multiple Domes

Plane Turbulent Wall Jets in Limited Tailwater Depth

Transitional Steps Zone in Steeply Sloping Stepped Spillways

Quantification of the Effects of Turbulence in Wind on the Flutter Stability of Suspension Bridges

Discharge Coefficient in Siphon Spillway with Different Cross Sections

MODELING OF CLIMATE CHANGE IMPACTS ON COASTAL STRUCTURES - CONTRIBUTION TO THEIR RE-DESIGN

A Study on the Distribution of the Peak Wind Pressure Coefficient for the Wind Resistant Design of Rooftop Hoardings in High-rise Buildings

A Comparative Study of Self-aerated Stepped Spillway and Smooth Invert Chute Flow: The effect of Step-induced Macro-roughness

Effects of wind incidence angle on wind pressure distribution on square plan tall buildings

Plan B Dam Breach Assessment

PhD student, January 2010-December 2013

Plane Turbulent Wall Jets in Limited Tailwater Depth

Broadly speaking, there are four different types of structures, each with its own particular function:

Wave Induced Flow around Submerged Sloping Plates

Wind Flow Model of Area Surrounding the Case Western Reserve University Wind Turbine

Rock Ramp Design Guidelines. David Mooney MS Chris Holmquist-Johnson MS Drew Baird Ph.D. P.E. Kent Collins P.E.

1. In most economical rectangular section of a channel, depth is kept equal to

Hydraulic Engineering

Stability of Concrete Macro-Roughness Linings for Overflow Protection of Earth Embankment Dams ( 1 ) - Discussion

Impact of anti-vortex blade position on discharge experimental study on the coefficient of morning glory spillway

Nearshore Placed Mound Physical Model Experiment

COST EFFECTIVE STORAGE CAPACITY INCREASE FOR ALUMINA TAILINGS DISPOSAL AREA THROUGH SPILLWAY OPTIMISATION

19.1 Problem: Maximum Discharge

IMPACT OF MAKING THE CREST OF WEIR MULTIFACETED ON DISCHARGE COEFFICIENT OF MORNING GLORY SPILLWAY

LARGE-SCALE MODEL TESTS ON SCOUR AROUND SLENDER MONOPILE UNDER LIVE-BED CONDITIONS

Evaluating Surge Potential in CSO Tunnels

DESIGN OF BELL-MOUTH SPILLWAY AT BARVI DAM

Morphological modelling of Keta Lagoon case

2. RAPIDLY-VARIED FLOW (RVF) AUTUMN 2018

EXPERIMENTAL STUDY ON THE DISCHARGE CHARACTERISTICS OF SLUICE FOR TIDAL POWER PLANT

Greenup Lock Filling and Emptying System Study

Pressure coefficient on flat roofs of rectangular buildings

Investigation of Suction Process of Scroll Compressors

ZIN Technologies PHi Engineering Support. PHi-RPT CFD Analysis of Large Bubble Mixing. June 26, 2006

International Journal of Technical Research and Applications e-issn: , Volume 4, Issue 3 (May-June, 2016), PP.

3D NUMERICAL MODELLING OF THE CAPACITY FOR A PARTIALLY PRESSURIZED SPILLWAY

Modelling of Pressurised Pipes within InfoWorks ICM and CS

Simulation of hydraulic regime and sediment transport in the Mekong delta coast

TRANSPORT OF NEARSHORE DREDGE MATERIAL BERMS

SUBMERGED VENTURI FLUME. Tom Gill 1 Robert Einhellig 2 ABSTRACT

THEORETICAL EVALUATION OF FLOW THROUGH CENTRIFUGAL COMPRESSOR STAGE

EXPERIMENTAL STUDY OF WIND PRESSURES ON IRREGULAR- PLAN SHAPE BUILDINGS

Critical Gust Pressures on Tall Building Frames-Review of Codal Provisions

Experimental Study of Water Hammer Pressure in a Commercial Pipe

Flow Characteristics and Energy Dissipation Over Traditional and Stepped Spillway with Semicircular Crest

IMAGE-BASED STUDY OF BREAKING AND BROKEN WAVE CHARACTERISTICS IN FRONT OF THE SEAWALL

A Review of the Bed Roughness Variable in MIKE 21 FLOW MODEL FM, Hydrodynamic (HD) and Sediment Transport (ST) modules

Using sea bed roughness as a wave energy dissipater

WMS 8.4 Tutorial Hydraulics and Floodplain Modeling HY-8 Modeling Wizard Learn how to model a culvert using HY-8 and WMS

Computational Analysis of Cavity Effect over Aircraft Wing

Experimental Investigation Of Flow Past A Rough Surfaced Cylinder

Tidal River Management (TRM)-An Innovative Scientific Approach for Sustainable Sediment Management

Experimental Analysis of Brink Depth in a Smooth Rectangular Channel

Wind tunnel effects on wingtip vortices

DEAD-ENDCHANNELFLUSHINGINHARBORS

AIR-WATER FLOW STRUCTURES AT AN ABRUPT DROP WITH SUPERCRITICAL FLOW

FL M FLUME FOR HYDRAULIC BENCH

CFD ANALYSIS AND COMPARISON USING ANSYS AND STAR-CCM+ OF MODEL AEROFOIL SELIG 1223

INTERACTION BETWEEN WIND-DRIVEN AND BUOYANCY-DRIVEN NATURAL VENTILATION Bo Wang, Foster and Partners, London, UK

Comparison of stepped and smooth spillway effects on stream reaeration

3. Observed initial growth of short waves from radar measurements in tanks (Larson and Wright, 1975). The dependence of the exponential amplification

The Hydraulic Design of an Arced Labyrinth Weir at Isabella Dam

An Experimental Investigation of Head loss through Trash Racks in Conveyance Systems

Undertow - Zonation of Flow in Broken Wave Bores

MODELLING ANCILLARIES: WEIR COEFFICIENTS

ICE INDENTATION TEST ON CONSOLIDATED LAYER MODEL MADE OF SALINE AND FRESHWATER ICE BLOCKS

Wave Breaking and Wave Setup of Artificial Reef with Inclined Crown Keisuke Murakami 1 and Daisuke Maki 2

AN EXPERIMENTAL INVESTIGATION OF SPILLING BREAKERS

LABORATORY STUDY ON TSUNAMI REDUCTION EFFECT OF TEIZAN CANAL

Shoreline Evolution Due to Oblique Waves in Presence of Submerged Breakwaters. Nima Zakeri (Corresponding Author), Mojtaba Tajziehchi

SEDIMENT BUDGET OF LIDO OF PELLESTRINA (VENICE) Written by Marcello Di Risio Under the supervision of Giorgio Bellotti and Leopoldo Franco

EXPERIMENTAL STUDY ON BEESWAX USING WATERJET DRILLING

Study on wind turbine arrangement for offshore wind farms

Numerical Simulation And Aerodynamic Performance Comparison Between Seagull Aerofoil and NACA 4412 Aerofoil under Low-Reynolds 1

Cross-shore sediment transports on a cut profile for large scale land reclamations

Transcription:

Experimental Investigation on Changes of Water Surface Profile with Gaussian Shaped Bottom and Side Md. Rafiue Islam a, Shariful Islam b*, Md. Abdul Qaiyum Talukder c, S. M. Rezwan Hossain d Abstract Bed roughness is a main factor on which the water surface profile depends. In this study the water surface profile over Gaussian shaped rough surface is experimentally investigated. The Gaussian shaped bed roughnesses were made of wood of different heights (2cm, 4cm and 6cm) with same base width of 12cm. Besides, these experiments were performed for same slope and varying Fraud number in the laboratory tilting flume. The multipurpose tilting flume is an artificial channel having a channel width of 8cm. Its maximum upstream water depth is 2cm. Discharge throughout the channel can be changed by rotating the key-wheel attached to the pump. Discharge throughout the channel varies from (.1-.43) m3/s. Its slope can also be changed by rotating the wheel at upstream of the channel. The pattern of water surface fluctuation over single Gaussian shaped roughness was found different than that of double Gaussian shaped roughness. For the same slope and same roughness, if the Fraud number is increased the maximum fluctuation was found to be increased. It is always found that the water depths at the upstream of roughness were constant which are gradually decreasing from the middle of roughness and at the downstream of roughness the water depth is the lowest. Index Terms Bed roughness, Water surface profile, Fraud number, Gaussian shaped roughness, Water depth fluctuation. formation of sand waves in river and seas are still one of the most interesting subjects. The water surface profile mainly depends on the bed surface and bed roughness of the channel. Nabi (28) reported that the water levels during river floods, and hence the risk of flooding, depend on the hydraulic roughness of the river. Based on the flow parameter, different types of bed forms are observed in the river bed. One of the examples of this type of periodic irregularities on the river bed is dunes. The development of dunes and the associated hydraulic roughness during a flood is complex. Initially dunes grow higher and make the river bed rougher, but in later stages the dunes grow longer with the opposite effect of making the river bed smoother. Fig. 1 shows the water surface profile over wavy bed surface. Nabi, A. M. (28) studied 3D Model of Detailed Hydrodynamics with Sediment Transport for Simulation of Subaqueous Dunes. The model aims at giving better insight into the development of the hydraulic roughness, and hence the flooding risk, during river floods. The good understanding thus obtained will allow the development of parameterized models for large spatial and temporal scales that can be used in operational model for flood early warning system the determination of design water levels. I. INTRODUCTION The interruption problem between the boundary and the water surface profile of open channel flow such as the Md. Rafiue Islam a, Department of Civil Engineering, Khulna University of Engineering & Technology, Khulna, Bangladesh, +881673922867, (E-mail: rafiue.anik3@gmail.com). Shariful Islam b*, Department of Civil Engineering, Khulna University of Engineering & Technology, Khulna, Bangladesh, +881739294, (E-mail: bappytalukder@gmail.com). Md. Abdul Qaiyum Talukder c, Department of Civil Engineering, Khulna University of Engineering & Technology, Khulna, Bangladesh, +8817672, (E-mail: qaiyumkuet@gmail.com). S. M. Rezwan Hossain d, Department of Civil Engineering, Khulna University of Engineering & Technology, Khulna, Bangladesh, +88167677226, (E-mail: smrh.rezwan@gmail.com). Fig. 1 Wavy bed surface Erpicum S. (26) studied the comparison of 2D turbulence models for steady flow computation in a macro-rough channel. In this study, experimental researches have been carried out at the laboratory of hydraulic construction (LCH-EPL) to assess the head losses generated by a succession of regular cavities at both bank of a straight channel. In this scope, numerical modeling of the experimental tests has been realized at the hydrology, applied hydrodynamics and hydraulic constructions research unit (HACH-ULg) with a 2D finite volume multiblock free surface flow solver. 1873

Ahmed, T. (211) studied the experimental Investigation on Changes of Water Surface Profile with Bottom for different slope, velocity and discharges in the laboratory tilting flume. In this study he found, I. For single roughness, the water surface was found to be abruptly falling down near the peak of the roughness and the water depth was found to be lowest just before the second trough. II. For double triangular roughness, the water depth is the lowest just before the crest of the triangular roughness and the water surface is peak at the end of the roughness. The pattern of water surface fluctuation over single triangular roughness was found different than that of double triangular roughness. III. If the fraud number is increase the maximum fluctuation was found to be increased. IV. For the series of triangular roughness (11 triangles), the pattern of wavy surface profile was found to be varied with flow condition. To study the fluctuation of water surface profile with single and double Gaussian Shaped roughness on channel bed and side wall for different flow condition and finally to find out how the water surface profile changes in the crest and near the trough of roughness if the bed surface is wavy or highly rough, is the main focus of this study. II. METHODOLOGY The tests were done in the laboratory multipurpose tilting flume, changing the following parameters. Size of roughness Discharges Weir tail gates are made by wood which were inserted at the downstream of the channel to control the depth of the flow. Fig. 3 shows the pictorial view of laboratory tilting flume in which experiments were performed to determine the water surface profile. Fig. 3 Pictorial view of laboratory tilting flume Table 1 shows the dimensions of the different types of Gaussian shaped roughness. Table 1 Dimensions of the different types of Gaussian shaped roughness Series No. Height of the roughness Z o(cm) Base length of roughness L o(cm) Length of the full roughness L R(cm) 1(Gaussian 2 12 12 2 2 12 24 3( 4 12 12 4 4 12 24 ( 6 12 12 6 6 12 24 The multipurpose tilting flume is an artificial channel having a channel width of 8cm. Its maximum upstream water depth is 2cm. Discharge throughout the channel can be changed by rotating the key-wheel attached to the pump. Discharge throughout the channel varies from (.1-.43) m3/s. Its slope can also be changed by rotating the wheel at upstream of the channel. In this project the data of water depths was measured from the crest of roughness up to 4cm right and 4cm left taking an interval of 1cm. Fig. 2 shows the pictorial view of Gaussian shaped roughness of 2cm, 4cm, 6cm height which is made by wood. III. RESULTS AND DISCUSSIONS For the different heights and base lengths of the Gaussian shaped roughness, experimental investigations on the changes of water surface profile in an open rectangular channel were performed. Now results are presented for different cases and discussions are made based on the comparison of the results. The experimental conditions (flow and roughness) are given in Table 2. The results and discussion are given below case by case. Base length of each roughness is 12cm and channel slope is 1/296. S1, S6, S8 for Single Gaussian Shaped S2, S3, S4, S, S7 and S9 for Double Gaussian Shaped. Fig. 2 Gaussian roughness of different height 1874

Water depth(cm) Water depth(cm) Serie s S1 S2 S3 S4 S S6 S7 Case s S1 S2 S3 S6 S7 S8 S9 Height of roughnes s Zo(cm) Length of the full roughness L R(cm) Table 2 Different Test Cases Discharg e Q (m 3 /s) u/s depth mm, bottom d/s depth e mm sid d/s depth mm side Velocity V (m/s) Fraud no. 6 12.33 141. 1 1.29.2 6 24.1 114 86 91.17.16 6 24.19 123 91 92.2.18 4 12 8.33 116 8 83.36.34 4 24.33 12 84 9.36.34 2 12.33 94 72 71.44.46 2 24.33 96 73 72.44.46 A. Flow profile for series S1 (Gaussian shaped roughness of 1 nos x 6cm x l2cm) Fig. 4 shows the comparison of Water surface profile over a Gaussian and Triangular shaped roughness placed at channel bottom and channel sidewall. It is found that the water depths at upstream for all the cases are constant that are gradually decreasing from the middle of roughness. At the downstream of the trough of roughness the water depth is the lowest. 2 2 1 Waterdepth(bottom roughness) Water depth(side roughness) 2 3 4 6 7 Fig. 4 Water surface profile over a Gaussian shaped roughness (Case S1) The water depth is increasing gradually towards downstream. The water level fluctuation for single Gaussian shaped roughness is found 4.cm when it is placed at bottom and 3.1cm when it is placed at side wall. Similarly for triangular roughness, it is found 3.cm for bottom and 3.1cm for sidewall. It is also found that the lowest point of water depth for Gaussian shaped roughness is formed earlier than the Triangular roughness. B. Flow profile for series S2 (Gaussian shaped roughness of 2 nos x 6cm x l2cm) Fig. shows the water surface profile over a Gaussian sidewall. It is found that the water depths at upstream for all the cases are constant that are starts to decrease gradually 2 2 1 Water depth( Gaussian shaped bottom roughness) Water depth( Gaussian shaped side roughness) 2 3 4 6 7 8 Fig. Water surface profile over a Gaussian shaped roughness (Case S2) Again the water surface is gradually increased and shows the peak at the end of roughness. Then it falls down towards the downstream next to 2nd roughness. Again it starts to increase gradually towards downstream. The water level fluctuation for double Gaussian shaped roughness is found 2.8cm, when it is placed at bottom and 2.3cm when it is placed at side wall. The positions of lowest points are formed very near to the trough of roughness. C. Flow profile for series S3 (Gaussian shaped roughness of 2 nos x 6cm x l2cm) Fig. 6 shows the water surface profile over a Gaussian sidewall. It is found that the water depths at upstream for all the cases are constant that are starts to decrease gradually from the 1st roughness. The water depth is the lowest before the crest of 2nd roughness. Again the water surface is gradually increased and shows the peak at the end of roughness. Then it falls down towards the downstream next to 2nd roughness. Again it starts to increase gradually towards downstream. The water level fluctuation for double Gaussian shaped roughness is found 3.2cm when it is placed at bottom and 3.1cm when it is placed at side wall. The positions of lowest points are formed very near to the trough of roughness. from the 1st roughness. The water depth is the lowest before the crest of 2nd roughness. 187

2 2 1 Water depth( Gausian shaped bottom roughness) 2 3 4 6 7 8 Fig. 6 Water surface profile over a Gaussian shaped roughness (Case S3) D. Flow profile for series S4 (Gaussian shaped roughness of 1nos x 4cm x l2cm) Fig. 7 shows the water surface profile over a Gaussian sidewall. It is found that the water depths at the upstream are constant those are gradually decreasing from the middle of roughness. At the downstream of the trough of roughness the water depth is the lowest. Then the water depth is increasing gradually towards downstream and flows at a constant water depth towards downstream. The water level fluctuation for single Gaussian shaped roughness is found 3.1cm, when it is placed at bottom and 3.3cm when it is placed at side wall. It is observed that the lower point of water depth is formed at a distance of 2cm downstream of trough for side roughness, it is found at a distance of 7cm downstream for the case of bottom roughness. 1 Water depth(gaussian shaped bottom roughness) upstream of the roughness are constant which are gradually decreasing from the middle of roughness. At the downstream of the trough of roughness the water depth is the lowest. Then the water depth is increasing gradually towards downstream and flows at a constant water depth towards downstream. The water level fluctuation for single Gaussian shaped roughness is found 2.2cm, when it is placed at bottom and 2.3cm when it is placed at side wall. It is observed that the lower point of water depth is formed at a distance of 4cm downstream of trough for both the side & bottom roughness. 1 Water depth(gausssian shaped bottom roughness) 2 3 4 6 7 8 9 Fig. 9 Water surface profile over a Gaussian shaped roughness (Case S6) G. Flow profile for series S7 (Gaussian shaped roughness of 2nos x 2cm x l2cm) Fig. shows the water surface profile over a double Gaussian shaped roughness placed at channel bottom and channel sidewall. It is found the water level fluctuation for single Gaussian shaped roughness is found 2.3cm, when it is placed at bottom and 2.4cm when it is placed at side wall. 2 3 4 6 7 8 9 Fig. 7 Water surface profile over a Gaussian shaped roughness (Case S4) E. Flow profile for series S (Gaussian shaped roughness of 2nos x 4cm x l2cm) Fig. 8 shows the water surface profile over a Gaussian sidewall. It is found that the water level fluctuation for double Gaussian shaped roughness is found 3.6cm, when it is placed at bottom and 3.cm when it is placed at side wall. 1 Water depth(gaussian shaped bottom roughness) 2 3 4 6 7 8 Fig. 8 Water surface profile over a Gaussian shaped roughness (Case S3) 1 Water depth(gausssian shaped bottom roughness) 2 3 4 6 7 8 Fig. Water surface profile over a Gaussian shaped roughness (Case S7C2) Comparison of Maximum Fluctuation In addition to experimental evaluation, analytical solution for lowering of water surface over a raised bottom is calculated analytically based on the specific energy equation using the concept of transition problem. The comparisons of experimentally and analytically computed values on water surface fluctuation are shown in Table 3. Table 3 shows the Water depth and maximum fluctuation of water surface for different cases. F. Flow profile for series S6 (Gaussian shaped roughness of 1nos x 2cm x l2cm) Fig. 9 shows the water surface profile over a single Gaussian shaped roughness placed at channel bottom and channel sidewall. It is found that the water depths at the 1876

Cases Table 3 Water depth and maximum fluctuation of water surface for different cases Size of the roughness Fraud no. Max. Fluct. (mm) Gaussian. bottom rough) Max. Fluct. (mm) (Gaussian side rough) Max. Fluct. (mm) Calc. 1. 1 nos 6 12.2 4. 31. 36. 2. 2 nos 6 12.16 28 23 29 3. 2 nos 6 12.18 32 31 33.6 4. 1 nos 4 12.34 31 33 34. 2 nos 4 12.34 36 3 34 6. 1 nos 2 12.46 22 23 2 7. 2 nos 2 12.46 23 24 21 From the Table 3, it was found that the fluctuation of water depth is greater for cases (1-3) when the roughness placed at bottom compared to its placement at side. However, for cases (4-7) when the roughness height reduced the fluctuation for two placements does not vary. This is probably due to the boundary effect of sidewall. It is also found that, with increase in fraud number the fluctuation of water depth is increased. The measured fluctuation is found very near to calculated amount. IV. CONCLUSION From the above results and discussion it can be concluded that The fluctuation of water depth is greater for greater height of roughness and lower for smaller height of roughness. For roughness height of 6 cm, the fluctuation of water depth is found greater when the roughness placed at bottom compared to its placement at side. However, when the roughness height reduced the fluctuation for two placements does not vary. This is probably due to the boundary effect of sidewall. guidance, suggestions, cooperation, and encouragement at every stage of this research work. Thanks are also due to the staff of the department and other officer of the Engineering Section for their informative helps or other type of co-operations. REFERENCES [1] Ahmed, T. (211). Experimental Investigation on Changes of Water Surface Profile with Bottom, UG Thesis, Dept. of Civil Engineering, KUET. Bangladesh. [2] Chow, V.T. (199). Open Channel Hydraulics, McGraw-Hill Company, Inc. New York, USA, pp. 39-61. [3] Erpicum, S. (26). Comparison of 2D turbulence Models for steady flow computation in a macro-rough channel, PhD Thesis. Liege, Unite HACH, University of Liege. Belgium. [4] Nabi, A.M. (28). A 3D Model of Detailed Hydrodynamics with Sediment Transport for Simulation of Subaqueous Dunes, 2nd International Symposium on Shallow Flows, HKUST, Hong Kong. pp. 83. Md. Rafiue Islam a B. Sc. in Civil Engineering, The measured fluctuation is found very near to calculated amount. Khulna University of Engineering & Technology. With increase in fraud no. fluctuation of water depth is increased. ACKNOWLEDGMENT The authors wish to express their gratitude and heartiest thanks to those persons without whom this work never be completed and specially indebtedness to Prof. Dr. Md. Shahjahan Ali, Department of Civil Engineering, Khulna University of Engineering & Technology, for his expert Shariful Islam b* B. Sc. in Civil Engineering, Khulna University of Engineering & Technology. 1877

Md. Abdul Qaiyum Talukder c B. Sc. in Civil Engineering, Khulna University of Engineering & Technology. S. M. Rezwan Hossain d B. Sc. in Civil Engineering Khulna University of Engineering & Technology 1878