6c64c6bb-bfb6-4dd4-885f-14030e2bb6ff 133 5 7 decimal

1. Example:Main Canal Design - To Find FSD, Vel, Froud Number etc, for Open Channel flow

]]>

Note:Input values highlighted. Default Units - SI

1. Input

]]>

Q.d 10 m 3^s / * :

Design discharge

Q.ol 12 m 3^s / * :

Overloaded discharge

Bed width

BW 4 m * :

Bed slope

S.0 0.00018 :

Side slope

Z 2 :

percent (#1) ％ 1100 / :

Manning's "n"

n 0.015 :

Free board

FB 0.75 m * :

2. General Equations used

]]>

Area y A y BW y Z * + (* :

Perimeter y P BW 2 y * 1 Z 2^+ sqrt * + :

Hydraulic radius y R y A y P / :

Velocity y V Q.d y A / :

Top width y T BW 2 Z * y * + :

Manning's Eqn y Vm 1 n / y R (23 /^* (S.0 sqrt * m 13 /^s / (* :

3. Find Normal Depth y₀ & Critical Depth y_c

It is advisable to plot variation of depth vs. Q to estimate the guess depth (y_g). We see that there

are two values as solutiins, ie. yy=-3.5m & yy=1.5 m approximately. Hence, take y_g = 1.5m.

]]>

Hence, assume y.g 1.5 m * :

Guess depths yy 4 - 2 - 1 - 012 2.5 71 mat m * :

j 1 yy rows range QQ j el yy j el A yy j el Vm * : for

XYPlot.XLimMin 1 el 5 - :

XYPlot.XLimMax 1 el 4 :

Variables for X-Y Plot

(Many thanks to "Mezentsev"

for the latest plug-in).

]]>

XYPlot.YLimMin 1 el 5 - :

XYPlot.YLimMax 1 el 25 :

XYPlot.XTick 1 el 1 :

XYPlot.YTick 1 el 2.5 :

yy QQ augment Q.d 2 1 sys

3.1. Find Normal depth

]]>

3.2. Find Critical depth

]]>

Critical Flow Eqn for 
Trapeziodal sections F1 Q.d 2^y.c T * g.e y.c A (3^* / 1 ≡ (:

Normal depth Eqn. to solve F Q.d y A y Vm * ≡ (:

Use symbolic optimization F Unknowns y 11 mat

F1 Unknowns y.c 11 mat

guess depth y.g 1.5 m *

guess depth y.g 1.5 m *

Try y_g=1m and we get

-3.68 m (-ve) which is

not admissible

]]>

Normal depth y.0 F y.g FindRoot : 1.559 m *

Ctitical depth y.c F1 y.g FindRoot : 0.754 m *

Normal velocity V.0 y.0 V : 0.901 m s / *

Critical velocity V.c y.c V : 2.409 m s / *

Check calculations Q.0 y.0 A V.0 * : 10 m 3^s / *

Check calculations Q.c y.c A V.c * : 10 m 3^s / *

4. Find flow conditions

]]>

