Thursday, January 24, 2008

Partial dambreak

To run a simple test of two-dimensional problem with SUIKOU-2D.FOR, there are seven input files required. Say, an example of simple dambreak problem. The computational domain is 200m x 200m and 5m x 5m mesh. The initial water level is maintained as 10m at the upstream dam and tailwater was 5m.

The input files to run the codes are:
geo.txt --> contains the geometry information on each node.

mesh.txt --> describe the mesh size and its numbers (dx,dy, nx,ny)
ctrl.txt
--> setting up parameters, and hydraulic constants.
elv.txt
--> contains the information of bed elevation data.
ho.txt
--> initial water depth.
uo.txt --> initial velocity in x-direction.
vo.txt
--> initial velocity in y-direction.

Screenshot of Running program :

The output file is ready formatted for GNUPLOT visualization.

  • Assumed it is in the directory of output files, then simply run the GNUPlot.
  • Type a command shown in the figure below (for advanced plotting, please refer GNUPLOT manual).


and it gave :



If it's need, the animations could be made.

The Water Surface animation

The vector velocity animation

(Reference: Open Channel-Flow, Chaudhry M.H, 1993. pp. 371-374)


No comments: