Hydrological Flow Calculator

with No Comments




Calculate the river height – Alpha

All default values are values of the “Schwarze Lütschine” – Switzerland

Created as part of the Accademic Curriculum of the University of Bern – Lecture: “Hydrological Modelling”

Kst – Gauckler–Manning coefficient of river bed:
I – Slope of the hydraulic grade line [m/m]:
X – Array of river cross section x-coordinates in WGS84 [x1, x2, x3,…,xn]:
Y – Array of river cross section y-coordinates in meters over sea level [y1, y2, y3,…,yn]:
Q*1 – Cross-sectional average velocity [m^3/s]:

dxH – Accuracy [m] (0.01=Centimeter accurate):

*1 To calculate the max velocity: leave empty. To calculate the river height at smaller velocities, please define cross-section velocity in m^3/s.

Information

This tool currently only supports cross-sections, where each point in the direction of the river middle is lower than its neighbour in direction of river side. The highest points are therefore the river sides

This tool is still in its Alpha-Phase. Tips are welcome. If someone is interested in a further development, feel free to contact me.

Possible further development

Automatic generation of a 3d visualisation
Allow rivers with complex forms


Changelog

----------------------------------------------------------
27.10.2014:
Dynamic D3 graphics generation (Alpha)
Replace ymin with D3 built in function
----------------------------------------------------------

----------------------------------------------------------
25.10.2014:
Tidy up script
Introduce functions
Introduce display number of iterations
Corrected bug in FOR-LOOP
----------------------------------------------------------

----------------------------------------------------------
18.10.2014:
Introduced option: if Q is empty, calculate highest possible flow
----------------------------------------------------------

----------------------------------------------------------
16.10.2014:
Corrected calculation of hydraulic radius
Automatic diagramm generation
----------------------------------------------------------

----------------------------------------------------------
15.10.2014:
Upload first version
----------------------------------------------------------

Notizen / Bugs

IF y[0] or y[len-1] are not highest point, check for highest point, make new y[0] and y[len-1]