CCHE2D 2.2 User's Manual
CCHE2D 2.2 User's Manual
CCHE2D 2.2 User's Manual
Yaoxin Zhang
School of Engineering
The University of Mississippi
University, MS 38677
May 2005
Yaoxin Zhang
Post-doctoral Research Associate
May 2005
Table of Contents
1 INTRODUCTION...................................................................................................................... 1
1.1 FAMILY OF CCHE2D MODEL ................................................................................................ 1
1.2 USING THIS MANUAL............................................................................................................. 2
1.3 RELATED DOCUMENTS........................................................................................................... 3
2 CCHE2D MODEL ..................................................................................................................... 4
2.1 NUMERICAL SIMULATION ...................................................................................................... 4
2.1.1 INTRODUCTION ............................................................................................................................ 4
2.1.2 GENERAL PROCEDURE OF USING CCHE2D MODEL ................................................................... 5
2.1.2.1 Mesh Generation.................................................................................................................. 5
2.1.2.2 Boundary Conditions ........................................................................................................... 6
2.1.2.3 Parameter Setting ................................................................................................................. 6
2.1.2.4 Simulation ............................................................................................................................ 6
2.1.2.5 Results Interpretation ........................................................................................................... 6
Table of Contents
ii
1 Introduction
Chapter 1 Introduction
CCHE2D-GUI
Graphical Users Interface
Input
Output
Chapter 1 Introduction
Chapter 4 describes in detail the procedures to run the CCHE2D model through the
GUI. You will learn how to edit grid, set initial conditions for flow, set model
parameters, set the initial bed composition file and sediment boundary conditions
files, set the initial bed material properties, set the monitor points, run simulations.
Chapter 5 shows you how to visualize simulation results for flow and sediment
transport, probe and extract data, set the flood properties and contour properties, edit
the texts, set legend properties, and set the background color and background image.
The appendix gives you a quick reference to the formats of files the users must
provide for the CCHE2D-GUI.
Readers of this manual may also need to read the manual of CCHE2D Mesh Generator, a
separate program that generates the numerical meshes required to start using the CCHE2DGUI.
2 CCHE2D Model
Numerical models approximate physical problems, it has however all the components to
represent the physics to be simulated. Mesh or grid is used to represent channels and
bathymetry; inflow and outflow are defined at inlet and outlet sections as boundary
conditions. All the above and model parameters have to be defined before a simulation can
start. To make setting up a simulation case efficiently, a mesh generator and a User Graphic
Interface have been developed.
Mesh generation
Parameters setting
Simulation
Inlet(s) and outlet(s) should be sufficiently far away from the zones of interest;
Mesh generation particularly for practical problems takes a lot of time, however, the time
shall be paid off if good quality is achieved. In many cases, the simulation code will run with
a low quality mesh but the results maybe less reliable.
2.1.2.4 Simulation
When the mesh is ready, boundary conditions have specified, one is ready to start a round of
simulation. Since the CCHE2D is a model using a time marching scheme, one should also reexam the initial condition. Because the initial condition for a water flow in natural condition
is unknown, cares have to be taken to make sure the guessed initial condition is reasonable,
particularly for unsteady problems.
Stability is a problem the modeler must aware. When the time step is too large for a
particular problem, the simulation will not continue or it may produce totally unreasonable
results. In these situations, one should reduce the time step used and retry, until the solutions
becomes stable. Although an implicit scheme is used in the CCHE2D model, time step can
not be set arbitrarily.
Numerical model is developed under certain assumptions, the results should be reliable if the
simulated flow satisfies these assumptions. However, the natural is always complicated, these
assumptions could never be satisfied one hundred percent. If the results are questioned, we
recommend validating them by comparing with physical model measurements or field data.
Once the model is validated with site specific data, it could be used to study the trend of the
flow or sediment transport processes by varying concerned parameters
(2.1)
Momentum Equations:
u
u
u
Z 1 ( h xx ) ( h xy ) bx
+u
+v
= g
+ [
+
]
+ f Cor v
t
x
y
x h
x
y
h
(2.2)
Z 1 (h yx ) (h yy ) by
v
v
v
+ [
+
]
f Cor u
= g
+u +v
h
y h
x
y
y
x
t
(2.3)
where u and v are the depth-integrated velocity components in the x and y directions
respectively; g is the gravitational acceleration; Z is the water surface elevation; is water
density; h is the local water depth; fCor is the Coriolis parameter; xx , xy , yx and yy are the
depth integrated Reynolds stresses; and bx and by are shear stresses on the bed surface.
xx = 2 t
u
x
(2.4a)
xy = yx = t (
yy = 2 t
u v
+ )
y x
v
y
(2.4b)
(2.4c)
t =
Axy
6
U * h
(2.5)
where Axy is an adjustable coefficient of eddy viscosity, is the von Karman constant, and
t = l 2 2(
U 2
u v
v
u 2
) + 2( ) 2 + ( + ) 2 + (
)
z
x x
x
x
(2.6a)
1
z
l = z (1 )dz = h 1 d 0.267h
h
h
0
(2.6b)
U
U*
= Cm
z
h
(2.6c)
where C m is a coefficient with a value of 2.34375 so that Equation (2.6) will cover Equation
(2.5) in the case of a uniform flow in which all the horizontal velocity gradients vanish.
k
k
k k
k
+u
+v [ t
] [ t
] = P + PkV
t
x
y x k x y k y
(2.7)
2
+u
+v
[
] [ t
] = c1 P c 2
+ PV
t
x
y x x y y
k
k
(2.8)
where
P = u ''i u 'j u i , j = t [2(
PkV = C k
U *3
h
u 2
v
u v
) + 2( ) 2 + ( + ) 2 ]
x
y
x y
PV = C
U *4
h2
(2.10)
U * = c f (u 2 + v 2 )
Ck =
(2.11)
C = 3 .6
cf
(2.9)
c 2
c 3f / 4
(2.12)
t =
c k 2
(2.13)
In the above equations, the following values are used for the empirical constants:
c = 0.09, c1 = 1.45, C 2 = 1.90, k = 1.0, = 1.3 .
10
Figure 2-1
11
3 Overview of CCHE2D-GUI
3.1 Introduction
This chapter explores the interface of the CCHE2D-GUI. The functions of the main window,
the toolbars, and the menus will be explained in details to help you get familiar with this
interactive graphical environment.
Figure 3-1
13
The main window provides an interactive platform for the users. As shown in Fig. 3-1, it
contains all the functional components or the accesses to these components. It is composed of
the title bar (displays the file path to the data set), the menu bar, the toolbars (main toolbar
and editing toolbars), the status bar, and the plot area.
3.3 Toolbars
Main Toolbar
Add/Remove XY Axis
Plot Points
Add/Remove Legend
About
Add/Remove Title
Simulation Wizards
Show/Hide
Run simulations
Add Texts
Select Pointer
Data Probe
Zoom In
Grid On/Off
Pan
Vector On/Off
Incremental Zoom In
Display variable
Background
Delete I Line
14
Add J Line
Delete J Line
Move J Line
Add I Line
Move I Line
Rotate Right
Downward Translation
Zoom In Image
Left Translation
Translation
Right Translation
Rotate Left
Upward Translation
3.4 Menus
3.4.1 File
This is most probably the first menu you will touch when using CCHE2D-GUI. It provides
almost all the accesses to files Input/Output, such as opening Geo and Layout files, saving
Layout files, importing and exporting Bitmap file, printing current image, and etc.
15
Figure 3-2
3.4.2 Simulation
After loading the Geo mesh file, you need to visit this menu to run simulations. According to
the functions, the menu items are divided into four groups, wizard group, run group,
flow group and sediment group (see Figure 3-3).
In wizard group, it provides you a quick guide on how to run simulations.
In run group, it provides you access to run simulations, i.e., Run CCHE2D Model.
In flow group, it provides you editing tools to prepare the initial conditions, boundary
conditions and the model parameters for flow.
In sediment group, it provides you editing tools to prepare the initial conditions, boundary
conditions and the model parameters for sediment.
16
Figure 3-3
3.4.3 View
This menu provides you necessary tools to view the current image more clearly and
accurately, such as Select, Zoom, and Pan, etc.
Figure 3-4
17
3.4.4 Settings
Visiting this menu, you can set up the visualization styles best fit your needs and personality,
such as plot type, colors, adding text, showing/hiding legend, title, grid, boundary,
simulation time, vectors, and etc. You can also set up all kinds of visualization
parameters/properties, such as properties of legend, title, text, flood, contour, vector,
background image, and etc.
Figure 3-5
18
3.4.5 Visualization
Through this menu, you can visualize all kinds of available result files, the intermediate file,
the final results file, and the history file.
Figure 3-6
3.4.6 Data
You can probe and extract data from the image using this menu.
Figure 3-7
3.4.7 Help
Have any difficulties? Found any bugs? Any comments? Contact us! Your suggestions and
critics are very important to us
Figure 3-8
4 Run Simulations
4.1 Introduction
One main function of the CCHE2D-GUI is to run simulations using CCHE2D numerical
model. Basically, there are five steps to run simulations. The first step in getting a case ready
for simulation is to generate a mesh for the domain of interest. A mesh file, a file with a
.geo extension, is required by both the CCHE2D-GUI and the CCHE2D numerical model.
You need the help of the CCHE2D Mesh Generator. The mesh file contains the geometry and
the initial conditions of the grid, such as the (x,y) coordinates, the initial bed elevation, the
initial water surface level, the bed roughness and the nodal type. In this manual, it is
assumed that you already have such a mesh file.
After loading the mesh file, you can set the initial conditions for flow, including editing the
grid (refining, coarsening, and modifying nodal type), setting the initial bed elevation, the
initial water surface level and the initial bed roughness.
After generating and editing the mesh, the next step is to set the model parameters. The
parameters are divided into two groups: flow parameters and sediment parameters. In flow
parameters, you need to set the simulation parameters, the bed roughness and the general
physical properties of water. In the sediment parameters, you need to define the sediment size
classes, set the transport mode, the transport equations or formulas, and the sediment
simulation options and parameters, etc.
Once the mesh generation process is complete and parameters for flow and sediment are set,
you need to specify proper initial and boundary conditions. If sediment transport option is
desired, you must specify initial composition of the bed. You can specify multiple inlet and
outlet sections for the domain, prescribe specific discharge, total discharge, or discharge
hydrograph at the inlet sections. At outlet sections, you can prescribe a constant water surface
level, a open boundary condition, a rating curve, or a stage hydrograph. In addition, you can
prescribe either steady or unsteady sediment boundary condition at the inlets.
After all the above steps are finished, the last step is to run simulation. You can run the
CCHE2D model locally, that is on the same computer you are running the CCHE2D-GUI.
21
Figure 4-1
To Add or Delete J Line,
Click
or
to enter the editing status. The I lines of the mesh will become gray
while the J lines will become highlighted. The cursor will become
or
.
Click the desired place (actually somewhere in between two existing J lines) to add a
J line.
Click the J line you want to delete.
22
Figure 4-2
To Add or Delete I Line,
Click
or
to enter the editing status. The J lines of the mesh will become gray
while the I lines will become highlighted. The cursor will become
or
.
Click the desired place (actually somewhere in between two existing I lines) to add an
I line.
Click a I line you want to delete.
23
Figure 4-3
. You cannot undo the changes after you saved the changes. If you want to save the
Note: adding or deleting mesh lines change the mesh size and the orderings of the mesh
nodes. As a result, the boundary conditions and the results associated with the old mesh will
be incompatible with the new mesh. And, moving mesh lines or nodes may degenerate the
mesh quality and hence the accuracy of the solutions. Therefore, you must be very careful
with the grid editing.
24
conditions must be specified before running the simulations. The initial conditions are very
important, since inappropriate initial conditions may slow down the convergence process or
even cause the simulation to fail. The initial conditions include the flow initial conditions and
the sediment initial conditions. This section will tell you how to set the flow initial
conditions.
In addition to the geometry information, the mesh file also contains the information of the
initial flow conditions that include the initial water surface level, bed roughness, and initial
bed elevation. To set the initial conditions for flow, you need to visit Set Flow Initial
Conditionsin Simulation menu which activates the Nodal Properties dialog window.
To Specify Initial Water Surface Level, Roughness and Bed Elevation,
Select the flow variable you want to specify by checking the radio button in Assign
WSL, Roughness or Bed Elevation.
If you click Whole Domain, the selected flow variable will be assigned a constant
value through the domain.
You can also give a spatial distribution to the selected flow variable by clicking SubRegion. After clicking this button, the GUI waits for you to define the region by
selecting two points which defines a rectangular area within the mesh.
The Assign Value dialog will appear right after clicking Whole Domain or after you
click the second point when defining Sub-Region. Enter the desired value and click
OK.
You can repeat the above steps until the initial conditions of all the flow variables are
set.
Note: Please make sure the initial water level at the inlets is higher than bed.
Tip: If the initial WSL is too low, there will be many dry nodes at the beginning
of the simulation. This may cause the simulation fail. To avoid it, specify the
initial WSL a little higher to reduce the number of dry nodes.
25
Figure 4-4
You can specify a linear distribution in either I direction or J direction of water
surface level in a region represented by a pair of diagonal points instead of a constant
distribution. To do the linear interpolation, check the option Interpolate initial water
surface level and then click two different points to define a region. Finally choose the
direction by checking Interpolate along constant I-lines or Interpolation along
constant J-lines and click OK.
Figure 4-5
26
Three nodal types, namely, internal nodes, boundary nodes, and external nodes, are
identified in CCHE2D model. You can change the nodal type for a specified region
represented by a pair of diagonal points. For example, if you want to construct an island in
the domain, follow the following steps:
Click the radio button Assign new ID to nodes and then select a region by clicking
on two different nodes.
Check the option Set as boundary and click OK. The region will be surrounded by
the blue lines.
Click two diagonal points to specify the region covering all the other nodes within the
blue lines.
Warning: You must be very careful when changing the nodal type. It may result
in an ill-defined mesh and the model may fail to run. Each time when you want
to construct an island, make sure follow the above correct steps. You should
always use the CCHE2D Mesh Generator to create high quality mesh.
27
Figure 4-6
28
o Simulation time (s): the total time period in seconds of the simulation. In
case of steady flow, this time should be sufficient so that a steady solution
could be achieved
o Time step (s): the step for time marching of the simulation.
In group Time Step for Output, you can set the step intervals to output results into
the Intermediate File, the History File, and for the Monitor Points, and set the step
intervals for screen outputs indicating Convergence process.
o For Parabolic Eddy Viscosity Model and Mixing Length Model, you can
set the Turbulence viscosity coefficient. This coefficient serves as a
multiplier, i.e., a value of 10 means that the turbulent viscosity is 10 times that
computed from the selected turbulence model. Although normally it equals 1,
this coefficient has been tested in a range of [0.1, 1000].
In group Unsteady Flow Computation, you can choose to Compute as quasisteady flow instead of the real unsteady flow. If the flow boundary condition at any
of the inlet section is unsteady, i.e., discharge hydrograph is specified, both the
sediment and flow time steps are the same and the flow and sediment transport
simulations are performed for each time step. However, for long-term simulation the
computation time may be quite long. To increase the computation efficiency, the user
has the choice to turn on the option that computes the unsteady flow as quasi-steady
flow. In this case, the discharge hydrograph provided by the user must be a step
function. The discharge during each step is considered as constant and the flow and
sediment simulations are performed as for steady flow. Suppose (Qi, Ti) and (Qi+1,
Ti+1) represents two consecutive hydrograph ordinates then for quasi-steady
computation the discharge Qi is assumed to be valid for the duration between Ti and
Ti+1. In case of unsteady flow computation, the discharge is interpolated for any time
between Ti and Ti+1. You can set the Time steps to reach steady state during each
step. Note that the time series in the hydrograph files of the inlets and outlets
must be the same. The quasi-steady simulation cannot start from rest and must
run from some base result.
In group Numerical,
o Wall slipness coefficient: It is used to specify the wall boundary condition at
no-flow boundaries. A value of 0.0 means no slip condition and a value of 1.0
means total slip, i.e., tangential velocity at no flow boundaries is allowed. A
value between 0.0 and 1.0 would mean partial slip. A value greater than 1.0
29
Figure 4-7
In page Bed Roughness, the parameters are divided into two groups: For Flow Simulation
Only and For Sediment Transport Simulation.
30
For flow simulation, you can choose to either Use Values in .geo File or Use Bed
Roughness Formula. The roughness value in Geo file can be either Mannings n or
Roughness height Ks. You must specify if the roughness values in the .geo file are
Mannings n or roughness height.
o Two bed roughness formulas are available for calculation: Wu and Wang
(1999); and van Rijn (1986). If you choose roughness formula for flow, the
selected formula will be also used in the bed roughness calculation for
sediment. You also need to specify the formula parameters, such as D16, D50,
D90 and Calibration Factor. The Calibration Factor is within the range of
[0.2, 5.0] and its default value is 1.0
For sediment transport simulation, there are five methods to calculate the bed
roughness. You can Use the value in *.geo file, or calculate bed roughness according
to the sediment diameter size of D90 or D50, or use Wu and Wang (1999)s formula
and van Rijn (1986)s formula.
The options Wu and Wang (1999) and van Rijn (1986) can be set in either group
For Flow Simulation Only or For Sediment Transport Simulation. Setting them
in one group equals to setting them in both groups at the same time.
31
Figure 4-8
In page Advanced, there are four parameters, namely, Coriolis force coefficient,
gravitational acceleration, von Karman constant, and kinematic viscosity of fluid, with
default values that suffice for most cases, however, if needed you can change the default
values.
32
Figure 4-9
Define the sediment size classes. All bed material and boundary condition data are
specified according to these classes.
Sediment parameters are specified through the Set Sediment Parameters in the
Simulations page. There are four pages, Sediment Size Classes, Sediment Transport,
Sediment, and Bed Roughness.
In page Sediment Size Classes,
33
Number of Bed Layers: is used to define the bed configurations in vertical direction.
The default value is 3. The first top layer is the mixing layer where the exchanges
between sediments in water and on bed occur.
Figure 4-10
In page Sediment Transport, four groups of parameters need to be specified.
34
In the first group, you can set the Transport Mode and the Transport Capacity
Formula. There are five transport modes and four capacity formulas available.
o Total Load as Bed Load Plus Suspended Load Model is for simulation of
both the bed load and the suspended load (total load).
o Total Load as Bed Load Model is for simulation of both the bed load and the
suspended load but with the bed load dominant.
o Total Load as Suspended Load Model is for simulation of both the bed load
and the suspended load but with the suspended load dominant.
o If you select Total Load as Bed Load Model or Total Load as Suspended
Load Model, you need to choose appropriate sediment transport capacity
formula.
In group Sediment Simulation Mode, two simulation modes are available. The first
option Slow Bed Change can be selected only if the flow is steady, while the option
Fast Bed Change is for unsteady flow. Note that if you choose to Compute as
quasi-steady flow in page Simulation Parameters of Set Flow Parameters (please
refer to section 4.4 for details), it implies the simulation mode of Fast Bed Change for
sediment transport and you will not be able to set other simulation mode.
In group Adaptation Length for Bedload, you can choose appropriate option to
calculate adaptation length for bed load for non-equilibrium sediment transport.
In group Adaptation Length Factor for Suspended Load, similarly, you can choose
appropriate option to calculate the adaptation length factor for suspended.
35
Figure 4-11
In page Sediment,
The Sediment specific gravity has a default value of 2.65 that suffices for most
applications.
The Curvature Effects can be included into the sediment transport simulation if the
domain has the curved parts. You need to set the Average channel width
accordingly.
36
Figure 4-12
The page Bed Roughness is the same as the page Bed Roughness in Set Flow Parameters.
Please refer to section 4.4 for details.
37
Figure 4-13
To edit these two files, select Set Sediment Boundary Conditionsin menu
Simulations. Please make sure you already defined the sediment size class, otherwise
an error message will appear.
In the File Editor, you can create (New), Import and Save the sediment boundary
conditions file selected from the File Selector.
38
Figure 4-14
Figure 4-15
To create a template for the selected file, click New. The sbc and bbc template is
created based on the Number of Size Classes and Number of Data Points. The
Number of Size Classes is a predefined value and you cannot set it here. The
Number of Data Points defines the number of time series points. The sbc and bbc
39
files have the same format. The unit of Sediment Discharge in sbc is kg / m 3 while
in bbc is kg / s .
To edit the cell value, click the active cell and then type the desired values. The cells
in gray color are inactive. You can Cut, Copy, Paste and Delete by right clicking the
selected cells to invoke the popup edit menu.
Figure 4-16
Note that once the Number of Data Points is set, you cannot add more data points by
directly inputting values into cells. The sum of the factions of the size classes at each
time series should be 1, other wise an error message will appear.
Figure 4-17
40
There are two kinds of boundary node strings, inlet boundary node string (denoted by an
arrow entering the domain), outlet boundary node string (denoted by an arrow going out of
the domain).
Figure 4-18
Figure 4-19
41
For Inlet Boundary Conditions, you need to set the flow inlet boundary conditions and the
sediment inlet boundary conditions.
In page Flow,
o You can specify either Total discharge or Discharge hydrograph with
Inflow angle (degree) at the inlet. If the Discharge hydrograph is selected,
you need to provide a discharge hydrograph file (*.dhg) whose file format can
be found in Appendix. Note that the hydrograph must start from t = 0.
o A shortcut of Set Flow Parametersis provided here. For details, please
refer to section 4.4.
o To detach a flow boundary condition, click Delete Flow BC.
Figure 4-20
In page Sediment, you need to provide a bed load boundary condition file (*.bbc) or
a suspended load boundary condition file (*.sbc) or both, which depends on the
option you selected in the Transport Mode of page Sediment Transport in Set
Sediment Parameters.
o To attach an existing sediment boundary condition, click Import BC File
in group BedLoad Sediment or Suspended Sediment. If no existing
sediment boundary condition is available, you can create a new sediment
boundary condition by clicking Create BC File.and the File Editor will
appear. For details on how to edit sediment boundary conditions using File
Editor, please refer to section 4.6.
42
Figure 4-21
o To detach a sediment boundary condition, click Delete Sediment BC.
For Outlet Boundary Conditions, four kinds of outlet flow boundary conditions are
available, Open boundary conditions, Water surface level, Rating curve and Stage
Hydrograph.
In case of Open Boundary Conditions, the model will estimate the water surface
level at the outlet boundary based on kinematic wave condition. This option is
valuable if you cannot assess the water level at the channel outlet.
If you select Rating curve, you need to provide a rating curve file (*.rcv) with stagedischarge relationship. Please see the Appendix for instructions.
If you select Stage hydrograph, you need to provide a stage hydrograph file (*.shg)
with time-stage relationship. Note that the stage hydrograph must start from t = 0.
43
Figure 4-22
44
If the bmt file exists, the defined bed samples will be displayed, otherwise a blank
sheet is displayed.
You need to set the Porosity and the fractions of each sample. You dont need to set
the Sample No.
Figure 4-23
To create a new sample, click Add Sample. A sample with default porosity and
equal-distributed factions will be added at the end of the records.
To edit a sample, click the corresponding cell and type the desired value. Please make
sure the sum of the factions of each sample be 1, otherwise an error message will
appear.
45
To save the samples, click OK. The defined bed samples will be saved into the bmt
file and a bed file which defines the initial bed material composition for each mesh
node will be created at the same time using the information from the bmt file. If the
bed file is already existing, you need to confirm if it is overwritten or not.
Figure 4-24
Figure 4-25
There are five properties divided into two groups for each mesh node, Layer Thickness, Bed
Sample Number, Erodibility, Maximum Erosion Thickness, and Maximum Deposition
Thickness.
46
Figure 4-26
In group Layer properties,
You must set the layer thickness and layer sample number for each layer. The Layer
Number is set in page Sediment Size Classes in Set Sediment Parameters (please
refer to section 4.5). The samples are set in Define Bed Material Samples (please
refer to section 4.8).
To set the layer thickness or layer sample number, type the desired value or select
desired sample, then click Whole Domain or Define Distribution. If you click
Whole Domain, the selected property will be assigned a constant value for each
node. You can also give a spatial distribution to the selected property by clicking
Define Distribution. After clicking this button, the GUI expects you to define the
region by selecting two points which are the diagonal points of a rectangular area
within the mesh.
47
Figure 4-27
Note that you need to specify the layer properties for EACH LAYER.
In group Nodal Erodibility and Thickness,
The erodibility helps differentiate the non-erodible and the erodible nodes. For
example, the nodes on the concrete hydraulic structure should be non-erodible.
The Maximum Erosion Thickness and Maximum Deposition Thickness are used
to limit the erosion and deposition process, respectively. By default, a uniform value
of 90 meter is used, that is, no limit is enforced.
For a selected property, you can set a constant distribution by clicking Whole
Domain. You can define the distribution by clicking Define Distribution. After
clicking this button, the GUI is waiting for you to define the region by selecting two
points which are the diagonal points of a rectangular area from the mesh.
To save your settings, click Save. The changes will be saved into the bed file.
To undo the changes of the previous step, click Undo.
48
To add monitor points through IJ Indices, input the I index and J index and then
click Add Point in group Monitor Point Indices.
You can also click Add Points from Grid and then click a grid point.
You can also click Select Monitor Point and then click the desired place on the
mesh. A grid point closest to this point will be added as a monitor point.
If you set the monitor points, you need to set time step for monitor points output on
page Simulation in Set Flow Parameters.
49
Figure 4-28
In Simulation Options, the current Run Number is displayed at the top of the window. This
Run Number will be used to name the corresponding results file. The flow result file has the
extension flw, while the sediment result file has the extension sdm. The name of the
result file has the following generic form:
CaseName_Run-Current Run Number (Start Run Number).flw (or sdm)
For example, if the CaseName is try, the Current Run Number is 4, and the Start Run
Number is 1, the result file will be: try_Run-4(1).flw or try_Run-4(1).sdm. The Start Run
Number represents the results from the corresponding run. If it is 0, it means the simulation
will begin from rest.
There are four simulation options.
50
Start Flow Simulation from Rest: It is also called cold start for flow simulation.
The flow simulation will begin with the initial water surface specified as in section
4.3 and the initial velocity field (velocities are set to zero by default)
Figure 4-29
Continue Flow Simulation from Flow Filed at Time: It is also called hot start for
flow simulation. The flow simulation will begin with the computed flow field at a
selected time. You can select the flow field from the corresponding flow filed
selector.
Figure 4-30
51
Start Sediment Transport using Flow Field at Time: It is also called cold start
for sediment transport simulation. The sediment transport will begin with the initial
bed defined in the mesh file (geo) and the computed flow field at a selected time. You
can select the flow field from the corresponding flow field selector.
Figure 4-31
Continue Sediment Transport from Results at Time: It is also called hot start for
sediment transport simulation. The sediment transport will begin with the computed
bed and flow field at a selected time. You can select it from the sediment results file
selector. Use this option to continue a sediment transport simulation.
Figure 4-32
52
To actually start the simulation, first select the simulation option and then click Start
Simulation. The GUI will check all the model parameters, and the initial and boundary
conditions. After the check, the simulation will be carried out in a console window.
Progress of the run is displayed in this window. You can close the CCHE2D-GUI while
the model is running. Closing the console window will interrupt the simulation. During
the simulation, you can visualize the intermediate results.
Figure 4-33
The File Status shows the status of files associated with the corresponding step.
Available denotes the associated files already exist, while Missing denotes
the non-existing files that should be created by the CCHE2D-GUI or provided by
the users.
The Action Status traces the users action on the corresponding step. Done
means the corresponding step has been successfully executed; Canceled means
the corresponding step has been executed but not completed; and Not Done
means the corresponding step has not been executed or has been skipped.
53
The Note shows the current properties of the steps. Required means the
corresponding step must be executed, while Optional means the corresponding
step may not necessarily be executed and can be skipped if not needed for the
desired simulation.
There are 9 steps, among which steps 3, 4, 5, 7, and 8 are related to the sediment
transport simulation. The properties of these steps are determined by the sediment
Transport Mode in page Sediment Transport of Set Sediment Parameters
(please refer to section 4.4 for details). If a certain transport mode is selected,
these steps become Required if you want to run sediment transport simulation.
For flow simulation, you can skip these steps (3, 4, 5, 7 and 8).
.Figure 4-34
The navigator buttons <<Back and Next>> are used to select a step. Once a step is selected,
click Execute to perform this step. Details on these 9 steps can be found in previous sections
in this chapter.
5 Visualize Results
5.1 Introduction
Visualization of modeling results is another main function of the CCHE2D-GUI. It provides
quite a few ways to visualize the results, such as flood color shading, contour lines, color
points, vector plotting, and etc.
In CCHE2D-GUI, there are three visualization modes, namely, Grid, Flow and Sediment. In
Grid mode, only the geometry file (*.geo) is loaded; in Flow mode, the geometry file and a
result file (*.mdw or *.flw or *.his) are loaded; and in Sediment mode, the geometry file and
a result file (*.mds or *.sdm or *.sed) are loaded.
55
To load it automatically, select Auto-check Flow Intermediate Result. You can set
the time interval to check the result.
Figure 5-2
For Flow Final Results File, you need to further select the computed flow field at a certain
time.
Figure 5-3
After a result file (mdw or flw) is loaded, you can select flow variables from variable
selector on the main toolbar.
Figure 5-4
56
When reloading the results files (*.mdw or *.str), you will be asked to set the Results
Settings. You can keep the flood scale settings, contour scale settings, or view coordinates
as the same as or different from the previous file loading.
Figure 5-5
Figure 5-6
57
In History File Editor, there are two groups. In the first group, you can set the plot options
of the min and max values for the current plot variable. If the third one is selected, you can
further set the properties of the current flood values.
Figure 5-7
The second group is the control panel. The buttons related to frame operation are listed in the
following table and others are explained as follows.
Flow time and Sediment time: To display the flow time and sediment time of the
current frame.
Extract Frame Data: To extract the data of the current frame and save it as the final
result file (*.flw).
Export Images: To export the frames as a series of Bitmap files with a range set in
From frame and To frame.
58
Create AVI: To create an AVI file consisting of the frames with a range set in
From frame and To frame.
Table 5-1
Controls
Notes
Controls
Notes
Current frame
Progress bar
Previous frame
Play
Next frame
Note: According to the tests, the AVI file would NOT be created successfully if
the screen resolution is higher than 1280 1024.
59
Figure 5-9
For Sediment Final Results File, you need to further select the computed sediment result at
a certain time.
Figure 5-10
After a result file (mds or sdm) is loaded, you can select sediment variables from variable
selector on the main toolbar.
Figure 5-11
If multiple sediment size classes are defined, you can select Sediment Size Class Display
in menu Visualization to view the results for each size.
60
Figure 5-12
If you load Sediment History File, you can get help from History File Editor. (For details
see section 5.2.2)
61
Figure 5-13
For Grid Probe, in addition to probe the nodal properties, you can also Modify properties
of the node.
You can select grid point either using I and J indices or using x and y coordinates.
Input the desired values in the x, y, z, roughness or wsl and then click Set.
You can Save you changes, Undo one step of previous change and Undo All
changes.
62
Figure 5-14
For Flow Probe, in addition to probe flow data, you can also extract flow data.
Click Extract Along I Line and the flow data will be extracted to a data file(*.dat)
along a constant I line set in I index and highlighted in purple color.
Click Extract Along J Line and the flow data will be extracted to a data file(*.dat)
along a constant J line set in J index and highlighted in purple color.
63
Figure 5-15
Similarly, for Sediment Probe, in addition to probe sediment data, you can also extract
sediment data.
Click Extract Along I Line and the sediment data will be extracted to a data
file(*.dat) along a constant I line set in I index and highlighted in purple color.
Click Extract Along J Line and the sediment data will be extracted to a data
file(*.dat) along a constant J line set in J index and highlighted in purple color.
64
Figure 5-16
Figure 5-17
To Extract Data From Polyline, go to menu Data and click Extract from Polyline, the
Polyline for Data Extraction will appear.
65
Figure 5-18
Click the desired place in the domain. For each segment, you need to specify the
number of points.
Figure 5-19
Figure 5-20
66
5.5 Settings
To personalize your view styles, go to menu Settings.
, Dry Area
, Title
, Frame, XY Axis
, Grid
, Boundary, Simulation
, and add
.
Text
You can select Color Scale, Gray Color or Reverse Color Map, set Background
Color and choose Plot Types (four types are available).
67
Figure 5-21
Figure 5-22
68
In Flood Properties,
Two flood options are available, Flood using color shading and Flood using color
contours. If you select the second option, you can specify the Number of contours.
For scales of color bar, you can edit the color bar by setting the new Minimum and
Maximum values and you can restore the Original Scale by clicking Reset Scale.
Figure 5-23
For scales of color bar, you can set the new Minimum and Maximum values and Reset
Scale to the Original Scale. You can also Set Contours Color
69
Figure 5-24
In Contour Labels, first select the contour line in Contour Selector and the current contour
will be displayed in Contour Value.
To add label, click Add Label and click the desired place on the corresponding
contour highlighted as purple color. You can add multiple labels for one contour line.
To set the rotation angle for current label, input the desired angle in Rotation Angle.
This angle will be applied to the currently selected contour value.
To remove label, click Remove Label and then click the label you want to delete.
To remove all labels of one contour line, click Remove All Labels of One Line.
To set font size and significant digit for all labels, input the desired values in Font
scaling factor and Significant digits and then click Apply in Label Properties (all
labels).
To set background for all labels, check or uncheck the Font background to make the
label show the background color or transparent.
Note: Please make sure to add labels onto the correct contour line when adding
labels.
70
Figure 5-25
71
Figure 5-5
Set Vector Scale Options: You can choose to Scale vectors to velocity or set the
Uniform vector size. If you check the Animated velocity field, an animation effect
will be shown on the velocity field.
Set Display Parameters: You can choose not to display the whole velocity field by
specifying the Number of I lines to be skipped and the Number of J lines to be
skipped.
Set Vector Parameters: You can set the Scale factor, Line Thickness, Arrow head
length, Arrow head angle.
72
Figure 5-6
5.5.4 Text
There are two statuses for texts, adding and editing. In adding status, the editing functions,
such as Replace Text, Delete Text, and Delete All, will be disabled. They will be enabled in
editing status.
Figure 5-28
73
Click
area.
In Add/Modify Text, input the desired text in Input Text and then click OK.
or Add Text in menu Settings and then click the desired place in the plot
To Edit Texts,
and then Double-click the text you want to edit. The current text will be
Click
highlighted and displayed in Input Text.
In Add/Modify Text,
o To Replace Text: Select the text in Text Selector and input the new text in Input
Text, then click Replace Text.
o To Delete Text: Select the text in Text Selector and then click Delete Text.
o To Delete All Texts: Click Delete All.
o You can set the current text coordinates dependent or independent of world
coordinates to make it attached or detached to the real coordinates system of the
mesh.
o To Set Font Size for Current Text: Input the desired Scale factor and then click
Apply.
o To Set Font Angle for Current Text: Input the desired Angle and then click
Apply.
o To Set Font Color for Current Text: Click Color...
o To Set Font for All Texts: Select the font in Font and then click Apply.
o To Move Text: Click
and then click the text you want to move and hold it.
Then drag it to the desired place and release the mouse.
Figure 5-29
74
5.5.5 Legend
To show/hide Legend, click
and then click the legend and hold it. Then drag it to the desired
To move legend, click
place and release the mouse.
To set legend properties, click
Figure 5-30
In Legend Properties
75
Figure 5-31
5.5.6 Background
You set background color and background image.
To set background color, click Background Colorin menu Settings and select the
background color.
To set background image,
Figure 5-32
76
Figure 5-33
77
Figure 5-34
5.5.7 Layout
You can save the current visualization settings as a layout file (*.lof).
Figure 5-35
File Name
Mesh File
(*.geo)
Discharge
Hydrograph
File
(*.dhg)
Stage
Hydrograph
Curve Files
(*.shg)
File Name
Imax, Jmax
Do j = 1, Jmax
Do i = 1, Imax
x (x or east coordinate),
y (y or north coordinate),
wsl (initial water surface),
bed (initial bed elevation),
lsl (nodal type),
roughness
Enddo
Enddo
Number of Ordinates
Do i =1, Number of Ordinates
Time (second), Discharge
(m3/s)
Enddo
Monitor
Points File
(*.mon)
Rating
Curve File
(*.rcv)
Number of Ordinates
Do i =1, Number of Ordinates
Stage (m), Discharge
(m3/s)
Enddo
Data
Extract
File
Before
Extraction
(*.fex)
Number of Points
Do i = 1, Number of Points
x, y
Enddo
Number of Ordinates
Do i =1, Number of Ordinates
Time (second), Stage(m)
Enddo
References
Jia, Y.F. and Wang, Sam S.Y. (2001). CCHE2D: Two-dimensional Hydrodynamic and
Sediment Transport Model for Unsteady Open Channel Flow Over Loose Bed. NCCHE
Technical Report. NCCHE-TR-2001-01, Aug. 2001
Wu, W.M. (2001). CCHE2D Sediment Transport Model version 2.1. NCCHE Technical
Report. NCCHE-TR-2001-03, Aug. 2001
Khan, A.A. (2002). CCHE2D-GUI: Graphical User Interface for CCHE2D Model. Users
ManualVersion 2.0. NCCHE Technical Report. NCCHE-TR-2002-04, 2002.
Zhang, Yaoxin (2005). CCHE2D Mesh Generator Users ManualVersion 2.50. NCCHE
Technical Report. NCCHE-TR-2005-01, Jan 2005.
Zhang, Yaoxin (2005). CCHE2D-GUI Quick Start Guide. NCCHE Technical Report.
NCCHE-TR-2005-04, June 2005.