Problem Statement
Open channel particularly high-speed channels are used for drainage in urban parts, since urban sprawl addition rainfall overflow due to altered land usage. Flood control channels are designed and built to safely pull off the awaited hydrologic burden. The desire is to minimise the H2O ‘s clip of abode in the urban country. The channels are designed to transport supercritical flow to cut down the H2O deepnesss and the needed path. Structures, such as decompression sicknesss and passages cause flow to choke and organize leaps. These hydraulic conditions by and large necessitate higher walls, Bridgess and other dearly-won containment constructions. A ill designed channel can do bank eroding, damaged equipment, increased operating disbursal, and decreased efficiency ( Berger et. al.1995 ) . Furthermore, crossings may be washed out, and the town may deluge.
Predicting the possible location of dazes and finding the lift of H2O surface in channel are necessary to measure and make up one’s mind the needed sidewall highs. Normally empirical equations are frequently used in the channel design due to its simple application. However, the presence of decompression sicknesss, contractions, passages, meetings, span wharfs and entree inclines can do the flow to choke or to bring forth a series of standing moving ridges and these all will perplex channel design.
In the yesteryear, applications of physical theoretical accounts are common for this H2O profile rating. Although physical theoretical account can reproduce a channel if decently conducted, but great attention must be taken in theoretical account dimension and graduated table. A major drawback of physical theoretical accounts is the job of scaling down a field state of affairs to the dimensions of a research lab theoretical account. Phenomena measured at the graduated table of a physical theoretical account are frequently different from conditions observed in the field. Though physical theoretical accounts can reproduce inside informations of existent hydraulic constructions, they are still subjected to the restriction of graduated table mold because sometimes it is impossible to reproduce the physical job to graduated table.
Changes to the physical theoretical account necessitate a “cut and try” technique that involves rupturing down the unwanted subdivisions of the channel and reconstructing them with the new coveted design. Due to the clip and cost restraints of physical theoretical accounts, it is non practical to analyze a broad scope of designs. This could ensue in hydraulic public presentation that is merely acceptable over a limited scope.
Mathematical theoretical accounts have been developed to get the better of the job mentioned supra. A mathematical theoretical account consists of a set of differential equations that are known to regulate the flow of surface H2O. The dependability of anticipations of theoretical accounts depends on how good the theoretical account approximates the field state of affairs. Inevitably, simplifying premises must be made because the field state of affairs is excessively complex to be simulated precisely. Normally, the premises necessary to work out a mathematical theoretical account analytically are reasonably restrictive. To cover with more realistic state of affairss, it is normally necessary to work out the mathematical theoretical account about utilizing numerical techniques. Therefore, an cheap and a readily available theoretical account for measuring these channels are needed. A numerical theoretical account is a logical attack.
An country of technology design that can profit the usage of numerical theoretical account is the design and alteration of high-velocity channels indispensable for the routing of floodwater through urban countries. The proper design of new channels and re-design of bing channels is required to avoid such things as bank eroding, damaged equipment, increased operating disbursals, implosion therapy, and higher building costs. By utilizing numerical theoretical account, a better channel design can be produced with minimal cost and clip.
Aim Of The Study
The primary intent of the research is to develop a methodological analysis and determine the effectivity of utilizing numerical theoretical account for unfastened channel mold. The challenges for this numerical theoretical account prevarication in stand foring supercritical passages and capturing the possible location and motion of the dazes. The specific aims of the survey are listed as followed:
1. To measure the practicality of utilizing planar numerical theoretical account to assistance in the design of a realistic unfastened channel.
2. To measure the public presentation of numerical theoretical account in managing daze capturing in assorted trial instances through comparing with published consequences, laboratory trials and analytical solutions.
Significance Of Research
In surface H2O modeling, the most ambitious portion is to observe the location and H2O lift of hydraulic leap or daze. The tallness of the leap is critical to the design of channel walls and Bridgess within high-speed channel. And through this anticipation besides, we can specify easy the critical location within bing channel so that betterment can be done rapidly before inundation happen in that location. A batch of flow theoretical accounts used late non able to execute this undertaking accurately. However, there is still some flow theoretical accounts were developed specially for this daze gaining control purpose but most of them in unidimensional ( 1D ) manner.
There was some concern to the adequateness of a unidimensional ( 1D ) analysis of the flow conditions such as contractions, enlargements, decompression sicknesss, hydraulic leaps and span wharfs which normally found in high-velocity channels. There was a inquiry as to whether calculating cross-sectional averaged flow variables provided a sufficiently accurate estimation of flow deepnesss and speeds within these boundary characteristics. Therefore, a planar ( 2D ) analysis was deemed necessary to measure these flow conditions which ever cause overhead problem in high- speed channels.
A numerical theoretical account HIVEL2D used to measure the design computationally before the building of the physical theoretical account Begins and to test options. Using a numerical theoretical account would speed up this design procedure and lead to an improved initial physical theoretical account therefore cut downing the clip spent on the physical theoretical account. This would let for geographic expedition of more design options in a shorter length of clip ensuing in a more cost-efficient solution.
Previous Study
A batch of research documents were published to demo the theoretical account simulation and confirmation of open-channel flows in assorted trial instances. Different techniques had been applied such as finite-difference method and finite-element method. MacCormack and Gabuti explicit finite-difference strategy were introduced by Fennema et. Al. ( 1990 ) to incorporate the equations depicting 2D, unsteady bit by bit varied flows, by presuming hydrostatic force per unit area distribution, little incline and unvarying speed distribution in perpendicular way.
The same finite-difference strategy ( MacCormack ) was used to imitate contraction instances ( Jimenez et. Al. 1988 ) . Here, the shallow H2O equation was used as a basic equation. For boundary status, Abbett process was applied. The basic thought of this process is to use the numerical strategy up to the wall utilizing nonreversible differences as a first measure. Then to implement the surface tangency demand, a simple moving ridge is superimposed on the solution to do the flow analogue to the wall.
The comparing between computed and measured consequences indicated that there are some instances for which the premise of hydrostatic force per unit area distribution is excessively restrictive. In these state of affairss, the usage of more general equations, e.g. , Boussinesq type equations that include perpendicular acceleration effects, becomes desirable. In that survey, computed consequences were compared with contraction trial instances which conducted by Ippen et. Al. ( 1951 ) . The fake H2O deepness increased four times within a short distance. The dissension between the experimental and computed consequences becomes big.
Gharangik et. Al. 1991 used unidimensional Boussinesq equations were used to work out hydraulic leap job in a horizontal rectangular channel. Again, MacCormack and two-four expressed finite-difference strategies were used for solution until a steady province was reached. Experiments with the Froude figure upstream of leap runing from 2.3 to 7.0 were conducted for theoretical account confirmation. The importance of the Boussinesq footings was investigated. Results show that the Boussinesq footings have small effects in finding the leap location. However, consequences from this survey will be used for theoretical account simulation in this survey, as discussed in the undermentioned subdivision.
In work outing open-channel flows job, shallow H2O equations are really frequently used by research workers together with finite-element method and Galerkin strategy. Schwanenberg et. Al. ( 2004 ) had developed a entire fluctuation decreasing Runge Kutta discontinuous Galerkin finite-element method for 2D depth-averaged shallow H2O equations. In his survey, the smooth parts utilizing the 2nd order strategy for additive elements and 3rd order for quadratic form maps both in clip and infinite. In that theoretical account, dazes were usually captured within two elements. 5 trial instances including the existent dike interruption of Malpasset, France, indicated a well public presentation of the strategy.
Hicks et. Al. ( 1997 ) proved that a 1D preparation besides can provides an first-class solution in patterning dam-break inundations in natural channels. St. Venant equations were used in the theoretical account, which solved with the characteristic dissipative Galerkin finite-element method ( CDG ) . The computational simulations were conducted utilizing both varied and unvarying spacial discretization. Confirmation was made by comparing dam break experiment from Bellos et. Al. ( 1992 ) , which was performed in a rectangular channel of changing breadths.
A discrepancy of the Galerkin strategy for preservation Torahs in 2D, about horizontal flow, which exhibits a singular shock-capturing ability, was presented ( Katopodes 1984 ) . The method was based on discontinuous burdening maps which introduce weather effects in the solution while keeping cardinal difference truth. However, the cardinal hypothesis concerns the perpendicular distribution of force per unit area is hydrostatic.
Katapodes presented comparison consequences between analytical solution, classical Galerkin solution and Pseudo-viscosity solution in a sudden H2O release trial instance.
The finite-element Galerkin was found really dissatisfactory, although non worse than non-dissipative finite-difference methods. In Galerkin solution, the jobs such as parasitic moving ridges behind the forepart and the spreading of discontinuity over elements were found. However, the survey demonstrated that better consequences can be obtained by a fluctuation of the Galerkin technique known as the Petrov-Galerkin preparation. The confirmation was made by comparison to analytical solutions for 4 trial instances ( 1D rush, rush through symmetric gradual bottleneck, rush through asymmetric disconnected bottleneck and enlargement ) .
The planar vertically averaged and minute equations theoretical account, developed by Ghamry et. Al. ( 2002 ) was used to analyze the consequence of using different distribution forms for speeds and force per unit area on the simulation of curving unfastened channels. Linear and quadratic distribution forms were assumed for the horizontal speed meanwhile a quadratic distribution form was considered for perpendicular speed. Linear hydrostatic and quadratic non-hydrostatic distribution forms were suggested for force per unit area. The finite component intercrossed Petrov-Galerkin and Bubnov-
Galerkin Schemes Were Used.
Comparisons of the theoretical account anticipations were made with the experimental consequences obtained in “S” form unfastened channel, “U” form of rectangular gulch and 270 grade curved rectangular gulch. Note that merely subcritical flows were simulated in all experiment with Fr & lt ; 5.0. In all comparing, merely the longitude speed distribution was focused. Consequences suggested that pre-assumed speed distribution forms are non really sensitive ; moreover the attained higher truth on using the non-hydrostatic premise theoretical account is undistinguished compared to linear hydrostatic theoretical account.
A entire fluctuation decreasing Runge Kutta Discontinuous Galerkin ( RKDG ) finite-element method for planar depth-averaged shallow H2O equations has been developed by D.Schwanenberg and M.Harms in twelvemonth 2004. The expressed clip integrating, together with the usage of extraneous form maps, makes it every bit efficient as comparable finite-volume strategies. The method was shown to hold 2nd or 3rd order of convergence in clip and infinite for additive and quadratic form maps in smooth parts of the solution and crisp representation of dazes. The trial indicate an first-class public presentation of the strategy and giving suggestion that advanced analysis utilizing full 3D Navier-Stokes equation is possible and can be conducted.
Models normally face trouble in managing leaps. One of the methods called “shock tracking” that track the leap location and enforce an internal boundary at that place. The shallow H2O equation so allows weak solutions in which a discontinuity represents the hydraulic leap. This is referred as “shock capturing” as originated by von Neumann and Richtmyer ( 1950 ) . Note that this might be non easy for research workers to track daze location accurately. Furthermore, great attention must be taken to guarantee that the mistakes merely local to the leap ( discontinuity location ) .
Normally a theoretical account with uninterrupted deepnesss will conserve mass and impulse through the leap but will besides bring forth oscillation at the shortest wavelengths to conserve energy. Energy dissipation which should look in leap does non be. In fact, when leaps happen, energy is being transferred into perpendicular gesture. And since perpendicular gesture is non included in shallow H2O equation, it causes some lost in theoretical account. Therefore, a strategy is needed to turn to this job will be dissipative and can fulfill the demand of daze capturing every bit good.
In 1995, a 2D finite-element theoretical account for the shallow-water equations was produced utilizing an extension of the streamline upwind Petrov-Galerkin ( SUPG ) construct. A mechanical was implemented which detects the presence of a leap by ciphering the mechanical energy fluctuation per component and so allows the theoretical account to increase the grade of up weaving in the daze locality while keeping more precise solutions in drum sander flow parts ( Berger et. Al. 1995 ) .
Consequences from Berger demonstrated the ability of theoretical account to reproduce the velocity and tallness of a traveling hydraulic leap and the ability of the shock-detection mechanism to follow the leap. This was a comparing with an analytical solution. A 2D illustration of a supercritical contraction was so demonstrated by comparing with flume consequences by Ippen and Dawson ( 1951 ) . Finally, the information from the survey of Margarita Channel was used for theoretical account confirmation excessively. Consequences showed that the theoretical account is equal to turn to hydraulic jobs affecting leaps and oblique dazes.
Previously, finite-element methods were found can non conserve mass locally. However, Berger et. Al. ( 2002 ) demonstrated that, by utilizing the flux inherent in the discrete, finite-element preservation statement, the amount of the fluxed around an component or group of elements exactly matches the internal mass alteration. These happening were supported by computations in one and three dimensional ( see Berger et. Al. 2002 ) .
A planar numerical flow theoretical account for trapezoidal high-velocity channels which holding spilling sidewall was developed by Stockstill et. Al. ( 1997 ) . This theoretical account was developed after bettering the theoretical account introduced by Berger et. Al. ( 1995 ) . When dainty with spilling sidewall where the deepness is unknown, an attack involve updating the traveling boundary supplanting merely one time each clip measure was applied. For interior nodes, big supplanting of the traveling boundary nodes can take to element form deformations. This job was solved by regridding the side slopes each clip measure as a map of the boundary nodal supplanting.
A trapezoidal gulch with horizontal curve was conducted in U.S. Army Engineer Waterways Experiment Station, Hydraulic Laboratory for theoretical account confirmation ( Stockstill et. Al. 1997 ) . The first trial status demonstrated that the theoretical account accurately solved the H2O lines through the passage where the flow accelerated from subcritical to supercritical. The experiment was so repeated by adding wharfs. The theoretical account was found unable to depict undular leaps which were formed in the trial, but accurately represented the clogged flow status and the maximal deepness. Overall consequences showed that this method is utile in subcritical flow but non so efficient in supercritical flow. However, it was proved to be stable at important Courant Numberss.
The numerical theoretical account which introduced by Berger and Stockstill was further extended its application on imitating flatboat drawdown and currents in channel and backwater countries ( Stockstill et. Al. 2001 ) . Vessel effects were modeled numerically by utilizing a moving force per unit area field to stand for the vas ‘s supplanting. Verification theoretical account included existent field informations such as Illinois State Water Survey, Mississippi River and Sundown Bay where located along Texas seashore between Aransas Bay and San Antonio Bay. The theoretical account was shown able to reproduce chief channel return currents in consecutive ranges of little channels ( Illinois Waterway ) and in the off-channel countries of broad rivers ( Mississippi River ) .
Another unpublished study from Army Engineer Waterways Experiment Station showed application of 2D numerical theoretical account which was introduced by Berger, in San Timoteo Creek which is tributary of Santa Ana River. The proposed design within the range studied includes a deposit basin, a concrete weir followed by a chute holding meeting sidewalls, a compound horizontal curve dwelling of spirals between a round curve and the upstream and downstream tangents with a banked invert, and a span wharf associated with the San Timoteo Canyon Road. The trial had been conducted utilizing two different discharge value, 19000 californium and 12000 californium. These series of trials demonstrated the application of the numerical theoretical account in site. The sensitiveness of fake consequences to the pick of dissipation coefficient and grid declaration was presented. The study concluded that the solution of flow field is non significantly influenced by the dissipation coefficient and grid polish.
Another trial parametric quantity was the roughness coefficient. Different Manning ‘s N ( n = 0.012 and 0.014 ) were applied in the trial and it was found that the maximal deepness was reduced and the moving ridge crests were located farther downstream big smaller Manning ‘s N.
The San Timoteo Creek study besides proved that hydrostatic premise is appropriate in the country of the oblique standing moving ridge initiated at the wharf nose. The perpendicular acceleration in this locality was calculated to be 0.4 comparative to gravitation. It proved that the hydrostatic premise is sensible even in parts where the flow is unsmooth.
In this research, the application of planar finite-element theoretical account for the shallow H2O equations derived by Berger ( 1995 ) , is demonstrated in assorted trial instances such as hydraulic leap, contraction, enlargement, open-channel junction, gradual contraction, span wharf and weir construction. The theoretical account is produced utilizing an extension of the Petrov-Galerkin strategy. A mechanism which detects the nowadays of dazes by ciphering the mechanical energy fluctuation per component is implemented. Model consequences will be compared with analytical solution and published research lab informations. A few research lab trials were carried out for theoretical account simulation. Datas from these experimental surveies will be presented and the general public presentation of flow under assorted trial instances will be described. Through this research, the public presentation of numerical theoretical account will be evaluated and the theoretical account can supply another alternate tool in planing open-channel construction.
Research Methodology
By and large research methodological analysis can be divided into two parts: experimental work and computing machine patterning utilizing an bing numerical theoretical account.
Experimental Plants
Three different hydraulic instances will be conducted in Hydraulic Laboratory. The trial instances consist of contraction & A ; 90 degree enlargement, hydraulic leap and span wharfs. These characteristics are normally found in high-velocity channels which will organize daze moving ridge in unfastened channel. In the experimental work, preliminary plants were conducted for puting up the trapezoidal channel ( perpendicular side incline ( rectangular and three side slopes 1:1, 1:1.5, 1:2 ) in the research lab and for the control trial. Depth measurings, speed measuring are conducted for all instances.
Numeric Model
Using numerical theoretical account to all instances and doing elaborate comparing with experimental informations in order to measure the pertinence of two dimensional theoretical accounts to high speed channel flow.
Mentions
-Berger, R. C. , and Howington, S. E. ( 2002 ) . Discrete Fluxes and Mass Balance in Finite Element. J. Hydr. Engrg. ASCE, 128 ( 1 ) , 87-92.
-Berger, R.C. , and Stockstill, R.L. ( 1995 ) . Finite Element for High-Velocity Channels. J. Hydr. Engrg. ASCE, 121 ( 10 ) , 710-716.
-Borghei S. M. , Jalili, M. R. , and Ghodsian, M. ( 1999 ) . Discharge Coefficient for Sharp-crested Side Weir in Subcritical Flow. J. Hydr. Engrg. ASCE, 125 ( 10 ) , 1051-1056.
-Chaudhry M. H. ( 1993 ) . Open-Channel Flow. Englewood Cliffs, N.J. : Prentice-Hall.
Fagherazzi, S. , Rasetarinera, P. , Hussaini, M. Y. , and Furbish, D. J. ( 2004 ) . Numeric Solution of the Dam-break Problem with a Discontinuous Gelerkin Method. J. Hydr. Engrg. ASCE, 130 ( 6 ) , 532-539.
-Ghamry, H. K. , and Steffler, P. M. ( 2002 ) . Consequence of Applying Different Distribution Shapes for Velocities and Pressure on Simulation of Curved Open Channels. J. Hydr. Engrg. ASCE, 128 ( 11 ) , 969-982.
-Gharangik, A. M. , and Chaudhry, M. H. ( 1991 ) . Numeric Simulation of Hydraulic Jump. J. Hydr. Engrg. ASCE, 117 ( 9 ) , 1195-1211.
-Henderson, F. M. ( 1966 ) . Open Channel Flow, MacMillian, New York, NY.
-Hicks, F. E. , and Steffler, P. M. ( 1992 ) . Characteristic Dissipative Galerkin Scheme for Open-Channel Flow. J. Hydr. Engrg. ASCE, 118 ( 2 ) , 337-352.
-Hicks, F. E. , Steffler, P. M. and Yasmin, N. ( 1997 ) . Unidimensional Dam-break Solutions for Variable Width Channels. J. Hydr. Engrg. ASCE, 123 ( 5 ) , 464-468.
-Hunt, John. , Brunner, G. W. , and Larock, B. E. ( 1999 ) . Flow Transitions in Bridge Backwater Analysis. J. Hydr. Engrg. , ASCE, 125 ( 9 ) , 981-983.
-Jimenez, O. , and Chaudhry, M. H. ( 1988 ) . Calculation of Supercritical Freesurface Flows. J. Hydr. Engrg. , ASCE, 114 ( 4 ) , 377-395.
-Katapodes, N. D. ( 1984 ) . A Dissipative Gelerkin Sheme for Open-channel Flow, Jour, Hydraulic Engineering, ASCE, 110, no. 4 ( April ) :450-66.
-Katopodes N. D. ( 1984 ) . Planar Surges and Shocks in Open Channels. J. Hydr. Engrg. , ASCE, 110 ( 6 ) , 794-812
-Lai Y.G. , Weber, L.J. , and Patal, V.C. ( 2003 ) . Nonhydrostatic 3-dimensional method for hydraulic flow simulation. I: Formuation and confirmation. J. Hydr. Engrg. , ASCE, 129 ( 3 ) , 196-205.
-Liu, M. , Rajaratnam, N. , and Zhu, D. Z. ( 2004 ) . Turbulence Structure of Hydraulic Jumps of Low Froude Numbers. J. Hydr. Engrg. , ASCE, 130 ( 69 ) , 511-520.
-Molinas, A. , Kheireldin, K. , and Wu, B. ( 1998 ) . Shear Stress Around Vertical Wall Abutments. J. Hydr. Engrg. , ASCE, 124 ( 8 ) , 822-830.
-Reinauer, R. , and Hager, W. H. ( 1998 ) . Supercritical Flow in Chute Contraction. J. Hydr. Engrg. , ASCE, 124 ( 1 ) , 55-64.
-Salaheldin T.M. , Imran J. , and Chaudhry, M. H. ( 2004 ) . Numeric Modelling of Three-Dimensional Flow Field Around Circular Piers. J. Hydr. Engrg. , ASCE, 130 ( 2 ) , 91-100.
-Schwanenberg, D. , and Harms, M. ( 2004 ) . Discontinuous Galerkin Finite-Element Method for Transcritical Two-Dimensional Shallow Water Flows. J. Hydr. Engrg. , ASCE, 130 ( 5 ) , 412-421.
-Seckin, G. , Yurtal, R. , and Haktanir, T. ( 1998 ) . Contraction and Expansion Losingss Through Bridge Constrictions. J. Hydr. Engrg. , ASCE, 124 ( 5 ) , 546-549.
-Paper HL-94-7, Waterways Experiment Station, Vicksburg, MS, Octorber. Stockstill, R.L. , Berger, R.C. , and Nece R.E. ( 1997 ) . Planar Flow Model for Trapezoidal High-Velocity Channels. J. Hydr. Engrg. , ASCE, 123 ( 10 ) , 844-852.
-Stockstill, R.L. , Berger, R.C. , and Nece R.E. ( 2001 ) . Imitating Barge Drawdown and Currents in Channel and Backwater Areas. J. Hydr. Engrg. , ASCE, 127 ( 5 ) , 290-298.
-Unami, K. , Kawachi, T. , Munir Babar, M. , and Itagaki, H. ( 1999 ) . Planar Numeric Model of Spillway Flow. J. Hydr. Engrg. , ASCE, 125 ( 4 ) , 369-375.
-Weber, L. J. , Schumate, E. D. , and Mawer, N. ( 2001 ) . Experiments on Flow at a 90° Open-channel Junction. J. Hydr. Engrg. , ASCE, 127 ( 5 ) , 340-350.
-Wu, B. , and Molinas, A. ( 2001 ) . Choked Flows Through Short Contractions. J. Hydr. Engrg. , ASCE, 127 ( 8 ) , 657-662.
-Wu, S. , Rajaratnam, N. ( 1996 ) . Passage from Hydraulic Jump to Open Channel Flow. J. Hydr. Engrg. , ASCE, 122 ( 9 ) , 526-528.