Froud Number at yo Fr.0 Q.d 2^y.0 T * g.e y.0 A (3^* / : 0.076

Froud Number at yc Fr.c Q.d 2^y.c T * g.e y.c A (3^* / : 1

Using Formatted label

Fr.0 1 < Sub-Critical Flow black light green 1 3 mat 1 Fr.0 ≤ 1.00001 < Critical Flow red light green 1 3 mat Super Critical red yellow 1 3 mat if if Fr.c 1 < Sub-Critical Flow black light green 1 3 mat 1 Fr.c ≤ 1.00001 < Critical Flow red light green 1 3 mat Super Critical red yellow 1 3 mat if if

Slope of given channel S.0 0.00018

Critical slope for given canal parameters S.c n Q.d * y.c A y.c R (23 /^* m 13 /^s / (* / (2^: 0.002807

Ratio of Sc/So S.c S.0 / 15.6

5. Advanced Method - suitable for programming (using "FindRoot")

5.1. General Equations used

]]>

Area y b z A1 y b y z * + (* :

Perimeter y b z P1 b 2 y * (1 z 2^+ sqrt * + :

Hydraulic radius y b z R1 y b z A1 y b z P1 / :

Manning's Eqn y b z SL V1.m 1 n / y b z R1 (23 /^* (SL sqrt * m 13 /^s / (* :

Top width y b z TT b 2 z * y * + :

y b z V1 Q.d y b z A1 / :

Eqn to find normal depth for 
trapezoidal section b z q SL F7 q y b z A1 y b z SL V1.m * ≡ (:

Eqn to find critical depth for 
Trapeziodal sections b z q FF1 q 2^yy.c b z TT * g.e yy.c b z A1 (3^* / 1 ≡ (:

b z q SL F7 Unknowns b q SL y z 51 mat

b z q FF1 Unknowns b q yy.c z 41 mat

Eqn 1: This is useful for

programmimg as shown below

]]>

guess depth

b z q SL yy yy b z q SL F7 yy FindRoot :

b z q y yy.c b z q FF1 y FindRoot :

Critical depth for Trapez. section for Qd yy.c0 BW Z Q.d 1 m * yy.c : 0.754 m *

Depth for critical slope. Before 
you use this, you must know Sc BW Z Q.d S.c y.g yy 0.754 m *

Normal depth for Trapezoidal section yy.0 BW Z Q.d S.0 y.g yy : 1.559 m *

Critical slope for given canal for Qd n Q.d * yy.c0 A yy.c0 R (23 /^* m 13 /^s / (* / (2^0.002807

Normal depth, for Rect. section BW 0 Q.d S.0 y.g yy 2.581 m *

Critical depth for Rect. section BW 0 Q.d 1 m * yy.c 0.861 m *

Froud Number at normal depth Fr1.0 Q.d 2^yy.0 BW Z TT * g.e yy.0 BW Z A1 (3^* / : 0.076

Froud Number at critical depth Fr1.c Q.d 2^yy.c0 BW Z TT * g.e yy.c0 BW Z A1 (3^* / : 1

Critical depth for Trapez. section for Qol yy2.c0 BW Z Q.ol 1 m * yy.c : 0.839 m *

Critical slope for given canal for Qol n Q.ol * yy2.c0 A yy2.c0 R (23 /^* m 13 /^s / (* / (2^0.002735

Program to find FSD, Vel, Froud Number for a range of bed widths b z q1 sl q2 Find_Y j 1 b rows range y_rect j el b j el 0 q1 sl 1 m * yy : FSD1 j el b j el z q1 sl y_rect j el yy : FSD2 j el b j el z q2 sl y_rect j el yy : AA1 j el FSD1 j el b j el z A1 : AA2 j el FSD2 j el b j el z A1 : YC1 j el b j el z q1 1 m * yy.c : YC2 j el b j el z q2 1 m * yy.c : SC1 j el n q1 * YC1 j el A YC1 j el R (23 /^* m 13 /^s / (* / (2^: SC2 j el n q2 * YC2 j el A YC2 j el R (23 /^* m 13 /^s / (* / (2^: SCR1 j el SC1 j el S.0 / : SCR2 j el SC2 j el S.0 / : VV1 j el q1 AA1 j el / : VV2 j el q2 AA2 j el / : T1 j el b j el 2 z * FSD1 j el * + : T2 j el b j el 2 z * FSD2 j el * + : fr1 j el VV1 j el g.e AA1 j el T1 j el / * sqrt / : fr2 j el VV2 j el g.e AA2 j el T2 j el / * sqrt / : sp j el - : 181 line for A b m / FSD1 m / YC1 m / VV1 m s / / SC1 SCR1 fr1 sp augment eval : A A FSD2 m / YC2 m / VV2 m s / / SC2 SCR2 fr2 augment eval : 31 line :

First find depth for rect. section using Eqn 1. Note we use z=0 for rect. section & guess depth = 1m

If we do not use this step, we may get -ve values for FSD. Then use y_rect as guess depth to find FSD.

]]>

---- FSD for Q_d, using Eqn 1

]]>

---- FSD for Q_ol, using Eqn 1

]]>

Areas for Q_d& for Q_ol

]]>

Critical depths for Q_d & for Q_ol

]]>

Critical slopes

]]>

Ratios of SC1/SCo of SC/S.0

for Q_d & for Q_ol

]]>

Velocities

]]>

Top widths

]]>

Froud numbers

]]>

---- Space

]]>

Define range of bed widths; 
This input can be modified 
suitably in ascending order bb 3 3.25 3.5 3.75 4 4.25 4.5 4.75 591 mat m * :

Ans bb Z Q.d S.0 Q.ol Find_Y :

title BW FSD1 YC1 V1 SC1 SC1/So Fr1 - FSD2 YC2 V2 SC2 SC2/So Fr2 augment :

title1 m m m m/s - - - - m m m/s - - - augment :

title title title1 stack :

Ans title Ans stack :

Q.d 10 m 3^s / *

Q.ol 12 m 3^s / *

Suffix 1: Canal parameters for Qd

Suffix 2: Canal parameters for Qol

]]>

Ans BW FSD1 YC1 V1 SC1 SC1/So Fr1 - FSD2 YC2 V2 SC2 SC2/So Fr2 m m m m/s - - - - m m m/s - - - 3 1.708 0.856 0.912 0.002 9.811 0.276 - 1.864 0.947 0.957 0.002 9.692 0.279 3.25 1.669 0.828 0.91 0.002 11.086 0.276 - 1.823 0.918 0.955 0.002 10.907 0.279 3.5 1.63 0.801 0.907 0.002 12.473 0.276 - 1.783 0.89 0.952 0.002 12.228 0.279 3.75 1.594 0.777 0.904 0.003 13.976 0.276 - 1.745 0.864 0.949 0.002 13.656 0.279 4 1.559 0.754 0.901 0.003 15.597 0.276 - 1.709 0.839 0.947 0.003 15.196 0.279 4.25 1.525 0.732 0.898 0.003 17.339 0.276 - 1.674 0.815 0.944 0.003 16.85 0.28 4.5 1.493 0.712 0.894 0.003 19.206 0.276 - 1.64 0.794 0.94 0.003 18.621 0.28 4.75 1.462 0.692 0.891 0.004 21.2 0.276 - 1.608 0.773 0.937 0.004 20.512 0.28 5 1.433 0.674 0.887 0.004 23.323 0.276 - 1.576 0.753 0.934 0.004 22.524 0.28 1114 mat

Get the full path of the current document MyDirectory DocumentDirectory : E:\All Smath Files in ASUS\

Export as Excel file to current document directory Ans Canal_Design.xlsx exportData.XLSX 1

The formatted Excel file is shown as an image below (for demonstration only).

]]>