SUIKOU for One-Dimensional Problems

For the particular case study of one-dimensional flow, SUIKOU-1D.FOR was developed. There are two input files need to run. First is GEO.TXT which contains the grid of computational area. Another one is a control file, CNTL.TXT which contains parameters and program controls.

Let’s take a simple example of 1D flow over rough-variable topography. There is a single straight channel consists with 4 reaches and different bed slopes (So1=0.03, Lo1=40m; So2=0.008, Lo2=40m; So3=0.03, Lo3=40m; So4=0.011, Lo4=40m). The unit discharge is 2.36m2/s at the upstream and the water depth at downstream was set to 0.8m. Manning roughness coefficient was assumed as 0.03. The grid size is 1m. The downstream end elevation is set to be 0m.

The input files is described as below .

GEO.TXT is containing the information of channel properties. They are 4-coloumns separated by tab. First column is a node number. Then it's followed by the channel width (B) at the second column. Third column is representing the datum elevation on each node and the fourth-column is representing the initial water level to start the simulation.

CNTL.TXT is containing simulation control and hydraulic constants. First row is representing the upstream and downstream boundary conditions respectively. The second row is n manning roughness, bank slope and the gravity constant which used in the simulation. Third row is describing the momentum coefficient (), optional switch of artificial viscosity and-coefficient . Lasty, the fourth column is containing time step (dt) and total simulation time in seconds.

In the 1D code, it writes only one output file namely as OUT.TXT. This is a text file containing six columns and separated by tab. First column is indicating the node numbers and the second one is datum elevation of nodes. The third column is the elevation of water surface. The fourth column is the discharge unit. Finally, the fifth column represents the flow velocity whilts the sixth column is the froude number. All units are described on the metric system (SI).

Here is the detail format of OUT.TXT.

To visualize the water level along the channel by GNUPLOT, the following simple commands can be used.

• Assume it's the directory contains the output files and GNU-Plot is running.
• Type a command as shown in the figure below (for advanced plotting, please refer GNUPLOT manual):

For the discharge along the channel :

Animations: