Mike FM HD 2D PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 108

MIKE 21 FLOW MODEL FM

Hydrodynamic Module
User Guide

MIKE by DHI 2011

Please Note
Copyright
This document refers to proprietary computer software which is protected
by copyright. All rights are reserved. Copying or other reproduction of
this manual or the related programs is prohibited without prior written
consent of DHI. For details please refer to your 'DHI Software Licence
Agreement'.
Limited Liability
The liability of DHI is limited as specified in Section III of your 'DHI
Software Licence Agreement':
'IN NO EVENT SHALL DHI OR ITS REPRESENTATIVES (AGENTS
AND SUPPLIERS) BE LIABLE FOR ANY DAMAGES WHATSOEVER INCLUDING, WITHOUT LIMITATION, SPECIAL, INDIRECT,
INCIDENTAL OR CONSEQUENTIAL DAMAGES OR DAMAGES
FOR LOSS OF BUSINESS PROFITS OR SAVINGS, BUSINESS
INTERRUPTION, LOSS OF BUSINESS INFORMATION OR OTHER
PECUNIARY LOSS ARISING OUT OF THE USE OF OR THE INABILITY TO USE THIS DHI SOFTWARE PRODUCT, EVEN IF DHI
HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
THIS LIMITATION SHALL APPLY TO CLAIMS OF PERSONAL
INJURY TO THE EXTENT PERMITTED BY LAW. SOME COUNTRIES OR STATES DO NOT ALLOW THE EXCLUSION OR LIMITATION OF LIABILITY FOR CONSEQUENTIAL, SPECIAL, INDIRECT,
INCIDENTAL DAMAGES AND, ACCORDINGLY, SOME PORTIONS
OF THESE LIMITATIONS MAY NOT APPLY TO YOU. BY YOUR
OPENING OF THIS SEALED PACKAGE OR INSTALLING OR
USING THE SOFTWARE, YOU HAVE ACCEPTED THAT THE
ABOVE LIMITATIONS OR THE MAXIMUM LEGALLY APPLICABLE SUBSET OF THESE LIMITATIONS APPLY TO YOUR PURCHASE OF THIS SOFTWARE.'
Printing History
June 2004
August 2005
November 2006
March 2007
October 2007
January 2009
September 2010

MIKE 21 Flow Model FM

CONTENTS

ABOUT THIS GUIDE . . . . . . . .


1.1 Purpose . . . . . . . . . . .
1.2 Assumed User Background
1.3 General Editor Layout . . .
1.3.1 Navigation tree . .
1.3.2 Editor window . . .
1.3.3 Validation window .
1.4 Online Help . . . . . . . . .

. . . . . . . . . . . . . . . . . . . . . . . . . . . 9
. . . . . . . . . . . . . . . . . . . . . . . . . . . 9
. . . . . . . . . . . . . . . . . . . . . . . . . . . 9
. . . . . . . . . . . . . . . . . . . . . . . . . . . 9
. . . . . . . . . . . . . . . . . . . . . . . . . . . 9
. . . . . . . . . . . . . . . . . . . . . . . . . . 10
. . . . . . . . . . . . . . . . . . . . . . . . . . 10
. . . . . . . . . . . . . . . . . . . . . . . . . .

10

INTRODUCTION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.1 Short Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2 Application Areas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

GETTING STARTED . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

EXAMPLES . . . . . . . . . . . . . . . . . . . . . . . .
4.1 General . . . . . . . . . . . . . . . . . . . . . . .
4.2 Wind Set-up in a Rectangular Lake . . . . . . .
4.2.1 Purpose of the Example . . . . . . . .
4.2.2 Defining the Problem . . . . . . . . . .
4.2.3 Analytical Solution . . . . . . . . . . . .
4.2.4 Presenting and Evaluating the Results
4.2.5 List of Data and Specification Files . .
4.3 Retention Basin by River . . . . . . . . . . . . .
4.3.1 Purpose of the Example . . . . . . . .
4.3.2 Defining the Hydrodynamic Model . .
4.3.3 Presenting and Evaluating the Results
4.3.4 List of Data and Specification Files . .
4.4 Flow Modelling in the Sound, Denmark . . . .
4.4.1 Purpose of the Example . . . . . . . .

BASIC PARAMETERS . . . . . . . .
5.1 Domain . . . . . . . . . . . . .
5.1.1 Mesh and bathymetry
5.1.2 Domain specification
5.1.3 Boundary names . .
5.1.4 GIS background . . .
5.2 Time . . . . . . . . . . . . . .
5.2.1 Remarks and hints .
5.3 Module Selection . . . . . . .

. . . . . . . . . . . . . . .

17

. . . . . . . . . . . . . . .

17

. . . . . . . . . . . . . . .

17

.
.
.
.
.

.
.
.
.
.

17
18
18
19
20

. . . . . . . . . . . . . . .

20

.
.
.
.

.
.
.
.

20
21
24
25

. . . . . . . . . . . . . . .

26

. . . . . . . . . . . . . . .

26

. . . . . . . . . . . . . . . . . . . . . . . . .

27

. . . . . . . . . . . . . . . . . . . . . . . . .

27

.
.
.
.

.
.
.
.

27
27
28
29

. . . . . . . . . . . . . . . . . . . . . . . . .

29

. . . . . . . . . . . . . . . . . . . . . . . . .

30

. . . . . . . . . . . . . . . . . . . . . . . . .

30

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.

HYDRODYNAMIC MODULE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
6.1 Solution Technique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
6.1.1 CFL number . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33

MIKE 21 Flow Model FM

6.2

6.3

6.4

6.5

6.6
6.7

6.8
6.9
6.10
6.11

6.12
6.13

6.14
6.15

6.1.2 Remarks and hints . . . .


Flood and Dry . . . . . . . . . . .
6.2.1 General description . . .
6.2.2 Recommended values .
6.2.3 Remarks and hints . . . .
Density . . . . . . . . . . . . . . .
6.3.1 Recommendation . . . .
6.3.2 Remarks and hints . . . .
Eddy Viscosity . . . . . . . . . . .
6.4.1 Horizontal Eddy Viscosity
6.4.2 General description . . .
6.4.3 Recommended values .
6.4.4 Remarks and hints . . . .
Bed Resistance . . . . . . . . . .
6.5.1 General description . . .
6.5.2 Recommended values .
6.5.3 Remarks and hints . . . .
Coriolis Forcing . . . . . . . . . .
Wind Forcing . . . . . . . . . . . .
6.7.1 Wind friction . . . . . . .
6.7.2 General description . . .
6.7.3 Remarks and hints . . . .
Ice coverage . . . . . . . . . . . .
Tidal Potential . . . . . . . . . . .
6.9.1 Remarks and Hints . . .
Precipitation-Evaporation . . . . .
Wave Radiation . . . . . . . . . .
6.11.1 Recommendations . . . .
6.11.2 Remarks and hints . . . .
Sources . . . . . . . . . . . . . . .
6.12.1 Source specification . . .
Structures . . . . . . . . . . . . .
6.13.1 Weirs . . . . . . . . . . .
6.13.2 Culverts . . . . . . . . . .
6.13.3 Gates . . . . . . . . . . .
6.13.4 Piers . . . . . . . . . . . .
6.13.5 Turbines . . . . . . . . .
6.13.6 Composite structures . .
6.13.7 Remarks and Hints . . .
Initial Conditions . . . . . . . . . .
Boundary Conditions . . . . . . .

. . . . . . . . . . . . . . . . . . . . . . . 34
. . . . . . . . . . . . . . . . . . . . . . . 35
. . . . . . . . . . . . . . . . . . . . . . . 35
. . . . . . . . . . . . . . . . . . . . . . . 36
. . . . . . . . . . . . . . . . . . . . . . . 36
. . . . . . . . . . . . . . . . . . . . . . . 36
. . . . . . . . . . . . . . . . . . . . . . . 37
. . . . . . . . . . . . . . . . . . . . . . . 37
. . . . . . . . . . . . . . . . . . . . . . . 37
.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

37
38
38
38

. . . . . . . . . . . . . . . . . . . . . . . 39
. . . . . . . . . . . . . . . . . . . . . . . 39
. . . . . . . . . . . . . . . . . . . . . . . 40
. . . . . . . . . . . . . . . . . . . . . . . 40
. . . . . . . . . . . . . . . . . . . . . . . 41
. . . . . . . . . . . . . . . . . . . . . . . 41
. . . . . . . . . . . . . . . . . . . . . . . 42
. . . . . . . . . . . . . . . . . . . . . . . 43
. . . . . . . . . . . . . . . . . . . . . . . 44
. . . . . . . . . . . . . . . . . . . . . . . 44
. . . . . . . . . . . . . . . . . . . . . . . 45
. . . . . . . . . . . . . . . . . . . . . . . 46
. . . . . . . . . . . . . . . . . . . . . . . 46
. . . . . . . . . . . . . . . . . . . . . . . 48
. . . . . . . . . . . . . . . . . . . . . . . 49
. . . . . . . . . . . . . . . . . . . . . . . 49
. . . . . . . . . . . . . . . . . . . . . . . 49
. . . . . . . . . . . . . . . . . . . . . . . 50
. . . . . . . . . . . . . . . . . . . . . . . 51
.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

51
57
62
64
66
68
70

. . . . . . . . . . . . . . . . . . . . . . . 73
. . . . . . . . . . . . . . . . . . . . . . . 74

6.15.1 Boundary specification


6.15.2 General description . .
6.16 Temperature/Salinity Module .
6.17 Decoupling . . . . . . . . . . . .
6.18 Outputs . . . . . . . . . . . . . .
6.18.1 Geographical View . .
6.18.2 Output specification . .
6.18.3 Output items . . . . . .
7

TEMPERATURE/SALINITY MODULE
7.1 Equation . . . . . . . . . . . . .
7.2 Solution Technique . . . . . . .
7.2.1 Remarks and hints . .
7.3 Dispersion . . . . . . . . . . . .
7.3.1 Horizontal dispersion .
7.3.2 Recommended values
7.4 Heat Exchange . . . . . . . . .
7.5 Precipitation-Evaporation . . .
7.6 Sources . . . . . . . . . . . . .
7.6.1 Source Specification .
7.6.2 Remarks and hints . .
7.7 Initial Conditions . . . . . . . . .
7.8 Boundary Conditions . . . . . .
7.8.1 Boundary Specification

. . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . .

75
80

. . . . . . . . . . . . . . . . . . . . . . . .

81

. . . . . . . . . . . . . . . . . . . . . . . .

81

. . . . . . . . . . . . . . . . . . . . . . . .

82

. . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . .

82
82
86

. . . . . . . . . . . . . . . . . . . . . . . .

89

. . . . . . . . . . . . . . . . . . . . . . . .

89

. . . . . . . . . . . . . . . . . . . . . . . .

89

. . . . . . . . . . . . . . . . . . . . . . . .

89

. . . . . . . . . . . . . . . . . . . . . . . .

90

. . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . .

90
91

. . . . . . . . . . . . . . . . . . . . . . . .

92

. . . . . . . . . . . . . . . . . . . . . . . .

93

. . . . . . . . . . . . . . . . . . . . . . . .

95

. . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . .

95
96

. . . . . . . . . . . . . . . . . . . . . . . .

96

. . . . . . . . . . . . . . . . . . . . . . . .

97

. . . . . . . . . . . . . . . . . . . . . . . .

97

TURBULENCE MODULE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101

LIST OF REFERENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103

Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105

MIKE 21 Flow Model FM

Purpose

ABOUT THIS GUIDE

1.1

Purpose
The main purpose of this User Guide is to enable you to use, MIKE 21
Flow Model FM, Hydrodynamic Module, for applications of hydraulic
phenomena in lakes, estuaries, bays, coastal areas and seas. The User
Guide is complemented by the Online Help.

1.2

Assumed User Background


Although the hydrodynamic module has been designed carefully with
emphasis on a logical and user-friendly interface, and although the User
Guide and Online Help contains modelling procedures and a large amount
of reference material, common sense is always needed in any practical
application.
In this case, common sense means a background in coastal hydraulics
and oceanography, which is sufficient for you to be able to check whether
the results are reasonable or not. This User Guide is not intended as a substitute for a basic knowledge of the area in which you are working: mathematical modelling of hydraulic phenomena.
It is assumed that you are familiar with the basic elements of MIKE Zero:
file types and file editors, the Plot Composer, the MIKE Zero Toolbox, the
Data Viewer and the Mesh Generator. An introduction to these can be
found in the MIKE Zero\manuals folder: .\MIKE_ZERO\MzGeneric.pdf.

1.3

General Editor Layout


The MIKE 21 Flow Model FM setup editor consists of three separate
panes.

1.3.1

Navigation tree
To the left is a navigation tree, that shows the structure of the model setup
file, and is used to navigate through the separate sections of the file. By
selecting an item in this tree, the corresponding editor is shown in the central pane of the setup editor.

About This Guide

1.3.2

Editor window
The edior for the selected section is shown in the central pane. The content
of this editor is specific for the selected section, and might contain several
property pages.
For sections containing spatial data - e.g. sources, boundaries and output a geographic view showing the location of the relevant items will be available. The current navigation mode is selected in the bottom of this view, it
can be zoom in, zoom out or recenter. A context menu is available from
which the user can select to show the bathymetry or the mesh, to show the
optional GIS background layer and to show the legend. From this context
menu it is also possible to navigate to the previous and next zoom extent
and to zoom to full extent. If the context menu is opend on an item - e.g. a
source - it is also possible to jump to this items editor.
Further options may be available in the context menu depending on the
section being edited.

1.3.3

Validation window
The bottom pane of the editor shows possible validation errors, and is
dynamically updated to reflect the current status of the setup specifications.
By double-clicking on an error in this window, the editor in which this
error occurs will be selected.

1.4

Online Help
The On-line Help can be activated in several ways, depending on the
user's requirement:
!

F1-key seeking help on a specific activated dialog:


To access the help associated with a specific dialog page, press the F1key on the keyboard after opening the editor and activating the specific
property page. See Figure 1.1.

Open the On-line Help system for browsing manually after a specific help page:
Open the On-line Help system by selecting Help Topics in the main
menu bar.

10

MIKE 21 Flow Model FM

Online Help

Figure 1.1

On-line Help system for MIKE 21 Flow Model FM

11

About This Guide

12

MIKE 21 Flow Model FM

INTRODUCTION
MIKE 21 Flow Model FM is a new modelling system based on a flexible
mesh approach. The modelling system has been developed for applications within oceanographic, coastal and estuarine environments.

Figure 2.1

Mariager Estuary, Denmark. Computation mesh used in MIKE 21


Flow Model FM for studying flow circulation due to combined tide,
wind and run-off.

MIKE 21 Flow Model FM is composed of following modules:


!

Hydrodynamic Module

Transport Module

ECO Lab/Oil Spill Module

Particle Tracking Module

Mud Transport Module

Sand Transport Module

The Hydrodynamic Module is the basic computational component of the


entire MIKE 21 Flow Model FM modelling system providing the hydrodynamic basis for the Transport Module, ECO Lab/Oil Spill Module, Mud
Transport Module, Particle Tracking Module and Sand Transport Module.

13

Introduction

2.1

Short Description
The Hydrodynamic Module is based on the numerical solution of the twodimensional shallow water equations - the depth-integrated incompressible Reynolds averaged Navier-Stokes equations. Thus, the model consists
of continuity, momentum, temperature, salinity and density equations. In
the horizontal domain both Cartesian and spherical coordinates can be
used.
The spatial discretization of the primitive equations is performed using a
cell-centered finite volume method. The spatial domain is discretized by
subdivision of the continuum into non-overlapping element/cells. In the
horizontal plane an unstructured grid is used comprising of triangles or
quadrilateral element. An approximate Riemann solver is used for computation of the convective fluxes, which makes it possible to handle discontinuous solutions.
For the time integration an explicit scheme is used.

2.2

Application Areas
The application areas are generally problems where flow and transport
phenomena are important with emphasis on coastal and marine applications, where the flexibility inherited in the unstructured meshes can be utilized.

14

MIKE 21 Flow Model FM

GETTING STARTED
A comprehensive step-by-step training example covering the MIKE 21
Flow Model FM is provided with the DHI Software installation. The
example is related to the construction of a fixed link across resund, the
sound between Denmark and Sweden, where one of the main constraints
were that the environment of the Baltic Sea should remained unaffected
by the link. The training guide (PDF-format) can be found in
C:\Program Files\DHI\2011\MIKE Zero\Manuals\MIKE_21\
FlowModel_FM\HD\ MIKE_FM_HD_Step_By_Step.pdf

Figure 3.1

Computation mesh used the MIKE 21 Flow Model FM step-by-step


training example

15

Getting Started

16

MIKE 21 Flow Model FM

General

EXAMPLES

4.1

General
One of the best ways of learning how to use a modelling system such as
MIKE 21 Flow Model FM is through practice. Therefore examples are
included which you can go through yourself and which you can modify, if
you like, in order to see what happens if one or other parameter is
changed.
The specification files for the examples are included with the installation
of MIKE Zero. A directory is provided for each example. The directory
names are as follows:
!

Wind set-up example:


Examples\MIKE_21\FlowModel_FM\HD\Lake

Retention Basin by River:


Examples\MIKE_21\FlowModel_FM\HD\Structure

The Sound (resund), Denmark. From scratch to calibrated model


example:
Examples\MIKE_21\FlowModel_FM\HD\Oresund

4.2

Wind Set-up in a Rectangular Lake

4.2.1

Purpose of the Example


This example has been chosen as a fairly simple one, so that it is possible
to check the results analytically.
The problem is to determine the wind set-up in the lake. The test conditions are:
!

The lake is rectangular in shape with a length of 4,9 km in the east/west


direction, and a length of 1.9 km in the north/south direction. The lake
has a uniform depth of 10 m.

The lake is connected to the sea by a 100 m wide channel in the middle
of the shore to the west. The sea is assumed to have a constant water
level of 0.0 m. Thus, there is one open boundary.

A westerly wind of 35 m/s is blowing.

17

Examples

4.2.2

Defining the Problem


The main condition defining the problem is:
!

An unstructured mesh with 1858 elements and 999 nodes is used. The
mesh has a minimum and maximum edge length of 71 m and 173 m,
respectively. The mesh is shown in Figure 4.1.

The time step of 2 seconds is selected and the duration time of the simulation is 4 hours (7200 time steps).

The horizontal eddy viscosity type has been chosen to Smagorinsky


formulation and a constant Smagorinsky coefficient of 0.28 is applied.

The bed resistance type has been chosen to Manning number and a
constant value of 32 m1/3/s is applied.

The wind is specified as constant with wind speed 35 m/s and wind
direction 270 N.

The length of soft start interval (warm-up period) for the wind has
been chosen to 2 hours (7200 seconds) to avoid chock effects.

A constant wind friction coefficient of 0.002425 is applied.

Figure 4.1

4.2.3

18

Wind set-up in a regular lake, model layout and mesh.

Analytical Solution
A steady wind blowing over an enclosed water body generates a circulation in the water body and a set-up of the water surface. Assuming a static
equilibrium between the wind set-up and the wind shear stress, the wind
set-up can be determined from a depth integrated balance of forces.

MIKE 21 Flow Model FM

Wind Set-up in a Rectangular Lake

The wind induced surface stress, w , is given by the following empirical


relation
w = a cd uw uw

(4.1)

where a is the density of air, cd is the empirical drag coefficient of air, and
u w = (u w,v w) is the wind speed 10 m above the ground.
Assuming that the water depths are large compared to the wind set-up, the
force balance can then be written as
h
g 0 d ------ = w
L

(4.2)

where h is the wind set-up, d is the still water depth, L is length of lake,
0 is density of water and g is the gravitational acceleration. Rearranging
yields an expression for the set-up
w L
h = -------- --g 0 d

4.2.4

(4.3)

Presenting and Evaluating the Results


Time series for point data can be specified directly as model output or
extracted from 2D data files using the MIKE Zero Data Viewer.
Figure 4.2 shows a time series of surface elevation at a position on the
middle of the shore to the east (4900 m, 950 m).
Entering the parameters into eq. (4.3) and assuming that the model has
reached a steady solution results in a set-up of 0.19 m. This corresponds
closely to the final result shown in Figure 4.2.

19

Examples

Figure 4.2

4.2.5

Wind set-up in a rectangular lake. Time series of total water depth at


a position on the middle of the shore to the east (4900 m, 950 m).

List of Data and Specification Files


The following data files (included in the \HD\Lake folder) are supplied
with MIKE 21 Flow Model HD FM:
File name: lake.mesh
Description: Mesh file including the mesh and bathymetry
File name: lake.m21fm
Description: MIKE 21 Flow Model FM specification file

4.3

Retention Basin by River

4.3.1

Purpose of the Example


This example has been chosen to describe the flow between a river and a
retention basin during a period of flooding using structures. The example
is divided into two parts:
1 Discover the benefit of using a weir structure to simulate overflow
2 How to combine weir and culvert to simulate a flooding scenario

Part 1
The problem is to simulate the overflow of a river bank and subsequent
filling of a retention basin. The test conditions are:

20

MIKE 21 Flow Model FM

Retention Basin by River

The river section is 210 m long and 25 m wide. The river bed level is
constant over the area at -4 m. The retention basin is 200 m long and
200 m wide with a uniform bed level at -2 m. Between the river and the
retention basin is a 5 m wide riverbank. The flow between the river and
the retention basin may take place along a 100 long section of the river
bank, where the bathymetry level is +1 m.

The natural flow conditions in the river correspond to an overall water


surface gradient of 0.0001, corresponding to a difference in water level
between the upstream boundary and the downstream boundary of 0.02
m, resulting in uni-directional flow.

Initially the water level is -1 m in the river and the basin. Due to flooding the water level upstream and downstream then increases linearly to
about +2 m during 3 hours whereafter the water level remains constant
at the boundaries for the next 6 hours.

Part 2
The problem is to simulate the filling and emptying of a retention basin.

4.3.2

Using the same conditions as above the time frame is extended as the
water level in the river decrease linearly from +2 m to -2 m during a 6
hour period followed by constant water level at the boundaries for the
next 9 hours.

The emptying of the retention basin is carried out using natural gravity
through culverts with valves that allow for flow from the retention
basin to the river but not vice versa. The culverts have a total cross-sectional area of 0.6 m2.

Defining the Hydrodynamic Model


The main condition defining the problem is:
!

A mesh with combined triangular and quadrangular mesh elements is


applied. The basin and the river is defined by triangular elements
whereas the river bank is defined by quadrangular elements. The mesh
is shown in Figure 4.3.

An overall time step of 60 seconds is selected.


The duration of simulation 1 is 9 hours (540 overall time steps).
The duration of simulation 2 is 24 hours (1440 overall time steps).

The horizontal eddy viscosity type has been chosen to a flux based
Smagorinsky formulation wih a constant value of 0.28.

The bed resistance type has been chosen to Manning number and a
constant value of 32 m1/3/s is applied.

21

Examples

22

The water level in the river is described by a time series shown in


Figure 4.4. The initial water level in both the river and the retention
basin is -1 m.

The length of soft start interval (warm-up period) has been chosen to
180 seconds.

Figure 4.3

Filling of a retention basin, mesh resolution

Figure 4.4

Water level in river section during simulations

MIKE 21 Flow Model FM

Retention Basin by River

Part 1
To model the filling of the retention basin two different approaches are
used:
!

Defining the river bank by specifying the bathymetry level to 1 m.


The bathymetry is shown in Figure 4.5.

Defining the river bank as a broad crested weir with crest level +1 m on
a bed level of -2 m. The weir is located along the line (30.0m, 52.5m)
to (30.0m, 152.5m). The headloss factors are defined by default values.
The bathymetry is shown in Figure 4.5.

Figure 4.5

Filling of a retention basin, model layout


Upper: Bathymetry, using modified bathymetry for river bank
Lower: Bathymetry, using weir structure for river bank

23

Examples

Part 2
To model the filling and following emptying of the retention basin the
river bank is modelled by a composite structure.

4.3.3

The river bank is defined as a broad crested weir with specifications as


above.

Three identical culverts are placed along the same line section as the
weir. The culverts are defined as closed, circular culverts with a diameter of 0.5 m and positioned at an invert level of -1.5 m both upstream
and downstream. The headloss factors are defined by default values.
The culverts have valves that only allows for positive flow (out into the
river).

Presenting and Evaluating the Results


Part 1: Filling the retention basin
The filling of the basin was simulated using two approaches. Figure 4.2
shows the surface elevation in position (100m, 100m) inside the basin
along with the upstream water level in the river.

Figure 4.6

Water level at position (100m, 100m) inside basin.


Red stippled line: upstream water level in river for reference
Black solid line: water level using modified bathymetry for river bank
Blue dotted line: water level using weir structure for river bank

Figure 4.2 shows that modelling the river bank directly by the bathymetry
level gives a faster filling of the basin than the use of the weir structure.
This is due to the fact that the weir computation includes contraction and
expansion loss for the free overflow thus slowing the flow whereas the
case with modified bathymetry does not include these energy losses
explicitly.

24

MIKE 21 Flow Model FM

Retention Basin by River

Part 2 Filling and emptying retention basin


Figure 4.7 shows the surface elevation in position (100m, 100m) inside
the basin along with the upstream water level in the river.

Figure 4.7

Water level at position (100m, 100m) inside basin


Red stippled line: upstream water level in river for reference
Black solid line: water level inside basin

Figure 4.7 shows that the emptying of the retention basin through the three
culverts are much slower than the time it takes to fill the basin.
This is obviously because the free flow over the weir is critical and
through a much larger cross-section whereas the flow out of the basin
using the culverts are sub-critical with a much smaller cross-section.

4.3.4

List of Data and Specification Files


The following data files (included in the \HD\Structure folder) are supplied with MIKE 21 Flow Model HD FM:
File name: Bathy0.mesh
Description: Mesh file including mesh and bathymetry, no river bank
File name: BathyB.mesh
Description: Mesh file including mesh and bathymetry, with river bank
File name: WLBoundary.dfs0
Description: Water levels at river boundaries
File name: Sim1_Weir.m21fm
Description: MIKE 21 Flow Model FM specification file

25

Examples

File name: Sim1_Bank.m21fm


Description: MIKE 21 Flow Model FM specification file
File name: Sim2_Weir.m21fm
Description: MIKE 21 Flow Model FM specification file

4.4

Flow Modelling in the Sound, Denmark

4.4.1

Purpose of the Example


The objectives of this example is to setup a MIKE 21 Flow Model FM for
resund, Denmark, from scratch and to calibrate the model to a satisfactory level.
An easy to follow Step-by-step guide is provided in the default installation
folder:
Manuals\MIKE_21\FlowModel_FM\HD

26

MIKE 21 Flow Model FM

Domain

BASIC PARAMETERS

5.1

Domain
Providing MIKE 21 Flow Model FM with a suitable mesh is essential for
obtaining reliable results from your model. Setting up the mesh includes
selection of the appropriate area to be modelled, adequate resolution of the
bathymetry, wave, wind and flow fields under consideration and definition
of codes for open and closed boundaries. Furthermore, the resolution in
the geographical space must also be selected with respect to stability considerations.
The MIKE 21 Flow Model FM is based on flexible mesh approach.

5.1.1

Mesh and bathymetry


You generate your mesh file in the MIKE Zero Mesh Generator, which is a
tool for the generation and handling of unstructured meshes, including the
definition and editing of boundaries.
The mesh file is an ASCII file including information of the geographical
position and bathymetry for each node point in the mesh. The file also
includes information of the node-connectivity in the mesh.

5.1.2

Domain specification
Map projection
If the mesh is generated by the MIKE Zero Mesh Generator, the map projection is defined in the mesh file and is only shown for reference in the
user interface. If the map projection information is not defined in the mesh
file, you have to select the correct map projection corresponding to the
data in the mesh file.
Minimum depth cutoff
If the bathymetry level in an element is above the minimum depth cutoff
value then the minimum depth cutoff value is replacing the actual bathymetry value in the computations. Please note that the minimum depth cutoff
value may be negative as the bathymetry levels is often so in the mesh file.
If you also apply a Datum shift - the depth cutoff is relative to the corrected depths.
For instance - you have a mesh file with values between +2 and -20
meters. You then shift these to a different datum with a shift of +1 meters.

27

BASIC PARAMETERS

Your corrected bathymetry now ranges between +1 and -21 m. You can
then cutoff all depths above -2m, leaving the bathymetry used in the
model to range between -2 and -21 m.
Datum shift
You can use any convenient datum for setting up the mesh of your model.
This can be Chart Datum (CD), Lowest Astronomical Tide (LAT) or Mean
Sea Level (MSL). The actual datum is unimportant.
What is important, however, is that for each simulation you must provide
the model with the correct height of the model reference level relative to
the datum used in the setup of your bathymetry. Specifying the datum shift
does this. In this way it is possible to carry out simulations using a range
of different water levels without having to alter the mesh file.
If you do not plan to apply different water levels in different simulations it
is recommended that you set up your bathymetry with the datum that you
plan use in the simulations, thus having a datum shift of 0 m.
A datum shift of e.g. 2 m (-2 m) means the water depth is increased
(decreased) by 2 m in all node points.
Mesh decomposition
To improve the performance of the numerical scheme it is possible to
include reordering of the mesh (renumbering of the element and node
numbers). This can significantly speed up the computational time by optimizing the memory access.
Note: When reordering is applied the numbering of the nodes and elements in the output files has been changed compared to the information in
the mesh file. The information in the log file corresponds to the new ordering.
To improve the performance of the numerical scheme a domain decomposition technique is applied. You can specify the number of subdomain. If
reordering is included the reordering is applied at subdomain level.

5.1.3

28

Boundary names
When you generate your mesh using the MIKE Zero Mesh Generator you
already have defined a code value for open water boundaries. Figure 5.1
shows the definition of codes in a simple application.

MIKE 21 Flow Model FM

Time

Figure 5.1

The definition of boundary codes in a mesh is made in the Mesh


Generator

In this case three open boundaries have been detected from the mesh file
specified in the domain parameters; code 2, code 3 and code 4. In the main
Boundary names dialog you can re-name the code values to more appropriate names, see Figure 5.2.

Figure 5.2

Change of default code names (from the mesh file) to more appropriate names

5.1.4

GIS background
If an ESRI ArcMap installation is found on the local machine, the user
will have the option to select a GIS layer file (.lyr) to use as geographic
background. This file can be any layer file created by ArcGIS. The layer
will be projected to the map projection chosen for the domain.

5.2

Time
The period to be covered by the simulation is specified in this dialog. You
have to specify the simulation start date, the overall number of time steps
and the overall time step interval (in seconds). The overall discrete time
steps specified on this page are used to determine the frequency for which
output can be obtained from the different modules and to synchronize the
coupling between the different modules.

29

BASIC PARAMETERS

The simulation always starts with time step number 0 and the simulation
start date is the historical data and time corresponding to time step 0. The
simulation end date is presented for reference.

5.2.1

Remarks and hints


The various modules in MIKE 21 Flow Model FM each use an internal
time step, see Figure 5.3.

Figure 5.3

Overall time step related to internal time step

The time steps for the hydrodynamic calculations, the advection-dispersion calculations and the spectral waves calculations are dynamic and
each determined to satisfy stability criterias.
All time steps within the simulation for the various modules are syncronized at the overall discrete time step. The time step for the hydrodynamic calculations is syncronized at the time step for the
advection/dispersion calculations. For the Sand Transport module and the
ECOLab module the time step for updating the process description can be
a multipla of the overall time step.

5.3

Module Selection
The MIKE 21 Flow Model FM comprises various modules, which you can
select to be included in the simulation.
You can make your selection to explicitly include the following modules

30

Transport

ECO Lab / Oil Spill

Mud Transport

MIKE 21 Flow Model FM

Module Selection

Particle Tracking

Sand Transport

by ticking one, or more, of the check boxes.


The Hydrodynamic model is obligatory.

31

BASIC PARAMETERS

32

MIKE 21 Flow Model FM

Solution Technique

HYDRODYNAMIC MODULE
The hydrodynamic module calculates the resulting flow and distributions
of salt, temperature, subject to a variety of forcing and boundary conditions. Baroclinic effect due to salt and temperature variations are considered as subordinated to the HD module and are set up here.

6.1

Solution Technique
The simulation time and accuracy can be controlled by specifying the
order of the numerical schemes which are used in the numerical calculations. Both the scheme for time integration and for space discretization
can be specified. You can select either a lower order scheme (first order)
or a higher order scheme. The lower order scheme is faster, but less accurate. For more details on the numerical solution techniques, see the Scientific documentation.
The time integration of the shallow water equations and the transport
(advection-dispersion) equations is performed using an explicit scheme.
Due to the stability restriction using an explicit scheme the time step interval must be selected so that the CFL number (p. 33) is less than 1. A variable time step interval is used in the calculation of both the shallow water
equations and the transport equations determined so that the CFL number
is less than a critical CFL number in all computational nodes. The CFL
number for the transport equations is usually less restrictive than the CFL
number for the shallow water equations. Hence, a larger time step interval
is normally used for the transport equations than the shallow water equations. To control the time step it is also possible for the user to specify a
minimum time step and a maximum time step. The time step interval for
the transport equations is synchronized to match the overall time step and
the time step interval for the shallow water equations is synchronized to
match the time step interval for the transport equations.

6.1.1

CFL number
For the shallow water equations in Cartesian coordinates the CourantFriedrich-Lvy (CFL) number is defined as
t
t
CFL HD = ( gh + u ) ------ + ( gh + v ) -----x
y

(6.1)

where h is the total water depth, u and v are the velocity components in the
x- and y-direction, g is the gravitational acceleration, x and y are a
characteristic length scale in the x- and y-direction for an element and t

33

HYDRODYNAMIC MODULE

is the time step interval. The characteristic length scale, x and y, is


approximated by the minimum edge length for each element and the water
depth and the velocity component is evaluated at element center.
For the transport equations in Cartesian coordinates the CFL number is
defined as
t
t
CFL AD = u ------ + v -----x
y

6.1.2

(6.2)

Remarks and hints


If the important processes are dominated by convection (flow), then
higher order space discretization should be chosen. If they are dominated
by diffusion, the lower order space discretization can be sufficiently accurate. In general, the time integration method and space discretization
method should be chosen alike.
Choosing the higher order scheme for time integration will increase the
computing time by a factor of 2 compared to the lower order scheme.
Choosing the higher order scheme for space discretization will increase
the computing time by a factor of 1 to 2. Choosing both as higher order
will increase the computing time by a factor of 3-4. However, the higher
order scheme will in general produce results that are significantly more
accurate than the lower order scheme.
The stability of the numerical scheme should be secure if the CFL number
less than 1. However, the calculation of the CFL number is only an estimate. Hence, stability problems can occur using this value. In these cases
you can reduce the critical CFL number. Therefore the default value of the
critical CFL number is set to 0.8. It must be in the range from 0 to 1. Alternatively, you can reduce the maximum time step interval. Note, that setting the minimum and maximum time step interval equal to the overall
time step interval specified on the Time dialog (p. 29), the time integration
will be performed with constant time step. In this case the time step interval should be selected so the the CFL number is smaller than 1.
For both the time integration of the shallow water equations and the transport equations, the total number of time steps in the calculation and the
maximum and minimum time interval during the calculation are printed in
the log-file for the simulation. The CFL number can be saved in an output
file.

34

MIKE 21 Flow Model FM

Flood and Dry

6.2

Flood and Dry


If your model is located in an area where flooding and drying occur, you
can enable the flood and dry facility. In this case you have to specify a drying water depth, a flooding water depth and a wetting depth.
When the water depth is less than the wetting depth the problem is reformulated and only if the water depth is less than the drying depth the element/cell is removed from the simulation. The flooding depth is used to
determine when an element is flooded (i.e. reentered into the calculation).
The reformulation is made by setting the momentum fluxes to zero and
only taking the mass fluxes into consideration.
If flooding and drying is not enabled, you should specify a minimum
depth cutoff (see Section 5.1.2, Domain specification), which is less than
zero. If the total water depth becomes less than zero a blow-up is detected
and the simulation halted.

6.2.1

General description
The approach for treatment of the moving boundaries (flooding and drying fronts) problem is based on the work by Zhao et al. (1994) and Sleigh
et al. (1998). When the depths are small the problem is reformulated and
only when the depths are very small the elements/cells are removed from
the calculation. The reformulation is made by setting the momentum
fluxes to zero and only taking the mass fluxes into consideration.
The depth in each element/cell is monitored and the elements are classed
as dry, partially dry or wet. Also the element faces are monitored to identify flooded boundaries.
!

An element face is defined as flooded if the following two criteria are


satisfied. Firstly, the water depth at one side of face must be less than a
tolerance depth, hdry , and the water depth at the other side of the face
larger than a tolerance depth, hflood . Secondly, the sum of the still water
depth at the side for which the water depth is less than hdry and the surface elevation at the other side must be larger than zero.

An element is dry if the water depth is less than a tolerance depth, hdry ,
and no of the element faces are flooded boundaries. The element is
removed from the calculation.

An element is partially dry if the water depth is larger than hdry and less
than a tolerance depth, hwet , or when the depth is less than the hdry and
one of the element faces is a flooded boundary. The momentum fluxes
are set to zero and only the mass fluxes are calculated.

35

HYDRODYNAMIC MODULE

An element is wet if the water depth is greater than hwet . Both the mass
fluxes and the momentum fluxes are calculated.

Note: When an element is removed from the calculation, water is removed


from the computational domain. However, the water depths at the elements, which are dried out, are saved and then reused when the element
becomes floodded again.

6.2.2

Recommended values
The default values are: drying depth hdry = 0.005m, flooding depth hflood =
0.05m and wetting depth hwet = 0.1m.
The wetting depth, hwet , must be larger than the drying depth, hdry , and
flooding depth, hflood , must satisfy
h dry < h flood < h wet

(6.3)

6.2.3

Remarks and hints


For very small values of the tolerance depth, hwet , unrealistic high flow
velocities can occur in the simulation and give cause to stability problems.

6.3

Density
The density is assumed to be a function of salinity and temperature. If you
select barotropic mode, both temperature and salinity (TS) will be constant and the density will not be updated during the simulation. When the
density is considered to be a function of temperature and/or salinity (baroclinic mode), the transport equation for the temperature and/or salinity
must be solved. The solution of these equations is automatically invoked
and the density is updated from the TS fields.
The specification of the setup for the solution of additional transport
equations is described in the Section 7 TEMPERATURE/SALINITY
MODULE.
The density is calculated using UNESCO's standard equation of state for
sea water. These relations are applicable for temperatures ranging from
2.1 C to 40.0 C and salinities in the range from 0 to 45 PSU (Practical
Salinity Unit).
The reference temperature and salinity can be used to improve the accuracy of the density calculation, as these are subtracted from the TS fields
before the density calculation. If density is a function of temperature

36

MIKE 21 Flow Model FM

Eddy Viscosity

alone, the density is calculated using the actual temperature and the reference salinity, and vice versa for salinity.

6.3.1

Recommendation
Density gradients can be included in the hydrodynamic simulation if the
water mass is well mixed, such that the temperature and salinity and thus
the density is almost constant over the depth. The density gradient is counteracted by the water-level gradients and in estuaries density gradients
may result in set-up of the water level of several cm. In most cases the
effect of the density gradients in 2D shallow water equations is small and
hence the density terms can be neglected. If the density variation over the
depth is significant you should apply the 3-dimensional flow model.

6.3.2

Remarks and hints


If the density is included, the CPU time for the simulation is increased significantly as one or two additional transport equations are solved.

6.4

Eddy Viscosity
The decomposition of the prognostic variables into a mean quantity and a
turbulent fluctuation leads to additional stress terms in the governing
equations to account for the non-resolved processes both in time and
space. By the adoption of the eddy viscosity concept these effects are
expressed through the eddy viscosity and the gradient of the mean quantity. Thus the effective shear stresses in the momentum equations contain
the laminar stresses and the Reynolds stresses (turbulence).

6.4.1

Horizontal Eddy Viscosity


The horizontal eddy viscosity can be specified in three different ways
!

No eddy

Constant eddy formulation

Smagorinsky formulation

Selecting the constant eddy formulation you must specify the eddy coefficient (in m2/s), and selecting Smagorinsky formulation you must specify
the Smagorinsky coefficient.
Data
The format of the eddy viscosity coefficient (in m2/s), or the Smagorinsky
coefficient, can be specified as

37

HYDRODYNAMIC MODULE

Constant (in domain)

Varying in domain

For the case with values varying in domain you have to prepare a data file
containing the eddy viscosity coefficient (in m2/s) or the Smagorinsky
coefficient before you set up the hydrodynamic simulation. The file must
be a 2D unstructured data file (dfsu) or a 2D grid data file (dfs2). The area
in the data file must cover the model area. If a dfsu-file is used piecewice
constant interpolation is used to map the data. If a dfs2-file is used bilinear
interpolation is used to map the data.
Eddy parameters
Using the Smagorinsky formulation you must specify a minimum and
maximum value for the eddy viscosity. The unit is (m2/s).

6.4.2

General description
In many applications a constant eddy viscosity can be used for the horizontal stress terms. Alternatively, Smagorinsky (1963) proposed to
express sub-grid scale transports by an effective eddy viscosity related to a
characteristic length scale. The sub-grid scale eddy viscosity is given by
A = c s2 l 2 2S ij S ij

(6.4)

where cs is a constant, l is a characteristic length and the deformation rate


is given by
1 u u
S ij = --- -------i + -------j ( i, j = 1,2 )
2 x j x i

(6.5)

For more details on this formulation, the reader is referred to Lilly (1967),
Leonard (1974), Aupoix (1984), and Horiuti (1987).

6.4.3

Recommended values
The Smagorinsky coefficient, CS, should be chosen in the interval of 0.25
to 1.0.

6.4.4

Remarks and hints


In the same way as for the bed resistance you can use the eddy coefficients
to damp out numerical instability (see section 6.5 Bed Resistance). You
should only use this as a last resort to your stability problem: the schematisation of the bathymetry and the boundary conditions are usually the primary causes for a blow-up.

38

MIKE 21 Flow Model FM

Bed Resistance

When you use the Smagorinsky formulation of the turbulence the CPU
time for a simulation is increased.

6.5

Bed Resistance
The bed resistance can be specifyed in three different ways
!

No bed resistance

Chezy number

Manning number

Selecting the one of the last two options you must specify the Checy
number (in m1/2/s) and the Manning number (in m1/3/s), respectively.
Data
The format of the Chezy number, C (m1/2/s), and the Manning number, M
(m1/3/s), can be specified in two ways
!

Constant (in domain)

Varying in domain

For the case with values varying in domain you have to prepare a data file
containing the Chezy number (in m1/2/s) or the Manning number (in
m1/3/s) before you set up the hydrodynamic simulation. The file must be a
2D unstructured data file (dfsu) or a 2D grid data file (dfs2). The area in
the data file must cover the model area. If a dfsu-file is used piecewice
constant interpolation is used to map the data. If a dfs2-file is used bilinear
interpolation is used to map the data.

6.5.1

General description
The bottom stress, b , is determined by a quadratic friction law

----b- = c f u b u b
0

(6.6)

where cf is the drag coefficient, u b is the flow velocity above the bottom
and 0 is the density of the water.

39

HYDRODYNAMIC MODULE

For two-dimensional calculations u b is the depth-average velocity and the


drag coefficient can be determined from the Chezy number, C, or Manning number, M
g
c f = -----2C

(6.7)

g
c f = --------------------( Mh 1 / 6 ) 2

(6.8)

where h is the total water depth and g is the gravitational acceleration. The
units of Chezy numbers and Manning numbers are (m1/2/s) and (m1/3/s),
respectively.
Please note that the relation between the Manning number and the bed
roughness length, ks , can be estimated using the following
25.4
M = --------k s1 / 6

(6.9)

Also note that the Manning number used here is the reciprocal value of the
Mannings n described in some textbooks.

6.5.2

Recommended values
If the relative variation of the water depth is considerable you should specify Manning numbers. Values in the range 20-40 m1/3/s are normally used
with a suggested value of 32 m1/3/s if no other information is available.
If you use Chezy numbers values in the range 30 - 50 m1/3/s are common
for coastal and marine applications.
Note that due to of the definition of the resistance numbers a smal resistance number corresponds to a high bed resistance and vice versa.

6.5.3

40

Remarks and hints


If the boundary conditions at one of your boundaries are inaccurate and
you therefore have stability problems (blow-ups) at this boundary, you can
specify a small area along the boundary with a very high resistance in
order to remove noise generated by the boundary. Manning numbers in the
range 5 -10 m1/3/s or roughness heights in the order of 1 m have been
applied successfully. However, this method should only be used if it is
impossible to improve the boundary conditions. Furthermore, the simula-

MIKE 21 Flow Model FM

Coriolis Forcing

tion results in the area which high bed resistance should be used with caution.

6.6

Coriolis Forcing
The effect of the coriolis force can be included in three different ways
!

No Coriolis force

Constant in domain

Varying in domain

If the constant in domain option is selected the Coriolis force will be calculated using a constant specified reference latitude (in Degree).
If the varying in domain option is selected the Coriolis force will be calculated based on the geografical information given in the mesh file.

6.7

Wind Forcing
It is possible to take into account the effect of the wind on the flow field.
The format of the wind data can be specified as
!

Constant. The wind is blowing from the same direction and with the
same magnitude for the whole simulation period and over the entire
model area.

Varying in time and constant in domain. The magnitude and direction


varies during the simulation period but is the same over the entire
model area.

Varrying in time and domain. The magnitude and direction varies during the simulation period and over the model area.

NOTE: the directions are given in degrees and measured clockwise from
true north to where the wind is blowing from.
Data
For the case with wind varying in time but constant in domain you have to
prepare a data file containing the wind speed (in m/s) and direction (in
degrees from true North) before you set up the hydrodynamic simulation.
The data file must be a time series file (dfs0). This can be done by entering
the data in an ASCII file using your normal editor and then reading this
file into the standard data file format using the MIKE Zero Time Series

41

HYDRODYNAMIC MODULE

Editor. The wind data must cover the complete simulation period. The
time step of the wind input data file does not, however, have to be the
same as the time step of the hydrodynamic simulation. A linear interpolation will be applied if the time steps differ.
For the case with wind varying both in time and domain you have to prepare a data file containing the wind speed components (in m/s) and air
pressure (in hPa) before you set up the hydrodynamic simulation. The file
must be a 2D unstructured data file (dfsu) or a 2D grid data file (dfs2).
This can be done by either using one of the two MIKE 21 wind generating
programs (cyclone generated wind and pressure, or wind generated on the
basis of digitised pressure fields). Or you can enter the data in an ASCII
file using your normal editor and then reading this file into the standard
data file format using the MIKE Zero Grid Series Editor. The area in the
data file must cover the model area. If a dfsu-file is used piecewice constant interpolation is used to map the data. If a dfs2-file is used bilinear
interpolation is used to map the data. The wind data must cover the complete simulation period. The time step of the wind input data file does not,
however, have to be the same as the time step of the hydrodynamic simulation. A linear interpolation will be applied if the time steps differ.
Neutral pressure
For the case with wind varying in both time and domain you have to specify a reference or neutral pressure level (in hPa).
Soft start interval
You can specify a soft start interval (in seconds) during which the wind
speed is increased from 0 to the specified wind speed in order to avoid
shock waves being generated in the model. By default the soft start interval is zero (no soft start). The wind direction is not altered by the soft start.

6.7.1

42

Wind friction
The wind friction can be specified either as a constant or varying with the
wind speed. In the latter case, the friction is linearly interpolated between
two values based on the wind speed and if the wind speed is below the
lower limit or above the upper limit the friction is given the value corresponding to that limit.

MIKE 21 Flow Model FM

Wind Forcing

6.7.2

General description
You can include the effects of a wind blowing over the model area in the
following way. The surface stress, s , is determined by the wind speed
above the water. The stress is given by the following empirical relation
s = a c d u w u w

(6.10)

where a is the density of air, cd is the empirical drag coefficient of air,


and u w = ( u w, v w ) is the wind speed 10 m above the sea surface.
Please notice that the direction of the wind is given in degrees blowing
from (relative to true north (see Figure 6.1).

Figure 6.1

Definition of wind direction

The drag coefficient can either be a constant value or depend on the wind
speed. The empirical formulae proposed by Wu (1980, 1994) is used for
the parametrization of the drag coefficient.
ca

cb c a
c d = c a + ----------------- ( w 10 w a )
w

b wa
cb

w 10 < w a
w a w 10 w b

(6.11)

w 10 > w b

where ca, cb, wa and wb are empirical factors and w10 is the wind speed 10
m above the sea surface.
The default values for the empirical factors are ca = 1.25510-3, cb =
2.42510-3, wa = 7 m/s and wb = 25 m/s. The value gives generally good
results for open sea applications. Field measurements of the drag coeffi-

43

HYDRODYNAMIC MODULE

cient collected over lakes indicate that the drag coefficient is larger than
open ocean data. For a detailed description of the drag coefficient see Geenaert and Plant (1990).

6.7.3

Remarks and hints


You can use the drag coefficient as a parameter in your model calibration.

6.8

Ice coverage
It is possible to take into account the effects of ice coverage on the flow
field. The ice coverage can be specified in four different ways
!

No coverage

Specified ice concentration

Specified ice thickness

Specified ice concentration and thickness

For the two cases where the ice concentration (the area coverage ratio of
the ice) is specified the effect of the ice is taken into account in areas
where the ice concentration becomes larger than a user-specified critical
concentration (default value 0.9). For the case where only the ice thickness is specified the effect of the ice is taken into account in areas where
the ice thickness is larger than zero.
If included in the simulation, maps of the ice concentration and/or thickness (in m) must be specified.
In areas where the sea is covered by ice the wind stress is excluded that is
the wind speed is set to zero. Additionally, the effect of the ice roughness
on the flow field can be included. For this case you must specify the
roughness height (in m).
Ice coverage data
You have to prepare a data file containing the ice concentration and/ot
thickness before you set up the hydrodynamic simulation. The file must be
a 2D unstructured data file (dfsu) or a 2D grid data file (dfs2). The data
file must cover the model area. If a dfsu-file is used piecewice constant
interpolation is used to map the data and if a dfs2-file is used bilinear
interpolation is used. The data must cover the complete simulation period.
The time step of the input data file does not, however, have to be the same
as the time step of the hydrodynamic simulation. A linear interpolation
will be applied if the time steps differ.

44

MIKE 21 Flow Model FM

Tidal Potential

Critical concentration
In areas where the ice concentration becomes larger than the specified
value of the critical concentration the effect on the flow field is taken into
account.
Roughness data
The format of the ice roughness height, ks (m), can be specified as
!

Constant (in domain)

Varying in domain

For the case with values varying in domain you have to prepare a data file
containing the drag coefficient or the roughness height (in m) before you
set up the hydrodynamic simulation. The file must be a 2D unstructured
data file (dfsu) or a 2D grid data file (dfs2). The data file must cover the
model area. If a dfsu-file is used piecewice constant interpolation is used
to map the data and if a dfs2-file is used bilinear interpolation is used.

6.9

Tidal Potential
The tidal potential is a force, generated by the variations in gravity due to
the relative motion of the earth, the moon and the sun. The forcing acts
throughout the computational domain. The forcing is considered as the
sum of several harmonic terms, each representing a specific constituent,
that are specified separately.
The format of the tidal potential data can be specified as
!

Specified from dialog

Specified from file

Data
The tidal potential is defined by the number of constituents that should be
included and each constituent is described by a number of parameters, as
shown in Table 6.1, where the parameters for the M2 tidal constituent are
shown as an example. The default is 11 constituents comprising M2, O1,
S2, K2, N2, K1, P1, Q1, Mm, Mf and Ssa. There is no limit on the number of

45

HYDRODYNAMIC MODULE

constituents that can be defined. The values for other constituents can be
found in standard tidal books as for example (Pugh, 1987).
Table 6.1

Specification of the M2 tidal constituent

Name

M2

Species

Constituent

Amplitude

0.242334

Earthtide

0.693

Period scaling

3600

Multiplies Period to obtain seconds

Period

12.42

Period of the constituent

Nodal no. 1

1.0

Nodal modulation factor

Nodal no. 2

-0.037

Nodal modulation factor

Nodal no. 3

-2.1

Nodal factor, angle

Argument 1

Doodson number, i1

Argument 2

Doodson number, i2

Argument 3

Doodson number, i3

Argument 4

Doodson number, i4

Argument 5

Doodson number, i5

Phase

Phase angle correction

Semi-diurnal, i0=2,
Diurnal, i0=1,
Seasonal, i0=0

Factor based on Love numbers

6.9.1

Remarks and Hints


The tidal potential is a weak force and is usually significant in deep
enclosed water bodies only as e.g. Mediterranean Sea or Gulf of Mexico
or in large scale models of e.g. the Pacific Ocean. In model setups with
open boundaries, these usually dominates the flow.

6.10

Precipitation-Evaporation
In applications where the rain fall is important for the flow, you can
include precipitation in your simulation. The precipitation can be included
in three different ways

46

No precipitation

Specified precipitation

Net precipitation

MIKE 21 Flow Model FM

Precipitation-Evaporation

The net precipitation is the precipitation minus the evaporation. Thus


evaporation corresponds to negative net precipitation. Selecting the No
precipitation or the Specified precipitation option evaporation can also be
included seperately in three different ways
!

No evaporation

Specified evaporation

Computed evaporation

For the case with computed evaporation the evaporation rate is calculated
as part of the latent heat flux. This option can only be selected if both density variations and heat exchange are included.
Selecting the specified precipitation (or net precipitation) option you must
specify the precipitation rate in mm/day. Selecting the specified evaporation option you must specify the evaporation rate in mm/day.
Data
The format of the precipitation rate (or net precipitation rate) and evaporation rate (in mm/day) can be specified as
!

Constant (in time and domain)

Varying in time, constant in domain

Varying in time and domain

If specified precipitation/evaporation is selected the precipitation/evaporation rates must be positive. If net precipitation is selected and the precipitation rate is negative it will act as evaporation.
For the case with precipitation rate varying in time but constant in domain
you have to prepare a data file containing the precipitation rate (in
mm/day) before you set up the hydrodynamic simulation. The data file
must be a time series file (dfs0). The data must cover the complete simulation period. The time step of the input data file does not, however, have to
be the same as the time step of the hydrodynamic simulation. A linear
interpolation will be applied if the time steps differ.
For the case with precipitation rate varying both in time and domain you
have to prepare a data file containing the precipitation rate (in mm/day)
before you set up the hydrodynamic simulation. The file must be a 2D
unstructured data file (dfsu) or a 2D grid data file (dfs2). The area in the
data file must cover the model area. If a dfsu-file is used piecewice constant interpolation is used to map the data. If a dfs2-file is used bilinear

47

HYDRODYNAMIC MODULE

interpolation is used to map the data. The data must cover the complete
simulation period. The time step of the input data file does not, however,
have to be the same as the time step of the hydrodynamic simulation. A
linear interpolation will be applied if the time steps differ.
For the case with evaporation rate varying in time but constant in domain
and for the case with evaporation rate varying both in time and domain the
data file needed must be prepared in the same way as described above for
the precipitation input.
Soft start interval
You can specify a soft start interval (in seconds) during which the precipitation/evaporation rate is increased linearly from 0 to the specified values
of the precipitation/evaporation rate. By default the soft start interval is
zero (no soft start).

6.11

Wave Radiation
The second order stresses due to breaking of short period waves can be
included in the simulation. The radiation stresses act as driving forces for
the mean flow and can be used to calculate wave induced flow.
If included in the simulation, a data file with the the three components of
the radiation stress divided by the density of the water, Sxx, Syy and Sxy
(m3/s2), must be specified.
The data file containing the wave radiation stresses can be generated by
the wave models MIKE 21 SW, MIKE 21 NSW or MIKE 21 PMS.
Data
You have to prepare a data file containing the three components of the
wave radiation stress divided by the density of water (in m3/s2) before you
set up the hydrodynamic simulation. The file must be a 2D unstructured
data file (dfsu) or a 2D grid data file (dfs2). The area in the data file must
cover the model area. If a dfsu-file is used piecewice constant interpolation is used to map the data. If a dfs2-file is used bilinear interpolation is
used to map the data. The data must cover the complete simulation period.
The time step of the input data file does not, however, have to be the same
as the time step of the hydrodynamic simulation. A linear interpolation
will be applied if the time steps differ.

48

MIKE 21 Flow Model FM

Sources

Soft start interval


You can specify a soft start interval (in seconds) during which the wave
radiation stresses are increased linearly from 0 to the specified values of
the radiation stresses in order to avoid shock waves being generated in the
model. By default the soft start interval is zero (no soft start).

6.11.1

Recommendations
In order not to generate a chock in the computations it is recommended
to apply these stresses in combination with a soft-start (warm-up) period.

6.11.2

Remarks and hints


Since the wave radiation stresses describe the average flow over one wave
period, the stresses are connected to a certain water depth. Application of
the wave radiation stresses in simulations with time varying water depths
(e.g. tide and/or storm surges) is possible, but the error introduced by the
changes in water depth should be considered.
If the Flood and Dry facility is applied, the user should be sure that the
wave radiation stresses are well defined in all grid points which will be
flooded during the simulation. In general, it is not recommended to use
wave radiation stresses together with the Flood and Dry facility.

6.12

Sources
The effects of rivers, intakes and outlets from power plants etc. can be
included in the simulation using point sources.
Depending on the choice of property page you can see a geographic view
or a list view of the sources.
There are two different methods for specification of the sources. In the list
view you can create a new source by clicking on the "New source" button.
By selecting a source in the Source list and clicking on the "Delete source"
you can remove this source. For each source you can specify the name of
the source and whether the source should be active or not. The specification of detailed information for each source is made subsequently. From
the list view page you can go to the dialog for specification by clicking on
the "Go to .." button or by selecting a source from the source list and clicking on the Edit source button.
In the geographic view it is also possible to create a new source by doubleclicking a position on the map or by selecting "New source" from the context menu. The name of the source must then be specified and the geo-

49

HYDRODYNAMIC MODULE

graphical position of the source can be edited. The additional information


for the source is made subsequently.

6.12.1

Source specification
Three types of sources can be selected
!

Simple source

Standard source

Connected source

For the simple source option only the source contribution to the continuity
equation is taken into account. For this option you have to specify only the
magnitude of the source (in m3/s). If the magnitude of the source is positive, water is discharged into the ambient water and if the magnitude is
negative water is discharged out of the ambient water.
For the standard source option both the source contribution to the continuity equation and the momentum equations is taken into account. For this
option you have to specify both the magnitude of the source (in m3/s) and
the velocity (in m/s) by which the water is discharged into the ambient
water. Note that the contribution to the momentum equation is only taken
into account when the magnitude of the source is positive (water is discharge into the ambient water).
For the connected source option both the source contribution to the continuity equation and the momentum equations is taken into account. For this
option you have to specify the number of the source which the source is
connected to .The magnitude of the source is then obtained as the magnitude of the source specified for the connected source, but with opposite
sign. You also have to specify the velocity (in m/s) by which the water is
discharged into the ambient water. Note that the contribution to the
momentum equation is only taken into account when the magnitude of the
source is positive (water is discharge into the ambient water).
Location
You must select the map projection (Longitude/Latitude, UTM etc.) in
which you want to specify the horizontal location of the source point.
You must specify the horizontal coordinates
Data
The format of the source information can be specified as
!

50

Constant in time

MIKE 21 Flow Model FM

Structures

Varying in time

For the case with source information varying in time you have to prepare a
data file containing the source information (discharge (in m3/s) and/or
velocity components (in m/s)) before you set up the hydrodynamic simulation. The data file must be a time series data file (dfs0). The data must
cover the complete simulation period. The time step of the input data file
does not, however, have to be the same as the time step of the hydrodynamic simulation. A linear interpolation will be applied if the time steps
differ.

6.13

Structures
The horizontal dimension of structures is usually much smaller than the
cell (element) sizes used in the computational grid. Therefore, the effect of
structures is modeled by a subgrid technique (the flow past a structure is
modelled by considering the upstream and downstream water level). Five
different types of structures can be included in the simulations:
!

Weirs

Culverts

Gates

Piers

Turbines

Furthermore you have the option to include Composite structures by combining two or more defined structures.
Note: A structure will only come into effect if placed in an area where
water flow would exist without the structure.

6.13.1

Weirs
Depending on the choice of property page you can see a Geographic View
or a List View of the weirs.
There are two different methods for specification of the weir. In the List
View you can create a new weir clicking on the "New weir" button. By
selecting a weir in the Weir list and clicking on the "Delete weir" button
you can remove this weir. For each weir you can specify the name of the
weir. The specification of detailed information for each weir is made subsequently. From the List View page you can go to the dialog for specifica-

51

HYDRODYNAMIC MODULE

tion by clicking on the "Go to .." button or by selecting a weir in the Weir
list and clicking on "Edit weir" button.
In the Geographic View it is also possible to create a new weir by doubleclicking a position on the map. This position will be the first point defining the cross-section for the weir or by selecting New weir. The name of
the weir must then be specified and the geographical position of the weir
can then be edited. The additional information for the weir is made subsequently.
Weir specification
The location and geometrical layout of the weirs must be specified.

Figure 6.2

Setup definition of contracted weir

Location
Weirs are defined in the domain as a cross section specified as a list of
points (a minimum of two points required). The section is composed of a
sequence of line segments. The line segments are straight lines between
two successive points. The cross section in the numerical calculations is
defined as a section of element faces. The face is included in the section
when the line between the two element centers of the faces crosses one of
the line segments. The geographical coordinates are taken from the dialog
or from an ASCII file. The file format is two space separated floats (real
numbers) for the x- and y-coordinate on separate lines for each of the
points. The faces defining the line section are listed in the log-file.
You must also select the map projection (Longitude/Latitude, UTM etc.)
in which you want to specify the location of the line section for the weir.
For more information please see Location of a structure (p. 70).

52

MIKE 21 Flow Model FM

Structures

Attributes
The required formula to be applied for the specific weir-structure is
defined in the attributes. A range of formulas are available:
!

Broad Crested Weir formula

Weir formula 1

Weir formula 2 (Honma formula)

Valve
Valve regulation of the structure flow can be defined as part of the structure definition.
Four different valve regulation types are available:
!

None: No valve regulation applies (flow is not regulated).

Only Positive Flow


Only flow in positive flow direction is allowed. Valve regulation does
not allow flow in negative flow direction and the flow through the
structure will be zero in this case.

Only Negative Flow


Only flow in negative flow direction is allowed. Valve regulation does
not allow flow in positive flow direction and the flow through the
structure will be zero in this case.

No Flow
No flow is allowed in the structure. Valve regulation closes completely
the structure.

For more information please see Flow direction (p. 72).


Alpha zero

Alpha zero is the water level difference at which the discharge calculation is described by a linear variation. If the water level difference is below this value the discharge gradients are suppressed.
The default setting is 0.01 meter. If a structure shows oscillatory
behavior it is recommended to increase this value slightly.
For more information please see Alpha zero value (p. 73).
Head Loss Factors
The factors determining the energy loss occurring for flow through
hydraulic structures.

53

HYDRODYNAMIC MODULE

Following head loss factors shall be defined (for positive and negative
Flow direction):
!

Inflow (contraction loss)

Outflow (expansion loss)

Free Overflow

Head Loss Factors are applied in structure flow calculation only for broad
crested weir and special weir types.
For more information please see Head Loss Factors (p. 56) and Flow
direction (p. 72).
Geometry
The weir geometry defines the geometrical shape of the active flow area
of the weir. Geometry is defined as a Level-width relationship, where the
Level/Width table defines the Weir shape as a set of corresponding set of
levels and flow widths. Values in the levels column must be continuous,
increasing values.
For more information please see Level-width relationship (p. 71).
Levels are defined relative to the datum (starting from the crest or sill
level and up). E.g. for a horisontal weir positioned at invert level -10m
(bed level) and extending 6 m above the bed, the weir could should be
defined by L0= -4m and L1= 0m.
Datum defines an offset which is added to the level column in the
level/width table during computation. This is typically used if the weir
geometry has been surveyed with respect to a local benchmark.
NOTE: Weir geometry shall only be defined for Broad Crested Weir.
General description
Broad Crested Weir
For a broad crested weir the user describes the shape of the "hole" through
a level/width table (see Figure 6.3). The datum value for the structure may
be used to shift the levels by a constant amount. The latter is typically
used if the weir geometry has been surveyed with respect to a local benchmark.

54

MIKE 21 Flow Model FM

Structures

Figure 6.3

Definition sketch of broad crested weir geometry

The standard formulations for flow over a broad crested weir are established automatically by the program on the basis of the weir geometry and
the user specified head loss and calibration coefficients. These formulations assume a hydrostatic pressure distribution on the weir crests. Different algorithms are used for drowned flow and free overflow, with an
automatic switching between the two.
Weir formula 1
For the Weir formula 1 description the parameters are given by Figure 6.4.
The width is perpendicular to the flow direction. Typically the invert level
coincides with the overall datum.
Weir formula 1 is based on a standard weir expression, reduced according
to the Villemonte formula:
H ds H w
k
q = WC ( H us H w ) 1 -------------------- H us H w

0,385

(6.12)

where Q is discharge through the structure, W is width, C is weir coefficient, k is the weir exponential coefficient, Hus is upstream water level, Hds
is downstream water level and Hw is weir level (see Figure 6.4). The invert
level is the lowest point in the inlet or outlet section respectively.

55

HYDRODYNAMIC MODULE

Figure 6.4

Definition sketch for Weir Flow

Weir formula 2
For the Weir formula 2 the geometry is given by a crest level and a width.
The crest level is taken with respect to the global datum. The width is perpendicular to the flow direction.
Weir formula 2 is the Honma formula:

C W ( H us H w ) ( H us H w )
Q = 1
C 2 W ( H ds H w ) ( H us H ds )

for

( H ds H w ) H us < 2 3

for

( H ds H w ) H us 2 3

(6.13)

where Q is the discharge through structure, W is the width, C1 is the weir


coefficient 1, and C2, the weir coefficient 2, is calculated according to
C 2 = 1,5 3C 1 , Hus is the upstream water level, Hds is the downstream
water level and Hw is weir level (see Figure 6.4).
Head Loss Factors
The flow description generally used for a structure is given by:
2

Vs
H = t -----2g

(6.14)

where H is the energy loss over the structure, t is the total head loss
coefficient and Vs is the mean cross sectional velocity at the structure.

56

MIKE 21 Flow Model FM

Structures

The total head loss coefficient (t) is composed of entrance (1) and exit
(2) coefficients. The coefficients are generally related to the input parameters for Inflow (in) and Outflow (out) and the changes in velocity (v) and
area (A):
v1
As 2
t = 1 + 2 = in ---- + out -----
vs
A 2

(6.15)

where suffix '1' and '2' represents velocity and Area on inflow and outflow
side of structure respectively, and 's' represents the velocity and Area in
the structure itself.
However, in the present implementation, upstream and downstream cross
sections are not extracted and accordingly, tabulated relations on cross
section areas as function of water levels are not known. Instead, upstream
and downstream areas are set to a large number resulting in a full loss contribution from the head loss factors defined. Viz,
t = 1 + 2 = in + out

(6.16)

Care must be taken when selecting loss coefficients, particularly in situations where both subcritical and supercritical flow conditions occur. When
flow conditions change from subcritical to supercritical (or the Froude
number FR becomes greater than 1), the loss coefficients in and out
(specified in the Head Loss Factors box) are modified:
!

If FR > 1 in upstream h-point then in = in / 2

If FR > 1 in downstream h-point then out = out / 2

The critical flows (and orifice flows for culverts as well) are multiplied by
the critical flow correction factor, c, specified as the Free Overflow Head
Loss Factor. Typically a value of 1.0 is used.

6.13.2

Culverts
Depending on the choice of property page you can see a Geographic View
or a List View of the culverts.
There are two different methods for specification of the culverts. In the
List View you can create a new culvert clicking on the "New culvert" button. By selecting a culvert in the Culvert list and clicking on the "Delete
culvert" button you can remove this culvert. For each culvert you can
specify the name of the culvert. The specification of detailed information

57

HYDRODYNAMIC MODULE

for each culvert is made subsequently. From the List View page you can
go to the dialog for specification by clicking on the "Go to .." button or by
selecting a culvert in the Culvert list and clicking on "Edit culvert" button.
In the Geographic View it is also possible to create a new culvert by double-clicking a position on the map or by selecting New culvert from the
context menu. This position will be the first point defining the cross-section for the culvert. The name of the culvert must then be specified and the
geographical position of the culvert can then be edited. The additional
information for the culvert is made subsequently.
Culvert specification
The location and geometrical layout of the culverts must be specified.

Figure 6.5

Setup definition of culvert

Location
Culverts are defined in the domain as a cross section specified as a list of
points (a minimum of two points required). The section is composed of a
sequence of line segments. The line segments are straight lines between
two successive points. The cross section in the numerical calculations is
defined as a section of element faces. The face is included in the section

58

MIKE 21 Flow Model FM

Structures

when the line between the two element centers of the faces crosses one of
the line segments. The geographical coordinates are taken from the dialog
or from an ASCII file. The file format is two space separated floats (real
numbers) for the x- and y-coordinate on separate lines for each of the
points. The faces defining the line section are listed in the log-file.
You must also select the map projection (Longitude/Latitude, UTM etc.)
in which you want to specify the location of the line section for the culvert.
For more information please see Location of a structure (p. 70).
Attributes
A number of parameters define the culvert characteristics:
!

Upstream Invert
Invert level at inflow location of the culvert.

Downstr. Invert
Invert level at outflow location of the culvert.

Length
Length of the culvert.

Manning's n
Manning's bed resistance number along the culvert (for friction loss
contribution). n = 1/M (Manning number)

No. of Culverts
'No. of Culverts' is a number identifying how many culverts exist at the
specific culvert location with identical geometrical definition.
An example; Five identical shaped draining pipes are placed just next
to each other in an earth dam, and in order not to make 5 individual
culvert definitions - one for each pipe - the 'No. of Culverts' in this case
can be defined as 5 and the simulation engine will recognize that 5 culverts of identical shape and size are located here and flow calculations
will take this into account accordingly.

Alpha zero
Alpha zero is the water level difference at which the discharge calculation is described by a linear variation. If the water level difference is
below this value the discharge gradients are suppressed.
The default setting is 0.01 meter. If a structure shows oscillatory
behavior it is recommended to increase this value slightly.
For more information please see Alpha zero value (p. 73).

59

HYDRODYNAMIC MODULE

Valve
Valve regulation of the structure flow can be defined as part of the structure definition.
Four different valve regulation types are available:
!

None: No valve regulation applies (flow is not regulated).

Only Positive Flow


Only flow in positive flow direction is allowed. Valve regulation does
not allow flow in negative flow direction and the flow through the
structure will be zero in this case.

Only Negative Flow


Only flow in negative flow direction is allowed. Valve regulation does
not allow flow in positive flow direction and the flow through the
structure will be zero in this case.

No Flow
No flow is allowed in the structure. Valve regulation closes completely
the structure.

For more information please see Flow direction (p. 72).


Section Type
A culvert structure can be defined as either a Closed or an Open section
type.
For more information please see Closed / Open Section switch (p. 61).
Head Loss Factors
The factors determining the energy loss occurring for flow through
hydraulic structures.
Following head loss factors shall be defined (for positive and negative
flow directions):
!

Inflow (contraction loss)

Outflow (expansion loss)

Free Overflow

Bends

For more information please see Head Loss Factors (p. 61).

60

MIKE 21 Flow Model FM

Structures

Geometry
The culvert geometry defines the geometrical shape of the active flow area
of the culvert. The cross sectional geometry of a culvert can be specified
as:
!

Rectangular
The width and height specify the geometry.

Circular
The geometry is specified by the diameter.

Irregular Level-Width Table


The geometry is specified using a level/width table. Values in the level
column must be increasing.

The Level/Width table defines the Weir shape as a set of corresponding


levels and flow widths. Values in the levels column must be continuous,
increasing values.
General description
Closed / Open Section switch
A culvert structure can be modelled as an open section if required by setting the Closed / Open switch. An example where this may be used is a
"long" weir where the friction along the length of the weir is of importance and/or the flow areas at the entrance and exit are significantly different.
If set to open the culvert will never run full or partially full, therefore only
those flow conditions which represent a free water surface are modelled.
When the water level is higher than the soffit the hydraulic parameters are
calculated based on a section extended vertically upwards with a width
equal to that at the soffit. For example, in the case of a rectangular section
the height value is essentially redundant as the cross-section will be modeled as an open section of constant width.
In the case of a circular section, this switch is invalid and will be set to
closed.
Head Loss Factors
The total head loss, H through a culvert is given by:
2
Q 1 f + b 2
- + --------------- + -----H loss = ------ -----2
2
2g A 2
A
A
s
s
s
1

(6.17)

61

HYDRODYNAMIC MODULE

where As is the mean cross section area along the length of the culvert and
Q is the discharge, 1 is the entrance or contraction loss, 2 is the outlet or
expansion loss, f is the Friction loss calculated using the Manning formula and b is the bend loss coefficient.
As
1 = in 1 -----

A 1

(6.18)

A 2
2 = out 1 -----s

A 2

(6.19)

Cross section areas, A1 and A2 are not processed and extracted in the
present implementation and hence, defined as an infinite value. Contraction and Expansion losses are therefore assumed to be applied in full using
the defined inflow and outflow loss coefficients, in and out.
2

2gLn
f = -------------43
R

(6.20)

where L is the culvert length, n is Manning's coefficient and R is the mean


Hydraulic Radius along the culvert. The Manning's n-value depends on
the interior surface of the culvert. Table values can be found in litterature.
For e.g. a concrete culvert n would typically range from 0.011 to 0.017 .
he bend loss coefficient, b, is provided for situations where head losses
other than from the above occur, for example bends, damaged culverts,
trapped debris. For straight culverts in good condition a value of zero
would apply.
The critical flows (and orifice flows for culverts as well) are multiplied by
the critical flow correction factor, c, specified as the Free Overflow Head
Loss Factor. Typically a value of 1.0 is used.

6.13.3

Gates
Depending on the choice of property page you can see a Geographic View
or a List View of the gates.
There are two different methods for specification of the gates. In the List
View you can create a new gate clicking on the "New gate" button. By
selecting a gate in the Gate list and clicking on the "Delete gate" button
you can remove this gate. For each gate you can specify the name of the
gate and whether the gate should be active or not. The specification of

62

MIKE 21 Flow Model FM

Structures

detailed information for each gate is made subsequently. From the List
View page you can navigate to the dialog for specification by clicking on
the "Go to .." button or by selecting a gate in the Gate list and clicking on
"Edit gate" button.
In the Geographic View it is also possible to create a new gate by doubleclicking on a position on the map or by selecting "New gate" from the
context menu. This position will be the first point defining the line-section
for the gate. The name of the gate must then be specified and the geographical position of the gate can then be edited. The additional information for the gate is entered subsequently.
Gate specification
The location of the gate and the control information for operation of the
gate is specified.
Location
Gates are defined in the domain as a line section specified as list of points
(a minimum of two points is required). The section is composed of the
sequence of line segments. The line segments are straight lines between
two successive points. The cross section in the numerical calculations is
defined as a section of the nearest element faces. The face is included in
the section when the line between the two element centers of the faces
crosses one of the line segments. The geographical coordinates are taken
from the dialog or from an ASCII file. The file format is two space separated floats (real numbers) for the x- and y-coordinate on separate lines for
each of the points. The faces defining the line section are listed in the logfile.
You must also select the map projection (Longitude/Latitude, UTM etc.)
in which you want to specify the location of the line section for the gate.
For more information please see Location of a structure (p. 70).
Gate Data
The operation of the gate is determined by a control factor. If the factor is
1 the gate is open and if the factor is 0 the gate is closed. When the factor
is in the interval between 0 and 1 the gate is partly closed. The format of
the control factor can be specified as
!

Constant in time

Varying in time

63

HYDRODYNAMIC MODULE

For the case with the control factor varying in time you have to prepare a
data file containing the control factor before you set up the hydrodynamic
simulation. The data file must be a time series data file (dfs0). The data
must cover the complete simulation period. The time step of the input data
file does not, however, have to be the same as the time step of the hydrodynamic simulation. A linear interpolation will be applied if the time steps
differ.

6.13.4

Piers
The effect of bridge piers is modeled as sub-grid structures using a simple
drag-law to capture the increasing resistance imposed by the piers as the
flow speed increases.
Depending on the choice of property page you can see a Geographic View
or a List View of the piers.
There are two different methods for specification of the piers. In the List
View you can create a new pier by clicking on the "New pier" button. By
selecting a pier in the Pier list and clicking on the "Delete pier" button you
can remove this pier. For each pier you can specify the name of the pier
and whether the pier should be active or not. The specification of detailed
information for each pier is made subsequently. From the List View page
you can navigate to the dialog for specification by clicking on the "Go to
.." button or by selecting a pier in the Pier list and clicking on the "Edit
pier" button.
In the Geographic View it is also possible to create a new pier by doubleclicking on a position on the map or by selecting "New pier" from the context menu. The name of the pier must then be specified and the geographical position of the pier can then be edited. The additional information for
the pier is made subsequently.
Pier specification
The location and geometrical layout of the piers must be specified.
Location
You specify the horizontal coordinates in which the pier should be placed.
You must also select the map projection (Longitude/Latitude, UTM etc.)
in which you want to specify the horizontal location of the pier point.
Pier Data
First you specify the angle of the pier, the streamline factor and the
number of vertical pier sections. Then you specify the geometry.

64

MIKE 21 Flow Model FM

Structures

The angle is measured from projection north to the alignment, positive


clockwise and in degrees, see Figure 6.6. Note that projection north is not
necessary the same as geographical north.

Figure 6.6

Definition of pier angle

The streamline factor is a factor that is multiplied on the total drag force to
take into account the increased flow velocity due to the blocking of piers.
A typical value is 1.02.
The number of sections means the number of pier segments, i.e. the
number of vertical sections with different geometrical layout.
Geometry
Additionally, you specify the geometry for each section. You have to specify the type of pier. The type can be one of three:
!

Circular

Rectangular

Elliptical

For each section you specify a number of parameters depending on the


type of section. For a circular section the height (in m) and the width
(diameter) (in m) must also be specified. For a rectangular section the
height (in m), the width (in m), the length (in m), and the radius of the corners (in m) must be specified. Finally, for an elliptical section the height
(in m), the width (in m) and the length (in m) must be specified. The sec-

65

HYDRODYNAMIC MODULE

tions are ordered so that the first section is at the top of the pier structure
and the last section is at the bottom of the pier structure (at the sea bed),
see Figure 6.7.

Figure 6.7

Definition of pier sections

General description
The effect on the flow due to the piers is modeled by calculating the current induced drag force on each individual pier.
The effective drag force, F, is determined from
1
2
F = --- w C D A e V
2

(6.21)

where w is the density of water, is the streamline factor, CD is the drag


coefficient, Ae is the effective area of pier exposed to current and V is the
current speed. The sign of F is such that a positive force acts against the
current direction.

6.13.5

Turbines
The effect of tidal turbines is modelled as sub-grid structures using a simple drag-law to capture the increasing resistance imposed by the turbine
blades as the flow speed increases. Turbines are assumed always to have
their axis aligned with the flow direction.
Depending on the choice of property page you can see a Geographic View
or a list view of the turbines. On the property page output it is possible to
specify if turbine information should be written to an output file.

66

MIKE 21 Flow Model FM

Structures

There are two different methods for specification of the turbines. In the
List View you can create a new turbine by clicking on the "New turbine"
button. By selecting a turbine in the Turbine list and clicking on the
"Delete turbine" you can remove this turbine. For each turbine you can
specify the name of the turbine and whether the turbine should be active or
not. The specification of detailed information for each turbine is made
subsequently. From the List View page you can navigate to the dialog for
specification by clicking on the "Go to .." button or by selecting a turbine
in the Turbine list and clicking on the "Edit turbine" button.
In the Geographic View it is also possible to create a new turbine by double-clicking a position on the map or by selecting "New turbine" from the
context menu. The name of the turbine must then be specified and the geographical position of the turbine can then be edited. The additional information for the turbine is made subsequently.
Turbine specification
The location of the turbine and the turbine data is specified.
Location
You specify the horizontal coordinates in which the turbine should be
placed.
You must also select the map projection (Longitude/Latitude, UTM etc.)
in which you want to specify the horizontal location of the turbine point.
Turbine Data
You must specify the diameter of the turbine (in m) and the z-coordinate
(in m) of the centroid of the turbine. It is also necessary to specify a constant drag coefficient and a correction factor. The correction factor is multiplied on the drag.
The format of the correction factor can be specified as
!

Constant in time

Varying in time

For the case with correction factor varying in time you have to prepare a
data file containing the correction factor before you set up the hydrodynamic simulation. The data file must be a time series data file (dfs0). The
data must cover the complete simulation period. The time step of the input
data file does not, however, have to be the same as the time step of the
hydrodynamic simulation. A linear interpolation will be applied if the time
steps differ.

67

HYDRODYNAMIC MODULE

Output
If output is included you must specify the name of the output file. The file
will, for each turbine, contain the drag force and the average current speed
used in the calculation of the drag force.
General description
The effect to the flow due to the turbines is modeled by calculating the
current induced drag force on each individual.
The effective drag force, F, is determined from
1
2
F = --- w C D A e V
2

(6.22)

where w is the density of water, is a correction factor, CD is the drag


coefficient, Ae is the effective area of turbine exposed to current and V is
the current speed. The sign of F is such that a positive force acts against
the current direction.
For three-dimensional calculations the current speed used to determine the
force is the average current speed over the water column covered by the
turbine and the shear stress is equally distributed over water column covered by the turbine.

6.13.6

Composite structures
Composite structures can be defined by combining two or more defined
structures. A composite structure can be composed of a combination of
Weirs and/or Culverts.
An example of a composite structure could be a bridge with multiple
waterways. Such a structure can be described by a number of culverts,
each defining an individual waterway. Additionally, for a potential bridge
deck overtopping a weir can be included to describe such overflow.
A set of structures forming a composite structure are recognized by the
program from the location definitions. Locations must be completely identical for all the structures forming the composite structure. That is, the
table of coordinates defining the structure locations must be exactly identical (number of coordinates and coordinate values) for all structures
defined.
Examples of composite structures are given below.

68

MIKE 21 Flow Model FM

Structures

Example 1 - A wide weir with a small opening


Consider a weir as illustrated in Figure 6.8.

Figure 6.8

Wide weir with a small opening

There are a number of possibilities when modeling this in MIKE 21:


1 One weir
2 Two weirs with width W1 and width (W2-W1), respectively
3 Three weirs with widths W1 and two with width (W2-W1)/2,
respectively
Using the first approach is only appropriate if the weir can be contained
within a single grid cell. The second approach may be used if the weir
spans multiple cells, keeping in mind that the flow over the highest crest
(L2) is uniformly distributed over all the affected cells. The third approach
will give the best representation of the flow. Note that the location needs
to be defined for each of the segments for case 2 and 3.
Example 2 - A wide weir with multiple culverts
Consider a structure as illustrated in Figure 6.9.

69

HYDRODYNAMIC MODULE

Figure 6.9

Wide weir with multive culverts

The composite structure should be implemented as four separate structures:


1 A weir with a constant crest level L1 and a location defined by the full
extent of the weir
2 A circular culvert
3 A rectangular culvert
4 A irregular culvert described by a level/width table
Note that the location needs to be defined for each of the four structure
components separately. The location line should correspond to the maximum width of the structure component while still obeying the minimum
requirement with respect to intersecting a line segments connecting cell
centers.

6.13.7

Remarks and Hints


Location of a structure
The location of a structure may be given by a number of geo-referenced
points which together make up a poly-line. The poly-line defines the width
of the structure perpendicular to the flow direction.

70

For a weir and a gate the poly-line will describe the location of the
crest in the horizontal plan view.

For a culvert the poly-line should be given as a line perpendicular to


the flow direction and further the line should intersect the mid point of
the culvert.

MIKE 21 Flow Model FM

Structures

Figure 6.10

The location of a weir. Note the affected cell faces.

Figure 6.11

The location of a culvert. Note the affected cell face.

A cell face is affected if the defining poly-line intersects the line segment
connecting the mid points of the two adjacent grid cells. The flow through
the structure is evenly distributed along the affected cell faces.
Level-width relationship
The geometry of a structure may be defined by a level-width relation,
where the Level/Width table defines the shape of the active flow area as a
set of corresponding levels and flow widths. Values in the levels column
must be continuous, increasing values.

71

HYDRODYNAMIC MODULE

Figure 6.12

Definition sketch of level-width relationship


Upper: broad crested weir
Lower: Irregular culvert

Flow direction
Positive and Negative flow directions as referred to in Valve and Head
loss coefficient definitions follows the definition as presented in
Figure 6.13. That is, positive flow direction is defined as the left perpendicular to the direction of a structure line schematization.

72

MIKE 21 Flow Model FM

Initial Conditions

Figure 6.13

Positive and Negative flow direction definition at structures

Alpha zero value


When the water level gradient across a structure is small the corresponding gradient of the discharge with respect to the water levels is large. This
in turn may result in a very rapid flow response to minor changes in the
water level upstream and downstream.
As a way of controlling this effect an Alpha zero value has been introduced. The Alpha zero value defines the water level difference below
which the discharge gradients are suppressed. The default setting is 0.01
meter. If a structure shows oscillatory behavior it is recommended to
increase this value slightly.

6.14

Initial Conditions
The initial values for the hydrodynamic variables can be specified in three
different ways
!

Constant.

Spatial varying surface elevation.

Spatially varying water depth and velocities

For the last case the initial conditions can be the result from a previous
simulation in which case the initial conditions effectively act as a hot start
of the flow field.

73

HYDRODYNAMIC MODULE

To avoid generation of shock waves, it is recommended that your initial


surface elevation roughly matches the boundary conditions at the start of
the simulation.
Data
For the case with spatially varying surface elevation you have to prepare a
data file containing the surface elevation (in m) before you set up the
hydrodynamic simulation. The file must be a 2D unstructured data file
(dfsu) or a 2D grid data file (dfs2). The area in the data file must cover the
model area. If a dfsu-file is used piecewice constant interpolation is used
to map the data. If a dfs2-file is used bilinear interpolation is used to map
the data. In case the input data file contains a single time step, this field is
used. In case the file contains several time steps, e.g. from the results of a
previous simulation, the actual starting time of the simulation is used to
interpolate the field in time. Therefore the starting time must be between
the start and end time of the file.
For the case with spatially varying surface elevation and velocities you
have to prepare one data file before you set up the hydrodynamic simulation. The file must contain the total water depth (in m) and the velocity
components in the x- and y-directions (in m/s). The file must be a 2D
unstructured data file (dfsu) or a 2D grid data file (dfs2). The area in the
data file must cover the model area. If a dfsu-file is used piecewice constant interpolation is used to map the data. If a dfs2-file is used bilinear
interpolation is used to map the data. In case the input data file contains a
single time step, this field is used. In case the file contains several time
steps, e.g. from the results of a previous simulation, the actual starting
time of the simulation is used to interpolate the field in time. Therefore the
starting time must be between the start and end time of the file.

6.15

Boundary Conditions
Initially, the set-up editor scans the mesh file for boundary codes (sections), and displays the recognized codes and suggest a default name for
each. You can re-name these names to more meningful names in the
Domain dialog (see Boundary names (p. 28)).
Depending on the choice of property page you can see a geographic view
or a list view of the boundaries.
The specification of boundary information for each code (section) is made
subsequently. From the list view you can go to the dialog for specification
by clicking on the "Go to .." button.

74

MIKE 21 Flow Model FM

Boundary Conditions

6.15.1

Boundary specification
For hydrodynamic boundaries there are seven types
!

Land (zero normal velocity)

Land (zero velocity)

Specified velocities

Specified fluxes

Specified level

Specified discharge

Flather condition

It should be noted that each boundary section must contain at least 2


nodes.
For stationary solid boundaries two types of boundary conditions can be
applied. For the "Land (normal velocity)" the full slip boundary conditions
is assumed to hold, that is the normal velocity component is zero, while
for the "Land (zero velocity)" the no slip condition is assumed to hold, that
is both the normal and tangential velocity components are zero. By default
the code value for land boundaries automatically identified by the Mesh
Generator is set to one. For these boundaries the full slip boundary condition are applied. If the no slip boundary conditions should be applied for
these boundaries the code values should change to a value larger than one.
If you select the specified discharge option, the discharge will be distributed as it would have been in a uniform flow field with the Manning
resistance law applied, i.e. is relative to h5/3, where h is the depth. This distribution is, in most cases, a good approximation.
The Flather (1976) condition is one of the most efficient open boundary
condition. It is very efficient in connection with downscaling coarse
model simulations to local areas (see Oddo and Pinardi (2007)). The instabilities which are often observed when imposing stratified density at a
water level boundary can be avoided using Flather conditions.
Data
Velocity boundary
If a velocity boundary is selected the velocities in the x- and y-direction
(in m/s) can be specified in three different ways
!

Constant (in time and along boundary)

75

HYDRODYNAMIC MODULE

Variable in time and constant along boundary

Variable in time and along boundary

For the case with velocities varying in time but constant along the boundary you have to prepare a data file containing the velocity components in
the x- and y-direction (in m/s) before you set up the hydrodynamic simulation. The data file must be a time series file (dfs0). The data must cover
the complete simulation period. The time step of the input data file does
not, however, have to be the same as the time step of the hydrodynamic
simulation. You can choose between different types of interpolation.
For the case with velocities varying both in time and along the boundary
you have to prepare a data file containing the boundary values before you
set up the hydrodynamic simulation. The file must be a profile file (dfs1).
The mapping from the input data file to the boundary section is described
in Interpolation type. The data must cover the complete simulation period.
The time step of the input data file does not, however, have to be the same
as the time step of the hydrodynamic simulation. You can choose between
different types of interpolation.
Flux boundary
If a flux boundary is selected the fluxes (the depth-integrated velocities) in
the x- and y-direction (in m2/s) can be specified in three different ways
!

Constant (in time and along boundary)

Variable in time and constant along boundary

Variable in time and along boundary

For the case with fluxes varying in time but constant along the boundary
you have to prepare a data file containing the flux components in the xand y-direction (in m2/s) before you set up the hydrodynamic simulation.
The data file must be a time series file (dfs0). The data must cover the
complete simulation period. The time step of the input data file does not,
however, have to be the same as the time step of the hydrodynamic simulation. You can choose between different type of interpolation.
For the case with fluxes varying both in time and along the boundary you
have to prepare a data file containing the boundary values before you set
up the hydrodynamic simulation. The file must be a profile file (dfs1). The
mapping from the input data file to the boundary section is described in
Interpolation type. The data must cover the complete simulation period.
The time step of the input data file does not, however, have to be the same

76

MIKE 21 Flow Model FM

Boundary Conditions

as the time step of the hydrodynamic simulation. You can choose between
different type of interpolation.
Level boundary
If a level boundary is selected the format of the water level (surface elevation, in m) can be specified as
!

Constant (in time and along boundary)

Variable in time and constant along boundary

Variable in time and along boundary

For the case with water level varying in time but constant along the
boundary you have to prepare a data file containing the water level (in m)
before you set up the hydrodynamic simulation. The data file must be a
time series file (dfs0). The data must cover the complete simulation
period. The time step of the input data file does not, however, have to be
the same as the time step of the hydrodynamic simulation. You can choose
between different type of interpolation.
For the case with water level varying both in time and along the boundary
you have to prepare a data file containing the boundary values before you
set up the hydrodynamic simulation. The file must be a profile file (dfs1).
The mapping from the input data file to the boundary section is described
in Interpolation type. The data must cover the complete simulation period.
The time step of the input data file does not, however, have to be the same
as the time step of the hydrodynamic simulation. You can choose between
different type of interpolation.
Discharge boundary
If a discharge boundary is selected the format of the total discharge (in
m3/s) can be specified as:
!

Constant (in time)

Variable in time

For the case with discharge varying in time you have to prepare a data file
containing the total discharge (in m3/s) before you set up the hydrodynamic simulation. The data file must be a time series file (dfs0). The data
must cover the complete simulation period. The time step of the input data
file does not, however, have to be the same as the time step of the hydrodynamic simulation. You can choose between different types of time interpolation.

77

HYDRODYNAMIC MODULE

Flather boundary
If "Flather condition" is selected the velocities in the x-and y-direction (in
m/s) and the water level (surface elevation, in m) must be specified. For
specification of these parameters see "Velocity boundary" and "Level
boundary". If discharge constrain is included the total discharge (in m3/s)
must also be specified. For specification of the total discharge see Discharge boundary.
Interpolation type
For the two cases with values varying in time two types of time interpolation can be selected:
!

linear

piecewise cubic

In the case with values varying along the boundary two methods of mapping from the input data file to the boundary section are available:
!

normal

reverse order

Using normal interpolation the first and last point of the line are mapped
to the first and the last node along the boundary section and the intermediate boundary values are found by linear interpolation. Using reverse order
interpolation the last and first point of the line are mapped to the first and
the last node along the boundary section and the intermediate boundary
values is found by linear interpolation.
Soft start
You can specify a soft start interval (in seconds) during which boundary
values are increased from a specified reference value to the specified
boundary value in order to avoid shock waves being generated in the
model. The increase can either be linear or follow a sinusoidal curve.
Discharge constrain
To ensure that the discharge through the boundary is correct when using a
Flather boundary the specified velocities can be corrected. The correction
is performed by adding a uniform correction to the normal velocity components.
Coriolis/wind correction
Using level boundary conditions in the presence of Coriolis and wind
force, may in some cases gives unrealistic results. If you specify a con-

78

MIKE 21 Flow Model FM

Boundary Conditions

stant water level variation (horizontal water level) you will get a large
inflow together with a large outflow at the same boundary, especially in a
steady state situation, as the water level should actually be tilted. Also if
you specify a variation of the water level along the boundary, the data very
often do not include the effect of the local wind, e.g. if the level variation
is obtained from tidal predictions. For these cases you can include tilting
to improve the boundary conditions.
If tilting is included a correction of the water level in each point along the
boundary is calculated based on the steady state Navier Stokes equations.
The correction takes into account the effect of Coriolis and wind force.
This non-linear approach provides the best estimate in cases where you
have a gentle slope on the seabed. If the seabed is very jaggered at the
boundary, you may have to smooth the bathymetry.
The difficulty in using the tilt facility lies in the specification of the tilting
point. Presently, the tilting point is selected as the deepest point along the
boundary. If the still water depth along the line is constant, the centre point
is selected.
If you have two adjacent boundaries you should also be careful not to create a conflicting situation in the corner.
By default Coriolis/wind correction is not included. Note that coriolis/wind correction is not possible for boundaries with less than 3 node
points.
Pressure correction
When level boundaries is selected in combination with spatially varying
wind, the level boundaries can optionally be adjusted to take into account
the variation of the pressure.
By default pressure correction is included.
Radiation stress correction
When level or discharge boundaries is selected in combination with wave
radiation the boundaries values can be optionally adjusted to take into
account the effect of wave radiation.
For level boundaries the correction is assumed to be zero at the deepest
point along the boundary.
By default radiation stress correction is not included.

79

HYDRODYNAMIC MODULE

6.15.2

General description
Convention of direction of boundary
When specifying a line series at the boundary it is important to know how
the model defines the first and last node of the boundary.
The rule is: follow the shoreline with the discretized domain on the left
hand side, see example in Figure 6.14. When a boundary is reached, this is
the first node of the boundary.

Figure 6.14

Convention of boundary start and end points.


Boundaries are defined by codes in the mesh file.

Convention of direction of discharge


By definition, discharge is positive for flow towards left when positioned
at the first point and looking forward along the boundary line. I.e. referring to the Convention of direction of boundary positive discharge values
will describe flow into the domain whereas negative discharge values will
desribe flow out of the domain.

80

MIKE 21 Flow Model FM

Temperature/Salinity Module

Correction of water level


The correction of the water level due to spatially varying air pressure is
adjusted according to
air neutral
= * ------------------------------g

(6.23)

where air is the air pressure, neutral is the neutral pressure, g is the gravitational acceleration, is the density of the water, is the specified surface elevation and is the corrected surface elevation applied as boundary
condition.

6.16

Temperature/Salinity Module
The Temperature/Salinity (TS) module is invoked from the specification
of the density, provided baroclinic density (density depends on temperature and/or salinity) is selected (see Section 6.3 Density).
A more detailed description of the Temperature/salinity module can be
found in section 7 TEMPERATURE/SALINITY MODULE.

6.17

Decoupling
In many applications a number of simulations are performed, where the
model parameters for the add-on modules are changed, but where the
model parameters for the basic hydrodynamic simulations are unchanged.
If the basic flow information from a hydrodynamic simulation is saved in
a number of special decoupling output files, it is possible to re-run the
simulation in a decoupled mode, where instead of calculating the flow
information, it is read from these files. Running a simulation in decoupled
mode can significantly reduce the computational time.
If you choose to save decoupling data you have to specify two output file
names. You must also specifed the frequency (time step factor) by which
data is going to be saved. Finally, you must give the name of the specification file for the decoupled simulation. Setting up and run a simulation
using this specification file you can not change the basic parameters in the
time and domain dialog and the parameters in the hydrodynamic model,
but you are only allowed to change the parameters in the add-on modules.

81

HYDRODYNAMIC MODULE

Note, that in decoupled mode the time step in not restricted by the CFL
stability condition for solution of flow equations, but only the stability
restriction for solution of the transport equations.

6.18

Outputs
Standard data files with computed results from the simulation can be specified here. Because result files tend to become large, it is normally not
possible to save the computed discrete data in the whole area and at all
time steps. In practice, sub areas and subsets must be selected.
In the main Outputs dialog you can add a new output file by clicking on
the "New output" button. By selecting a file in the Output list and clicking
on the "Delete output" you can remove this file. For each output file you
can specify the name (title) of the file and whether the output file should
be included or not. The specification of the individual output files is made
subsequently. You can go to the dialog for specification by clicking on the
"Go to .." button. Finally, you can view the results using the relevant
MIKE Zero viewing/editing tool by clicking on the "View" button during
and after the simulation.

6.18.1

Geographical View
This dialog shows the geographical position of the output data.

6.18.2

Output specification
For each selected output file the field type, the output format, the treatment of flood and dry, the output file (name and location) and time step
must be specified. Depending on the output format the geographical
extend of the output data must also be specified.
Field type
For a 2D simulation 2D field parameters can be selected. The mass budget
for a domain and the discharge through a cross section can also be
selected.
Output format
The possible choice of output format depends on the specified field type.
For 2D field variables the following formats can be selected

82

Point series. Selected field data in geographical defined points.

Lines series. Selected field data along geographical defined lines.

MIKE 21 Flow Model FM

Outputs

Area series. Selected field data in geographical defined areas.

If mass budget is selected for the field type you have to specify the domain
for which the mass budget should be calculated. The file type will be a
dfs0 file.
If discharge is selected for the field type you have to specify the cross section through which the discharge should be calculated. The file type will
be a dfs0 file.

Table 6.2

List of tools for viewing, editing and plotting results

Output format

File type Viewing/editing tools Plotting tools

Point series

dfs0

Time Series Editor

Plot Composer

Line series

dfs1

Profile Series Editor

Plot Composer

Area series

dfsu

Data Viewer

Data Viewer

Output file
A name and location of the output file must be specified along with the
type of data (file type).
Treatment of flood and dry
For 2D and 3D field parameters the flood and dry can be treated in three
different ways
!

Whole area

Only wet area

Only real wet area

Selecting the Only wet area option the output file will contain delete values for land points. The land points are defined as the points where the
water depth is less than a drying depth. Selecting the Only real wet area
option the output file will contain delete values for points where the water
depth is less than the wetting depth. The drying depth and the wetting
depth are specified on the Flood and Dry dialog. If flooding and drying is
not included both the flooding depth and the wetting depth are set to zero.

83

HYDRODYNAMIC MODULE

Time step
The temporal range refers to the time steps specified under Simulation
Period in the Time dialog.
Point series
You must specify the type of interpolation. You can select discrete values
or interpolated values.
The geographical coordinates of the points are either taken from the dialog
or from a file. The file format is an ascii file with four space separated
items for each point on separate lines. The first two items must be floats
(real numbers) for the x- and y-coordinate. For 2D field data the third item
is unused (but must be specified). The last item (the remaining of the line)
is the name specification for each point.
You must also select the map projection (Long/Lat, UTM-32 etc.) in
which you want to specify the horizontal location of the points.
If "discrete values" is selected for the type of interpolation, the point values are the discrete values for the elements in which the points are located.
The element number and the coordinates of the center of the element are
listed in the log-file.
If "interpolated values" is selected for the type of interpolation, the point
values are determined by 2nd order interpolation. The element in which
the point is located is determined and the point value is obtained by linear
interpolation using the vertex (node) values for the actual element. The
vertex values are calculated using the pseudo-Laplacian procedure proposed by Holmes and Connell (1989). The element number and the coordinates are listed in the log-file.
Line series
You must specify the first and the last point on the line and the number of
discrete points on the line. The geographical coordinates are taken from
the dialog or from a file. The file format is an ascii file with three space
separated items for each of the two points on separate lines. The first two
items must be floats (real numbers) for the x- and y-coordinate. For 2D
field data the third item is unused (but must be specified). If the file contains information for more than two points (more than two lines) the information for the first two points will be used.
You must also select the map projection (Long/Lat, UTM-32 etc.) in
which you want to specify the horizontal location of the points.

84

MIKE 21 Flow Model FM

Outputs

The values for the points on the line are determined by 2nd order interpolation. The element in which the point is located is determined and the
point value is obtained by linear interpolation using the vertex (node) values for the actual element. The vertex values are calculated using the
pseudo-Laplacian procedure proposed by Holmes and Connell (1989).
The element number and the coordinates are listed in the log-file.
Note: If spherical coordinates (map projection LONG/LAT) is used for a
3D model simulation the line must be either a horizontal or a vertical line.
Area series
The discrete field data within a polygon can be selected. The closed region
is bounded by a number of line segments. You must specify the coordinates of the vertex points of the polygon. Two successive points are the
endpoints of a line that is a side of the polygon. The first and final point is
joined by a line segment that closes the polygon. The geographical coordinates of the polygon points are taken from the dialog or from a file. The
file format is an ascii file with three space separated items for each of the
two points on separate lines. The first two items must be floats (real numbers) for the x- and y-coordinate. For 2D field data the third item is unused
(but must be specified).
You must also select the map projection (Long/Lat, UTM-32 etc.) in
which you want to specify the horizontal location of the points.
Cross section series
You must specify the first and the last point between which the cross section is defined. The cross section is defined as a section of element faces.
The face is included in the section when the line between the two element
centers of the faces crosses the line between the specified first and last
point. The geographical coordinates are taken from the dialog or from a
file. The file format is an ascii file with three space separated items for
each of the two points on separate lines. The first two items must be floats
(real numbers) for the x- and y-coordinate. The third item is unused (but
must be specified). If the file contains information for more than two
points (more than two lines) the information for the first two points will be
used. The faces defining the cross section are listed in the log-file.
You must also select the map projection (Long/Lat, UTM-32 etc.) in
which you want to specify the horizontal location of the points.
Domain series
The domain for which mass budget should be calculated is specified as a
polygon in the horizontal domain. The closed region is bounded by a

85

HYDRODYNAMIC MODULE

number of line segments. You must specify the coordinates of the vertex
points of the polygon. Two successive points are the endpoints of a line
that is a side of the polygon. The first and final point is joined by a line
segment that closes the polygon. The geographical coordinates of the polygon points are taken from the dialog or from a file. The file format is an
ascii file with three space separated items for each of the two points on
separate lines. The first two items must be floats (real numbers) for the xand y-coordinate. The third item is unused (but must be specified).
You must also select the map projection (Long/Lat, UTM-32 etc.) in
which you want to specify the horizontal location of the points.

6.18.3

Output items
Field variables
All output is optional, i.e. the user is free to select among the variables.
Flow directions in horizontal domain are given in degrees positive clockwise from true North (going against). Flow directions in vertical domain
are given in degrees positive clockwise from the upward pointing z-axis.
The convergence angle is the angle from true North to projection North
(positive clockwise).
Mass Budget
You can select the mass budget calculation to be included for the flow and
for the temperature and the salinity if calculation of the temperature/salinity is included. For each selected component the following items are
included in the output file

86

Total area - total volume/energy/mass within polygon

Wet area - volume/energy/mass in the area within polygon for which


the water depth is larger than the drying depth.

Real wet area - volume/energy/mass in the area within polygon for


which the water depth is larger than the wetting depth

Dry area - volume/energy/mass in the area within polygon for which


the water depth is less than the drying depth

Transport - accumulated volume/energy/mass transported over lateral


limits of polygon

Source - accumulated volume/energy/mass added/removed by sources


within polygon

MIKE 21 Flow Model FM

Outputs

Process - accumulated volume/energy/mass added/removed by processes within polygon

Error - accumulated volume/energy/mass error within polygon determined as the difference between the total mass change and the accumulated mass due to transport, sources and processes

The accumulated volume/energy/mass error contains the contribution due


to correction of the transported component when the values become larger
than the specified maximum value or lower than the specified minimum
value. For the water volume the minimum value is 0, while there is no
upper limit. For temperature and salinity the maximum and minimum values are specified on the Equation dialog.
For temperature the relative energy (relative to the energy at 0 C), E (J),
in a volume of water V (m3) is calculated by
E = c water water T V

(6.24)

where T (C) is the temperature, cwater = 4217.0 J/(K.kg) is the specific


heat for water and water = 1010.0 kg/m3 is the density of water.
For salinity the mass of salt, M (kg), in a volume of water V (m3) is calculated by
M = water S V 1000

(6.25)

where S (PSU) is the salinity and water = 1010.0 kg/m3 is the density of
water.
Discharge
You can select the discharge calculation to be included for the flow and for
the temperature and the salinity if calculation of the temperature/salinity is
included. For each selected component the following items are included in
the output file
!

Discharge - volume/energy/mass flux through the cross section

Acc. discharge - accumulated volume/energy/mass flux through the


cross section

By definition, discharge is positive for flow towards left when positioned


at the first point and looking forward along the cross-section line. The
transports are always integrated over the entire water depth.

87

HYDRODYNAMIC MODULE

For temperature the relative energy (relative to the energy at 0 C), E (J),
in a volume of water V (m3) is calculated by
E = c water water T V

(6.26)

where T (C) is the temperature, cwater = 4217.0 J/(K.kg) is the specific


heat for water and water = 1010.0 kg/m3 is the density of water.
For salinity the mass of salt, M (kg), in a volume of water V (m3) is calculated by
M = water S V 1000

(6.27)

where S (PSU) is the salinity and water = 1010.0 kg/m3 is the density of
water.

88

MIKE 21 Flow Model FM

Equation

TEMPERATURE/SALINITY MODULE
The temperature/salinity (TS) module is invoked from the specification of
the density, provided baroclinic density (density depends on temperature
and/or salinity) is selected (see Density (p. 36)).
The TS module sets up additional transport equations for temperature and
salinity. Additionally the calculated temperature and salinity are feed-back
to the hydrodynamic equations through bouyancy forcing induced by density gradients.

7.1

Equation
Here you can specify maximum and minimum values for the temperature
(in C) and the salinity (in PSU).

7.2

Solution Technique
The simulation time and accuracy can be controlled by specifying the
order of the numerical schemes which are used in the numerical calculations. Both the scheme for time integration and for space discretization
can be specified. You can select either a lower order scheme (first order)
or a higher order scheme. The lower order scheme is faster, but less accurate. For more details on the numerical solution techniques, see the Scientific documentation.
The time integration of the transport (advection-dispersion) equations is
performed using an explicit scheme. Due to the stability restriction using
an explicit scheme the time step interval must be selected so that the CFL
number is less than 1. A variable time step interval is used in the calculation and it is determined so that the CFL number is less than a critical CFL
number in all computational nodes. To control the time step it is also possible for the user to specify a minimum time step and a maximum time
step. The time step interval for the transport equations is synchronized to
match the overall time step specified on the Time dialog.
The minimum and maximum time step interval and the critical CFL
number is specified in the Solution Technique dialog in the HYDRODYNAMIC MODULE.

7.2.1

Remarks and hints


If the important processes are dominated by convection (flow), then
higher order space discretization should be chosen. If they are dominated

89

TEMPERATURE/SALINITY MODULE

by diffusion, the lower order space discretization can be sufficiently accurate. In general, the time integration method and space discretization
method should be chosen alike.
Choosing the higher order scheme for time integration will increase the
computing time by a factor of 2 compared to the lower order scheme.
Choosing the higher order scheme for space discretization will increase
the computing time by a factor of 1 to 2. Choosing both as higher order
will increase the computing time by a factor of 3-4. However, the higher
order scheme will in general produce results that are significantly more
accurate than the lower order scheme.
The default value for the critical CFL number is 1, which should secure
stability. However, the calculation of the CFL number is only an estimate.
Hence, stability problems can occur using this value. In these cases you
can reduce the critical CFL number. It must be in the range from 0 to 1.
Alternatively, you can reduce the maximum time step interval. Note, that
setting the minimum and maximum time step interval equal to the overall
time step interval specified on the Time dialog, the time integration will be
performed with constant time step. In this case the time step interval
should be selected so the the CFL number is smaller than 1.
The total number of time steps in the calculation and the maximum and
minimum time interval during the calculation are printed in the log-file for
the simulation. The CFL number can be saved in an output file.
The higher order scheme can exhibit under and over shoots in regions with
steep gradients. Hence, when the higher order scheme is used in combination with a limitation on the minimum and maximum value of the concentration mass conservation cannot be guarenteed.

7.3

Dispersion
In 2D models the dispersion usually describes transport due to nonresolved processes. In coastal areas it can be transport due to non-resolved
turbulence or eddies. Especially in the horizontal directions the effects of
non-resolved processes can be significant, in which case the dispersion
coefficient formally should depend on the resolution.

7.3.1

Horizontal dispersion
The horizontal dispersion can be formulated in two different ways
!

90

No dispersion

MIKE 21 Flow Model FM

Dispersion

Scaled eddy viscosity formulation

Dispersion coefficient formulation

Using the scaled eddy viscosity formulation the dispersion coefficient is


calculated as the eddy viscosity used in solution of the flow equations
multiplied by at scaling factor. For specification of the eddy viscosity, see
Eddy Viscosity (p. 37).
Selecting the dispersion coefficient formulation you must specify the dispersion coefficient (in m2/s).
Data
Selecting Scaled eddy viscosity option the format of the scaling factor can
be specified as
!

Constant

Varying in domain

For the case with values varying in domain you have to prepare a data file
containing the scaling factor before you set up the hydrodynamic simulation. The file must be a 2D unstructured data file (dfsu) or a 2D data grid
file (dfs2). The area in the data file must cover the model area. If a dfsufile is used piecewice constant interpolation is used to map the data. If a
dfs2-file is used bilinear interpolation is used to map the data.
Selecting Dispersion coefficient option the format of the dispersion coefficent (in m2/s) can be specified as
!

Constant

Varying in domain

For the case with values varying in domain you have to prepare a data file
containing the dispersion coefficient (in m2/s) before you set up the hydrodynamic simulation. The file must be a 2D unstructured data file (dfsu) or
a 2D data grid file (dfs2). The area in the data file must cover the model
area. If a dfsu-file is used piecewice constant interpolation is used to map
the data. If a dfs2-file is used bilinear interpolation is used to map the data.

7.3.2

Recommended values
When more sophisticated eddy viscosity models, such as the Smagorinsky
or k- models, are used the scaled eddy formulation it is recommended to
use the scaled viscosity formulation.

91

TEMPERATURE/SALINITY MODULE

The scaling factor can be estimated by 1/, where is the Prandtl


number. The default value here for the Prandtl number is 0.9 corresponding to a scaling factor of 1.1.
The dispersion coefficient is usually one of the key calibration parameters
for the Temperature/Salinity Module. It is therefore difficult to device
generally applicable values for the dispersion coefficient. However, using
Reynolds analogy, the dispersion coefficient can be written as the product
of a length scale and a velocity scale. In shallow waters the length scale
can often be taken as the water depth, while the velocity scale can be given
as a typical current speed.
Values in the order of 1 are usually recommended for the scaling factor.

7.4

Heat Exchange
The heat in the water can interact with the athmosphere through heat
exchange. The heat exchange is calculated on basis of the four physical
processes
!

the long wave radiation

the sensible heat flux (convection)

the short wave radiation

the latent heat flux (evaporation)

The heat exchange affects the temperature (if varying).


If you include heat exchange you have to specify the the air temperature
(in C), the relative humidity (in %) and clearness coefficient (in %). A
clearness coefficient of 100% specifies a clear sky and 0% specifies
cloudy weather.
As many as possible of the parameters in the heat exchange description
have been hardcoded. However, some need to be specified as part of your
simulation set-up.
For the latent heat flux you need to specify the constants a1 and b1, called
Constant in Dalton's law and Wind coefficient in Dalton's law, respectively. Furthermore you should specify the constants a2 and b2 known as
Sun constants in ngstrm's law, the time displacement (if summer time
is effective) and the local standard meridian for the time zone. Heating by
short wave radiation in the water column is included. The short wave penetration is dependent on the visibility, which is specified on the menu as a

92

MIKE 21 Flow Model FM

Precipitation-Evaporation

Light extinction coefficient. You can also specify an exchange coefficient, which on the menu is known as Beta in Beer's Law.
For more details see the scientific documentation.
Data
The format of the air temperature (in C), the relative humidity (in %) and
clearness coefficient (in %) can be specified as
!

Constant (in both time and domain)

Varying in time and constant in domain

Varying in time and domain

For the case with temperature varying in time but constant in domain you
have to prepare a data file containing the relevant variable (the air temperature (in C), the relative humidity (in %) or clearness coefficient. (in %)),
before you set up the hydrodynamic simulation. The data file must be a
time series file (dfs0). The data must cover the complete simulation
period. The time step of the input data file does not, however, have to be
the same as the time step of the hydrodynamic simulation. A linear interpolation will be applied if the time steps differ.
For the case with temperature varying both in time and domain you have
to prepare a data file containing the relevant variable (the air temperature
(in C), the relative humidity (in %) or clearness coefficient. (in %)) before
you set up the hydrodynamic simulation. The file must be a 2D unstructured data file (dfsu) or a 2D grid data file (dfs2). The area in the data file
must cover the model area. If a dfsu-file is used piecewice constant interpolation is used to map the data. If a dfs2-file is used bilinear interpolation
is used to map the data. The data must cover the complete simulation
period. The time step of the input data file does not, however, have to be
the same as the time step of the hydrodynamic simulation. A linear interpolation will be applied if the time steps differ.

7.5

Precipitation-Evaporation
If your simulation include precipitation and/or evaporation, you need to
specify the temperature of the precipated and evaporated water mass. The
precipitation and evaporation can be included in one of two ways:
!

Ambient water temperature


The temperature of the precipated/eavaporated water mass is set equal
to the concentration of the ambient sea water.

93

TEMPERATURE/SALINITY MODULE

Specified temperature
The temperature (in C) of the precipated/evaporated water mass is
specified explictly.

If you have chosen the net precipitation option in the Precipitation-Evaporation dialog in the HYDRODYNAMIC MODULE the precipitation temperature will be used when the specified net-precipation is positive. When
the net-precipitation is negative, the evaporation temperature is used.
The salinity of the precipated/evaporated water mass is assumed to be
zero.
Data
Selecting the specified temperature option the format of the temperature
(in C) can be specified as
!

Constant (in both time and domain)

Varying in time and constant in domain

Varying in time and domain

For the case with temperature varying in time but constant in domain you
have to prepare a data file containing the temperature (in C) before you
set up the hydrodynamic simulation. The data file must be a time series
file (dfs0). The data must cover the complete simulation period. The time
step of the input data file does not, however, have to be the same as the
time step of the hydrodynamic simulation. A linear interpolation will be
applied if the time steps differ.
For the case with temperature varying both in time and domain you have
to prepare a data file containing the temperature (in C) before you set up
the hydrodynamic simulation. The file must be a 2D unstructured data file
(dfsu) or a 2D grid data file (dfs2). The area in the data file must cover the
model area. If a dfsu-file is used piecewice constant interpolation is used
to map the data. If a dfs2-file is used bilinear interpolation is used to map
the data. The data must cover the complete simulation period. The time
step of the input data file does not, however, have to be the same as the
time step of the hydrodynamic simulation. A linear interpolation will be
applied if the time steps differ.
Soft start interval
You can specify a soft start interval (in seconds) during which the precipitation/evaporation temperature is increased linearly from 0 to the specified
values of the precipitation/evaporation temperature. By default the soft
start interval is zero (no soft start).

94

MIKE 21 Flow Model FM

Sources

7.6

Sources
Point sources of temperature (heat) and salinity are important in many
applications as e.g. outlets from rivers, intakes and outlets from cooling
water or desalination plants.
In the Temperature/Salinity Module the temperature and/or salinity in
every sources point can be specified. The number of sources, their generic
names, location and discharges magnitude are specified in the Sources
dialog in the HYDRODYNAMIC MODULE.
Depending on the choice of property pge you can see a geographic view or
a list view of the sources.
The source information is specified individually for each source and for
temperature and salinity. From the list view you can go to the dialog for
specification by clicking on the Go to.. button.

7.6.1

Source Specification
The type of sources can be specified as
!

Specified temperature/salinity

Excess temperature/salinity

The source flux is calculated as Qsource.Csource where Qsource is the magnitude of the source and Csource is the temperature/salinity of the source. The
magnitude of the source is specified in the Sources dialog in the HYDRODYNAMIC MODULE.
Selecting the temperature/salinity option the source temperature/salinity is
the specified temperature/salinity if the magnitude of the source is positive
(water is discharge into the ambient water). The source temperature/salinity is the temperature/salinity at the source point if the magnitude of the
source is negative (water is discharge out the ambient water). This option
is pertinent to e.g. river outlets or other sources where the concentration is
independent of the surrounding water.
Selecting the excess temperature/salinity option the source temperature/salinity is the sum of the excess temperature/salinity and temperature/salinity at a point in the model if the magnitude of the source is
positive (water is discharge into the ambient water). If it is an isolated
source the point is the location of the source. If it is a connected source the
point is the location where water is discharged out of the water. The source
temperature/salinity is the temperature/salinity at the source point if the

95

TEMPERATURE/SALINITY MODULE

magnitude of the source is negative (water is discharge out the ambient


water). This type can be used to describe e.g. a heat exchange or other
processes where the temperature (heat) or salinity is added to the water by
a diffusion process.
Data
The format of the source information can be specified as
!

Constant (in time)

Varying in time

For the case with source concentration varying in time you have to prepare
a data file containing the concentration (temperature (in C) or salinity (in
PSU)) of the source before you set up the hydrodynamic simulation. The
data file must be a time series file (dfs0). The data must cover the complete simulation period. The time step of the input data file does not, however, have to be the same as the time step of the hydrodynamic simulation.
A linear interpolation will be applied if the time steps differ.

7.6.2

Remarks and hints


Point sources are entered into elements, such that the inflowing mass of
the component initially is distributed over the element where the source
resides. Therefore the temperature and salinity seen in the results from the
simulation usually is lower than the source specified temperature and
salinity.

7.7

Initial Conditions
The initial conditions are the spatial distribution of temperature and salinity throughout the computational domain at the beginning of the simulation. The initial conditions can be the result from a previous simulation in
which case the initial conditions effectively act as a hot start of the temperature and salinity field.
The initial conditions are specified individually for temperature and salinity.
Data
The format of the initial conditions for temperature (in C) and salinity (in
PSU) can be specified as
!

96

Constant (in domain)

MIKE 21 Flow Model FM

Boundary Conditions

Varying in domain

For the case with varying in domain you have to prepare a data file containing the initial conditions (temperature (in C) or salinity (in PSU))
before you set up the hydrodynamic simulation. The file must be a 2D
unstructured data file (dfsu) or a 2D grid data file (dfs2). The area in the
data file must cover the model area. If a dfsu-file is used piecewice constant interpolation is used to map the data. If a dfs2-file is used bilinear
interpolation is used to map the data. In case the input data file contains a
single time step, this field is used. In case the file contains several time
steps, e.g. from the results of a previous simulation, the actual starting
time of the simulation is used to interpolate the field. Therefore the starting time must be within the start and end time of the file.

7.8

Boundary Conditions
Initially, the set-up editor scans the mesh file for boundary codes (sections), and displays the recognized codes and suggest a default name for
each. You can re-name these names to more meningful names in the
Domain dialog (see Boundary names (p. 28)).
Depending on the choice of property page you can get a geographic view
or a list view of the boundaries.
The boundary conditions are specified individually for each code (section)
and for the temperature and the salinity. From the list view you can go to
the dialog for specification by clicking on the "Go to .." button.

7.8.1

Boundary Specification
You can choose between the following three boundary types
!

Land

Specified values (Dirichlet boundary condition)

Zero gradient (Neumann boundary condition).

Data
If specified values (Dirichlet boundary condition) is selected the format
of the temperature (in C) and/or salinity (in PSU) at the boundary can be
specified as
!

Constant (in time and along boundary)

Varying in time and constant along boundary

97

TEMPERATURE/SALINITY MODULE

Varying in time and along boundary

For the case with boundary data varying in time but constant along the
boundary you have to prepare a data file containing the boundary values
(temperature (C) or Salinity (PSU)) before you set up the hydrodynamic
simulation. The data file must be a time series file (dfs0). The data must
cover the complete simulation period. The time step of the input data file
does not, however, have to be the same as the time step of the hydrodynamic simulation. You can choose between different types of interpolation
(see Interpolation type).
For the case with boundary data varying both in time and along the boundary you have to prepare a data file containing the boundary values (temperature (C) or Salinity (PSU)) before you set up the hydrodynamic
simulation. The data file must be a profile file (dfs1). The data must cover
the complete simulation period. The time step of the input data file does
not, however, have to be the same as the time step of the hydrodynamic
simulation. You can choose between different types of time interpolation
(see Interpolation type).
Interpolation type
For the two cases with values varying in time two types of time interpolation can be selected:
!

linear

piecewise cubic.

In the case with values varying along the boundary two methods of mapping from the input data file to the boundary section are available:
!

normal

reverse order

Using normal interpolation the first and last point of the line are mapped
to the first and the last node along the boundary section and the intermediate boundary values are found by linear interpolation. Using reverse order
interpolation the last and first point of the line are mapped to the first and
the last node along the boundary section and the intermediate boundary
values is found by linear interpolation.
Soft start interval
You can specify a soft start interval (in seconds) during which boundary
values are increased from a specified reference value to the specified

98

MIKE 21 Flow Model FM

Boundary Conditions

boundary value in order to avoid shock waves being generated in the


model. The increase can either be linear or follow a sinusoidal curve.

99

TEMPERATURE/SALINITY MODULE

100

MIKE 21 Flow Model FM

TURBULENCE MODULE
The Turbulence module is available for 3D simulations only.
Please see User Guide for MIKE 3 Flow Model FM.

101

TURBULENCE MODULE

102

MIKE 21 Flow Model FM

LIST OF REFERENCES
Aupoix, B., (1984), Eddy Viscosity Subgrid Scale Models for Homogeneous Turbulence, in Macroscopic Modelling of Turbulent Flow, Lecture
Notes in Physics, Proc. Sophie-Antipolis, France.
Delhez, E. J. M., Deelersijder, E., Maouchet, A., Beckers, J. M.(2003), On
the age of radioactive tracers, J. Marine Systems, 38, pp. 277-286.
Flather, R. A. (1976), A tidal model of the northwest European continental
shelf, Memories de la Societe Royale des Sciences de Liege 6 (10), 141164.
Geernaert G.L. and Plant W.L (1990), Surface Waves and fluxes, Volume 1
- Current therory, Kluwer Academic Publishers, The Netherlands.
Holmes, D. G. and Connell, S. D. (1989), Solution of the 2D NavierStokes on unstructured adaptive grids, AIAA Pap. 89-1932 in Proc.
AIAA 9th CFD Conference.
Horiuti, K. (1987), Comparison of Conservative and Rotational Forms in
Large Eddy Simulation of Turbulent Channel Flow, Journal of Computational Physics, 71, pp 343-370.
Leonard, A. 1974), Energy Cascades in Large-Eddy Simulations of Turbulent Fluid Flows, Advances in Geophysics, 18, pp 237-247.
Lilly, D.K. (1966), On the Application of the Eddy Viscosity Concept in
the Inertial Subrange of Turbulence, NCAR Manuscript No. 123, National
Center for Atmospheric Research, Boulder, Colorado.
Munk, W., Anderson, E. (1948), Notes on the theory of the thermocline,
Journal of Marine Research, 7, 276-295.
Oddo P. and N. Pinardi (2007), Lateral open boundary conditions for
nested limited area models: A scale selective approach, Ocean Modelling
20 (2008) 134-156.
Rodi, W. (1980), Turbulence Models and Their Application in Hydraulics
- A State of the Art Review, Special IAHR Publication.
Smagorinsky (1963), J. General Circulation Experiment with the Primitive Equations, Monthly Weather Review, 91, No. 3, pp 99-164.

103

LIST OF REFERENCES

Sleigh, P.A., Gaskell, P.H., Bersins, M. and Wright, N.G. (1998), An


unstructured finite-volume algorithm for predicting flow in rivers and
estuaries, Computers & Fluids, Vol. 27, No. 4, 479-508.
Song, Y and Haidvogel D., (1994) A semi-implicit ocean circulation
model using a generalized topography-following coordinate system, Journal of Comp. Physics, 115, pp. 228-244.
Wu, Jin (1994), The sea surface is aerodynamically rough even under light
winds, Boundary layer Metheorology, 69, 149-158.
Wu, Jin (1980), Wind-stress Coefficients over sea surface and near neutral
conditions - A revisit, J. Phys. Oceanogr., 10, 727-740.
Zhao, D.H., Shen, H.W., Tabios, G.Q., Tan, W.Y. and Lai, J.S. (1994),
Finite-volume 2-dimensional unsteady-flow model for river basins, Journal of Hydraulic Engineering, ASCE, 1994, 120, No. 7, 863-833.

104

MIKE 21 Flow Model FM

INDEX

105

Index

About this guide . . . . . . . . . . . . .9


Area series . . . . . . . . . . . . . . . 85

Heat exchange . . . . . . . . . . . . . 92
Horizontal dispersion . . . . . . . . . . 90

Bed resistance . . . . . . . . . . . . . 39
Boundary conditions . . . . . . . . 74, 97
Bridge . . . . . . . . . . . . . . . . . . 68

Ice coverage . . . . . . . . . . . . . . . 44
Initial conditions . . . . . . . . . . . 73, 96

L
C
CFL number . . . . . . .
Composite structures .
Connected source . . .
Convention of boundary
Convergence angle . .
Coordinate type . . . . .
Coriolis forcing . . . . .
Courant-Friedrich-Lvy
Culvert . . . . . . . . . .

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

33
68
50
80
86
27
41
33
57

Latent heat flux . . . . . . . . . . . . . 92


Line series . . . . . . . . . . . . . . . . 84
Long wave radiation . . . . . . . . . . 92

M
Mesh file . . . . . . . . . . . . . . . . . 27
MIKE Zero Mesh . . . . . . . . . . . . 27
Minimum depth cutoff . . . . . . . . . 27

N
Neutral pressure

. . . . . . . . . . . . 42

D
28
28
36
80
80
90
29
43
35

Eddy viscosity . . . . . . . . . . . . . 37
Evaporation . . . . . . . . . . . . . 46, 93
Excess temperature . . . . . . . . . . 95

Datum shift . . . . . .
Decomposition . . . .
Density . . . . . . . . .
Direction of boundary
Direction of discharge
Dispersion . . . . . . .
Domain parameters .
Drag coefficient . . . .
Drying depth . . . . .

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

On-line help . . . . . . . . . . . . . . . 11
Outputs . . . . . . . . . . . . . . . . . . 82

P
Pier . . . . . . . . .
Pier section . . . .
Point series . . . .
Precipitation . . . .
Pressure correction

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

. . . 64
. . . 66
. . . 84
. 46, 93
. . . 79

E
Radiation stress . . . . . . . . . . . . . 48

S
F
Flood and dry . . . . . . . . . . .
Flooding depth . . . . . . . . . .
Flow direction . . . . . . . . . . .
Free Overflow Head Loss Factor

.
.
.
.

. .
. .
72,
. .

35
35
86
57

Salinity . . . . . . . .
Sensible heat flux . .
Short wave radiation
Simple source . . . .
Smagorinsky . . . .
Sources . . . . . . .
Standard source . .

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

. . . 89
. . . 92
. . . 92
. . . 50
. . . 38
. 49, 95
. . . 50

Gate . . . . . . . . . . . . . . . . . . . 62

106

MIKE 21 Flow Model FM

Index

T
Temperature . . .
Time parameters
Turbine . . . . . .
Turbulence . . .

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

. 89
. 29
. 66
101

U
User background . . . . . . . . . . . . 9

W
Wave radiation . .
Weir . . . . . . . .
Wetting depth . . .
Wind forcing . . . .

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

. 48
. 51
. 35
. 41

107

Index

108

MIKE 21 Flow Model FM

You might also like