Telemac3d_user_manual_v7p1_ap_verif_JF Telemac3d User Manual V7p1
telemac3d_user_manual_v7p1
User Manual: Pdf
Open the PDF directly: View PDF 
.
Page Count: 108
| Download | |
| Open PDF In Browser | View PDF |
TELEMAC MODELLING SYSTEM
3D hydrodynamics
TELEMAC-3D Software
Release 7.1
OPERATING MANUAL
JULY 2016
The information given in this manual is subject to revision without notice. EDF R&D
disclaims any responsibility for or in relation to the contents hereof.
The TELEMAC system is the property of EDF R&D.
© Copyright 2016 EDF R&D
************************
BLUE KENUE is the property of the Canadian Hydraulics Centre, Ottawa, Ontario,
Canada
Copyright ©1998-2012 Canadian Hydraulics Centre, National Research Council
http://www.nrc-cnrc.gc.ca/fra/idp/chc/logiciels/kenue/blue-kenue.html
FUDAA-PREPRO is the property of the CEREMA, Compiègne, France
JANET is the property of the Smile Consult GmbH, Hannover, Deutschland
DELWAQ is the property of the DELTARES, Delft, The Netherlands
TECPLOT is the property of Tecplot Inc., USA
EVOLUTIONS OF THE
DOCUMENT
DATE
EDITOR
EVOLUTION
01/2013
Jonathan DESOMBRE
jonathan.desombre@ingerop.com
General update for release version 6.2
11/2014
Chi-Tuân PHAM, Cédric GOEURY, Antoine
JOLY
General update for release version 7.0
01/2016
Chi-Tuân PHAM, Antoine JOLY
General update for release version 7.1
Typing conventions used in this manual
The computational items (variable names, file names, etc.) are written in courier
font.
The keywords are written in UPPER CASE ITALICS.
The literature references are given between brackets [ ].
FOREWORD
This operating manual has been drafted for the 7.1 software release.
This manual does not deal with the sedimentology subroutines nor with the coupling to
DELWAQ, which are included in TELEMAC-3D.
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page a
TABLE OF CONTENT
1.
2.
INTRODUCTION ......................................................................................................... 7
1.1.
POSITION OF THE TELEMAC-3D CODE WITHIN THE TELEMAC MODELLING SYSTEM......8
1.2.
SOFTWARE ENVIRONMENT ...................................................................................................8
1.3.
USER PROGRAMMING............................................................................................................8
THEORETICAL ASPECTS ....................................................................................... 10
2.1.
NOTATIONS ........................................................................................................................... 10
2.2.
EQUATIONS ........................................................................................................................... 10
2.2.1. EQUATION WITH THE HYDROSTATIC PRESSURE HYPOTHESIS ............................ 10
2.2.2. NON-HYDROSTATIC NAVIER-STOKES EQUATIONS ................................................. 13
2.2.3. THE LAW OF STATE .................................................................................................... 13
2.2.4. K-ε MODEL ................................................................................................................... 14
2.2.5. EQUATIONS OF TRACERS.......................................................................................... 14
2.3.
MESH DESCRIPTION............................................................................................................. 15
2.3.1. THE MESH STRUCTURE ............................................................................................. 15
2.3.2. THE TWO-DIMENSIONAL MESH ................................................................................. 15
2.3.3. THE THREE-DIMENSIONAL MESH.............................................................................. 16
3.
THE INPUTS / OUTPUTS.......................................................................................... 17
3.1.
PRELIMINARY COMMENTS .................................................................................................. 17
3.1.1. BINARY FILES FORMAT .............................................................................................. 18
3.2.
THE STEERING FILE ............................................................................................................. 19
3.3.
THE GEOMETRY FILE ........................................................................................................... 21
3.4.
THE BOUNDARY CONDITIONS FILE..................................................................................... 21
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page b
3.5.
THE FORTRAN FILE .............................................................................................................. 21
3.6.
THE LIQUID BOUNDARIES FILE............................................................................................ 22
3.7.
THE PREVIOUS COMPUTATION FILE .................................................................................. 22
3.8.
THE REFERENCE FILE.......................................................................................................... 22
3.9.
THE STAGE-DISCHARGE CURVES FILE .............................................................................. 22
3.10. THE SOURCES FILE .............................................................................................................. 23
3.11. THE 3D RESULT FILE ............................................................................................................ 23
3.12. THE 2D RESULT FILE ............................................................................................................ 24
3.13. THE OUTPUT LISTING........................................................................................................... 25
3.14. THE FILE FOR SCOPE........................................................................................................... 26
3.15. THE AUXILIARY FILES........................................................................................................... 26
3.15.1.
THE DICTIONNARY FILE ................................................................................. 27
3.16. TOPOGRAPHIC AND BATHYMETRIC DATA ......................................................................... 27
4.
GENERAL SETUP OF THE HYDRODYNAMIC COMPUTATION (NAVIER-STOKES
EQUATIONS) ............................................................................................................ 28
4.1.
MESH DEFINITION................................................................................................................. 28
4.2.
PRESCRIBING THE INITIAL CONDITIONS ............................................................................ 31
4.2.1. PRESCRIPTION THROUGH KEYWORDS.................................................................... 31
4.2.2. PRESCRIBING PARTICULAR INITIAL CONDITIONS (PROGRAMMING THE CONDIM
SUBROUTINE) ............................................................................................................. 32
4.2.3. RESUMING THE COMPUTATION ................................................................................ 32
4.3.
PRESCRIBING THE BOUNDARY CONDITIONS .................................................................... 33
4.3.1. THE BOUNDARIES IN TELEMAC-3D ........................................................................... 34
4.3.2. THE BOUNDARY-RELATED TYPES ............................................................................ 35
4.3.3. DESCRIPTION OF THE VARIOUS TYPES ................................................................... 36
4.3.4. THE BOUNDARY CONDITIONS FILE ........................................................................... 36
4.3.5. PROGRAMMING THE BOUNDARY CONDITIONS TYPE ............................................. 39
4.3.6. PRESCRIBING VALUES THROUGH KEYWORDS ....................................................... 39
4.3.7. BOUNDARY CONDITION ON THE BOTTOM ............................................................... 40
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page c
4.3.8. USING THE LIQUID BOUNDARIES FILE ...................................................................... 42
4.3.9. PRESCRIBING VALUES THROUGH PROGRAMMING ................................................ 43
5.
4.3.10.
STAGE-DISCHARGE CURVES ........................................................................ 43
4.3.11.
PRESCRIBING COMPLEX VALUES................................................................. 44
4.3.12.
PRESCRIBING A PROFILE .............................................................................. 44
4.3.13.
THOMPSON CONDITIONS .............................................................................. 45
4.3.14.
TIDAL HARMONIC CONSTITUENTS DATABASES.......................................... 46
PHYSICAL SETUP OF THE HYDRODYNAMIC COMPUTATION ............................ 49
5.1.
HYDROSTATIC PRESSURE HYPOTHESIS ........................................................................... 49
5.2.
MODELLING TURBULENCE .................................................................................................. 49
5.2.1. CONSTANT VISCOSITY ............................................................................................... 50
5.2.2. MIXING LENGTH (VERTICAL MODEL) ........................................................................ 50
5.2.3. SMAGORINSKY............................................................................................................ 52
ε
5.2.4. K- ............................................................................................................................... 52
5.3.
SETTING UP THE FRICTION ................................................................................................. 53
5.3.1. BOTTOM FRICTION ..................................................................................................... 54
5.3.2. SIDEWALL FRICTION .................................................................................................. 54
5.4.
PUNCTUAL SOURCE TERMS................................................................................................ 55
5.5.
SETTING UP THE WATER-ATMOSPHERE EXCHANGES..................................................... 56
5.5.1. THE WIND .................................................................................................................... 56
5.5.2. THE TEMPERATURE ................................................................................................... 57
5.5.3. THE PRESSURE .......................................................................................................... 57
5.5.4. RAIN OR EVAPORATION ............................................................................................. 57
5.5.5. ATMOSPHERE-WATER EXCHANGE MODELS ........................................................... 57
5.6.
ASTRAL POTENTIAL ............................................................................................................. 59
5.7.
CONSIDERATION OF WAVE DRIVEN CURRENTS ............................................................... 59
5.8.
OTHER PHYSICAL PARAMETERS ........................................................................................ 59
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
6.
Page d
NUMERICAL SETUP OF THE COMPUTATION ....................................................... 61
6.1.
GENERAL SETUP .................................................................................................................. 61
6.1.1. ADVECTION STEP ....................................................................................................... 61
6.1.2. DIFFUSION STEP......................................................................................................... 61
6.1.3. PROPAGATION STEP .................................................................................................. 62
6.2.
THE ADVECTION SCHEME ................................................................................................... 62
6.2.1. ADVECTION OF THE THREE-DIMENSIONAL VARIABLES ......................................... 63
6.2.2. CONFIGURATION OF THE SUPG SCHEME ................................................................ 64
6.2.3. CONFIGURATION OF THE WEAK CHARACTERISTICS.............................................. 64
6.3.
SPECIFIC PARAMETERS IN THE NON-HYDROSTATIC VERSION....................................... 64
6.4.
IMPLICITATION ...................................................................................................................... 65
6.5.
SOLUTION OF THE LINEAR SYSTEMS................................................................................. 66
6.5.1. SOLVERS ..................................................................................................................... 66
6.5.2. ACCURACIES............................................................................................................... 67
6.5.3. PRECONDITIONINGS .................................................................................................. 67
6.6.
TIDAL FLATS.......................................................................................................................... 68
6.7.
HYDROSTATIC INCONSISTENCIES...................................................................................... 69
6.8.
OTHER PARAMETERS .......................................................................................................... 70
6.8.1. MASS-LUMPING........................................................................................................... 70
6.8.2. CONVERGENCE AID ................................................................................................... 70
6.8.3. MATRIX STORAGE ...................................................................................................... 70
6.8.4. VELOCITIES PROJECTION.......................................................................................... 70
7.
TRACER TRANSPORT ............................................................................................ 72
7.1.
GENERAL SETUP .................................................................................................................. 72
7.1.1. PRESCRIBING THE INITIAL CONDITIONS .................................................................. 72
7.1.2. PRESCRIBING THE BOUNDARY CONDITIONS .......................................................... 72
7.2.
PHYSICAL SETUP.................................................................................................................. 73
7.2.1. ACTIVE TRACERS ....................................................................................................... 73
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page e
7.2.2. PUNCTUAL SOURCE TERMS...................................................................................... 74
7.2.3. GENERAL SOURCE TERMS ........................................................................................ 75
7.3.
8.
9.
NUMERICAL SETUP .............................................................................................................. 75
DROGUES ................................................................................................................ 76
8.1.
CONFIGURATION OF SIMULATION ...................................................................................... 76
8.2.
VISUALISATION OF RESULTS .............................................................................................. 76
OIL SPILL MODELLING ........................................................................................... 78
9.1.
INPUT FILES .......................................................................................................................... 78
9.2.
STEERING FILE ..................................................................................................................... 78
9.3.
OIL SPILL STEERING FILE .................................................................................................... 79
9.4.
THE OIL_FLOT SUBROUTINE ............................................................................................... 80
9.5.
OUTPUT FILES ...................................................................................................................... 81
9.5.1. THE 3D RESULT FILE .................................................................................................. 81
9.5.2. THE OUTPUT DROGUES FILE .................................................................................... 81
10. OTHERS CONFIGURATIONS .................................................................................. 83
10.1. MODIFICATION OF BOTTOM TOPOGRAPHY (CORFON)..................................................... 83
10.2. MODIFYING COORDINATES (CORRXY) ............................................................................... 83
10.3. SPHERICAL COORDINATES (LATITU) .................................................................................. 83
10.4. ADDING NEW VARIABLES .................................................................................................... 84
10.5. ARRAY MODIFICATION OR INITIALIZATION ........................................................................ 85
10.6. VALIDATING A COMPUTATION (VALIDA) ............................................................................. 86
10.7. COUPLING ............................................................................................................................. 86
10.8. CHECKING THE MESH (CHECKMESH)................................................................................. 87
11. PARALLELISM ......................................................................................................... 88
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page f
LIST OF FIGURES
FIGURE 1: A THREE-DIMENSIONAL MESH. ....................................................... 17
FIGURE 2: EFFECT OF THE MESH TRANSFORMATION KEYWORD – VALUE 1:
SIGMA.
29
FIGURE 3: EFFECT OF THE MESH TRANSFORMATION KEYWORD – VALUE 2:
ZSTAR.
29
FIGURE 4: EFFECT OF THE MESH TRANSFORMATION KEYWORD – VALUE 3:
USER DEFINED. 30
FIGURE 5: THE VARIOUS BOUNDARIES IN TELEMAC-3D (BRIDGE PIERS CASE
35
STUDY).
FIGURE 6: MIXING LENGTHS VERSUS DEPTH. .................................................. 51
FIGURE 7: MUNK AND ANDERSON DAMPING FUNCTION. ................................... 52
LIST OF APPENDIX
APPENDIX N° 1. LAUNCHING THE COMPUTATION ............................................ 89
APPENDIX N° 2. LIST OF USER SUBROUTINES ................................................. 90
APPENDIX N° 3. DESCRIPTION OF THE SERAFIN FORMAT ................................. 91
APPENDIX N° 4. GENERATING OUTPUT FILES FOR DELWAQ .......................... 93
APPENDIX N° 5. POSTEL-3D ........................................................................ 94
APPENDIX N° 6. REFERENCES ...................................................................... 98
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
1.
Page 7
INTRODUCTION
The TELEMAC-3D code solves such three-dimensional equations as the free surface flow
equations (with or without the hydrostatic pressure hypothesis) and the transport-diffusion
equations of intrinsic quantities (temperature, salinity, concentration). Its main results, at each point
in the resolution mesh in 3D, are the velocity in all three directions and the concentrations of
transported quantities. Water depth is the major result as regards the 2D surface mesh. The
TELEMAC-3D’s prominent applications can be found in free surface flow, in both seas and rivers;
the software can take the following processes into account:
•
Influence of temperature and/or salinity on density,
•
Bottom friction,
•
Influence of the Coriolis force,
•
Influence of weather elements: air pressure, rain or evaporation and wind,
•
Consideration of the thermal exchanges with the atmosphere,
•
Sources and sinks for fluid moment within the flow domain,
•
Simple or complex turbulence models (k-ε) taking the effects of the Archimedean force
(buoyancy) into account,
•
Dry areas in the computational domain: tidal flats,
•
Current drift and diffusion of a tracer, with generation or disappearance terms,
•
Oil spill modelling.
The code is applicable to many fields. The main ones are related to the marine environment
through the investigations of currents being induced either by tides or density gradients, with or
without the influence of such an external force as the wind or the air pressure. It can be applied
either to large extent areas (on a sea scale) or to smaller domains (coasts and estuaries) for the
impact of sewer effluents, the study of thermal plumes or even sedimentary transport. As regards
the continental waters, the study of thermal plumes in rivers, the hydrodynamic behaviour or
natural or man-made lakes can be mentioned as well.
TELEMAC-3D is developed by the LNHE (Laboratoire National d’Hydraulique et Environnement) of
the Research and Development Division of EDF (EDF-R&D). As for previous versions, the 7.1
release of the code complies with the Quality Assurance procedures of scientific and technical
softwares of EDF-R&D. It is a process of construction and verification of the product quality in the
different phases of his life. In particular, a software following the Quality Assurance procedures
comes with a Validation Folder that describes the intended use of the software and a set of test
cases. This document allows you to judge the performance and limitations of the software, situating
the field of application. These tests are also used in development of the software and are checked
at every new release.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
1.1.
Page 8
POSITION OF THE TELEMAC-3D CODE WITHIN THE TELEMAC
MODELLING SYSTEM
The TELEMAC-3D software is part of the TELEMAC modelling system developed by the LNHE of
EDF R&D. TELEMAC is a set of modelling tools allowing to treat every aspects of natural free
surface hydraulics: currents, waves, transport of tracers and sedimentology.
The pre-processing and post-processing of simulations can be done either directly within the
TELEMAC system or with different software that present an interface of communication with the
system. We can particularly mention the following tools:
• The FUDAA-PREPRO software, developed from the FUDAA platform by the CEREMA’s
Recherche, Informatique et Modélisation Department, covers all the pre-processing tasks
involved by the achievement of a numerical hydraulic study, as well as a graphical postprocessing tool,
• The BLUE KENUE software, developed the Hydraulic Canadian Center, proposes a
powerful mesh generation tool and a user-friendly post-processing tool,
• The JANET software, developed by Smile Consult GmbH, which offers among others, a
mesh generation tool,
• The PARAVIEW software, developed by Sandia National Laboratories, Los Alamos
National Laboratory and Kitware, which enables to visualise 3D results, big data in
particular and is open source,
• The QGIS software, which is an open source Geographic Information System.
1.2.
SOFTWARE ENVIRONMENT
All the simulation modules are written in FORTRAN 90, with no use of the specific language
extensions in a given machine. They can be run on all the PCs (or PC "clusters") under Windows
and Linux operating systems as well as on the workstations under the Unix operating system.
1.3.
USER PROGRAMMING
When using a simulation module from the TELEMAC system, the user may have to program
specific subroutines which are not in the code’s standard release. In particular, that is made
through a number of so-called « user » subroutines. These subroutines are written so that they can
be modified, provided that the user has a basic knowledge in FORTRAN language, with the help of
the « Guide for programming in the Telemac system » [5].
The procedure to be carried out in that case comprises the steps of:
RELEASE 7.1
•
Recovering the standard version of the user subroutine(s) as supplied in the distribution
and copying it into the current directory,
•
Amending the subroutine(s) according to the model to be constructed,
•
Concatenating the whole set of subroutines into a single FORTRAN file which will be
compiled during the TELEMAC-3D launching process.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 9
During that programming stage, the user can gain access to the various variables of the software
through the FORTRAN 90 structures.
All the data structures are gathered within FORTRAN files, which are known as modules. For
TELEMAC-3D, the file name is DECLARATION_TELEMAC3D. To gain access to the TELEMAC-3D
data, just insert the command USE DECLARATIONS_TELEMAC3D into the beginning of the
subroutine. Adding the command USE BIEF may also be necessary in order to reach the
structure in the BIEF library.
Nearly all the arrays which are used by TELEMAC-3D are declared in the form of a structure. For
example, the access to the water depth array will be in the form H%R, %R meaning it is a real-typed
pointer. In case of an integer-typed pointer, the %R is replaced by a %I. However, in order to avoid
having to handle too many %R and %I, a number of aliases are defined, such as, the NPOIN3,
NELEM3 and NPTFR2 variables. For further details, the user can refer to the programming guide in
TELEMAC [5].
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
2.
THEORETICAL ASPECTS
2.1.
NOTATIONS
Page 10
TELEMAC-3D is a three-dimensional computational code describing the 3D velocity field (U, V, W)
and the water depth h (and, from the bottom depth, the free surface S) at each time step. Besides,
it solves the transport of several tracers which can be grouped into two categories, namely the socalled “active” tracers (primarily temperature and salinity1), which change the water density and act
on flow through gravity), and the so-called “passive” tracers which do not affect the flow and are
merely transported.
2.2.
EQUATIONS
The reader will refer to the J.-M. Hervouet’s book [1] for a detailed statement of the theoretical
aspects which TELEMAC-3D is based on.
2.2.1.
EQUATION WITH THE HYDROSTATIC PRESSURE HYPOTHESIS
In its basic release, the code solves the three-dimensional hydrodynamic equations with the
following assumptions:
1
•
Three-dimensional Navier-Stokes equations with a free surface changing in time,
•
Negligible variation of density in the conservation of mass equation (incompressible fluid),
•
Hydrostatic pressure hypothesis (that hypothesis results in that the pressure at a given
depth is the sum of the air pressure at the fluid surface plus the weight of the overlying
water body),
•
Boussinesq approximation for the momentum (the density variations are only taken into
account as buoyant forces).
Sediment transport with TELEMAC-3D is not described in this manual.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 11
Due to these assumptions, the three-dimensional equations being solved are:
∂U ∂V ∂W
+
+
=0
∂x ∂y
∂z
∂Z
∂U
∂U
∂U
∂U
+U
+V
+W
= − g S + υ ∆ (U ) + Fx
∂t
∂x
∂y
∂z
∂x
∂Z
∂V
∂V
∂V
∂V
+U
+V
+W
= − g S + υ ∆ (V ) + F y
∂t
∂x
∂y
∂z
∂y
Z ∆ρ
p = p atm + ρ 0 g (Z s − z ) + ρ 0 g ∫z
dz '
S
ρ0
∂T
∂T
∂T
∂T
+U
+V
+W
= div (υ grad T ) + Q
∂t
∂x
∂y
∂z
wherein:
h
(m)
water depth,
ZS
(m)
free surface elevation,
U, V, W (m/s)
three-dimensional components of velocity,
T
(°C, g/L…)
passive or active (acting on density) tracer,
p
(X)
pressure,
patm
(X)
atmospheric pressure,
g
(m/s2) acceleration due to gravity,
υ
(m2/s) cinematic viscosity or tracer diffusion coefficients,
Zf
(m)
bottom depth,
ρ0
(X)
reference density,
∆ρ
(X)
variation of density around the reference density,
t
(s)
time,
x, y
(m)
horizontal space components,
z
(m)
vertical space component,
Fx, Fy
(m/s2) source terms,
Q
(tracer unit)
tracer source of sink.
h, U, V, W and T are the unknown quantities, also known as computational variables.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 12
Fx and Fy are source terms denoting the wind, the Coriolis force and the bottom friction (or any
other process being modelled by similar formulas). Several tracers can be taken into account
simultaneously. They can be of two different kinds, either active, i.e. influencing the flow by
changing the density, or passive, without any effect on density and then on flow.
The TELEMAC-3D basic algorithm can be split up in three computational steps (three fractional
steps).
The first step consists in finding out the advected velocity components by only solving the
advection terms in the momentum equations.
The second step computes, from the advected velocities, the new velocity components taking into
account the diffusion terms and the source terms in the momentum equations. These two solutions
enable to obtain an intermediate velocity field.
The third step is provided for computing the water depth from the vertical integration of the
continuity equation and the momentum equations only including the pressure-continuity terms (all
the other terms have already been taken into account in the earlier two steps). The resulting twodimensional equations (analogous to the Saint-Venant equations without diffusion, advection and
source terms) are written as:
∂h ∂(uh) ∂ (vh )
+
+
=0
∂t
∂x
∂y
∂Z
∂u
= −g S
∂t
∂x
∂Z
∂v
= −g S
∂t
∂y
The u and v in lower case denote the two-dimensional variables of the vertically integrated
velocity.
These two-dimensional equations are solved by the libraries in the TELEMAC-2D code and enable
to obtain the vertically averaged velocity and the water depth.
The water depth makes it possible to recompute the elevations of the various mesh points and then
those of the free surface.
Lastly, the computation of the U and V velocities is simply achieved through a combination of the
equations linking the velocities. Finally, the vertical velocity W is computed from the continuity
equation.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
2.2.2.
NON-HYDROSTATIC
Page 13
NAVIER-STOKES EQUATIONS
The following system (with an equation for
solved:
W which is similar to those for U and V ) is then to be
∂U ∂V ∂W
+
+
=0
∂x
∂y
∂z
1 ∂p
∂U
∂U
∂U
∂U
+U
+V
+W
=−
+ υ ∆(U ) + Fx
∂t
∂x
∂y
∂z
ρ ∂x
1 ∂p
∂V
∂V
∂V
∂V
+U
+V
+W
=−
+ υ ∆(V ) + Fy
∂t
∂x
∂y
∂z
ρ ∂y
∂W
∂W
∂W
∂W
1 ∂p
+U
+V
+W
=−
− g + υ ∆(W ) + Fz
∂t
∂x
∂y
∂z
ρ ∂z
In order to share a common core as much as possible with the solution of the equations with the
hydrostatic pressure hypothesis, the pressure is split up into a hydrostatic pressure and a
"dynamic" pressure term.
ZS
∆ρ
z
ρ0
p = p atm + ρ 0 g (Z S − z ) + ρ 0 g ∫
dz + p d
The TELEMAC-3D algorithm solves a hydrostatic step which is the same as in the previous
paragraph, the only differences lying in the continuity step ("projection" step in which the dynamic
pressure gradient changes the velocity field in order to provide the required zero divergence of
velocity) and the computation of the free surface.
2.2.3.
THE LAW OF STATE
Two laws of state can be used by default through TELEMAC-3D.
In most of the simulations, salinity and temperature make it possible to compute the variations of
density. The first law expresses the variation of density from only these two parameters. The
second law is more general and enables to construct all the variations of density with the active
tracers being taken into account in the computation.
The first law is written as:
[
ρ = ρ ref 1 − (T (T − Tref ) − 750S )10−6
With
2
]
Tref as a reference temperature of 4°C and ρ ref as a reference density at that temperature
when the salinity S is zero, then
ρ ref
= 999.972 kg/m3. That law remains valid for 0°C < T < 40°C
and 0 g/L < S < 42 g/L.
The second law is written as:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 14
ρ = ρ ref 1 − ∑ β i (Ti − Ti 0 )i
ρ ref ,
the reference density can be modified by the user together with the volumetric expansion
coefficients
2.2.4.
i
βi
related to the tracers
Ti .
K-ε MODEL
The turbulent viscosity can be given by the user, as determined either from a mixing length model
or from a k-ε model the equations of which are:
∂k
∂k
∂k
∂k
+U
+V
+W
=
∂t
∂x
∂y
∂z
∂ε
∂ε
∂ε
∂ε
+U
+V
+W
=
∂t
∂x
∂y
∂z
∂ νt
∂x σ k
∂ νt
∂x σ ε
+ C lε
wherein:
k=
∂k ∂ ν t
+
∂x ∂y σ k
∂ε ∂ ν t
+
∂x ∂y σ ε
ε
∂k ∂ ν t
+
∂y ∂z σ k
∂ε ∂ ν t
+
∂y ∂z σ ε
[P + (1 − C3ε )G ] − C 2ε ε
k
∂k
+ P−G −ε
∂z
∂ε
∂z
2
k
1 ' '
u i ui denotes the turbulent kinetic energy of the fluid,
2
ui' = U i − u i denotes the ith component of the fluctuation of the velocity U(U, V, W),
∂u i' ∂ui'
ε =υ
is the dissipation of turbulent kinetic energy,
∂x j ∂x j
P is a turbulent energy production term,
G is a source term due to the gravitational forces,
∂U
∂U j
P =ν t i +
∂x j
∂x i
and
νt
∂U i
∂x j
verifies the equality: ν t = C µ
k2
ε
G=−
ν t g ∂ρ
Prt ρ ∂z
,
,
C µ , Prt , C1ε , C 2ε , C3ε , σ k , σ ε are constants in the k-ε model.
2.2.5.
EQUATIONS OF TRACERS
The tracer can be either active (it affects hydrodynamics) or passive in TELEMAC-3D.
Temperature, salinity and in some cases a sediment are active tracers. The tracer evolution
equation is formulated as:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 15
∂T
∂T
∂T
∂T
∂ ∂T ∂ ∂T ∂ ∂T
+ ν T
+U
+V
+W
= ν T
+ ν T
+Q
∂t
∂x
∂y
∂z ∂x ∂x ∂y ∂y ∂z ∂z
with:
•
T
•
νT
(m2/s) tracer diffusion coefficients,
•
t
(s)
time,
•
x, y, z
(m)
space components,
•
Q
(tracer unit)
tracer either passive or affecting the density,
(tracer(s) unit) tracer source or sink.
2.3.
MESH DESCRIPTION
2.3.1.
THE MESH STRUCTURE
The TELEMAC-3D mesh structure is made of prisms (possibly split in tetrahedrons). In order to
prepare that mesh of the 3D flow domain, a two-dimensional mesh comprising triangles which
covers the computational domain (the bottom) in a plane is first constructed, as for TELEMAC-2D.
The second step consists in duplicating that mesh along the vertical direction in a number of curved
surfaces known as "planes". Between two such planes, the links between the meshed triangles
make up prisms.
The computational variables are defined at each point of the three-dimensional mesh, inclusive of
bottom and surface. Thus, they are "three-dimensional variables” except, however, for the water
depth and the bottom depth which are obviously defined only once along a vertical. Thus, they are
"two-dimensional variables". Some TELEMAC-3D actions are then shared with TELEMAC-2D and
use the same libraries, such as the water depth computation library. Therefore, it is well understood
that TELEMAC-3D should manage a couple of mesh structures: the first one is two-dimensional
and is the same as that used by TELEMAC-2D, and the second one is three-dimensional. That
implies managing two different numberings a detailed account of which is given below.
2.3.2.
THE TWO-DIMENSIONAL MESH
The two-dimensional mesh, which is made of triangles, can be prepared using a mesh generator
software compatible with the TELEMAC system (MATISSE, BLUE KENUE, JANET…).
Using a mesh generator that does not belong to the TELEMAC chain involves converting the
resulting file, through the STBTEL interface for instance, to the SERAFIN format, which can be
read by TELEMAC as well as by the RUBENS, BLUE KENUE or FUDAA-PREPRO postprocessors. In addition, STBTEL checks such things as the proper orientation of the local
numbering of the mesh elements.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 16
The two-dimensional mesh (included in the GEOMETRY FILE) consists of NELEM2 elements and
NPOIN2 vertices of elements which are known through their X, Y, and Zf co-ordinates (the
BOTTOM variable). Each element is identified by a code known as IELM2 and includes NDP nodes
(3 for a triangle with a linear interpolation). The nodes on an element are identified by a local
number ranging from 1 to NDP. The link between that element-wise numbering (local numbering)
and the mesh node numbering ranging from 1 to NPOIN2 (global numbering) is made through the
connectivity table IKLE2. The global number of the IDP local-number node in the IELEM2 element
is IKLE2(IELEM2,IDP).
2.3.3.
THE THREE-DIMENSIONAL MESH
The three-dimensional mesh, which is made of prisms, is automatically constructed by TELEMAC3D from the previous mesh. The data in the three-dimensional mesh of finite elements are as
follows:
•
NPOIN3: the number of points in the mesh (NPOIN3 = NPOIN2 × NPLAN),
•
NELEM3: the number of elements in the mesh,
•
NPLAN: the number of planes in the mesh,
•
X, Y, Z: NPOIN3 dimensional arrays. X and Y are obtained by merely duplicating the abovedescribed arrays of the two-dimensional mesh. Dimension Z obviously depends on the
mesh construction being selected (keyword MESH TRANSFORMATION),
•
IKLE3: dimensional arrays (NELEM3,6). IKLE3(IELEM3,IDP) provides with the global
number of the IDP point in the IELEM3 element. IKLE3 defines a numbering of the 3D
elements and a local numbering of the points in each element, it provides for the transition
from that local numbering to the global numbering.
From these data, TELEMAC-3D constructs other arrays, such as the edge point global address
array.
Figure 1 herein below illustrates a TELEMAC-3D three-dimensional mesh.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 17
Figure 1: A three-dimensional mesh.
3.
THE INPUTS / OUTPUTS
3.1.
PRELIMINARY COMMENTS
In a computation, the TELEMAC-3D code uses a number of input and output files, some of which
are optional. Most of these files are similar or identical to their counterparts in TELEMAC-2D.
The input files are:
• The steering file (mandatory), which contains the "configuration" of the simulation,
• The geometry file (mandatory), which contains the information regarding the mesh,
• The boundary conditions file (mandatory), which contains the description of the type of
each boundary,
• The FORTRAN file, which contains the specific subroutines of the simulation (modified
TELEMAC subroutine or specifically created),
• The bottom topography file, which contains the description of the bottom topography.
Generally, the topographic data are already contained in the geometry file and the bottom
topography file is generally not used,
• The liquid boundaries file, which contains the information on the imposed values at liquid
boundaries,
• The previous computation file, which provides the initial state of the calculation in the case
of a restart calculation,
• The reference file, which contains the calculation of "reference" used for the validation
process,
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
Page 18
The stage-discharge curves file, which contains the information on the imposed values at
liquid boundaries in case of height / flow rate law,
The sources file, which contains the information regarding the sources.
The output files are:
• The 3D result file, which contains the graphical results associated to the 3D mesh,
• The 2D result file, which contains the graphical results associated to the 2D mesh,
• The output listing, which is simulation report. In case of difficulty in performing a
calculation, the user can request the printing of additional information by activating the
logical keyword DEBUGGER. DEBUGGER = 1 provides the calling sequences of
subroutines from the main program telemac3d.f. This technique is useful in case of
critical computation crash to identify the responsible subroutine,
• The file for Scope, which is an additional test file available for the user,
• The restart file, allowing to perform a restart computation without information loss.
In addition, the user may have to manage additional files:
• The binary data files 1 and 2, in input,
• The formatted data files 1 and 2, in input,
• The binary results file, in output,
• The formatted results file, in output.
3.1.1.
BINARY FILES FORMAT
The binary files used within the TELEMAC system can have various formats. The most common
format is the SERAFIN format (also known as SELAFIN, for misidentified historical reasons) which
is the standard internal TELEMAC system format (described in Appendix N° 3). This SERAFIN
format can be configured so as to store the real-typed values in single or double precision. The
other available format is the MED format which is used by the SALOME platform jointly developed
by EDF and the French atomic energy commission (CEA). In the current TELEMAC release, this
format is restricted to EDF internal use.
Depending on the specified format, the binary files may be treated by different software. However,
in the current release of TELEMAC, only the SERAFIN single precision format can be read by the
post-processing tools like RUBENS, FUDAA-PREPRO or BLUE KENUE.
The selection of the format of a file is done by the corresponding keyword. Thereby, the keyword
GEOMETRY FILE FORMAT specifies the format of the geometry file. Those keywords may take 3
different values (as an 8 character string): the value ‘SERAFIN ‘ corresponds to the standard single
precision SERAFIN format, which is the default and recommended value (do not forget the space
in last position). The ‘SERAFIND‘ value corresponds to the double precision SERAFIN format
which allows to increase the precision of the results, especially in case of a restart or reference file.
Finally, the ‘MED‘ value corresponds to the hdf5 MED format (internal EDF use), double precision.
The keywords involved are:
• 2D RESULT FILE FORMAT,
• 3D RESULT FILE FORMAT,
• BINARY DATA FILE 1 FORMAT,
• FILE FOR 2D CONTINUATION FORMAT,
• GEOMETRY FILE FORMAT,
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
3.2.
Page 19
PREVIOUS COMPUTATION FILE FORMAT,
REFERENCE FILE FORMAT,
RESTART FILE FORMAT.
THE STEERING FILE
The steering file contains all the data about the selection of computational options (physical,
numerical, etc.). It is an ASCII file which can be generated either through the FUDAA-PREPRO
software or directly using a text editor. In a way, it serves as the computation dashboard. It includes
a set of keywords to which values are assigned. If a keyword does not occur in that file, then
TELEMAC-3D will assign it the default value as defined in the dictionary file. If such a default value
is not defined in the dictionary, then the computation will be interrupted with an error message. For
instance, the command TIME STEP = 10.0 specifies that the computation time step value is 10
seconds.
TELEMAC-3D reads the steering file at the beginning of the computation.
Both dictionary file and steering file are read using the DAMOCLES library which is included in the
TELEMAC chain. It is therefore necessary, when generating the steering file, to observe the
DAMOCLES syntax rules (what is performed automatically if the file is generated using FUDAAPREPRO).
The writing rules are as follows:
•
The keywords can be of the Integer, Real, Logical or Character type,
•
The sequence order in the steering file is of no importance,
•
Several keywords can be on the same line,
•
Each line cannot exceed 72 characters. However, one can start a new line as many times
as one wishes provided that the keyword name is not astride two lines,
•
For the array-types keywords, the character separating successive values is the
semicolon. For example:
PRESCRIBED FLOWRATES = 10.0; 20.0
•
The symbols ":" and "=" can both be used in order to separate a keyword name and its
value. They can be either preceded or followed by any number of blanks. The value itself
can occur on the following line. For example:
TIME STEP = 10.
or
TIME STEP : 10.
Or else
TIME STEP =
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 20
10
•
The characters occurring between two "/" on one line are taken as comments. Likewise,
those characters occurring between a "/" and a line ending are considered as comments.
For example:
VERTICAL TURBULENCE MODEL = 3
•
/ k-epsilon model
A line beginning with a "/" in a first column is wholly considered as a comment, even
though there is another ¨/¨ on the line. For example:
/ The geometry file is ./mesh/geo
•
Writing the integers: do not exceed the maximum size allowed by the machine (in a
machine with a 32 bit architecture, the extreme values range from -2 147 483 647 to
+2 147 483 648. Do not insert any blank between the sign (optional for the sign +) and the
number. A point after the end of the number is tolerated,
•
Writing the reals: the point and the comma are accepted as a decimal separator, as well as
the FORTRAN formats E and D. ( 1.E-3 0.001 0,001 1.D-3 represent the same value),
•
Writing the logical values: the values 1, OUI, YES, .TRUE., TRUE, or VRAI on the one
hand, and 0, NON, NO, .FALSE., FALSE, or FAUX on the other hand, are accepted,
•
Writing the character strings: the strings including blanks or reserved symbols ("/",":", "=",
"&") should be inserted between quotes ('). The value of a character keyword may contain
up to 144 characters. As in FORTRAN, the quotes included in a string should be doubled.
A string may not begin or end with a blank. For example:
TITLE = 'COASTAL ENVIRONMENT STUDY'
In addition to the keywords, a number of pseudo instructions or metacommands which are
interpreted during the sequential reading of the steering file can also be used:
RELEASE 7.1
•
The &FIN command specifies the file end (even though the file is not over). Thus, some
keywords can be disabled simply by placing them behind that command so that they can
easily be reactivated later on. However, the computation keeps running,
•
The &ETA command prints the list of keywords and the values which are assigned to them
when DAMOCLES meets that command. That display will take place at the head of the
output listing,
•
The &LIS command prints the list of keywords. That display will take place at the head of
the output listing,
•
The &IND command prints the detailed list of keywords. That display will take place at the
head of the output listing,
•
The &STO command causes the program to be halted, whereas the computation does not
keep running.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
3.3.
Page 21
THE GEOMETRY FILE
It is the same file as that used by TELEMAC-2D. It is a SERAFIN-formatted binary file which can
then be read by RUBENS or FUDAA-PREPRO and is generated either by MATISSE, BLUE
KENUE or JANET or by the STBTEL module (from the file(s) originating from of mesh generator).
The SERAFIN format structure is described in Appendix N° 3.
That file contains all the data about the two-dimensional mesh (see in paragraph 2.3.2). It includes
the number of points in the mesh (NPOIN2 variable), the number of elements (NELEM2 variable),
the number of vertices per element (NDP variable), the X and Y arrays containing the co-ordinates
of all the points and, lastly, the IKLE array containing the connectivity table.
That file may also contain bathymetric and bottom friction data for each point in the mesh.
NOTE: TELEMAC-3D retrieves the geometry data at the beginning of the 2D result file. Any
computational 2D result file can therefore be used as a geometry file when a further simulation on
the same mesh is desirable.
That file name is provided using the keyword: GEOMETRY FILE and its format is specified by the
keyword: GEOMETRY FILE FORMAT (default value: ‘SERAFIN ’), see paragraph 3.1.1.
3.4.
THE BOUNDARY CONDITIONS FILE
It is the same file as that used by TELEMAC-2D. It is a formatted file which is generated
automatically by MATISSE or STBTEL and can be amended using FUDAA-PREPRO or a text
editor. Each line in that file is dedicated to a point at the 2D mesh boundary. The edge point
numbering is that of the file lines; it first describes the domain outline in the counter clockwise
direction from the bottom left-hand point (point where the sum of horizontal co-ordinates of which is
minimum), then the islands in the clockwise direction.
For a thorough description of that file, refer to the specific paragraph 4.3.4.
That name of file is provided using the keyword: BOUNDARY CONDITIONS FILE.
3.5.
THE FORTRAN FILE
The FORTRAN file may include a number of subroutines (so-called "user" subroutines) available
under the TELEMAC-3D tree structure which the user can modify as well as those subroutines
which have been specifically developed for the computation.
The user subroutines from the various libraries used by TELEMAC-3D are listed in Appendix N° 2.
Since TELEMAC-3D is available in Open Source, every subroutine can be freely used. Every user
subroutine copied into the user FORTRAN file is automatically substituted for the same named
subroutine occurring in the TELEMAC-3D compiled libraries.
Upon the creation and every amendment of the FORTRAN file, a new executable program is
generated (compilation and link) for the simulation.
That file name is provided using the keyword: FORTRAN FILE.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
3.6.
Page 22
THE LIQUID BOUNDARIES FILE
It is an ASCII file enabling the user to specify time-varying boundary conditions values (flow rate,
depth, velocity, tracer concentration) at all the liquid boundaries. That file can be generated under
the FUDAA-PREPRO software interface.
For a thorough description of that file, refer to the specific paragraph 4.3.8.
That name of file is provided using the keyword: LIQUID BOUNDARIES FILE.
3.7.
THE PREVIOUS COMPUTATION FILE
It is a TELEMAC-3D result file which is used for initializing a new computation. In order to activate
the optional use of that file, the keyword COMPUTATION CONTINUED should be activated. In
order to specify the previous computation file, its name should be stated through the keyword:
PREVIOUS COMPUTATION FILE. The initial conditions of the new computation are defined by the
last recorded time step in the previous computation file. The whole set of data from the steering file
is read and it makes possible to redefine or amend the variables (time step, turbulence model,
addition or deletion of a tracer…).
A computation can also be initialized from a TELEMAC-2D result. In order to activate that option,
the 2D CONTINUATION keyword should be validated. The TELEMAC-2D result file should then be
associated with the FILE FOR 2D CONTINUATION keyword (default value: ‘SERAFIN ’), see
paragraph 3.1.1.
3.8.
THE REFERENCE FILE
During the validation step of a calculation, this file contains the reference result. At the end of the
calculation, the result of the simulation is compared to the last time step stored in this file. The
result of the comparison is given in the control printout in the form of a maximum difference in
elevation and the three components of velocity.
The name of this file is given by the keyword: REFERENCE FILE and its format is specified by the
keyword: REFERENCE FILE FORMAT (default value: ‘SERAFIN ’), see paragraph 3.1.1.
3.9.
THE STAGE-DISCHARGE CURVES FILE
This text file enables the user to configure the evolution of the prescribed value on specific open
boundaries. This file is used when the prescribed elevation is determined by a discharge-elevation
law. The descriptions of the appropriate laws are given through this file.
See paragraph 4.3.10 for a complete description of this file.
The name of this file is specified with the keyword: STAGE-DISCHARGE CURVES FILE.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 23
3.10. THE SOURCES FILE
This text file enables the user to specify values for time-dependent conditions for sources
(discharge, tracers concentration).
See section 5.4 for a complete description of this file.
The file name is specified with the keyword: SOURCES FILE.
3.11. THE 3D RESULT FILE
It is the file into which TELEMAC-3D stores the information during the computation. It is a
SERAFIN-formatted file (refer to Appendix N° 3) or a MED-formatted file. It first contains all the
data about mesh geometry, then the names of the stored variables. It also contains, for each
graphic printout and for each mesh point, the values of the various recorded variables.
Its content varies according to the values of the following keywords:
NUMBER OF FIRST TIME STEP FOR GRAPHIC PRINTOUTS: provided to set from which time
step onwards the information will be stored, in order to prevent too large files, especially when the
computation begins with an uninteresting transient stage related to the definition of unrealistic initial
conditions (e.g. invariably zero currents). The default value is 0 (writing of the graphic printouts at
the beginning of the simulation).
GRAPHIC PRINTOUT PERIOD: sets the period (in number of time steps) of printouts in order to
prevent a too large file. Default value is 1 (writing at every time step). For instance:
TIME STEP
= 60.0
GRAPHIC PRINTOUT PERIOD
= 30
*The results will be backed up every 1,800th second, i.e. 30th minute.
VARIABLES FOR 3D GRAPHIC PRINTOUTS: keyword specifying the list of variables which will be
stored in the result file. Each variable is identified by means of a name from the list below.
• Z
elevation (m),
• U
velocity along x axis (m/s),
• V
velocity along y axis (m/s),
• W
velocity along z axis (m/s),
• TA concentrations for tracers (TA1 for the 1st one, TA2 for the 2nd one…),
• NUX viscosity for U and V along x axis (m2/s),
• NUY viscosity for U and V along y axis (m2/s),
• NUZ viscosity for U and V along z axis (m2/s),
• NAX viscosity for tracers along x axis (m2/s),
• NAY viscosity for tracers along y axis (m2/s),
• NAZ viscosity for tracers along z axis (m2/s),
• RI
Richardson number in case of a mixing length model,
• K
turbulent energy for k-ε model (J/kg),
• EPS dissipation of turbulent energy (W/kg),
• DP dynamic pressure (multiplied by DT/RHO),
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
•
•
•
PH
RHO
P1
P2
P3
P4
Page 24
hydrostatic pressure (Pa),
relative density,
private variable 1,
private variable 2,
private variable 3,
private variable 4.
That file name is provided using the keyword: 3D RESULT FILE and its format using 3D RESULT
FILE FORMAT (default value: ‘SERAFIN ’). Stored variables by default are ‘Z, U, V, W’.
Using the private arrays requires updating the keyword NUMBER OF PRIVATE ARRAYS in order
to perform the required memory allocations (see section 10.4).
3.12. THE 2D RESULT FILE
It is the file into which TELEMAC-3D stores the specifically two-dimensional data during the
computation (such as the free surface, the depth-averaged horizontal components of velocity and
the depth-averaged tracers). It has a SERAFIN format. The free surface and the horizontal
components of velocity will then physically correspond to the same data as those being supplied by
TELEMAC-2D. The obtained values, however, may be different from an analogue computation
being directly made with TELEMAC-2D when the flow is specifically three-dimensional.
Its content varies according to the values of the following keywords:
NUMBER OF FIRST TIME STEP FOR GRAPHIC PRINTOUTS: same keyword as that described
in section 3.11.
GRAPHIC PRINTOUT PERIOD: same keyword as that described in section 3.11.
VARIABLES FOR 2D GRAPHIC PRINTOUTS: keyword specifying the list of variables which will be
stored in the result file. Each variable is identified by means of a name from the list below.
• U
average velocity along x axis (m/s),
• V
average velocity along y axis (m/s),
• C
celerity (m/s),
• H
water depth (m),
• S
free surface elevation (m),
• B
bottom elevation (m),
• TA averaged concentrations for tracers (TA1 for the 1st one, TA2 for the 2nd one…),
• F
Froude number,
• Q
scalar discharge (m2/s),
• I
discharge along x (m2/s),
• J
discharge along y (m2/s),
• M
norm of velocity (m/s),
• X
wind along x axis (m/s),
• Y
wind along y axis (m/s),
• P
air pressure (Pa),
• W
friction coefficient,
• RB non erodible bottom elevation (m),
• HD thickness of the sediment bed layer (m),
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
•
•
•
•
•
•
•
•
•
•
Page 25
EF erosion rate (kg/m2/s),
DF deposition flux (kg/m2/s),
DZF bed evolution,
PRIVE1
private array PRIVE 1,
PRIVE2
private array PRIVE 2,
PRIVE3
private array PRIVE 3,
PRIVE4
private array PRIVE 4,
QS
solid discharge (m2/s),
QSX
solid discharge along x (m2/s),
QSY
solid discharge along y (m2/s),
US
friction velocity (m/s),
MAXZ
maximum value of the water elevation during the computation (m),
TMXZ
time corresponding to this maximum elevation (s).
That file name is provided by means of the keyword: 2D RESULT FILE and its format using 2D
RESULT FILE FORMAT (default value: ‘SERAFIN ’). Stored variables by default are ‘U, V, H, B’.
Using the private arrays requires updating the keyword NUMBER OF 2D PRIVATE ARRAYS in
order to perform the required memory allocations (see section 10.4).
3.13. THE OUTPUT LISTING
It is a formatted file which can be created by TELEMAC-3D during the computation (program
launched with the –s option). It contains the report of a TELEMAC-3D running. Its contents vary
according to the values of the following keywords:
NUMBER OF FIRST TIME STEP FOR LISTING PRINTOUTS: provided to set at which time step it
is desired to begin printing the data, in order to prevent too large files.
LISTING PRINTOUT PERIOD: sets the period between two time step printings. The value is given
as a time step number. For instance, the following sequence:
TIME STEP = 30.0
LISTING PRINTOUT PERIOD = 2
Prints the output listing every minute of simulation.
MASS-BALANCE: if it is requested, the user will get information about the mass flow (or rather the
volumes) and the errors (primarily linked to the precision achieved by the solvers) of that
computation in the domain. This is not done by default.
INFORMATION ABOUT MASS-BALANCE FOR EACH LISTING PRINTOUT: if they are requested
(that is done by default), the user will get, at each time step, information about the flows within the
domain.
The file name is directly managed by the TELEMAC-3D launching procedure. Generally, the file
has a name which is created from the name in the steering file and the number of the processors
used in the computation, followed by the suffix ".sortie".
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 26
3.14. THE FILE FOR SCOPE
TELEMAC-3D gives the user an opportunity to retrieve computation variables along a 1D profile
during a computation. These profiles will then be stored in a dedicated file the name of which is
defined by the keyword FILE FOR SCOPE. For that purpose, the user should program the SCOPE
subroutine.
That subroutine is provided to define profiles of computed variables, or other variables, as created
by a user along a segment with (X1, Y1, Z1) and (X2, Y2, Z2) co-ordinates. The user also decides
the number of points as distributed along that segment. The data is automatically saved as per the
SCOPE format at all the time steps (for a description of this format, see Appendix of RUBENS
Version 4.0. Manuel Utilisateur, O. Quiquempoix, EDF report HE-45/95/031/B).
3.15. THE AUXILIARY FILES
Other files can be used by TELEMAC-3D:
• One or two binary data files, as specified by the keywords BINARY DATA FILE 1 and
BINARY DATA FILE 2. These files can be used to provide data to the program, the user
has to handle their reading within the FORTRAN FILE. The data from these files shall be
read using the T3DBI1 logic unit for binary data file 1 and the T3DBI2 logic unit for binary
data file 2,
• One or two formatted data files, as specified by the keywords FORMATTED DATA FILE 1
and FORMATTED DATA FILE 2. These files can be used to provide data to the program,
the user has to handle their reading within the FORTRAN FILE. The data from these files
shall be read using the T3DFO1 logic unit for formatted data file 1 and the T3DFO2 logic
unit for formatted data file 2,
• A binary results file, as specified by the keyword BINARY RESULTS FILE. This file can be
used to store additional results, the user has to handle their writing within the FORTRAN
program using the T3DRBI logic unit.
• A formatted results file, as specified by the keyword FORMATTED RESULTS FILE. This
file can be used to store additional results (for instance, results readable by an external 1D
simulation code used for a coupling with TELEMAC), the user has to handle their writing
within the FORTRAN program using the T3DRFO logic unit.
The read or write operations from/into these files should be thoroughly managed by the user (the
files are opened and closed by the program). That management can be performed from any point
which the user can gain access to. The logic unit numbers are stated in the
DECLARATIONS_TELEMAC3D module and the user can access them through a USE
DECLARATIONS_TELEMAC3D command at the beginning of a subroutine. For instance, using a file
for providing the initial conditions will result in managing that file within the CONDIM subroutine.
Likewise, using a file for inserting boundary conditions will be possible at the BORD3D subroutine. In
case of conflicting statements, one can use, for example: DECLARATIONS_TELEMAC3D, ONLY :
T3DBI1.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 27
3.15.1. THE DICTIONNARY FILE
This is a file containing all the information on the keywords (French name, English name, default
value, type, keywords documentation). This file can be read by the user, but it should not be
modified in any case.
3.16. TOPOGRAPHIC AND BATHYMETRIC DATA
Topographical and bathymetric data may be supplied to TELEMAC-3D at three levels:
• Either directly in the geometry file by a topographical or bathymetric value associated with
each mesh node. In this case, the data are processed while the mesh is being built using
MATISSE, JANET or BLUE KENUE, or when the STBTEL module is run before
TELEMAC-3D is started. STBTEL reads the information in one or more bottom topography
files (5 at most) and interpolates at each point in the domain,
• Or in the form of a cluster of points with elevations that have no relation with the mesh
nodes, during the TELEMAC-3D computation. TELEMAC-3D then makes the interpolation
directly with the same algorithm as STBTEL. The file name is provided by the keyword
BOTTOM TOPOGRAPHY FILE. Unlike STBTEL, TELEMAC-3D only manages one single
bottom topography file. This may be in SINUSX format or more simply a file consisting of
three columns X, Y, Z. The SINUSX format is described in the RUBENS user manual,
• Or using the T3D_CORFON subroutine. This is usually used for schematic test cases.
In all cases, TELEMAC-3D offers the possibility of smoothing the bottom topography in order to
obtain a more regular geometry. The smoothing algorithm can be iterated several times depending
on the degree of smoothing required. The keyword NUMBER OF BOTTOM SMOOTHINGS then
defines the number of iterations carried out in the T3D_CORFON subroutine. The default value of
this keyword is 0.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
4.
Page 28
GENERAL SETUP OF THE HYDRODYNAMIC
COMPUTATION (NAVIER-STOKES EQUATIONS)
The general setup of the computation is only performed at the steering file level.
The time data is provided by the two keywords TIME STEP (real, set at 1. by default) and
NUMBER OF TIME STEPS. The first keyword sets the period of time between two consecutive
computational moments (but not necessarily two outputs in the result file). The global duration of
the computation is provided through a number of time steps (keyword NUMBER OF TIME STEPS,
set at 1 by default) or a duration in seconds (keyword DURATION, set at 0. by default). If both are
given, TELEMAC-3D follows the instruction leading to the longest computation.
Both date and hour corresponding to the initial time of the computations can be specified using the
keywords ORIGINAL DATE OF TIME (AAAA, MM, JJ format; default value = 1900; 1; 1) and
ORIGINAL HOUR OF TIME (HH, MM, SS format, default value = 0; 0; 0). These two data are
mandatory if using the tidal data bases. They can be taken into account in programming by means
of the MARDAT and MARTIM variables.
The computation title is specified by the keyword TITLE.
4.1.
MESH DEFINITION
The three-dimensional mesh, consisting of prisms possibly cut into tetrahedrons, is automatically
constructed by TELEMAC-3D from the two-dimensional mesh. This construction is done in the
subroutine CALCOT from the information given by the subroutine defining the initial conditions
CONDIM.
The number of prisms is specified in the data file by means of the keyword NUMBER OF
HORIZONTAL LEVELS (default value = 2). That number of levels is equivalent to the number of
stacked prisms plus 1. Its minimum value is 2 (1 prism in the vertical direction).
TELEMAC-3D uses a change of variables in order to freeze the mesh on a time step (without such
a change, the mesh dimensions z would vary in accordance with the free surface evolution). The
frequently adopted change of variables is the sigma transform which consists in shifting from the z
(x,y,t) co-ordinate to the z* (x,y) co-ordinate. The user should enter the z* co-ordinates in the
CONDIM subroutine. The normalized co-ordinates will then range from 0 (the bottom) to 1 (the
surface).
The numbering of levels is made according to upward vertical. Level 1 follows the bottom and level
N corresponds to the free surface (N being specified by the keyword NUMBER OF HORIZONTAL
LEVELS).
The vertical mesh definition is based on the TRANSF_PLANE table which allows defining the
behaviour of each level.
The keyword MESH TRANSFORMATION sets the kind of level distribution along the vertical. The
value 0 corresponds to a distribution directly defined by the user in the CALCOT subroutine.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 29
Figure 2: Effect of the MESH TRANSFORMATION keyword – Value 1: sigma.
The default value is 1 (Figure 2) and results in a homogeneous distribution of levels in the vertical
direction (classical sigma transformation). The height of the levels varies depending on the water
depth, all planes can move (except for the bottom). In this case, no programming in the CONDIM
subroutine is required.
A value of 2 (Figure 3) will allow the user to define the distribution of levels (e.g. refinement near
surface) while maintaining the levels mobility (sigma transformation with given proportions). The
latter choice implies that the user will program his/her distribution in the CONDIM subroutine to
define the array ZSTAR that describes the distribution of levels along the vertical as a percentage of
the water depth. Changes to make are:
• Specifying the variable TRANSF_PLANE with a value of 2 for every level,
• Specifying the level distribution along the vertical through the array ZSTAR which describes
the distribution along the vertical as a percentage of the water depth (the values are
between 0. and 1.).
For example (Figure 3):
DO IPLAN = 1,NPLAN
TRANSF_PLANE%I(IPLAN)=2
ENDDO
ZSTAR%R(1)=0.D0
ZSTAR%R(2)=0.02D0
ZSTAR%R(3)=0.1D0
ZSTAR%R(4)=0.4D0
ZSTAR%R(5)=0.8D0
Figure 3: Effect of the MESH TRANSFORMATION keyword – Value 2: zstar.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 30
In order to better represent the densimetric stratification areas (thermoclines, halocline and/or
outfall), prescribing a maximum number of horizontal "levels" (particularly those where the
gradients are the highest) is sometimes suitable. For that purpose, the user can select the value 3
for the keyword MESH TRANSFORMATION. In this configuration, the user can freely use the 3
types of level definition available in TELEMAC-3D to create the mesh along the vertical:
• Fixed levels at a given altitude (correspond to value 3 of the TRANSF_PLANE variable, the
altitude is specified by the ZPLANE variable),
• Irregularly distributed movable levels between two fixed levels (correspond to value 2 of
the TRANSF_PLANE variable, the distribution is specified by the ZSTAR variable),
• Evenly distributed movable levels between two fixed levels (correspond to value 1 of the
TRANSF_PLANE variable).
This latter choice requires the user to make changes in the CONDIM subroutine:
• Specifying the variable TRANSF_PLANE at value 1, 2 or 3 for each level,
• Specifying the level distribution along the vertical through the array ZSTAR for levels of
type 2,
• Specifying the level altitude through the array ZPLANE for levels of type 3.
For example (Figure 4):
DO IPLAN = 1,5
TRANSF_PLANE%I(IPLAN)=2
ENDDO
ZSTAR%R(1)=0.D0
ZSTAR%R(2)=0.2D0
ZSTAR%R(3)=0.5D0
ZSTAR%R(4)=0.7D0
ZSTAR%R(5)=0.8D0
DO IPLAN = 7,NPLAN
TRANSF_PLANE%I(IPLAN)=1
ENDDO
TRANSF_PLANE%I(6)=3
ZPLANE%R(6)=0.D0
Figure 4: Effect of the MESH TRANSFORMATION keyword – Value 3: user defined.
Various programming examples are provided as comments in the CONDIM subroutine.
The triangular based prismatic elements can optionally be split into tetrahedrons. This option is
enabled using the ELEMENT keyword can take the value 'PRISM' (default value) or
'TETRAHEDRON'.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
4.2.
Page 31
PRESCRIBING THE INITIAL CONDITIONS
The initial conditions aim at defining the model condition at the beginning of the simulation.
In case of a continuing computation, the initial conditions are provided at one time step in the result
file of the previous computation (refer to section 3.7). The mandatory variables (at least the velocity
components) when resuming the computation should then have been stored into the file being
used as PREVIOUS COMPUTATION FILE.
Otherwise, the default initial condition is defined as follows:
•
Free surface set to an elevation equal to 0,
•
Zero velocities,
•
Steady zero active and passive tracers.
If that initial condition is not suitable for a computation, then it should be changed using keywords
in the simple cases or through a programming as described in the subsequent paragraphs.
4.2.1.
PRESCRIPTION THROUGH KEYWORDS
In all the cases, the kind of initial conditions is set by the keyword INITIAL CONDITIONS. That
keyword can have one of the following six values:
RELEASE 7.1
•
'ZERO ELEVATION': Initializes the free surface elevation to 0 (default value). The initial
water depths are then computed from the bottom depth,
•
'CONSTANT ELEVATION': Initializes the free surface elevation to the values as provided
by the keyword INITIAL ELEVATION. The initial water depths are then computed by
getting the difference between the free surface elevation and the bottom depth. In those
areas where the bottom depth exceeds the initial elevation, the initial water depth is zero,
•
'ZERO DEPTH': All the water depths are initialized with a zero value (free surface
coinciding with bottom). In other words, the whole domain is "dry" at the beginning of the
computation,
•
'CONSTANT DEPTH': Initializes the water depths to the value as provided by the keyword
INITIAL DEPTH,
•
‘TPXO SATELLITE ALTIMETRY’: The initial conditions are set using information provided
by the OSU harmonic constants database (TPXO for instance) in the case of the use of
this database for the imposition of maritime boundary conditions (see paragraph 4.3.14),
•
'PARTICULAR' or ‘SPECIAL’: The initial conditions are defined as programmed by the user
in the CONDIM subroutine (refer to the next paragraph). That procedure should be used
whenever the initial conditions of the model do not correspond to one of above four cases.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
4.2.2.
Page 32
PRESCRIBING PARTICULAR INITIAL CONDITIONS (PROGRAMMING THE CONDIM
SUBROUTINE)
The CONDIM subroutine should be programmed whenever the initial conditions programmed by
default are to be modified.
By default, the standard version of the CONDIM subroutine interrupts the computation if the
keyword INITIAL CONDITIONS is set to 'PARTICULAR' without any actual amendment of the
subroutine.
The CONDIM subroutine successively initializes the two-dimensional variables, then the threedimensional variables:
•
The water depth,
•
The 3D component of velocities,
•
The active and passive tracers.
The user can quite freely fill that subroutine. For instance, he/she can retrieve information in a
formatted or binary file, using the corresponding keywords.
4.2.3.
RESUMING THE COMPUTATION
TELEMAC-3D enables the user to resume a computation by taking as the initial condition the last
time step of a computation which was previously computed on the same mesh, or possibly with a
different number of levels. Thus, some computational parameters such as the time step, some
boundary conditions, the turbulence model can be modified, or else a computation can be initiated
once a steady state is achieved.
The file to be retrieved shall then inevitably contain all the data required for TELEMAC-3D, i.e. not
only the co-ordinates of the X, Y and Z computational points which it necessarily contains, but also
the 3D velocities and the tracers.
If some variables do not appear in the PREVIOUS COMPUTATION FILE, then they are
automatically set to zero values. A usual application consists in using the result of a hydrodynamic
computation in order to perform a tracer transport computation. Generally, the PREVIOUS
COMPUTATION FILE does not include any result for the tracer.
To resume a computation, it is required to use two keywords into the steering file.
The keyword COMPUTATION CONTINUED should be set to the YES value.
The keyword PREVIOUS COMPUTATION FILE should provide the name of the file which will
provide the initial conditions (default value = 0).
Optionally, the keyword RECORD NUMBER FOR RESTART can be used to define the record
number to read if it is not the last one (defined by a default value set to 0).
WARNING: the two-dimensional mesh on which the useful results have been computed should be
strictly identical to the mesh of the case to be handled.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 33
Resuming the computation usually leads to small differences in results compared to the same
calculation without interruption. This difference is mainly due to the fact that the velocity advection
is not treated properly at the first time step, because this operation requires information from the
previous time step. To correct this, the user has a specific recovery procedure to improve the
accuracy of calculations, using double precision format SERAFIN files:
• In the first computation, the keyword RESTART MODE is set to YES, which generates a
specific file containing the full information only of the last time step of the simulation (in
particular information on the advection field of the last time step). The name of this file is
specified using the keyword RESTART FILE,
• In the second computation, this specific file must be used as PREVIOUS COMPUTATION
FILE specifying the PREVIOUS COMPUTATION FILE FORMAT is 'SERAFIND' (SERAFIN
double precision). In this case, the keyword RECORD NUMBER FOR RESTART cannot
be used because the recovery file only contains the last time step of the simulation.
However, it has to be mentioned that even if it is not advisable, the creation of specific restart file
can be done not only SERAFIND format, but also with any other available format in the TELEMAC
system, especially in single precision. In this case, the keyword RESTART FILE FORMAT (by
default set at 'SERAFIND') must be set to the proper value.
A particular aspect of the resuming technique of computation is the value of the start time of the
second simulation. By default, the start time of the second calculation is equal to the value of the
last time step of restart file. This can be changed by using the logical keyword INITIAL TIME SET
TO ZERO if the user wants to start from zero (default value is NO).
It is also possible to resume a computation from a 2D results file. This is generally useful in river
hydraulic offering the possibility to initialise the model in 2D before shifting to 3D simulation. In this
case, the horizontal velocities are considered as constant on the vertical and equal to the 2D
velocities and the vertical velocities are initialised to zero. This possibility is activated using the 2D
CONTINUATION logical keyword (default value = NO). The 2D results file must be given using
FILE FOR 2D CONTINUATION. The keyword FILE FOR 2D CONTINUATION FORMAT gives the
format of the file and can takes the following values: ‘SERAFIN ’ (default value), ‘SERAFIND’ and
‘MED’.
4.3.
PRESCRIBING THE BOUNDARY CONDITIONS
The boundary conditions are handled through types of conditions which are related to the
computational variables. The combination of these types (from a list of possible choices) describes
whether the boundary is liquid or solid and how it should be processed.
In TELEMAC-3D, the water depth H, the horizontal velocities U and V and the tracers are the only
variables which type of boundary conditions has to be defined. Those types of boundary conditions
applicable to the vertical velocity and the k and ε functions are managed by TELEMAC-3D by the
user directly in the FORTRAN source files of TELEMAC-3D and therefore do not have a type. If the
computation takes tracers into account, then a single type (common to all the tracers) should also
be defined for a given boundary.
Once all the types of the boundary are defined, the user should enter the related values for the
computational variables (at least H, U and V).
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 34
For example, the user may want to set the sea level and leave the velocity field free (e.g. the tide
case). The type of boundary will be: "prescribed depth and free velocity". The values required for
that type are only water depth at every instant at that boundary. The values of velocities (if they are
entered) are not taken into account for that boundary.
Thus, for each TELEMAC-3D boundary, the computational variables (at least H, U and V) are
necessarily associated with one type and each type may be associated with one value (either used
or not).
The maximum number of boundaries is set to 30 by default but it can be changed by the user with
the keyword MAXIMUM NUMBER OF BOUNDARIES. This avoids changing the previously
hardcoded values (until version 7.0), which required recompiling the whole package.
After such a description of what is a boundary in TELEMAC-3D, we will describe the types, then
the related values.
4.3.1.
THE BOUNDARIES IN TELEMAC-3D
Water depth is the only two-dimensional variable computed. Its processing at the boundaries is like
that being performed by TELEMAC-2D. The boundary points to be handled are those of the twodimensional mesh (refer to Figure 5).
For the other variables (velocities and tracers), the boundary conditions should be handled over all
the boundaries of the three-dimensional mesh which includes:
RELEASE 7.1
•
The lateral boundary points (vertical column points linked to the boundaries of the twodimensional mesh), whether it is a liquid or solid boundary,
•
The points belonging either to the free surface or the bottom.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 35
Figure 5: The various boundaries in TELEMAC-3D (bridge piers case study).
By default, TELEMAC-3D automatically handles all the surface and bottom points which do not
belong to the side walls. The user, however, can modify them. This can be done by modifying the
FORTRAN sources.
All the remaining points (on the lateral boundaries) are linked to the two-dimensional mesh
boundaries at each horizontal level. Thus, they will be processed in a similar fashion as in
TELEMAC-2D. The number of required data, however, increases so much that a full external
handling (either through the steering file or the boundary conditions file) would become excessively
complex. That is why the range of options offered to the user for dealing with these boundary
conditions is narrower than in TELEMAC-2D and definitely implies programming the sources of the
user-available software. The next following paragraphs describe the way the boundary nodes are
handled.
4.3.2.
THE BOUNDARY-RELATED TYPES
The boundary condition type for H, U, V and T of the edge points is read in the boundary conditions
file. It can be either modified or directly defined by the user in the LIMI3D subroutine.
The various types of boundary conditions can be combined in order to prescribe the conditions of
different physical kinds (liquid inflow or outflow in supercritical conditions, open sea, wall, etc.).
Some combinations, however, are not physical (refer to paragraph 4.3.3. hereinafter).
Some boundary conditions are applicable to such facts as friction at the walls or wall
impermeability. However, the wall definition is ambiguous if one only retains a definition of point
wise boundary conditions. The following convention is then observed in order to determine the
nature of a segment lying between two different kinds of points: a liquid segment is a segment
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 36
linking two liquid-types points. Thus, under that convention, the connecting point between the shore
and the marine boundary (or between the river and the bank) is preferably of the liquid type.
Therefore, liquid + solid = solid
Any sequential arrangement of the boundary types may exist along an outline (for instance, one
may have a liquid boundary with a prescribed depth followed by a liquid boundary with a prescribed
velocity). The only condition to be met is that a boundary should consist of at least two points (it is
a computational requirement, a number of at least four points being highly advisable from a
physical point of view).
4.3.3.
DESCRIPTION OF THE VARIOUS TYPES
The type of boundary condition at a given point is provided, in the boundary conditions file, in the
form of four integers which are referred to as LIHBOR, LIUBOR, LIVBOR and LITBOR, with
values which can range from 0 to 6.
The available options are as follows:
•
Depth condition:
• Prescribed depth liquid boundary: LIHBOR=5,
• Free depth liquid boundary: LIHBOR=4,
• Solid boundary (wall): LIHBOR=2.
It is noteworthy that a depth/rate law is considered as a prescribed depth condition. The flow rate
value should then explicitly be computed, according to the water depth, by programming the Q3
subroutine.
4.3.4.
•
Rate or velocity condition:
• Prescribed flow rate liquid boundary: LIUBOR/LIVBOR=5,
• Prescribed velocity liquid boundary: LIUBOR/LIVBOR=6,
• Free velocity liquid boundary: LIUBOR/LIVBOR=4,
• Solid boundary with sliding or friction: LIUBOR/LIVBOR=2,
• Solid boundary with one or two zero velocity components: LIUBOR and/or
LIVBOR=0.
•
Tracer condition:
• Prescribed tracer liquid boundary: LITBOR=5,
• Free tracer liquid boundary: LITBOR=4,
• Solid boundary (wall): LITBOR=2.
THE BOUNDARY CONDITIONS FILE
That file is provided as standard by MATISSE, JANET, BLUE KENUE or STBTEL, but it can be
created or amended by means of FUDAA-PREPRO or a text editor. Each line of that file is
dedicated to one point at the two-dimensional mesh boundary. The boundary point numbering is
the same as that of the file lines, it first describes the domain outline in the counter clockwise
direction, then the islands in the opposite direction.
The convention being observed in TELEMAC implies that the first liquid boundary is that which is
defined, within the boundary conditions file, by the first two liquid-typed consecutive numbers. In
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 37
the example below (channel case study), the first liquid boundary is defined by the nodes 42-47
(edge numbering) and corresponds to a prescribed depth (codes 5 4 4 at the beginning of lines).
The second boundary begins at number 76 and ends at number 1 and corresponds to a prescribed
rate (codes 4 5 5).
4 5 5
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
1
1
2 2 2
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
5
2
2 2 2
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
6
3
2 2 2
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
44
41
5 4 4
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
2
42
5 4 4
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
45
43
5 4 4
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
86
44
5 4 4
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
3
47
2 2 2
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
87
48
2 2 2
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
73
74
2 2 2
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
74
75
4 5 5
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
4
76
4 5 5
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
88
77
4 5 5
0.000
0.000
0.000
0.0
2
0.000
0.000
0.000
89
85
…
…
…
…
For each point, and each line in the boundary conditions file, the following values are entered:
LIHBOR, LIUBOR, LIVBOR, HBOR, UBOR, VBOR, AUBOR, LITBOR, TBOR,
ATBOR, BTBOR, N, K
•
LIHBOR, LIUBOR, LIVBOR and LITBOR are boundary-typed integers for each of the
variables.
• HBOR (real) denotes the prescribed depth value when the LIHBOR value is set to 5,
• UBOR (real) denotes the prescribed U velocity value when the LIUBOR value is set to 6,
• VBOR (real) denotes the prescribed V velocity value when the LIVBOR value is set to 6,
• AUBOR denotes the value of the boundary friction coefficient when the LIUBOR or LIVBOR
value is set to 2. The friction law is then written as:
υT
dU
= AUBOR × U
dn
and/or υ T
dV
= AUBOR × V
dn
The AUBOR coefficient is applicable to the segment included between the edge point being
considered and the next point (in the counter clockwise direction for the outside outline and in the
clockwise direction for the islands). By default, the AUBOR value is 0. A friction corresponds to a
negative value. With the k-ε model, the value of AUBOR is automatically computed by TELEMAC3D, the indications in the boundary conditions file will then be ignored.
• TBOR (real) denotes the prescribed tracer value when the LITBOR value is set to 5,
• ATBOR et BTBOR denote the coefficient values of the flux law which is written as:
υT
RELEASE 7.1
dT
= ATBOR × T + BTBOR
dn
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 38
The ATBOR and BTBOR coefficients are applicable to the segment included between the edge point
considered and the next point (in the counter clockwise direction for the outside outline and in the
clockwise direction for the islands).
• N denotes the edge point global number,
• K denotes the point number in the edge point numbering. This number also represents a
node colour (as an integer). This number named BOUNDARY_COLOR, can be used in
parallel simulations to simplify the implementation of specific cases. Without particular
modification, this value is the rank of the border point in the global numbering. For
example, a test like IF
(I.EQ.144)
THEN can be replaced by IF
(BOUNDARY_COLOUR%I(I).EQ.144) THEN which is compatible with the parallel mode.
However, this only concerns the 2D mesh (Table BOUNDARY_COLOUR is only given for
level 1). Be careful not to modify the last column of the boundary conditions file that
contains this BOUNDARY_COLOUR table, when using tidal harmonic constants databases
(cf. [6]).
As regards the horizontal velocities, all the points in one water column will have the same type of
boundary condition defined by LIUBOR or LIVBOR. That principle is intrinsic to the TELEMAC-3D
formulation. Prescribing a different type of boundary condition in the vertical direction (e.g. for a
subterranean stream) may, indeed, induce severe inconsistencies with the hydrostaticity
hypothesis and generate, for instance, unrealistic vertical velocities. It is then advisable that the
user will follow that principle. Nonetheless, the boundary condition type in the vertical direction can
be altered through direct programming in the LIMI3D subroutine.
The so-called LIHBOR, LIUBOR and LIVBOR integers (which define the boundary type) can
assume a value ranging from 0 to 6. The available options are as follows:
•
Depth-related condition:
• Prescribed depth liquid boundary: LIHBOR=5,
• Free depth liquid boundary: LIHBOR=4,
• Solid boundary (wall): LIHBOR=2,
•
Velocity-related condition:
• Prescribed velocity liquid boundary: LIUBOR/LIVBOR=6,
• Prescribed rate liquid boundary: LIUBOR/LIVBOR=5,
• Free velocity liquid boundary: LIUBOR/LIVBOR=4,
• Solid boundary with sliding or friction: LIUBOR/LIVBOR=2,
• Solid boundary with one or two zero velocity components: LIUBOR
LIVBOR=0.
and/or
The boundary conditions of physical nature are defined by the relationship among the types of
variables. In most cases, the boundary type can be set by a mesh generator (e.g. MATISSE or
JANET) in the TELEMAC chain. The table below summarizes the physical relationship among the
boundary types.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
4.3.5.
Page 39
LIHBOR
LIUBOR
LIVBOR
LITBOR
2
2
2
2
2
0
2
2
2
2
0
2
2
0
0
2
4
4
4
4
5
4
4
4
5
4
0
4
Prescribed H, free U, zero V, free T.
5
0
4
4
Prescribed H, zero U, free V, free T.
1
1
1
4
Incident wave, free tracer.
4
5
5
5
Free H, prescribed Q, prescribed T.
4
5
0
5
Free H, prescribed Q with zero V, prescribed
T.
4
0
5
5
Free H, prescribed Q with zero U, prescribed
T.
4
6
6
5
5
5
5
5
5
6
6
5
Solid wall.
Solid wall with zero U.
Solid wall with zero V.
Solid wall with zero U and V.
Free H, free velocities, free T.
Prescribed H, free velocities, free T.
Free H, prescribed velocities, prescribed T.
Prescribed H and Q, prescribed T.
Prescribed H and velocities, prescribed T.
PROGRAMMING THE BOUNDARY CONDITIONS TYPE
The LIMI3D subroutine can be programmed to handle specific boundary conditions, for the edge
points as well as the surface and bottom points.
That subroutine is called upon each time step. Therefore, it can be used to change the boundary
condition type in time, if required.
4.3.6.
PRESCRIBING VALUES THROUGH KEYWORDS
In most simple cases, the boundary conditions will be prescribed by means of keywords. However,
if the values to be prescribed are variable in time, it is necessary to program the adequate functions
or to use the liquid boundaries file.
The appropriate keywords to prescribe the boundary values are as follows:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
Page 40
PRESCRIBED ELEVATIONS: provided to set the elevation value of a prescribed height
liquid boundary. It is an array which can contain up to 100 reals, and therefore up to 100
boundaries of that kind can be handled. The values defined by that keyword overwrite the
depth values read from the boundary conditions file,
Warning: the free surface level is set with this keyword, whereas the water depth is
set in the boundary conditions file.
PRESCRIBED FLOWRATES: provided to set the flow rate value of a prescribed flow at a
liquid boundary. It is an array which can contain up to 100 reals, and therefore up 100
boundaries of that kind can be handled. A positive value corresponds to a domain inflow
rate. The values defined by that keyword overwrite the velocity values read from the
boundary conditions file,
PRESCRIBED VELOCITIES: provided to set the velocity value of a prescribed velocity
liquid boundary. The scalar value is the wall normal velocity. A positive value corresponds
to a domain inflow. It is an array which can contain up to 100 reals, and therefore up 100
boundaries of that kind can be handled. The values defined by that keyword overwrite the
values read from the boundary conditions file.
In addition, several simple rules should be observed:
The boundary type as specified in the boundary conditions files should obviously be in accordance
with the keywords in the steering file (do not insert the keyword PRESCRIBED FLOWRATES if
there are no boundary points the LIUBOR and LIVBOR of which are set to 5). The keyword,
however, is ignored if no type matches it.
For each keyword, the number of specified values should be equal to the whole number of liquid
boundaries, whatever their types may be. When a boundary is inconsistent with the keyword, then
the specified value is ignored (a 0.0 value or, on the contrary, a very high value such as 999.0 may
be systematically inserted).
For example, in the channel test case, the first boundary (downstream boundary) is of prescribed
level type whereas the second one (upstream boundary) is of prescribed flow rate type. The
steering file contains a sequence of the following type:
PRESCRIBED ELEVATIONS = 0.5, 0.0
PRESCRIBED FLOWRATES = 0.0, 50.0
4.3.7.
BOUNDARY CONDITION ON THE BOTTOM
By default, the boundary condition on the bottom is an impermeable slip boundary (Neumann
condition of the same type as vertical conditions).
However, bottom velocities can be set to zero by using value 2 of the keyword BOUNDARY
CONDITION ON THE BOTTOM (Default value of 1 correspond to a slip condition). This option is
valid only if the vertical mesh is refined at bottom level.
Since version 7.1, it is possible to prescribe a flux on the bed in TELEMAC-3D (e.g.: a flow rate on
several liquid boundaries placed on the bed). To do so, it is necessary to define the imposed flow
rates using the keywords OPEN BOUNDARY CONDITIONS ON THE BED set to YES (default
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 41
value = NO) and PRESCRIBED FLOWRATES ON THE BED with values following the same
structure as for other prescribed flow rates in the TELEMAC-MASCARET system. It should be a list
of numbers separated by a semi-colon, one number per liquid boundary on the bed must be given.
The maximum number of boundaries on the bed is set to 30 by default but it can be changed by the
user with the keyword MAXIMUM NUMBER OF BOUNDARIES ON THE BED.
At the moment, the BOUNDARY CONDITIONS FILE only deals with horizontal boundaries,
therefore the user has to define the liquid boundary on the bed by hand. This can be done by
modifying the subroutine LIMI3D in the FORTRAN FILE. For example to add a circular boundary
of radius 5 m centred around coordinate (2 000, 2 000) m, the following modifications can be done:
...
!
BOUNDARY CONDITIONS ON VELOCITIES
!
*********************************
!
!
BOTTOM
!
======
!
!
DEFAULT: IMPERMEABILITY AND LOG LAW (SEE ALSO BORD3D)
!
IF(BC_BOTTOM.EQ.1) THEN
!
DO IPOIN2 = 1,NPOIN2
LIUBOF%I(IPOIN2) = KLOG
LIVBOF%I(IPOIN2) = KLOG
LIWBOF%I(IPOIN2) = KLOG
!
USEFUL ? SHOULD NOT BE USED ANYWAY
UBORF%R(IPOIN2) = 0.D0
VBORF%R(IPOIN2) = 0.D0
WBORF%R(IPOIN2) = 0.D0
IF(SQRT((X(IPOIN2)-2000.D0)**2
&
+(Y(IPOIN2)-2000.D0)**2)
&
.LE.50.D0) THEN
!5: IMPOSED FLOW RATE
LIUBOF%I(IPOIN2) = 5
LIVBOF%I(IPOIN2) = 5
LIWBOF%I(IPOIN2) = 5
NLIQBED%I(IPOIN2) = 1
WRITE(LU,*) '========================'
WRITE(LU,*) 'FOR POINT ',IPOIN2
WRITE(LU,*) 'BEDFLO',BEDFLO(1)
ENDIF
ENDDO
!
...
In this example, it should be noted that NLIQBED%I(IPOIN2) = 1 defines the position for the first
liquid boundary defined in the steering file.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 42
This is all that needs to be defined by the user to deal with fluxes on the bed. However, in version
7.1, only constant velocity profile is available. It should also be noted that it has not been possible
to prescribe a tracer or turbulence yet.
4.3.8.
USING THE LIQUID BOUNDARIES FILE
In case of time variable values, which are nonetheless constant in space along the relevant liquid
boundary, the prescription can be made using the liquid boundaries file (as an alternative to
programming).
It is a user-edited text file the name of which should be given by the keyword LIQUID
BOUNDARIES FILE. That file has the following format:
•
The optional line(s) begin(s) with the sign # (1st character on the line) will be treated as
comments,
•
It should contain a header line beginning with T for identifying the supplied time dependent
value(s) within that file. The identification is performed through mnemonic means which are
identical to the variable names: Q for the flow rate, SL for the level, U and V for the
velocities and TR for the tracer. These characters are directly followed by an integer in
between brackets which is used to specify the current boundary. That line is necessarily
followed by another line indicating the unit of the variables (lines of comments can be
inserted, but the line of units should be present). The units are given for information only
and TELEMAC-3D does not handle the conversion of units (thus, the user has to enter the
values using the standard unit),
•
The values to be prescribed are provided through a sequence of lines the format of which
should be consistent with the identification line. The time value should be increasing and
the last time value supplied should be higher than or equal to the value of the last time
step, otherwise the computation is suddenly interrupted.
Upon the retrieval of that file, TELEMAC-3D performs a linear interpolation in order to compute the
value prescribed at a particular time step. The value which is actually prescribed by the code is
printed on the check listing.
An example of a liquid boundaries file is given below.
# Example of liquid boundaries file
# 2 boundaries are managed
#
T
Q(1) SL(2)
s
m3/s m
0.
0.
135.0
25.
15.
135.2
100. 20.
136.
500. 20.
136.
In that example, the flow rate is prescribed at the first boundary and the free surface is prescribed
at the second boundary.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
4.3.9.
Page 43
PRESCRIBING VALUES THROUGH PROGRAMMING
Still in the case of time variable values that are constant in time along the liquid boundary
processed, the prescription can be made simply by programming particular functions:
• VIT3 function for prescribing a velocity,
• Q3 function for prescribing a flow rate,
• SL3 function for prescribing an elevation.
Functions Q3, VIT3 and SL3 are similarly programmed. In each case, the user knows the time,
the boundary rank (e.g. to determine whether the first or the second prescribed flow rate boundary
is processed). By default, the functions prescribe values that are read from the boundary conditions
file or provided by the keywords.
For instance, the body of function Q3 used to prescribe a flow rate ramp for the first 1,000 seconds
from 0 to 400 m3/s can take such a form as:
IF (AT.LT.1000.D0) THEN
Q3 = 400.D0 * AT/1000.D0
ELSE
Q3 = 400.D0
ENDIF
Or
Q3 = 400.D0 * MIN(1.D0,AT/1000.D0)
4.3.10. STAGE-DISCHARGE CURVES
TELEMAC-3D allows managing liquid boundaries for which the prescribed water elevation value is
a function of local flow rate. This situation is encountered particularly in river hydraulics.
First, it is necessary to specify which boundaries are concerned with the keyword STAGEDISCHARGE CURVES. This keyword provides an integer value for each boundary. This value can
be:
• 0: no stage-discharge curves (default),
• 1: elevation as a function of local flow rate.
The keyword STAGE-DISCHARGE CURVES FILE provides the name of the text file containing
information about the curves. An example is shown below:
#
# STAGE-DISCHARGE CURVE BOUNDARY 1
#
Q(1)
Z(1)
m3/s
m
61.
0.
62.
0.1
63.
0.2
#
# STAGE-DISCHARGE CURVE BOUNDARY 2
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
#
Z(2)
m
10.
20.
30.
40.
50.
Page 44
Q(2)
m3/s
1.
2.
3.
4.
5.
The order of curves has no significance. The column order can be reversed, as is the case for the
second boundary in the example. Lines beginning with # are comments. The lines defining the
units are mandatory but the units are not checked. The number of points of each curve is
completely free and need not be the same for each curve.
Warning: at initial conditions, the flow at the exit can be null. The initial level must correspond to
that of the calibration curve otherwise a sudden change is imposed. To avoid extreme situations,
the curve should be limited to a certain level of flow rate. In the example of boundary 1 above, the
flow rates below 61 m3/s generate a water elevation of 0 m above the flow of 63 m3/s to produce an
elevation equal to 0.2 m.
4.3.11. PRESCRIBING COMPLEX VALUES
If the values to be prescribed vary in both space and time, then a programming in the BORD3D
subroutine becomes necessary, since that subroutine can be used to prescribe the values in a
node wise way.
That subroutine describes all the liquid boundaries (loop on NPTFR2). For each boundary point, it
determines the boundary type in order to prescribe the adequate value (velocity, elevation or flow
rate). BORD3D programming to prescribe a flow rate, however, hardly makes any sense, since the
flow rate value is generally known for the whole boundary rather than on each boundary segment.
If a prescribed flow rate inlet is surrounded by walls with an adherence, then the corner velocities
are cancelled.
Note that, the BORD3D subroutine also makes it possible to prescribe the complex boundary values
of the tracers.
4.3.12. PRESCRIBING A PROFILE
4.3.12.1. HORIZONTAL PROFILE
When processing a prescribed flow rate or prescribed velocity boundary, the user has the keyword
VELOCITY PROFILES to specify which "horizontal" velocity profile should be prescribed by
TELEMAC-3D. The following options are possible:
• 1: The profile is normal and homogeneous along the boundary (default option),
• 2: The values of U and V are read from the boundary conditions file (UBOR and VBOR
values). In case of a prescribed flow rate, these values are multiplied by a constant in order
to get the desirable flow rate,
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
Page 45
3: The velocity vector is normal to the boundary and its norm is read from the boundary
conditions file as the UBOR value. In case of a prescribed flow rate, that value is multiplied
by a constant in order to get the desirable flow rate,
4: The velocity vector is normal to the boundary and its norm is proportional to the square
root of the water depth,
5: The velocity vector is normal to the boundary and its norm is proportional to the square
root of a virtual water depth computed from lowest point of the free surface on the
boundary.
4.3.12.2. VERTICAL PROFILE
When processing a prescribed flow rate or prescribed velocity boundary, the user has the keyword
VELOCITY VERTICAL PROFILES to specify which "vertical" velocity profile should be prescribed
by TELEMAC-3D. The options for that keyword are:
• 0: programmed by the user,
• 1: constant (default value for all the liquid boundaries),
• 2: logarithmic.
The user programming is done within the VEL_PROF_Z subroutine.
Activating the keyword DYNAMIC BOUNDARY CONDITION (default value = FALSE) also enables
to prescribe a velocity at the free surface coherent with the dynamic boundary condition.
4.3.13. THOMPSON CONDITIONS
In some cases, not all the necessary information concerning the boundary conditions is available.
This is usual for coastal domains where only the values of the sea level on several points are
known. This kind of model is referred to as an “under-constrained” model.
To solve this problem, the Thompson method uses the characteristics method to calculate the
missing values. For example, TELEMAC-3D will compute the velocity at the boundary in the case
of a prescribed elevation.
This method can also be used for “over-constrained” models. In this case, the user specifies too
much information at the boundary. If the velocity information and the level information are not
consistent, too little or too much energy is going into the model. For this, the Thompson method
computes a new value for the velocity and performs small adjustments to cancel the
inconsistencies in the information.
For this, the user can use the keyword OPTION FOR LIQUID BOUNDARIES, which offers two
values (the user must specify 1 value for every open boundary):
• 1: strong setting (default value for all boundaries),
• 2: Thompson method.
Taking a simplified view, it may be said that, in the case of the first option, the values are
“imposed”, in the case of the second option, the values are “suggested”.
However it is important to note that, given the two-dimensional aspect, the Thompson method can
only be used in the case of a zero velocity gradient imposed on the vertical (uniform velocity along
the vertical).
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 46
4.3.14. TIDAL HARMONIC CONSTITUENTS DATABASES
4.3.14.1. GENERAL PARAMETERS
To prescribe the boundary conditions of a coastal boundary subject to tidal evolution, it is generally
necessary to have the information characterizing this phenomenon (harmonic constants). One of
the most common cases is to use the information provided by large scale models.
4 databases of harmonic constants are interfaced with TELEMAC-3D:
• The JMJ database resulting from the LNH Atlantic coast TELEMAC model by Jean-Marc
JANIN,
• The global TPXO database and its regional and local variants from the Oregon State
University (OSU),
• The regional North-East Atlantic atlas (NEA) and the global atlas FES (e.g. FES2004 or
FES2012...) coming from the works of Laboratoire d’Etudes en Géophysique et
Océanographie Spatiales (LEGOS),
• The PREVIMER atlases.
However it is important to note that, in the current version of the code, the latter 2 databases are
not completely interfaced with TELEMAC-3D and their use is recommended only for advanced
users.
The keyword OPTION FOR TIDAL BOUNDARY CONDITIONS activates the use of one of the
available database when set to a value different from 0 (the default value 0 means that this function
is not activated). Since version 7.1 this keyword is an array of integers separated by semicolons
(one per liquid boundary) so that the user can describe whether tidal boundary conditions should
be computed or not (e.g. a weir) on a liquid boundary. When this keyword is activated, every tidal
boundary is treated using the prescribed algorithms for the boundaries with prescribed water
depths or velocities, with the same option for tidal boundary conditions (the values not equal to 0
have to be the same). The databases provide only a single value of the depth-averaged velocity,
thus TELEMAC-3D prescribes the same value of the velocity at each point along the vertical (for
more information, see Méthodologie pour la simulation de la marée avec la version 6.2 de
TELEMAC-2D et TELEMAC-3D, C.-T. Pham et al., EDF report H-P74-2012-02534-FR).
The database used is specified using the keyword TIDAL DATA BASE which can take the values:
• 1: JMJ,
• 2: TPXO,
• 3: MISCELLANEOUS (LEGOS-NEA, FES20XX, PREVIMER...).
Depending on the database used, some keywords have to be specified.
• If using the JMJ database, the name of the database (typically bdd_jmj) is given by the
keyword ASCII DATABASE FOR TIDE and the corresponding mesh file is specified using
the keyword TIDAL MODEL FILE,
• If using the TPXO database, the name of the water level database is given by the keyword
BINARY DATABASE 1 FOR TIDE (for example h_tpxo7.2) and the name of the velocity
database is given by the keyword BINARY DATABASE 2 FOR TIDE (for example
u_tpxo7.2). Moreover, it is possible to activate an interpolation algorithm of minor
constituents from data read in the database using the logical keyword MINOR
CONSTITUENTS INFERENCE, activation not done by default.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 47
The keyword OPTION FOR TIDAL BOUNDARY CONDITIONS specifies the type of tide to
prescribe. The default value 0 means no prescribed tide or that the tide is not treated by standard
algorithms. Value 1 corresponds to prescribing a real tide considering the time calibration given by
the keywords ORIGINAL DATE OF TIME (YYYY ; MM ; DD format) and ORIGINAL HOUR OF
TIME (HH ; MM ; SS format). Other options are the following, available for every tidal database
(JMJ, TPXO-type from OSU, LEGOS-NEA, FES, PREVIMER…). We call them “schematic tide” for
values from 2 to 6:
• 2: exceptional spring tide (French tidal coefficient approximately equal 110),
• 3: mean spring tide (French tidal coefficient approximately equal 95),
• 4: mean tide (French tidal coefficient approximately equal 70),
• 5: mean neap tide (French tidal coefficient approximately equal to 45),
• 6: exceptional neap tide (French tidal coefficient approximately equal to 30),
• 7: real tide (before 2010 methodology).
In the case of options 2 to 6 (schematic tides), the boundary conditions are imposed so that the
reference tide is approximately respected. In order to shift the phases of the waves of the tidal
constituents so that the computation starts close to a High Water, two keywords are available. If
using a TPXO-type tidal database from Oregon State University, the keyword GLOBAL NUMBER
OF THE POINT TO CALIBRATE HIGH WATER has to be filled with the global number of the point
with respect to which the phases are shifted (mandatory, otherwise the computation stops). If using
one of the other tidal databases (JMJ, NEA/FES, PREVIMER) the keyword LOCAL NUMBER OF
THE POINT TO CALIBRATE HIGH WATER should be filled in with the local number between 1
and the number of tidal boundary points of the HARMONIC CONSTANTS FILE; If not filled in
(default value = 0), a value is then automatically calculated. However, it is usually necessary to wait
for the second or third modelled tide in order to overcome the transitional phase of start-up of the
model. It is also necessary to warn the user that the French tidal coefficients shown are
approximate.
During a simulation, data contained in the tidal database are interpolated on boundary points.
When using of the JMJ database, this spatial interpolation can be time consuming if the number of
boundary points is important, and is not yet available in case of parallel computing. It is therefore
possible to generate a file containing harmonic constituents specific to the model treated. The
principle is at a first step, to perform a calculation on a single time step whose only goal is to
extract the necessary information and to generate a file containing for each boundary point of the
model, the harmonic decomposition of the tidal signal. Subsequent calculations directly use that
specific file rather than directly addressing to the global database. The harmonic constants specific
file is specified using the keyword HARMONIC CONSTANTS FILE, this file is an output file in the
first calculation, and an input file in subsequent calculations.
4.3.14.2. HORIZONTAL SPATIAL CALIBRATION
In order to perform the spatial interpolation of the tidal data, it is imperative to provide to
TELEMAC-3D information on the spatial positioning of the mesh model relative to the grid of the
tidal database. To do this, the user has two keywords:
The first keyword specifies the geographic system used to establish the coordinates of the 2D
mesh of TELEMAC-3D. This keyword GEOGRAPHIC SYSTEM, which has no default value, may
take the following values:
• 0: User Defined,
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
•
•
Page 48
1: WGS84 longitude/latitude in real degrees,
2: WGS84 UTM north,
3: WGS84 UTM south,
4: Lambert,
5: Mercator projection.
The second keyword is used to specify the area of the geographic system used to establish the
coordinates of the 2D mesh of TELEMAC-3D. This keyword ZONE NUMBER IN GEOGRAPHIC
SYSTEM which has no default value, may take the following values:
• 1: Lambert 1 north,
• 2: Lambert 2 center,
• 3: Lambert 3 south,
• 4: Lambert 4 Corsica,
• 22: Lambert 2 extended,
• X: UTM zone value of the WGS84 (X is the number of the zone).
4.3.14.3. CALIBRATION OF THE INFORMATION
The transfer of information between a large scale model and the boundaries of a more local model
generally requires calibration.
To do this, the user has three keywords:
• The keyword COEFFICIENT TO CALIBRATE SEA LEVEL (default real value 0.0) is used
to calibrate the mean tide level (the harmonic decomposition of information provided by the
various databases are used to generate the tidal signal oscillating around mean tide level).
The calibration of the mean tide level must obviously be made depending on the altimetric
reference used in the model,
• The keyword COEFFICIENT TO CALIBRATE TIDAL RANGE (default real value 1.0) is
used to specify a calibration coefficient applied on the amplitude of the tidal wave. This
coefficient is applied to the amplitude of the overall signal, and not on the amplitude of
each of the elementary waves,
• The keyword COEFFICIENT TO CALIBRATE TIDAL VELOCITIES (default real value
999,999.0) is used to specify the coefficient applied on velocities. The default value (999,
999.0) means that the square root of the value specified by the keyword COEFFICIENT
TO CALIBRATE TIDAL RANGE tidal is used.
For more information, the reader may refer to methodological guide for tide simulation with version
6.2 [6] available in the “Training material and tutorials” section of the opentelemac website.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
5.
Page 49
PHYSICAL SETUP OF THE HYDRODYNAMIC
COMPUTATION
In addition to the general setup given in the steering file, a number of physical parameters may or
should be specified upon a simulation.
5.1.
HYDROSTATIC PRESSURE HYPOTHESIS
Firstly, it should be specified whether one wants to use the hydrostatic pressure hypothesis or not.
That choice is made using the keyword NON-HYDROSTATIC VERSION which, by default, is set to
NO. As a reminder, the hydrostatic pressure hypothesis consists in simplifying the W vertical
velocity, ignoring the diffusion, advection and other terms. Therefore, the pressure at a point is only
related to the weight of the overlying water column and to the atmospheric pressure at the surface.
Without the hydrostatic pressure hypothesis (NON-HYDROSTATIC VERSION = YES), TELEMAC3D solves a W vertical velocity equation which is similar to the U and V equations, with the
additional gravity term.
5.2.
MODELLING TURBULENCE
The Reynolds numbers (Re = U.L/υ) reached for tidal flows or in an estuary are excessively high
and illustrate basically turbulent flows (L, the scale of eddies, for example, assumes the value of
water depth h for a vertically homogeneous flow). For such a kind of flow, the turbulence-induced
momentum prevails (in relation to the molecular diffusion). That diffusion is strictly defined by a
tensor which could be anisotropic.
The concept of eddy scale, however, is spatially constrained by the horizontal and vertical scales of
the modelled domain. At sea, for example, a one kilometre long cape can generates eddies the
dimensions of which are horizontally related to that scale. Vertically, however, the eddy size is
constrained by the water depth or even more by possible effects of stratifications. Synthetically, a
common practice consists in separating the vertical and horizontal turbulence scales which are not
relevant to the same dynamics for the standard applications of TELEMAC-3D. That involves
defining horizontal as well as vertical viscosities rather than a single viscosity. On the open sea, for
instance, the horizontal and vertical viscosities differ by several orders of magnitude.
Thus, the implementation of TELEMAC-3D requires defining two models of horizontal and vertical
turbulence (HORIZONTAL TURBULENCE MODEL, VERTICAL TURBULENCE MODEL).
Turbulence modelling is an awkward task and TELEMAC-3D offers the user several approach
options which are different, but also increasingly complex, and are applicable to velocities as well
as to active and passive tracers.
The Von Karman constant and Prandtl number (ratio between eddy viscosity and eddy diffusivity)
can be changed with the keywords KARMAN CONSTANT (default value = 0.4) and PRANDTL
NUMBER (default value = 1.).
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
5.2.1.
Page 50
CONSTANT VISCOSITY
The simplest turbulence model consists in using a constant viscosity coefficient (option for
parameters: 1="CONSTANT VISCOSITY", default value). In that case, the latter includes the
effects of molecular viscosity and dispersion (refer to Theoretical note [1]). The horizontal and
vertical turbulent viscosities are then constant throughout the domain. The global (molecular +
turbulent) viscosity coefficients are provided by the user by means of the keywords COEFFICIENT
FOR HORIZONTAL DIFFUSION OF VELOCITIES and COEFFICIENT FOR VERTICAL
DIFFUSION OF VELOCITIES, set by default to 10-6.
The value of that coefficient has a definite effect on both size and shape of the recirculations and
eddies. A low value will tend to only dissipate the small-sized eddies, a high value will tend to
dissipates large-sized recirculations. The user shall then carefully select that value according to the
case studied. Usually, that value becomes a model calibration data by comparison with
measurements. Besides, it is worth mentioning that a value bringing about the dissipation of
recirculations of a smaller than two mesh cell extent has nearly no influence on the computation
(i.e. there is a threshold beyond which the viscosity or turbulence value has substantially no effect).
TELEMAC-3D makes it possible to get a space- and time-variable coefficient. The VISCOS
subroutine will necessarily be programmed. Within that subroutine, geometrical information, basic
hydrodynamic information (water depth, velocity components) and time are made available to the
user.
That option theoretically aims at enabling the user to define the turbulent viscosity by programming
the VISCOS subroutine.
5.2.2.
MIXING LENGTH (VERTICAL MODEL)
The user also has the opportunity to use a vertical mixing length model (VERTICAL TURBULENCE
MODEL: 2="MIXING LENGTH"). The vertical diffusivity of velocities is then automatically computed
by TELEMAC-3D by means of the selected mixing length model taking or not taking the effects of
density into account. The mixing length model expresses the turbulent viscosity (or diffusion
coefficient) as a function of the mean velocity gradient and the mixing length (Prandtl’s theory):
ν = L2m 2 Dij Dij
, where Dij =
1 ∂U i ∂U j
+
2 ∂x j
∂xi
Due to that choice, the user should enter the following options for the mixing length model (MIXING
LENGTH MODEL):
• 1: Prandtl (default). Standard Prandtl’s model. That formulation suits such flows with a
strong barotropic component as the tidal flows,
• 3: Nezu and Nakagawa. Nezu and Nakagawa model,
• 5: Quetin. Better representation of wind drift. In windy weather, a surface boundary layer is
formed and viscosity decreases,
• 6: Tsanis. Better representation of wind drift.
The graph below shows the variations of the mixing length for the various models.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 51
Figure 6: Mixing lengths versus depth.
In the presence of a vertical density gradient, the environment stability (respectively the instability)
hinders (enhances) the vertical exchanges of mass and momentum.
In order to quantify the effects of the gravity terms in the turbulent power balance, the
dimensionless Richardson number is commonly used. It is a local number the value of which can
obviously be different at each flow point.
In order to take the mixture reduction into a stable stratified flow into account, a damping law is
introduced into the turbulence model according to the Richardson number. The user can set the
damping function through the keyword DAMPING FUNCTION. The available options are:
•
0: nothing (default value),
•
1: user-performed; Law programming in the DRIUTI subroutine,
•
2: Viollet,
•
3: Munk and Anderson.
The graph below illustrates the variation of the Munk and Anderson damping function according to
the Richardson number for velocity and salinity. In the case of a stable stratification, the pressure
fluctuations more readily transmit a momentum flux than a mass flux and the diffusion coefficient
becomes higher for the velocities than for the mass.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 52
Figure 7: Munk and Anderson damping function.
When using the mixing length model, the user can also configure the calculation of the vertical
derivative of the velocities with the keyword VERTICAL VELOCITY DERIVATIVES. Default value 1
corresponds to a vertical derivative that is linear. Value 2 activates a logarithmic computation
between the bottom and 0.2 times the water depth. This allows getting better results when
modelling the velocity profile near the bottom. This algorithm is implemented in the subroutine
VISCLM.
5.2.3.
SMAGORINSKY
That option is activated by setting the horizontal or vertical turbulence model to 4 (Smagorinsky).
The Smagorinsky scheme is recommended, in particular, in the presence of a highly non-linear
flow.
5.2.4.
K-
ε
TELEMAC-3D gives an opportunity to use the so-called k-ε model as proposed by Rodi and
Launder for solving the turbulence equations. That model is activated by setting the keywords of
turbulence models (HORIZONTAL TURBULENCE MODEL and VERTICAL TURBULENCE
MODEL) to the value 3.
The k-ε model is defined through a couple of equations solving the balance equations for k
(turbulent energy) and ε (turbulent dissipation). Applying the k-ε model often requires using a finer
two-dimensional mesh than the constant viscosity model and then increases the computation
times.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 53
For detailed information about the formulation of the mixing length and k-ε models, the user can
refer to the TELEMAC-3D Theoretical Note.
Strictly speaking, except for the constant viscosity model, the diffusion coefficient should be equal
to the molecular diffusion of water:
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES = 1.D-6
COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES = 1.D-6
The default value of that viscosity coefficient may have to be increased to ensure a minimum
diffusion, especially during the first time steps of the computation.
There are two options to compute the lateral boundary conditions of k and ε (in subroutine KEPCL3)
with the keyword OPTION FOR THE BOUNDARY CONDITIONS OF K-EPSILON:
5.3.
•
1 = No turbulence (k and ε takes the minimum values KMIN and EMIN defined in
subroutine CSTKEP), which is the default value,
•
2 = Hans and Burchard’s formula (introduced in version 7.0).
SETTING UP THE FRICTION
The bottom or sidewall friction reflects the continuity of the constraint at the fluid-solid interface.
Knowing the constraint involves knowing the flow in the vicinity of the bottom. The turbulence
models provide a modelling for that flow.
The constraint can be written in several forms:
r
r
∂U
1
= − ρC f U 2 + V 2 U
τ =µ
∂n
2
= − ρU *2
r
r
where U * denotes the friction, C f - a dimensionless friction, U - the velocity of current recorded
far enough from the wall.
The friction condition is then provided:
•
Either by a turbulence model which indicates the constraint through a friction velocity
formulation,
•
Or by the knowledge of the friction coefficient C f and the related velocity U (here, the
r
vertically-averaged velocity). That approach will then use the Chézy, Strickler, Manning…
laws.
The same approach is adopted for both sidewalls and bottom.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
5.3.1.
Page 54
BOTTOM FRICTION
The friction law used for the bottom friction modelling is set by the keyword LAW OF BOTTOM
FRICTION which can assume the following values:
• 0: No friction,
• 1: Haaland law,
• 2: Chézy law (default value),
• 3: Strickler law,
• 4: Manning law,
• 5: Nikuradse law.
As regards the 1-5 values, the value of the friction coefficient corresponds to the selected law, and
shall be given by means of the keyword FRICTION COEFFICIENT FOR THE BOTTOM. Obviously,
that only holds true if the friction is constant in both space and time. The default value for that
parameter is 60.
The computation of turbulent constraint at the bottom depends on the velocity profile above the
bottom (within the boundary layer). The profile depends on the ratio of the wall asperity size to the
viscous sub-layer thickness (for further details, refer to the Theoretical Manual). When the
asperities are larger than the viscous sub-layer thickness, the latter cannot be established and the
friction regime is rough. On the other hand, when there is a viscous sub-layer, the friction regime is
smooth.
The computation of the turbulent constraint depends on the keyword TURBULENCE REGIME FOR
THE BOTTOM. The available options are:
• 1: smooth regime,
• 2: rough (default value).
In smooth friction regime conditions, the friction law is not used and the constraint is computed
from Reichard law of velocity profile (a law giving the friction velocity valueU * ).
In rough friction regime conditions and for the bottom friction laws 0, 2, 3 and 4, the constraint is
computed from the friction velocity U * and its relation to the C f coefficient. For law 5, the friction
velocity is computed from the velocity profile within the logarithmic layer and from the asperity size
k S (FRICTION COEFFICIENT FOR THE BOTTOM).
5.3.2.
SIDEWALL FRICTION
The friction law used to model the sidewall friction is set by the keyword LAW OF FRICTION ON
LATERAL BOUNDARIES which can take the following values:
• 0: No friction (default value),
• 5: Nikuradse law.
The size of asperities (which is used in the Nikuradse law) is given by the FRICTION
COEFFICIENT FOR LATERAL SOLID BOUNDARIES (default value 60!).
The friction is activated using the keyword TURBULENCE REGIME FOR LATERAL SOLID
BOUNDARIES. The available options are:
• 1: smooth regime,
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
Page 55
2: rough (default value).
That option changes the formulation of the velocity profile and consequently, the friction velocity.
See previous section for more information.
5.4.
PUNCTUAL SOURCE TERMS
TELEMAC-3D offers an opportunity to place momentum sources or sinks in any point of the
domain.
The user places horizontally the various sources using the keywords ABSCISSAE OF SOURCES
and ORDINATES OF SOURCES. They are arrays of reals giving the co-ordinates of sources, in
meters. Actually, TELEMAC-3D will position a source at the closest mesh point to the point as
specified by these keywords. The software will determine itself the number of source according to
the number of values given to each keyword.
The vertical positioning of the sources is done using the keyword ELEVATIONS OF SOURCES.
TELEMAC-3D places the sources on the nearest mesh level. In this case, it is recommended to
use fixed levels at sources elevations in order to avoid unwanted vertical movement of the sources
during the simulation. Also note that the sources cannot be placed on the bottom level (level
number 1). This is to ensure an impermeable bottom boundary.
At each source, the user should specify the liquid flow rate. That liquid flow rate is given (in m3/s)
using the keyword WATER DISCHARGE OF SOURCES.
In case of sources with variable characteristics, the user can then either use specific programming
in the T3D_DEBSCE subroutine (and T3D_TRSCE in the presence of tracer), or use the source file
whose name is given by the keyword SOURCES FILE. This file has exactly the same structure as
the liquid boundary file. An example is shown below with two sources and two tracers. Between
two specified times, the information used by TELEMAC-3D at sources is obtained by linear
interpolation.
#
#
#
#
#
#
#
#
#
#
#
#
#
#
T
s
0.
RELEASE 7.1
FLOW RATES ANS TRACERS CONCENTRATIONS AT SOURCES 1 ET 2
T IS THE TIME
Q(1) IS FLOW RATE AT SOURCE 1
Q(2) IS FLOW RATE AT SOURCE 2
TR(1,1)
TR(1,2)
TR(2,1)
TR(2,2)
Q(1)
m3/s
0.
IS
IS
IS
IS
TRACER
TRACER
TRACER
TRACER
1
2
1
2
TR(1,1)
°C
99.
CONCENTRATIONS
CONCENTRATIONS
CONCENTRATIONS
CONCENTRATIONS
TR(1,2)
°C
20.
AT
AT
AT
AT
SOURCE
SOURCE
SOURCE
SOURCE
Q(2)
m3/s
0.
1
1
2
2
TR(2,1)
°C
30.
TR(2,2)
°C
40.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
2.
4.
1.
2.
Page 56
50.
25.
20.
80.
2.
4.
30.
30.
20.
20.
Besides, TELEMAC-3D is capable of taking into account an injection velocity (in m/s) at the
sources in the dynamic equations. By default, the injection takes places without any momentum
input. The user may prescribe a particular velocity. If the latter is constant throughout the
simulation, then its value can be given using the keywords VELOCITIES OF THE SOURCES
ALONG X, VELOCITIES OF THE SOURCES ALONG Y and VELOCITIES OF THE SOURCES
ALONG Z. Otherwise, the user should program the SOURCE subroutine in order to amend USCE (for
the velocity at sources along X) and VSCE (for the velocity at sources along Y). The user can use
the time and all the parameters of the sources within that subroutine.
From a theoretical point of view, complete mass conservation can only be ensured if the source is
treated as a Dirac function and not as a linear function. The type of treatment is indicated by the
user with the keyword TYPE OF SOURCES, which can have a value of 1 (linear function, default
value) or 2 (Dirac function). It should be noted that in the second case, the solutions are of course
less smoothed. It is the same implementation as in TELEMAC-2D and the Dirac option is
recommended with a big number of sources.
The maximum number of sources is set to 20 by default but it can be changed by the user with the
keyword MAXIMUM NUMBER OF SOURCES. This avoids changing the previously hardcoded
values (until version 7.0), which required recompiling the whole package.
5.5.
SETTING UP THE WATER-ATMOSPHERE EXCHANGES
5.5.1.
THE WIND
TELEMAC-3D can be used to simulate flow while taking into account the influence of a wind
blowing on the water surface. The logical keyword WIND is used first of all for determining whether
this influence is taken into account and if so, the coefficient is then provided with the keyword
COEFFICIENT OF WIND INFLUENCE (see next paragraph). Lastly, if the wind is constant in time
and space, wind speed in directions X and Y is supplied with the keywords WIND VELOCITY
ALONG X and WIND VELOCITY ALONG Y. Default values for these three coefficients are 0.
The coefficient of wind influence hides complex phenomena. In fact, the influence of the wind
depends on the smoothness (or, lack of it) of the free surface and the distance over which it acts
(called the “fetch”). The coefficient value can be obtained from many different formulas.
This is the formula used by the Institute of Oceanographic Sciences (United Kingdom):
If U vent < 5 m/s
If 5 < U vent < 19.22 m/s
If
U vent
> 19.22 m/s
avent
= 0.565 10-3
avent
= (- 0.12 + 0.137 U vent ) 10-3
avent
= 2.513 10-3
The parameter COEFFICIENT OF WIND INFLUENCE asked for by TELEMAC-3D is: avent (ρair /
ρ) and not only avent. ρair is approximately 1.2 kg/m3 and ρ is approximately 1,000 kg/m3. Thus it
is necessary to divide the value of avent by 1,000 to obtain the value of the TELEMAC-3D keyword.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 57
The whole formulation used to consider the wind effects, through the keyword COEFFICIENT OF
WIND INFLUENCE (refer to the Theoretical Note for the definition of that coefficient), on the
surface flows is fully stated in the BORD3D subroutine.
If the wind velocity is space- or time-variable, the user should modify the METEO subroutine.
WARNING: the METEO subroutine is provided to define the wind velocity and direction even though
they are space- or time-variable. The BORD3D subroutine describes the law of wind-induced drift of
water bodies.
If there are tidal flats or dry zones in the domain, the wind may trigger unphysical velocities as it
becomes the only driving term in the equations. To avoid this, the influence of the wind is cancelled
below a threshold value of depth, with the key-word THRESHOLD DEPTH FOR WIND (default
value at 1 m).
5.5.2.
THE TEMPERATURE
Air temperature may be specified using the keyword AIR TEMPERATURE if it is constant in time
and space.
5.5.3.
THE PRESSURE
The influence of air pressure is taken into account from the moment when the keyword AIR
PRESSURE is set to YES (the default value is NO). The value of that pressure is directly set in the
METEO subroutine. By default, the latter initializes a pressure of 105 Pa (≈ 1 atm) over the whole
domain.
5.5.4.
RAIN OR EVAPORATION
The modelling of the influence of precipitation or evaporation is activated with the logical keyword
RAIN OR EVAPORATION. The value of the contribution or the loss of water at the surface is
specified using the keyword RAIN OR EVAPORATION IN MM PER DAY which default value is 0
(a negative value reflects an evaporation).
In case of calculation with consideration of tracers, it is possible to specify the contribution related
to the rain with the keyword VALUES OF TRACERS IN THE RAIN (default value is 0.). It is
important to note that, in the case of evaporation, no tracer is taken into account in the water loss,
which is incorrect if the tracer is the temperature.
5.5.5.
ATMOSPHERE-WATER EXCHANGE MODELS
Before version 7.0, heat exchange between water and atmosphere could have been done with a
linearised formula of the balance of heat exchange fluxes at the free surface. An example of an
exchange with a constant atmosphere temperature and a constant sea salinity was given as
standard (as comments) through a direct programming in the BORD3D subroutine.
A much more elaborated model has been introduced in TELEMAC-3D. This module calculates the
complete balance of exchanged fluxes involved:
• The solar radiation,
• The atmospheric radiation,
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
•
•
Page 58
The water radiation,
The latent heat due to evaporation,
The sensitive heat of conductive origin.
It takes into account the solar radiation penetration in the water column.
The choice of the heat exchange model can be done with the keyword ATMOSPHERE-WATER
EXCHANGE MODEL (default value = 0: no exchange model). Value 1 will use with the linearised
formula at the free surface, whereas value 2 will use with the model with complete balance.
These calculations require additional data (wind magnitude and direction, air temperature,
atmospheric pressure, relative humidity, nebulosity and rainfall, all these variables may vary in
time) in a standard format, see the example “heat_exchange”. The format may be changed but the
user has to change the implementation of the reading and the interpolation of the meteorological
data. When using the complete module, evaporation is calculated by TELEMAC-3D, but the user
has to provide rainfall data with units homogeneous with length over time.
The
main
developments
of
EXCHANGE_WITH_ATMOSPHERE.
this
module
are
implemented
in
the
module
Some physical parameters have been hard-coded in the module, often imposed as the mean
value: the type of sky related to the luminosity of the site (very pure sky, mean pure sky or
industrial zone sky), the type of cloud (cirrus, cirro stratus, alto cumulus, alto stratus, stratus), but
these values can be changed in the module EXCHANGE_WITH_ATMOSPHERE).
Because the site of a study may not be equipped with local wind measurements and these kinds of
data are available at a different location, possibly far from the studied site, a wind function is used.
This is a linear function with a single coefficient of calibration b: f(U2) = b(1+U2) where U2 is the
wind velocity at 2 m high.
To get the wind velocity at 2 m high from classical wind data at 10 m high, a roughness length of
z0 = 0.0002 m has been chosen in the code, that leads to U2 ≈ 0.85U10. This value of 0.85 (or the
roughness length) may be changed by the user if needed.
Examples of solar radiation penetration are given in comments in the SOURCE_TRAC subroutines.
Two laws are suggested: the first one uses the in situ measurements of Secchi length and is
recommended if available; the second one uses two exponential laws that may be difficult to
calibrate and require an estimation of the type of water from turbidity.
Except for the coefficient to model the penetration of solar radiation in the water column, the
parameter b that appears in the wind function is the single calibration parameter of this module. Its
value is given by the keyword COEFFICIENT TO CALIBRATE THE ATMOSPHERE-WATER
EXCHANGE MODEL. (default value = 0.0025 but recommended values are between 0.0017 and
0.0035). This keyword is both used for the linearised formula at the free surface and the model with
complete balance (values 1 and 2 for the keyword ATMOSPHERE-WATER EXCHANGE MODEL).
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
5.6.
Page 59
ASTRAL POTENTIAL
When modelling large maritime areas, it is sometimes necessary to take into account the effect of
astral forces generating tide inside the area. For this, the user has several keywords at his
disposal.
First of all, the logical keyword TIDE GENERATING FORCE (default value NO) allows these
phenomena to be taken into account.
The keyword LONGITUDE OF ORIGIN POINT must be positioned at the right value.
Lastly, the two keywords ORIGINAL DATE OF TIME (format YYYY;MM;DD) and ORIGINAL HOUR
OF TIME (format HH;MM;SS) must be used to give the beginning time of the simulation. This
information is necessary for TELEMAC-3D to compute the respective position of the moon and the
sun.
5.7.
CONSIDERATION OF WAVE DRIVEN CURRENTS
It is possible to take into account the wave driven currents by retrieving information calculated by a
wave propagation module of the TELEMAC modelling system (mainly TOMAWAC but also
ARTEMIS). In TELEMAC-3D, these wave driven currents appear as a source term (wave
stresses). Two ways are possible: taking into account a steady state with the help of a file
previously computed by a wave propagation module or using a coupling with TOMAWAC. The first
procedure is as follows:
• Perform a calculation of wave propagation on the same mesh as the TELEMAC-3D
calculation requesting the storage of the driving forces. In the case of TOMAWAC, the
variables FX and FY,
• Get the wave results file and provide its name through the keyword BINARY DATA FILE 1,
• Activate the keyword WAVE DRIVEN CURRENTS,
• Fill the keyword RECORD NUMBER IN WAVE FILE. This value corresponds to the
iteration number stored in the wave file which must be taken into account by TELEMAC3D. Usually, this is the last iteration stored, by default this number is set to 1. The name of
the variables to read is "FORCE FX" and "FORCE FY", but this can be changed within the
TRISOU subroutine.
In the current version of TELEMAC-3D, this driving force is considered as being constant over the
vertical.
If the user wishes to take into account several results of the wave propagation module (e.g. to take
into account changes in the level of the sea), FORTRAN programming or coupling TELEMAC-3D
with TOMAWAC (see 10.7) are required.
5.8.
OTHER PHYSICAL PARAMETERS
When modelling large domains, the influence of the Coriolis force of inertia has to be taken into
account. This is done by activating the logical keyword CORIOLIS (which is set to NO by default).
In such a case, the value of the Coriolis coefficient (refer to the Theoretical Note) is defined by the
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 60
keyword CORIOLIS COEFFICENT (default value is 0.). The latter should be computed according
to latitude λ through the formula:
•
2 ω sin (λ ) where ω is the Earth’s rotational velocity of 7.27 x 10-5 rad/s and λ is the
average latitude of the model.
In the case of very large domains such as portions of oceans, it is necessary to carry out a
simulation with spherical coordinates, in which case the Coriolis coefficient is adjusted
automatically at each point of the domain (see 10.3), by activating the keyword SPHERICAL
COORDINATES.
TELEMAC-3D additionally offers the opportunity to set the acceleration due to gravity (keyword
GRAVITY ACCELERATION which, by default, is set to 9.81 m/s²).
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
6.
Page 61
NUMERICAL SETUP OF THE COMPUTATION
The numerical setup is comparatively common to a hydrodynamic computation alone or with a
tracer. Thus, in the following sections of this chapter, the numerical parameters as applied to the
solution of a tracer equation are integrated into the hydrodynamic parameters.
6.1.
GENERAL SETUP
TELEMAC-3D solves the Navier-Stokes equations in several stages, possibly through the three
stages of the fractional step method (see the theoretical note). The first stage consists in finding out
the advected velocity components by only solving the convection terms of the momentum
equations. The second stage computes, from the advected velocities, the new velocity components
by taking into account both diffusion and source terms of the momentum equations. These two
solutions enable to get an intermediate velocity field. The third stage computes the water depth
from the vertical integration of the continuity equation and momentum equations only including the
pressure-continuity terms. This step is called the propagation step.
The user can activate or deactivate, either globally or individually, some of these stages.
6.1.1.
ADVECTION STEP
Whether the convection terms will be considered or not will be determined by means of the logical
keyword ADVECTION STEP (default value YES). However, even though that keyword is set to
YES, then some advection terms can be deactivated by means of the following complete keywords
(value 0 = "NO ADVECTION "):
•
SCHEME FOR ADVECTION OF VELOCITIES: for the advection of velocities,
•
SCHEME FOR ADVECTION OF DEPTH: for taking the advection of depth into account,
•
SCHEME FOR ADVECTION OF K-EPSILON: for the advection of power and turbulent
dissipation,
•
SCHEME FOR ADVECTION OF TRACERS: for the advection of tracers (one value per
tracer).
See section 6.2 for more information on the possible choices.
6.1.2.
DIFFUSION STEP
Whether the diffusion terms are taken into account or not is established by means of the logical
keyword DIFFUSION STEP (default value YES). However, even though that keyword is set to
YES, then some diffusion terms can be deactivated by means of the following complete keywords
by the following keywords:
RELEASE 7.1
•
SCHEME FOR DIFFUSION OF VELOCITIES,
•
SCHEME FOR DIFFUSION OF TRACERS (one value per tracer),
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
•
Page 62
SCHEME FOR DIFFUSION OF K-EPSILON.
The 0 value at each keyword cancels the diffusion, whereas the 1 value (default value) leads to the
implicit calculation of diffusion.
For the treatment of the diffusion, two choices are possible for the keyword OPTION FOR THE
DIFFUSION. The default value 1 means an implicit treatment whereas the value 2 means an
uncoupled treatment between the horizontal diffusion and the vertical diffusion.
6.1.3.
PROPAGATION STEP
The velocity and water depth propagation processes are taken into account by the logical keyword
PROPAGATION STEP (default value YES).
Since the version 6.0 of TELEMAC, the propagation step can only be treated with TELEMAC-2D
"wave equation” option.
The propagation step can be linearized by activating the keyword LINEARIZED PROPAGATION
(default value = NO) particularly when one conducts a case study for which an analytical solution in
the linearized case is available. It is then necessary to set the water depth around which the
linearization is performed by means of the keyword MEAN DEPTH FOR LINEARIZATION (default
value 0).
With the « wave equation » option, the keyword FREE SURFACE COMPATIBILITY GRADIENT
can be used. A value lower than 1. (which is the default value) makes it possible to delete the
spurious oscillations of the free surface, but slightly alters the consistency between the water depth
and the velocities in the continuity equation.
When using the non-hydrostatic version, it is possible to enable the dynamic pressure term in
treatment of the wave equation. This option is enabled by the logical keyword DYNAMIC
PRESSURE IN WAVE EQUATION but it leads a dilemma:
• If the keyword is set to YES, the dynamic pressure gradient is taken into account when
calculating the evolution of the water depth (advantage), but the evolution of the water
depth is not known when calculating the dynamic pressure (disadvantage),
• If the keyword is set to NO (which is the default value), the dynamic pressure is calculated
taking into account the evolution of the water depth (advantage) but it is not taken into
account in the calculation of the water depth (disadvantage).
Therefore, when the effect of the dynamic pressure is important (nonlinear waves), it is
recommended to set this keyword to YES in combination with setting sub-iterations for
nonlinearities (see test case 310_NonLinearWave).
6.2.
THE ADVECTION SCHEME
The procedure for taking the advection terms into account is individualized for each of the variables
liable to be processed. It has previously been explained that the zero option corresponds to a
deactivation of the term.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
6.2.1.
Page 63
ADVECTION OF THE THREE-DIMENSIONAL VARIABLES
The advection schemes of the three-dimensional variables (i.e. all the variables but the water
depth) are (for further details about these options, the reader shall refer to the Theoretical Note):
•
•
•
•
•
•
•
•
0: Deactivation,
1: Method of characteristics. That method involves mutually independent advection and
diffusion steps. The method consists in writing that the value of the advected variable is
equal to the value of the same variable in the previous instant traced back on the path
travelled during the time step,
2: Explicit scheme + SUPG (Streamline Upwind Petrov Galerkin). That method uses test
functions which are deformed in the direction of the current for the variational method,
3: Explicit Leo Postma scheme,
4: Explicit scheme + MURD (Multidimensional Upwind Residual Distribution) N scheme,
5: Explicit scheme + MURD PSI scheme,
13: Explicit Leo Postma scheme for tidal flats,
14: Explicit scheme + MURD (Multidimensional Upwind Residual Distribution) N scheme
for tidal flats.
The latter five schemes are primarily recommended for the tracers, since they are advantageous in
being conservative and monotonic, i.e. they do not generate any numerical oscillation. On the other
hand, they are more diffusive than SUPG. In that respect, scheme 5 is an improvement to scheme
4, being less diffusive perpendicularly to the flow but quite obviously somewhat more computation
time-consuming. Nevertheless, it is still globally less time-consuming than SUPG.
Distributive schemes (options 3, 4, 5, 13 or 14) are schemes whose stability is conditioned by a
Courant number less than 1. When using one of these schemes, at each time step, TELEMAC-3D
performs a test for checking the Courant number point by point. In case it exceeds the value 1,
TELEMAC-3D will automatically execute sub-iterations to satisfy the stability criterion. However, if
the number of sub-iterations exceeds 100, TELEMAC-3D considers that the treatment of the
advection term is problematic and the calculation is interrupted, printing an error message in the list
control.
The default value for both velocity and k-ε is 1. That value is advisable, since it is satisfactory in
many instances and is by far the fastest. On the contrary, the default value for the tracers is 5. It is
the most reliable scheme, since the "mass" conservation of the active tracers is a frequently
essential point in TELEMAC-3D.
Since version 6.0, the value concerning the water depth advection scheme is ignored by
TELEMAC-3D. The optimum advection scheme is automatically selected by the software
(conservative scheme).
According to the schemes used, the mass conservation can be improved by performing subiterations. That consists in a updating, for one given time step, both advective field and propagative
field during several sub-iterations. Upon the first sub-iteration, the field of velocities is yielded by
the results achieved during the previous time steps. Through that procedure, the non-linearities can
be taken into account in a better way and the mass conservation can be significantly improved in
the cases of schemes 2 and 3. The number of sub-iterations is set by the keyword NUMBER OF
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 64
SUB ITERATIONS FOR NON LINEARITIES the default value of which is 1 (also refer to the
Theoretical Note).
The SUPG scheme can be configured using specific keywords (see § 6.2.2).
6.2.2.
CONFIGURATION OF THE SUPG SCHEME
When using the SUPG method, the user must determine the type of upwinding desired using the
keyword SUPG OPTION which is an array of 4 integers related, to the velocity, the water depth,
tracers and k-ε model respectively.
The possible values are:
•
•
•
0: no upwinding,
1: upwinding with the conventional SUPG method, that is to say, the upwinding is 1,
2: upwinding with the modified SUPG method, that is to say, the upwinding is equal to the
Courant number.
In theory, option 2 is more accurate when the Courant number is less than 1, but should not be
used when the latter is important. Thus, option 2 should be used in models for which the Courant
number remains low. If the Courant number cannot be estimated, it is strongly recommended to
use option 1 (which can be considered as the most universal).
6.2.3.
CONFIGURATION OF THE WEAK CHARACTERISTICS
When choosing the method of characteristics, two forms can be used with the keyword OPTION
FOR CHARACTERISTICS:
•
•
1: the strong form (by default),
2: the weak form.
None of them are recommended for the advection of tracers because they are not mass
conservative. The weak form will decrease the diffusion. If the keyword MASS-LUMPING FOR
WEAK CHARACTERISTICS = 1. (default value = 0. i.e. no mass-lumping), monotonicity of the
scheme appears. This weak form should be more conservative than the strong form. The NUMBER
OF GAUSS POINTS FOR WEAK CHARACTERISTICS defines the number of Gauss points used
to compute the weak characteristics. Possible choices are 1, 3 (default value) and 6. The bigger
the number is, the more conservative the scheme is, but the higher the computational costs are.
6.3.
SPECIFIC PARAMETERS IN THE NON-HYDROSTATIC VERSION
The application of the software NON-HYDROSTATIC VERSION (default = NO) requires that
complementary keywords be defined.
An equation for vertical velocity is initially solved in a same way as the U and V components. That
equation is only written with the hydrostatic pressure which, under that hypothesis, is cancelled out
with the gravity term. The convection scheme in that equation is identical to that chosen for U and
V (keyword SCHEME FOR ADVECTION OF VELOCITIES). The solution of the linear system (for
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 65
the diffusion step integrating the advection terms or not) is managed by the various following
keywords:
•
•
•
•
SOLVER FOR VERTICAL VELOCITY refer to paragraph 6.5.1 below,
MAXIMUM NUMBER OF ITERATIONS FOR VERTICAL VELOCITY refer to paragraph
6.5.2 below,
ACCURACY FOR VERTICAL VELOCITY refer to paragraph 6.5.2 below,
PRECONDITIONING FOR VERTICAL VELOCITY refer to paragraph 6.5.3 below.
Afterwards, TELEMAC-3D solves a Poisson equation for the dynamic pressure. The dynamic
pressure gradient plays the part of a correction providing the required zero divergence on velocity.
The solution of the linear system in that equation is managed by the following keywords:
•
•
•
•
•
SOLVER FOR PPE refer to paragraph 6.5.1 below,
MAXIMUM NUMBER OF ITERATIONS FOR PPE refer to paragraph 6.5.2 below,
ACCURACY FOR PPE refer to paragraph 6.5.2 below,
OPTION OF SOLVER FOR PPE refer to paragraph 6.5.1 below,
PRECONDITIONING FOR PPE refer to paragraph 6.5.3 below.
Once that pressure is computed, all the velocity components are updated with the dynamic
pressure gradient which will ensure the zero divergence condition. Updating that "solenoidal"
velocity does not require that a linear system be solved.
6.4.
IMPLICITATION
Apart from the terms of the time derivative, the unknowns f
components) can be considered in both extreme instants t
n
(the velocity and water depth
(the equation is then referred to as
n +1
explicit) or t
(the equation is then referred to as implicit). Strictly speaking, and for a 2-order
solution in time, an approach consists in considering the terms in the intermediate
instant ( f
n
+f
n +1
) / 2 . Practically, the latter approach is unstable and it becomes necessary to
define an implicitation coefficient for which the unknowns are actually discretized in time in the
following form:
θf
n +1
+ (1 − θ ) f n .
The implicitation coefficients are theoretically always higher than 0.5 (0.55 or 0.6 will generally yield
good results).
The user can use the keyword IMPLICITATION FOR VELOCITIES (default value: 1.) which defines
the value of the θ u coefficient for the velocity components. The keyword IMPLICITATION FOR
DEPTH (default value: 0.55) is provided to set the value of "propagation height" multiplying
coefficient θ h . Lastly, in order to make the construction of the various numerical schemes more
versatile, a diffusion-specific coefficient
θ ud
is provided (keyword IMPLICITATION FOR
DIFFUSION (default value: 1.) and can be different from θ u ).
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
6.5.
Page 66
SOLUTION OF THE LINEAR SYSTEMS
Both discretization and variational formulation of the equations lead to a linear system which now
has to be solved. The direct solution methods are often not suitable and are overly expensive as
soon as there are many unknowns. Thus, the main ways developed in TELEMAC-3D consist in
solving the linear systems by means of iterative solvers. Nevertheless for specific applications, a
direct solver can be used.
6.5.1.
SOLVERS
According to the relevant numerical parameters, various linear systems are liable to be solved. The
solver used for solving one of these systems can be selected by the user through the following
keywords:
•
•
•
•
•
•
SOLVER FOR DIFFUSION OF VELOCITIES (default value: 1),
SOLVER FOR PROPAGATION (default value: 1),
SOLVER FOR PPE (default value: 1),
SOLVER FOR VERTICAL VELOCITY (default value: 1),
SOLVER FOR DIFFUSION OF TRACERS (default value: 1, one value for each tracer),
SOLVER FOR DIFFUSION OF K-EPSILON (default value: 1).
Each of these keywords can assume a value ranging from 1 to 8, which values correspond to the
following possibilities:
•
•
•
•
•
•
•
•
1: conjugate gradient method (when the matrix of the system to solve is symmetric),
2: conjugate residual method,
3: conjugate gradient method on a normal equation,
4: minimum error method,
5: square conjugate gradient method,
6: CGSTAB (stabilized conjugate gradient) method,
7: GMRES (Generalised Minimum RESidual) method,
8: direct solver.
The GMRES method is well suited for improperly conditioned systems. This method requires that
the dimension of the Krylov space be defined. That parameter is set by means of the following
keywords:
•
•
•
•
•
OPTION OF SOLVER FOR DIFFUSION OF VELOCITIES,
OPTION OF SOLVER FOR PROPAGATION,
OPTION OF SOLVER FOR PPE,
OPTION OF SOLVER FOR DIFFUSION OF TRACERS,
OPTION OF SOLVER FOR DIFFUSION OF K-EPSILON.
The default values are set to 3. The larger that parameter is, the higher the memory requirements
and the number of matrix-vector products per iteration (and consequently the computational time
as well) are, but the better the convergence is.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
6.5.2.
Page 67
ACCURACIES
The principle of iterative methods consists in getting gradually closer to the solution of the problem
during the iterations. The systems to be solved imply a relative accuracy within a range from 10-4 to
10-10 with a restricted number of iterations. Both accuracy and maximum number of iterations
should be set for each system.
Accuracy is specified by the following keywords:
•
•
•
•
•
•
ACCURACY FOR DIFFUSION OF VELOCITIES (default value: 1.E-5),
ACCURACY FOR PROPAGATION (default value: 1.E-6),
ACCURACY FOR PPE (default value: 1.E-4),
ACCURACY FOR VERTICAL VELOCITY (default value: 1.E-6),
ACCURACY FOR DIFFUSION OF TRACERS (default value: 1.E-6),
ACCURACY FOR DIFFUSION OF K-EPSILON (default value: 1.E-6).
The maximum number of iterations is specified by the following keywords:
•
•
•
•
•
•
•
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF VELOCITIES (default value:
60),
MAXIMUM NUMBER OF ITERATIONS FOR PROPAGATION (default value: 200),
MAXIMUM NUMBER OF ITERATIONS FOR PPE (default value: 100),
MAXIMUM NUMBER OF ITERATIONS FOR VERTICAL VELOCITY (default value: 100),
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF TRACERS (default value:
60),
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF K-EPSILON (default value:
200),
MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION SCHEMES (default value: 10)
only for schemes 13 and 14.
The user automatically gets information about the solvers upon each listing printout. The
information provided in the listing can be of two types:
6.5.3.
•
Either the treatment converged before reaching the maximum allowable number of
iterations, and then TELEMAC-3D will provide the number of actually performed iterations,
as well as the achieved accuracy,
•
Or the treatment did not converge early enough. TELEMAC-3D will then provide the
message "MAXIMUM NUMBER OF ITERATIONS IS REACHED" and give the actually
achieved accuracy. In some cases, and if the maximum number of iterations is already set
to a high value (for instance over 100), convergence can then be improved by reducing the
time step or, quite often, by improving the quality of the mesh.
PRECONDITIONINGS
The iterative methods are sensitive to the "conditioning" of matrices, so that a complementary
preconditioning is necessary in order to reduce the number of iterations for getting a prescribed
accuracy.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 68
TELEMAC-3D offers several opportunities for preconditioning. The selection is made by means of
the following keywords:
•
•
•
•
•
•
PRECONDITIONING FOR DIFFUSION OF VELOCITIES,
PRECONDITIONING FOR PROPAGATION,
PRECONDITIONING FOR PPE,
PRECONDITIONING FOR VERTICAL VELOCITY,
PRECONDITIONING FOR DIFFUSION OF TRACERS (one value for each tracer),
PRECONDITIONING FOR DIFFUSION OF K-EPSILON.
The available options are:
•
•
•
•
•
•
•
•
•
•
•
0:
no preconditioning,
2:
diagonal preconditioning,
3:
diagonal preconditioning with the condensed matrix,
5:
diagonal preconditioning with absolute values,
7:
Crout preconditioning per element (downgraded in parallel),
11: Gauss-Seidel preconditioning per element (downgraded in parallel),
13: preconditioning matrix is provided by the user,
14: cumulated diagonal preconditioning and Crout preconditioning per element,
17: preconditioning through direct solution along each vertical direction,
21: cumulated diagonal preconditioning with the condensed matrix and Crout
preconditioning per element,
34: cumulated diagonal preconditioning with direct solution along each vertical direction.
The default value is 2 for all the preconditionings. Some preconditionings can be cumulated,
namely the diagonal preconditionings with other ones. Since the basic values are prime numbers, a
couple of preconditionings are cumulated by giving to the keyword the value of the product of the
two preconditionings which one wants to cumulate (e.g. numbers 14 and 21).
6.6.
TIDAL FLATS
TELEMAC-3D offers several treatment options as regards the tidal areas.
First, if the user has ascertained that his/her model has no tidal area throughout the simulation, the
processing of such areas can be deactivated by setting the keywords TIDAL FLATS to NO (the
default value is YES). That option makes it possible to save computational time (by deleting the
tidal flat testing).
The tidal flats can be treated in two different ways:
RELEASE 7.1
•
In the first case, the equations are treated all over the domain and in a thorough way. The
tidal areas are detected and such terms as the free surface gradient (in the absence of
water, the free surface gradient becomes the bottom gradient and generates spurious
motive terms) are corrected in them,
•
In the second case, the tidal areas are withdrawn from the computation. The exposed
elements are always part of the mesh, but all their contributions to the computations are
cancelled by a so-called "masking" array. Thus, the data structure and the computations
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 69
remain formally unchanged, to within the masking coefficient. That method, however,
raises issues as regards the mass conservation and the exposure and coverage dynamics.
The treatment will be selected by means of the keyword OPTION FOR THE TREATMENT OF
TIDAL FLATS which can be set either to 1 or 2, the default value being 1.
The treatment of negative depths can be specified using the keyword TREATMENT OF NEGATIVE
DEPTHS. A value of 1 (default), consists in a conservative smoothing of negative depths. The
second option is to limit the flux between the elements to ensure strictly positive water depths. This
second option should be used with advection schemes consistent with tidal flats (+ mass lumping
option for height at value 1.). The value 0 means that no special treatment is performed.
The keyword MINIMAL VALUE FOR DEPTH the default value of which is -1,000 enables to set the
threshold below which the smoothing is. For example, MINIMAL VALUE FOR DEPTH set to 0.01
means the minimum depth is 1 cm.
The following three keywords are for setting, after coverage, the value of the variable which has
been masked:
•
•
•
TREATMENT ON TIDAL FLATS FOR VELOCITIES,
TREATMENT ON TIDAL FLATS FOR TRACERS,
TREATMENT ON TIDAL FLATS FOR K-EPSILON.
The available options for these keywords are:
•
•
0: that option corresponds to a setting to zero of the variable on the element (default
value),
1: that option sets to its prior-to-masking value.
The keyword THRESHOLD FOR VISCOSITY CORRECTION ON TIDAL FLATS (default value is
0.2) allows to specify the minimum water depth from which the viscosity is gradually reduced (see
programming within the VISCLIP subroutine).
When the three-dimensional mesh has crushed levels (null water depth or fixed level "hitting" the
bottom), it is recommended to activate a specific treatment that prevents the transfer of very small
amounts of water at the calculation points which have no volume (this situation also tends to
degrade the mass conservation when using distributives PSI and N schemes). This algorithm is
activated with the keyword logical BYPASS VOID VOLUMES. When using PSI and N schemes
compatible with tidal flats, the option is automatically enabled, even if the keyword is set to NO.
6.7.
HYDROSTATIC INCONSISTENCIES
Hydrostatic inconsistencies (linked to the truncature errors in the computation of the buoyancy
terms) are liable to occur on the nearly-zero volume prisms. The keyword HYDROSTATIC
INCONSISTENCY FILTER (default = NO) is provided for the forces caused by the spurious
horizontal pressure gradients and the various diffusion coefficients on those prisms where at least
one of the lower base nodes has a higher elevation than one of the upper base nodes.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 70
6.8.
OTHER PARAMETERS
6.8.1.
MASS-LUMPING
Upon the solution of the linearized system, TELEMAC-3D makes it possible to perform a masslumping on the mass matrices. That procedure consists in partly or wholly aggregating the mass
matrix on its diagonal and enables to substantially shorten the computational times. The resulting
solutions, however, become smoothed, except for in steady flow conditions in which they are
unchanged. The mass-lumping rate is set by means of the keywords MASS-LUMPING FOR
DEPTH, MASS-LUMPING FOR VELOCITIES, MASS-LUMPING FOR DIFFUSION and MASSLUMPING FOR WEAK CHARACTERISTICS. Value 1. means maximum mass-lumping (the mass
matrices are diagonal), value 0. (default value) corresponds to the normal treatment without any
mass-lumping. For further details, the reader shall refer to the TELEMAC-3D Theoretical Note.
6.8.2.
CONVERGENCE AID
Another way to speed up the system convergence when solving the propagation step consists in
acting upon the initial solution rather than the matrix proper. To that purpose, the initial value being
set for h (actually, the unknown h
is replaced by the incrementationδh = h
− h ) is modified
at the beginning of computation, the user can take action at the keyword INITIAL GUESS FOR
DEPTH which can assume the following values:
n +1
6.8.3.
n +1
n
•
0: the initial value of δh = h
•
1: the initial value of
value),
•
2: δh = 2 h − δh
where δh is the value of δh at the previous time step, and δh
the value of δh two time steps earlier. It is actually an extrapolation.
n
n +1
− h n is zero,
δh is equal to the value of δh at the previous time step (default
n −1
n
n −1
is
MATRIX STORAGE
TELEMAC-3D provides a couple of procedures for storing the various matrices it has to handle,
namely the conventional EBE (Element By Element) method and the segment-wise storage. The
latter procedure is faster (by approx 20%) in most cases.
The choice among these two types of storage is made by means of the keyword MATRIX
STORAGE qui can assume the following values:
•
•
6.8.4.
1: classical EBE method,
3: edge-based storage (default value).
VELOCITIES PROJECTION
At the end of the time loop, it is possible to do a treatment which aim is to cancel the component of
the normal velocity at the bottom or the normal velocity on the solid lateral walls. This check is
activated with the logical keywords VELOCITY PROJECTED ON SOLID LATERAL BOUNDARIES
and VELOCITY PROJECTED ON BOTTOM.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 71
These two options are activated by default.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
7.
Page 72
TRACER TRANSPORT
The TELEMAC-3D software makes it possible to take into account the transport of passive or
active tracers (active tracers affect the hydrodynamics), being either conservative or not.
This chapter discloses the tracer transport features.
The maximum number of tracers is set to 20 by default but it can be changed by the user with the
keyword MAXIMUM NUMBER OF TRACERS. This avoids changing the previously hardcoded
values (until version 7.0), which required recompiling the whole package.
7.1.
GENERAL SETUP
The number of tracers is defined by the keywords NUMBER OF TRACERS. If that number is set to
zero (default value), then the tracers will not be taken into account by TELEMAC-3D.
In addition to the number of tracers, the user should enter the NAMES OF TRACERS. A tracer
name should be written with 32 characters (16 for the name and 16 for the unit).
NUMBER OF TRACERS: 2
NAMES OF TRACERS: 'TEMPERATURE
7.1.1.
°C
', 'SALINITY
G/L'
PRESCRIBING THE INITIAL CONDITIONS
If the initial values of tracers are constant all over the domain, just insert, into the steering file, the
keyword INITIAL VALUES OF TRACERS with the desired value(s) separated with ; if more than
one.
In more complex cases, an action shall be taken directly at the CONDIM subroutine, in the same
fashion as described in paragraph 4.2.2 dealing with the initial hydrodynamic conditions.
When resuming a computation, the initial condition of tracers corresponds to the condition of the
last time step which was stored into the restart file. The tracer management sequence order during
the previous computation needs to be well known and the same sequence order shall be followed
in the computational suite in order to prevent on confusion. If the restart file includes no information
about the tracer, TELEMAC-3D will use the value as set by the keyword INITIAL VALUES OF
TRACERS.
7.1.2.
PRESCRIBING THE BOUNDARY CONDITIONS
The tracer boundary conditions are prescribed according to the same principle as the
hydrodynamic boundary conditions (see section 4.3).
The boundary condition type will be yielded by the value of LITBOR in the boundary conditions file.
In case of an entering liquid boundary with a prescribed tracer (the value of LITBOR of 5), then the
tracer value can be yielded in various ways:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 73
•
If that value is constant both along the boundary and in time, then it is provided in the
steering file by means of the keyword PRESCRIBED TRACERS VALUES. The writing
convention is as follows: value of tracer 1 at boundary 1, value of tracer 2 at boundary 1,
…, value of tracer N at boundary 1, value of tracer 1 at boundary 2, value of tracer 2 at
boundary, …, value of tracer N at boundary 2, etc. The boundary order is the same as in
the case of hydrodynamic boundary conditions.
•
If the value is constant in time but varies along the boundary, it will be set directly by the
TBOR variable in the BOUNDARY CONDITIONS FILE,
•
If the value is constant along the boundary but varies in time, the user may either use the
LIQUID BOUNDARIES FILE or take action at the function TR3. The latter will be
programmed somewhat like the functions VIT3, Q3 and SL3. Note that the liquid
boundaries file will not be taken into account if the keywords PRESCRIBED ELEVATIONS
and PRESCRIBED FLOWRATES do not appear in the steering file.
T
s
0
1040400
SL(1)
m
0.47
0.57
TR(2,1)
°C
24.7
28.0
TR(2,2)
°C
38.0
36.7
In the above example of a liquid boundaries file, for the indices of tracers values, the first value is
the number of the liquid boundary and the second the number of the tracer.
•
If the value varies in both time and space, the user should then modify the BORD3D
subroutine, at the part regarding the tracer.
The keyword TREATMENT OF FLUXES AT THE BOUNDARIES enables, during the convection
step (with the SUPG, PSI and N schemes), to set a priority among the tracer flux across the
boundary and tracer value at that wall. Option 2 ("Priority to fluxes") will then induce a change in
the tracer prescribed value, but will bring about a good assessment of the "mass" of tracer passing
across the boundary. On the other hand, option 1 ("Priority to prescribed values", default value)
sets the tracer value without checking the fluxes.
Finally, in the case of an input boundary, it is possible to specify a concentration profile in the
vertical using the keyword TRACERS VERTICAL PROFILES. The options are:
•
•
•
•
0: user programming (in BORD3D),
1: constant profile (default),
2: constant profile (tracer diluted) or Rouse profile (sediment),
3: Rouse (normalised) and imposed concentration.
7.2.
PHYSICAL SETUP
7.2.1.
ACTIVE TRACERS
The active tracers affect the flow through the hydrostatic pressure gradient term. As a rule, indeed,
the pressure is written as:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 74
Zs
p = p h + p d = ρg ( Z s − z ) + ρ 0 g ∫
z
∆ρ
ρ0
dz + pd
where ph and pd denote the hydrostatic pressure and the dynamic pressure, respectively.
Thus, two elements should be defined, namely:
The term
ρ0
ρ0
and
∆ρ
.
ρ0
is defined by the keyword AVERAGE WATER DENSITY. The default value is 1,025;
it corresponds to sea (ocean) water.
The second term, which operates in the buoyancy source terms, directly depends on the values of
the active tracers and is defined by the keyword DENSITY LAW.
The available values for that keyword DENSITY LAW are:
•
•
•
•
•
0: no interaction with the tracers (default value),
1: variation of density according to temperature,
2: variation of density according to salinity,
3: variation of density according to temperature and salinity,
4: variation as a function of the spatial expansion coefficients.
With the 1-3 options, the variations are given by the law as defined in TELEMAC-3D. In such a
case, the name of the salinity tracer (expressed in kg/m3) shall necessarily begin with SALINI and
the name of the temperature tracer (in °C) shall begin with TEMPER.
With option 4, the term
∆ρ
is described by a linear function of the Ti tracer of the type:
ρ0
∆ρ
= −∑ β i (Ti − Ti 0 ) i
ρ0
i
The
βi
coefficients (spatial expansion coefficients) are set by the values of the keyword BETA
EXPANSION COEFFICIENT FOR TRACERS (default = 0.). They can be either positive
0
(temperature) or negative (salinity, suspended sediment). The values Ti are defined by the values
of the keyword STANDARD VALUES FOR TRACERS (default = 0.).
The user shall enter the expansion coefficients, the standard values and the tracer names in the
same sequence order to ensure that each tracer will have the correct parameters.
7.2.2.
PUNCTUAL SOURCE TERMS
For each source, the user shall enter the tracer value at the sources by means of the keyword
VALUE OF THE TRACERS AT THE SOURCES. Thus, it is an array of reals specifying the
concentration of tracers at the source. The writing convention is as follows: source value 1 of tracer
1; source value 2 of tracer 1; …, source value n of tracer 1; source value 1 of tracer 2; …, source
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 75
value n of tracer 2 etc. In case of time dependent value, it is necessary to use the source file or a
specific FORTRAN programming (function T3D_TRSCE). The user is invited to read the subsection
5.4 to deal with the hydrodynamic part of punctual source terms (flow rates, velocities).
7.2.3.
GENERAL SOURCE TERMS
If one wants to take the tracer generation or disappearance source terms into account, then this
has to be implemented within the SOURCE_TRAC subroutine.
7.3.
NUMERICAL SETUP
As with hydrodynamics, the advection schemes SCHEME FOR ADVECTION OF TRACERS (refer
to paragraph 6.1.1) and the diffusion schemes SCHEME FOR DIFFUSION OF TRACERS (refer to
paragraph 6.1.2) can be modified.
Otherwise, the horizontal and vertical diffusion of the tracers can be set with the keywords
COEFFICIENT FOR HORIZONTAL DIFFUSION OF TRACERS and COEFFICIENT FOR
VERTICAL DIFFUSION OF TRACERS. Their default values are 10-6 m2/s. These two keywords are
arrays since version 7.1, with one value per tracer, separated by semicolons, so that different
values can be given for different tracers.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
8.
Page 76
DROGUES
During a hydrodynamic simulation, TELEMAC-3D offers an opportunity to follow the paths of a
number of particles (drogues) which are released into the fluid from discharge points. The result is
provided in the form of an ASCII file (directly readable by TECPLOT) which contains the various
positions of the drogues.
8.1.
CONFIGURATION OF SIMULATION
Three parameters should be entered into the steering file. First, the user should mention the
number of drogues by means of the keywords NUMBER OF DROGUES the default value of which
is 0. Secondly, the user should enter the name of the file into which TELEMAC-3D will store the
successive positions of the drogues. This is defined through the keyword DROGUES FILE. Lastly,
the user can configure the printout period within that file by means of the keyword PRINTOUT
PERIOD FOR DROGUES (default value = 1). That value is expressed in a number of time steps
and is quite independent from of the printout period of the other results in TELEMAC-3D.
The subroutine ADD_PARTICLE is called within the FLOT3D subroutine to set the initial values of
variables XFLOT, YFLOT, ZFLOT and TAGFLO, which are the three-dimensional coordinates of the
release point and an identifier of the particle for each drogue. The release time step is also to be
amended. A commented example can be found within the subroutine FLOT3D or in the example
“particles”. This subroutine is to be inserted in the FORTRAN file.
The example below illustrates the programming of the steering file in the case of two drogues being
release at different times.
NUMBER OF DROGUES
DROGUES FILE
=2
= './drogues'
PRINTOUT PERIOD FOR DROGUES
8.2.
= 10
VISUALISATION OF RESULTS
The results are stored in the file specified by the keyword DROGUES FILE. This is an ASCII file
written in a format compatible with the TECPLOT software. If this tool is not available, it is quite
easy to get the coordinates of the different drogue positions to export them to another viewer. It is
also possible to develop a new drogue output format, but this must be done in the subroutine
DERIVE.
Within TECPLOT, in order to add the TECPLOT DROGUES FILE to TELEMAC result data already
loaded, select “Add to current data” set in the Load Data File Warning dialogue (cf. Figure
8). The Load Data File Warning dialogue will appear after the user has selected the file and zones
and/or variables to load.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 77
Figure 8: Load Data File Warning dialogue in TECPLOT
Once you have loaded your data with TECPLOT, the drogue positions will be considered as
“Scatter plots” by TECPLOT Software. Scatter plots are plots of symbols centered at the
data points in a field. To add a scatter layer to your plot, activate the “Scatter” toggle in the
Sidebar. To be visible in your plot, the Scatter layer which contains the TECPLOT DROGUES FILE
must be turned on and the Scatter layer containing the 3D RESULT FILE data must be turned off.
This is done by selecting “Yes” or “No” from the [Scat Show] button drop-down menu on the
Scatter page of the Zone Style dialogue (cf. Figure 9).
Figure 9: Zone Style dialogue in TECPLOT
Then, you can modify your Scatter plot using the Scatter page of the Zone Style dialogue and the
Scatter submenu of the Plot menu. You can control any of the following attributes for a zone or
group of zones from the Scatter page of the Zone Style dialogue. For complementary information
on the TECPLOT plot procedure, the reader may refer to the TECPLOT user guide.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
9.
Page 78
OIL SPILL MODELLING
During a hydrodynamic simulation, TELEMAC-3D offers an opportunity to follow the paths of an oil
spill which is released into the fluid from a discharge point.
The oil spill model introduced here combines an Eulerian and a Lagrangian approach. The
Lagrangian model simulates the transport of an oil spill near the surface. The oil slick is
represented by a large set of hydrocarbon particles. Each particle is considered as a mixture of
discrete non-interacting hydrocarbon components. Particles are therefore represented by
component categories (soluble and unsoluble components), and the fate of each component is
tracked separately. Each particle has associated to it, amongst other properties, an area, a mass,
its barycentric coordinates within the element it is located in, and the physico-chemical properties
of each of its components. The model accounts for the main processes that act on the spilled oil:
advection, effect of wind, diffusion, evaporation and dissolution. Though generally considered to be
a minor process, dissolution is important from the point of view of toxicity. To simulate soluble oil
component dissolution in water, an Eulerian advection-diffusion model is used. The fraction of each
dissolved component is represented by a tracer whose mass directly depends on the dissolved
mass of oil particles. The hydrodynamic data required for either Lagrangian and Eulerian transport
approach are provided by the TELEMAC-3D hydrodynamic model. The oil spill theoretical
background is explained in [7].
The result of the oil spill modelling is provided in the form of a TECPLOT formatted file which
contains the various positions of the oil particles and a TELEMAC-3D result file (SERAFIN format)
storing the oil dissolved components in water column during the computation.
9.1.
INPUT FILES
In addition to the minimum set of input files necessary to run a TELEMAC-3D case, an oil spill
computation also needs an oil spill steering file. Furthermore, to run an oil spill model the
subroutine OIL_FLOT needs to be modified in the FORTRAN file.
9.2.
STEERING FILE
In addition to the necessary information for running the TELEMAC-3D hydrodynamic model the
following essential information must be specified in the TELEMAC-3D steering file to run an oil spill
propagation model:
• The use of the oil spill model must be declared: OIL SPILL MODEL (= YES, default = NO),
• The name of the oil spill steering file which contains the oil characteristics: OIL SPILL
STEERING FILE (= name chosen by the user),
• The number of oil releases during oil spill: NUMBER OF DROGUES (= number chosen by
the user),
• The frequency of the drogues printout period: PRINTOUT PERIOD FOR DROGUES (=
number chosen by the user),
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 79
• The name of the TECPLOT oil file containing the oil displacement: DROGUES FILE (= name
chosen by the user, default = 1).
With the oil spill module, it is possible to take into account the transport of soluble oil components
in water (whose presence has no effect on the hydrodynamics). These may or may not be diffused
within the flow but their characteristics have to be defined in the OIL SPILL STEERING FILE. If
these components are allowed to diffuse in the flow, they are then treated with the tracer transport
computations of TELEMAC-3D. This implies that the NUMBER OF TRACERS must be set to the
number of the oil soluble components. In addition the TRACER keywords, described in chapter 7,
can be specified.
9.3.
OIL SPILL STEERING FILE
As seen previously, the OIL SPILL STEERING FILE name is given by the user in the TELEMAC
steering file. This file contains all the informations for an oil spill calculations based on the
composition considered by the user, i.e.:
• The number of unsoluble components in oil,
• The parameters of these components such as the mass fraction (%) and boiling point of
each component (K),
• The number of soluble components in oil,
• The parameters of these components such as the mass fraction (%), boiling point of each
component (K), solubility (kg.m-3) and the mass transfer coefficient of the dissolution and
volatilization phenomena (m.s-1)
• The oil density,
• The oil viscosity (m2.s-1),
• The volume of the spilled oil (m3),
• The water surface temperature (K),
• The spreading model chosen by the user:
1. Fay’s model,
2. Migr’Hycar model,
3. Constant area model.
WARNING:
• The parameters of soluble (or unsoluble) components need to be informed only if the
number of these components is not null
• If the sum of all mass fraction components is not equal to 1, the run is interrupted and the
following error message is displayed:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 80
”WARNING::THE SUM OF EACH COMPONENT MASS FRACTION IS NOT EQUAL TO 1.”
“PLEASE, MODIFY THE INPUT STEERING FILE”
An example of the oil spill steering file is given.
NUMBER OF UNSOLUBLE COMPONENTS IN OIL
6
UNSOLUBLE COMPONENTS PARAMETERS (FRAC MASS, TEB)
5.1D-02
,402.32D0
9.2D-02
,428.37D0
3.16D-01
,458.37D0
3.5156D-01
,503.37D0
8.5D-02
,543.37D0
9.4D-02
,628.37D0
NUMBER OF SOLUBLE COMPONENTS IN OIL
4
SOLUBLE COMPONENTS PARAMETERS (FRAC MASS, TEB, SOL, KDISS, KVOL)
1.D-02
,497.05D0, 0.018D0
, 1.25D-05 ,5.0D-05
3.2D-02 ,551.52D0, 0.00176D0 , 5.63D-06 ,1.51D-05
1.D-04
,674.68D0, 2.0D-04
, 2.D-06
,4.085D-07
2.D-05
,728.15D0, 1.33D-06 , 1.33D-06 ,1.20D-07
OIL DENSITY
830.D0
OIL VISCOSITY
4.2D-06
OIL SPILL VOLUME
2.02D-05
WATER TEMPERATURE
292.05D0
SPREADING MODEL (1=FAY’S MODEL, 2=MIGR’HYCAR MODEL, 3=CONSTANT AREA)
2
If in the oil spill steering file, the SPREADING MODEL is set to 3, two lines must be added to the
previous example:
CONSTANT AREA VALUE CHOSEN BY THE USER FOR EACH OIL PARTICLE
1 (example if the user wants area particle equal to 1 m2)
9.4.
THE OIL_FLOT SUBROUTINE
After inserting the OIL_FLOT subroutine in the FORTRAN file, it must be modified it in order to
indicate the release time step, together with the coordinates of the release point. If the release
point coordinates are outside the domain, the run is interrupted and an error message is displayed.
In addition, if a particle leaves the domain during the simulation, it is of course no longer monitored
but its previous track remains in the results file for consultation.
An example of modifications in the OIL_FLOT subroutine is given.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 81
The release time step in the first condition statement and the coordinates of the release point must
be changed:
...
IF(LT.EQ.10000)THEN
NUM_GLO=0
NUM_MAX=0
NUM_LOC=0
COORD_X=0.D0
COORD_Y=0.D0
NUM_MAX=INT(SQRT(REAL(NFLOT_MAX)))
DO K=1,NUM_MAX
DO J=1,NUM_MAX
COORD_X=336000.D0+REAL(J)
COORD_Y=371000.D0+REAL(K)
NUM_GLO=NUM_GLO+1
NFLOT_OIL=0
CALL ADD_PARTICLE(COORD_X,COORD_Y,0.D0,NUM_GLO,NFLOT_OIL,
&
1,XFLOT,YFLOT,YFLOT,TAGFLO,
&
SHPFLO,SHPFLO,ELTFLO,ELTFLO,MESH,1,
&
0.D0,0.D0,0.D0,0.D0,0,0)
...
END DO
END DO
END IF
9.5.
OUTPUT FILES
During an oil spill computation, the TELEMAC-3D software produces at least two output files:
• The 3D RESULT FILE,
• The output DROGUES FILE.
9.5.1.
THE 3D RESULT FILE
This is the file in which TELEMAC-3D stores information during the computation. It is normally in
SERAFIN format. First of all, it contains information on the mesh geometry, then the names of the
stored variables. It then contains the time for each time step and the values of the different
variables for all mesh points. For complementary information on the 3D RESULT FILE, the reader
may refer to 3.11.
9.5.2.
THE OUTPUT DROGUES FILE
This is an ASCII file created by TELEMAC-3D during the computation. It stores drogue positions in
TECPLOT format. To visualize the drogue positions with TECPLOT software, the user must:
• Use the File>Load Data File(s) command to load the 3D RESULT FILE
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 82
• Use the File>Load Data File(s) command to load the TECPLOT drogue file
See section 8.2 for more information as it is the same file as for drogues.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 83
10. OTHERS CONFIGURATIONS
10.1. MODIFICATION OF BOTTOM TOPOGRAPHY (CORFON)
Bottom topography may be introduced at various levels, as stated in section 3.16.
TELEMAC-3D offers the possibility of modifying the bottom topography at the beginning of a
computation using the T3D_CORFON subroutine. This is called up once at the beginning of the
computation and enables the value of variable ZF to be modified at each point of the mesh. To do
this, a number of variables such as the point coordinates, the element surface value, connectivity
table, etc, are made available to the user.
By default, the T3D_CORFON subroutine carries out a number of bottom smoothings equal to
LISFON, i.e. equal to the number specified by the keyword NUMBER OF BOTTOM SMOOTHINGS
for which the default value is 0 (no smoothing).
The T3D_CORFON subroutine is not called up if a computation is continued. This avoids having to
carry out several bottom smoothings or modifications to the bottom topography during the
computation.
10.2. MODIFYING COORDINATES (CORRXY)
TELEMAC-3D also offers the possibility of modifying the mesh point coordinates at the beginning
of a computation. This means, for example, that it is possible to change the scale (from that of a
reduced-scale model to that of the real object), rotate or translate the object.
The modification is done in the CORRXY subroutine (BIEF library), which is called up at the
beginning of the computation. This subroutine is empty by default and gives an example of
programming a change of scale and origin, within commented statements.
It is also possible to specify the coordinates of the origin point of the mesh. This is done using the
keyword ORIGIN COORDINATES which specify 2 integers. These 2 integers will be transmitted to
the results file in the SERAFIN format, for a use by post-processors for superimposition of results
with digital maps (coordinates in meshes may be reduced to avoid dealing with large real
numbers). These 2 integers may also be used in subroutines under the names I_ORIG and
J_ORIG. Otherwise they do not have a use yet.
10.3. SPHERICAL COORDINATES (LATITU)
If a simulation is performed over a large domain, TELEMAC-3D offers the possibility of running the
computation with spherical coordinates.
This option is activated when the keyword SPHERICAL COORDINATES is set to YES (default
value is NO). In this case, TELEMAC-3D calls a subroutine named LATITU through the subroutine
INBIEF at the beginning of the computation. This calculates the cosinus and sinus of the latitude
of each point. To do this, it uses the Cartesian coordinates of each point provided in the geometry
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 84
file, and the latitude of the point of origin of the mesh provided by the user in the steering file with
the keyword LATITUDE OF ORIGIN POINT.
The spatial projection type used for the mesh is then specified with the keyword SPATIAL
PROJECTION TYPE. That can take the following values:
•
1: Lambert Cartesian not geo referenced (default value – cannot be used in spherical)
•
2: Mercator;
•
3: Latitude/longitude (in degrees).
In this case of option 3, the coordinates of the mesh nodes should be expressed with latitude and
longitude in degrees. TELEMAC-3D then converts with the information with the help of the
Mercator’s projection.
The LATITU subroutine (BIEF library) may be modified by the user to introduce any other latitudedependent computation.
10.4. ADDING NEW VARIABLES
A standard feature of TELEMAC-3D is the storage of certain computed variables. In certain cases,
the user may wish to compute other variables and store them in the results file (the number of
variables is currently limited to four).
Since TELEMAC-3D uses 2D and 3D variables, the treatments linked to these variables may differ
and call three subroutines:
•
•
•
NOMVAR_2D_IN_3D: to manage 2D variables names,
NOMVAR_TELEMAC3D: to manage 3D variables names,
PRERES_TELEMAC3D: to compute new variables (2D and 3D).
TELEMAC-3D has a numbering system in which, for example, the array containing the Froude
number has the number 7. The new variables created by the user may have the numbers 25, 26,
27 and 28 (for 3D variables) and 27, 28, 29 and 30 (for 2D variables).
In the same way, each variable is identified by a letter in the keywords VARIABLES FOR 2D
GRAPHIC PRINTOUTS and VARIABLES FOR 3D GRAPHIC PRINTOUTS. The new variables are
identified by the strings PRIVE1, PRIVE2, PRIVE3 and PRIVE4 for 2D variables and P1, P2,
P3 and P4 for 3D variables.
At the end of the NOMVAR_XX subroutine (XX = 2D or 3D), it is possible to change the
abbreviations (mnemonics) used for the keywords VARIABLES FOR GRAPHIC 2D PRINTOUTS
and VARIABLES FOR GRAPHIC 3D. Sequences of 8 letters may be used. Consequently, the
variables must be separated by spaces, commas or semi-colons in the keywords, e.g.:
VARIABLES FOR GRAPHIC PRINTOUTS : 'U, V, H, B'
In the software data structure, these four variables correspond to the tables
PRIVE%ADR(1)%P%R(X),
PRIVE%ADR(2)%P%R(X),
PRIVE%ADR(3)%P%R(X)
and
PRIVE%ADR(4)%P%R(X) (in which X is the number of nodes in the mesh). These may be used in
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 85
several places in the programming, like all TELEMAC variables. For example, they may be used in
the subroutines CORRXY, CORSTR, BORD3D, etc. If a PRIVE table is used to program a case, it is
essential to check the value of the keyword NUMBER OF PRIVATE ARRAYS. This value fixes the
number of tables used (0, 1, 2, 3 or 4) and then determines the amount of memory space required.
The user can also access the tables via the aliases PRIVE1, PRIVE2, PRIVE3 and PRIVE4.
An example of programming using the second PRIVE table is given below. It is initialised with the
value 10.
DO I=1,NPOIN2
PRIVE%ADR(2)%P%R(I) = 10.D0
ENDDO
New variables are programmed in two stages:
•
•
Firstly, it is necessary to define the name of these new variables by filling in the
NOMVAR_TELEMAC3D (or NOMVAR_2D_IN_3D) subroutine. This consists of two equivalent
structures, one for English and the other for French. Each structure defines the name of
the variables in the results file that is to be generated and then the name of the variables to
be read from the previous computation if this is a restart. This subroutine may also be
modified when, for example, a file generated with the English version of TELEMAC-3D is to
be continued with the French version. In this case, the TEXTPR table of the French part of
the subroutine must contain the English names of the variables.
Secondly, it is necessary to modify the PRERES_TELEMAC3D subroutine in order to
introduce the computation of the new variable(s). The variables LEO, SORG2D and
SORG3D are also used to determine whether the variable is to be printed in the printout file
or in the results file at the time step in question.
User arrays can be handled to store extra variables in 2D with the help of two keywords to define
the number and the name of the extra variables in the 2D private arrays: NUMBER OF 2D
PRIVATE ARRAYS (up to 4, default value = 0) and NAMES OF 2D PRIVATE VARIABLES. It is the
names of the user arrays PRIVE%ADR(1)%P, PRIVE%ADR(2)%P… up to 4, that will be seen in
the results files. The great advantage is that these variables will be read if present in the
GEOMETRY FILE.
10.5. ARRAY MODIFICATION OR INITIALIZATION
When programming TELEMAC-3D subroutines, it is sometimes necessary to initialize a table or
memory space to a particular value. To do that, the BIEF library furnishes a subroutine called
FILPOL that lets the user modify or initialize tables in particular mesh areas.
A call of the type CALL FILPOL (F, C, XSOM, YSOM, NSOM, MESH) fills table F with the C
value in the convex polygon defined by NSOM nodes (coordinates XSOM, YSOM). The variable
MESH is needed for the FILPOL subroutine but have no meaning for the user.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 86
10.6. VALIDATING A COMPUTATION (VALIDA)
The structure of the TELEMAC-3D software offers an entry point for validating a computation, in
the form of a subroutine named VALIDA which has to be filled by the user in accordance with each
particular case. Validation may be carried out either with respect to a reference file (which is
therefore a file of results from the same computation that is taken as reference, the name of which
is supplied by the keyword REFERENCE FILE), or with respect to an analytical solution that must
then be programmed entirely by the user.
When using a reference file, the keyword REFERENCE FILE FORMAT specifies the format of this
binary file (SERAFIN by default).
The VALIDA subroutine is called at each time step when the keyword VALIDATION has the value
YES, enabling a comparison to be made with the validation solution at each time step. By default,
the VALIDA subroutine only does a comparison with the last time step. The results of this
comparison are given in the output listing.
10.7. COUPLING
The principle of coupling two (or in theory more) simulation modules involves running the two
calculations simultaneously and exchanging the various results at each time step. For example, the
following principle is used to couple a hydrodynamic module and a sediment transport module:
•
The two codes perform the calculation at the initial instant with the same information (in
particular the mesh and bottom topography),
•
The hydrodynamic code runs a time step and calculates the water depth and velocity
components. It provides this information to the sediment transport code,
•
The sediment transport code uses this information to run the solid transport calculation
over a time step and thus calculates a change in the bottom,
•
The new bottom value is then taken into account by the hydrodynamic module at the next
time step, and so on.
Two modules can be coupled in the current version of the code: the sedimentary transport module
SISYPHE and the sea state computational module TOMAWAC. The time step used for the two
calculations is not necessarily the same and is managed automatically by the coupling algorithms
and the keyword COUPLING PERIOD FOR SISYPHE and COUPLING PERIOD FOR TOMAWAC
with default values 1 (coupling at every iteration).
This functionality requires two keywords. The keyword COUPLING WITH indicates which
simulation code is to be coupled with TELEMAC-3D. The values of this keyword can be:
•
•
•
RELEASE 7.1
COUPLING WITH= ‘SISYPHE’ for coupling with the SISYPHE module,
COUPLING WITH= ‘TOMAWAC’ for coupling with the TOMAWAC module,
COUPLING WITH= ‘SISYPHE, TOMAWAC’ for coupling with both.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 87
Depending on the module(s) used, the keywords SISYPHE STEERING FILE and TOMAWAC
STEERING FILE indicate the names of the steering files of the coupled modules.
The keyword COUPLING WITH is also used if the computation has to generate the appropriate
files necessary to run a water quality simulation with DELWAQ. In that case, it is necessary to
specify COUPLING WITH= ‘DELWAQ’. Please refer to Appendix N° 4 for all informations
concerning communications with DELWAQ.
10.8. CHECKING THE MESH (CHECKMESH)
The subroutine CHECKMESH of the BIEF library is available to look for errors in the mesh, e.g.
superimposed points… The keyword CHECKING THE MESH (default value = NO) should be
activated to YES to call this subroutine.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 88
11. PARALLELISM
For simulations requiring a high computational power, it can be advisable to run the computations
in multi-processor machines, or in clusters of workstations. TELEMAC-3D is available in a parallel
version in order to take advantage of that kind of computational architecture.
The TELEMAC-3D parallel version uses the MPI library which has to be installed beforehand to be
implementable. The interface between TELEMAC-3D and that MPI library is achieved through the
parallel library which is common to all the TELEMAC system modules.
Lots of pieces of information concerning the implementation of the parallel version can be found in
the system’s installation literature.
The user shall initially specify the number of processors used by means of the keyword PARALLEL
PROCESSORS. That integer type keyword can take the following values:
•
0: Use of the conventional TELEMAC-3D version (default),
•
1: Use of the parallel TELEMAC-3D version on one processor,
•
2 ...: Use of the parallel TELEMAC-3D version using the specified number of processors.
Different partitioning tools may be used if installed. They can be selected by the keyword
PARTITIONING TOOL (default value is METIS). The possible choices are:
RELEASE 7.1
•
METIS,
•
SCOTCH,
•
PARMETIS,
•
PTSCOTCH…
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 89
APPENDIX N° 1. LAUNCHING THE COMPUTATION
A computation is launched through the telemac3d.py command. That command activates the
execution of a script which is common to all the computation modules in the TELEMAC system.
Depending on the platform, some options may be unavailable.
The syntaxes in that command are as follows:
telemac3d.py [cas]
[--options]
• cas: name of the steering file,
• --ncsize=NCSIZE: specifies the number of processors forced in parallel mode, default = the
number defined in the steering file,
• -c CONFIGNAME or –configname=CONFIGNAME: specifies the configuration name,
default is randomly found in the configuration file,
• -f CONFIGFILE, --configfile=CONFIGFILE: specifies the configuration file, default =
systel.cfg,
• -s, --sortiefile: specifies whether there is a sortie file, default is no,
• - t or --tmpdirectory: the temporary work directory is not destroyed on completion of
computation.
By default, the procedure runs the computation in an interactive mode and displays the control
listing on the monitor.
The operations performed by that script are as follows:
•
•
•
•
•
•
•
•
Creation of a temporary directory (name_cas_YYYY-MM-DD_HHhMMminSSs),
Duplication of the dictionary and the input files into that directory,
Execution of the DAMOCLES software in order to determine the work file names,
Creation of the computation launching script,
Allocation of the files,
Compilation of the FORTRAN file and link editing (as required),
Launching of the computation,
Retrieval of the results files and destruction of the temporary directory.
The procedure takes place with slight differences according to the selected options.
The detailed description of that procedure can be obtained through the help command:
telemac3d.py --help
RELEASE 7.1
or
telemac3d.py -h
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 90
APPENDIX N° 2. LIST OF USER SUBROUTINES
Even though all subroutines can be modified by the user, some subroutines have been specifically
designed to define complex simulation parameters. They are listed below:
BORD3D
CALCOT
CONDIN
CONDIS
CORFON
CORRXY
CORSTR
DECLARATIONS_TELEMAC3D
DRIUTI
DRSURR
FLOT3D
LIMI3D
NOMVAR_TELEMAC3D
PRERES_TELEMAC3D
Q3
SCOPE
SL3
SOURCE
SOURCE_TRAC
TR3
TRISOU
UTIMP
VEL_PROF_Z
VISCLM
VISCOS
VIT3
RELEASE 7.1
Management of the boundary conditions
Preparation of the array of mesh elevations between the bottom and
the free surface
Management of the initial conditions
Initialization of the arrays of the physical sedimentological quantities
Modification of bottoms
Modification of the mesh node co-ordinates
Correction of the bottom friction coefficient when it is time-variable
Statement of the TELEMAC-3D structures
User damping function
Computation of density (equation of state)
Initial conditions of the drogues
Management of the boundary conditions
Management of the names of variables for the graphic printouts
Computation of the output variables (free surface, flow rate…)
Management of the flow rates in a boundary condition
Creation of 1D sections
Management of the open surface in a boundary condition
User source term in the hydrodynamic equation
User source term in the equations of tracers
Management of the tracers in a boundary condition
Source terms for the velocity components
Additional variable writing
Definition of the vertical velocity profile
Computation of the viscosity in the mixing length models
Computation and initialization of the constant viscosity
Management of the velocities in a boundary condition
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 91
APPENDIX N° 3. DESCRIPTION OF THE SERAFIN FORMAT
It is a binary file.
The list of records is as follows:
•
1 record containing the title of the study (72 characters) and 8 characters indicating the
format (SERAFIN or SERAFIND),
•
1 record containing the couple of integer values NBV(1) and NBV(2) (number of linear
and quadratic discretization variables, NBV(2) being 0),
•
NBV(1) records containing both name and unit of each variable (over 32 characters –
normally 16 for the variable’s name and 16 for the unit),
•
A record containing the IPARAM array consisting of 10 integers,
o
o
o
o
o
o
o
RELEASE 7.1
If IPARAM(3) ≠ 0: the value corresponds to the x-coordinate of the origin of the
mesh,
If IPARAM(4) ≠ 0: the value corresponds to the y-coordinate of the origin of the
mesh,
If IPARAM(7) ≠ 0: the value corresponds to the number of planes on the vertical
(3D computation),
If IPARAM(8) ≠ 0: the value corresponds to the number of boundary points (in
parallel),
If IPARAM(9) ≠ 0: the value corresponds to the number of interface points (in
parallel),
If IPARAM(8) or IPARAM(9) ≠ 0: the array IPOBO below is replaced by the array
KNOLG (total initial number of points). All the other numbers are local to the subdomain, including IKLE,
If IPARAM(10) = 1 : the file contains the record of both date and time of the
computation start (6 integers) which take the values of the keywords ORIGINAL
DATE OF TIME and ORIGINAL HOUR OF TIME from the steering file.
•
A record containing the integers NELEM3, NPOIN3, NDP, 1 (number of elements, number of
points, number of points per element and the value 1),
•
A record containing the IKLE3 integer array ((NDP, NELEM3-dimensioned array), the
connectivity table,
•
A record containing the IPOBO integer array (NPOIN3-dimensioned array). The value of an
element is 0 for an inner point and yields the edge point numbers for the others),
•
A record containing the X real array (NPOIN3-dimensioned array of the node abscissae),
•
A record containing the Y real array (NPOIN3-dimensioned array of node ordinates),
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 92
Afterwards, the following record can be found for each time step:
RELEASE 7.1
•
A record containing the AT time (real),
•
NBV(1)+NBV(2) records containing the result arrays for each variable at the AT time.
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 93
APPENDIX N° 4. GENERATING OUTPUT FILES FOR DELWAQ
The TELEMAC-3D software is able to generate the appropriate files necessary to run a DELWAQ
water quality simulation. This generation is managed only through the following keyword:
BOTTOM SURFACES DELWAQ FILE
DELWAQ PRINTOUT PERIOD
DELWAQ STEERING FILE
DIFFUSION FOR DELWAQ
DIFFUSIVITY DELWAQ FILE
EXCHANGE AREAS DELWAQ FILE
EXCHANGES BETWEEN NODES DELWAQ FILE
NODES DISTANCES DELWAQ FILE
SALINITY DELWAQ FILE
SALINITY FOR DELWAQ
TEMPERATURE DELWAQ FILE
TEMPERATURE FOR DELWAQ
VELOCITY DELWAQ FILE
VELOCITY FOR DELWAQ
VERTICAL FLUXES DELWAQ FILE
VOLUMES DELWAQ FILE
More information about these keywords can be found in the TELEMAC-3D reference manual. For
more information, please refer to the DELWAQ user documentation.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 94
APPENDIX N° 5. POSTEL-3D
POSTEL-3D allows extracting 2D vertical or horizontal cross sections from the 3D result file. The
resulting cross section files are readable with the different post-processing tools of the TELEMAC
system.
The files generated by POSTEL-3D are in SERAFIN single or double precision format or MED.
A computation is launched through the command postel3d.py [cas] (cas: the name of the
steering file).
POSTEL-3D is implemented like all the codes in the TELEMAC treatment chain.
Its execution is structured around the STEERING FILE, which is a priori the single file which the
user will have to consult and amend. It gathers the names of all the files which define the
computation to be carried out.
The input files are:
•
The 3D RESULT FILE as provided by TELEMAC-3D, in the TEL3D format,
•
The 3D RESULT FILE FORMAT to specify the file format; by default it is SERAFIN single
precision,
•
The ASCII FORTRAN FILE containing the amended subroutines. That file is optional.
The output files are:
•
The HORIZONTAL CROSS SECTION FILES, in the SERAFIN format (can be run under
any post-processing tool that can read SERAFIN files),
•
The VERTICAL CROSS SECTION FILES, in the SERAFIN format (can be run under any
post-processing tool that can read SERAFIN files).
THE HORIZONTAL CROSS SECTIONS
The horizontal cross sections are stored into binaries files in the SERAFIN format, on the basis of
one file per cross section. The name of each such file consists in a common radical, which is given
by the keyword HORIZONTAL CROSS SECTION FILE, followed by an extension specifying the
cross section number.
A horizontal cross section is not necessarily horizontal, i.e. parallel to the level Z = 0. It may also be
in the form of a 2D level of the TELEMAC-3D computation, with a possible vertical offset.
The mesh on which each of these cross sections is based is the 2D mesh which was used for
conducting the TELEMAC-3D computation.
The horizontal cross sections are defined by means of the keywords:
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 95
•
NUMBER OF HORIZONTAL CROSS SECTIONS specifies the number of cross sections to
be made,
•
REFERENCE LEVEL FOR EACH HORIZONTAL CROSS SECTION specifies the
TELEMAC-3D vertical level from which the cross section shape is defined,
•
ELEVATION FROM REFERENCE LEVEL specifies the vertical distance from the
reference level at which the cross section has to be made.
The keyword REFERENCE LEVEL assumes a value ranging from 0 and NPLAN, NPLAN is the
number of levels selected for the TELEMAC-3D computation which made it possible to prepare the
3D result file read. If the value is between 1 and NPLAN, the reference level is the relevant mesh
level which is liable to move in time; if the value is zero, then the reference level is the level z = 0.
The cross section level is then inferred from the reference level through a mere vertical translation
the amount of which is set up by the keyword ELEVATION FROM REFERENCE LEVEL.
In doing so, there is a slight difficulty because of the choice the user has to define a cross section is
through two parameters: REFERENCE LEVEL FOR EACH HORIZONTAL CROSS SECTION and
ELEVATION FROM REFERENCE LEVEL. That only defines the cross section level, and in such a
case, two problematic scenarios can occur:
•
Either that level locally occurs above the free surface,
•
Or it occurs locally under the bottom.
We have therefore adopted another approach which consists in performing a linear extrapolation at
these points from the closest 2 values, always occurring vertically above that point. For those
points located above the surface, that extrapolation is then done from the values computed at the
levels NPLAN-1 and NPLAN whereas the points located below the bottom, the extrapolation is
computed from the values at levels 1 and 2.
Though the result of that extrapolation generally gives a value which is realistic, since it is not far
from the values found in the domain, the user shall be made aware of the fact that such a result
has no physical meaning and that it does not occur on the result planes.
In order to indicate the locations of these nodes, an INDICATEUR_DOM variable is provided for the
user in the cross section file. When that variable is negative, that means the points are outside the
domain. By inserting a coloured (e.g. white) surface of the INDICATEUR_DOM variable with a ]-∞, 0]
threshold, the user can then mask the out-of-domain areas that have meaningless values.
As regards the VITESSE_U and VITESSE_V variables, we have preferred to preset these
variables to 0 at the points located outside the domain. That treatment is more suitable for a
vectorial plot.
THE VERTICAL CROSS SECTIONS
A vertical cross section can be defined in a 2D mesh as a sequence of linked points making up a
pecked line consisting of segments. That pecked line is vertically extended from the surface down
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 96
to the bottom. The minimum point number is 2 (1 segment) and the maximum number is 9 (8
segments).
The vertical cross sections are defined by means of the keywords:
•
NUMBER OF VERTICAL CROSS SECTIONS specifying the number of cross sections to
be made. That number cannot be in excess of 9. If over 9 vertical cross sections are
desired, then several software executions should be planned,
•
NUMBER OF NODES FOR VERTICAL CROSS SECTION DISCRETIZATION setting the
number of interpolation points in the horizontal direction of the cross section. The points
are evenly spaced. That number should be in excess of 2,
•
ABSCISSAE OF THE VERTICES OF CROSS SECTION X specifying the abscissa of each
point along the pecked line making up the vertical cross section. That number should be
higher than or equal to 2 and lower than or equal to 9. X specifies the cross section
number ranging from 1 to 9. WARNING: if X is higher than the NUMBER OF VERTICAL
CROSS SECTIONS, the information in that keyword will merely not be treated,
•
ORDINATES OF THE VERTICES OF CROSS SECTION X specifying the ordinate of each
point along the pecked line making up the vertical cross section. That number should be
higher than or equal to 2 and lower than or equal to 9. X specifies the cross section
number ranging from 1 to 9. WARNING: if X is higher than the NUMBER OF VERTICAL
CROSS SECTIONS, the information in that keyword will merely not be treated.
The vertical cross sections are stored into binary files in the SERAFIN format, on the basis of one
file per cross section and per recorded time step. The name of each such file consists of a common
radical, as given by the keyword VERTICAL CROSS SECTION FILE, followed by an extension
specifying the number of the cross section, then by an extension specifying the number of the
recorded time step. That increased number of files is necessary because the meshes are distorted
in time, and so are these cross sections, due to the free surface motions.
The horizontal velocity components are given in a cross section-related co-ordinate system
provided for directly drawing the projection of the velocity vector onto the cross section level. The
components in the new co-ordinate system are known as:
•
VITESSE_UT: tangential component,
•
VITESSE_UN: normal component.
The TELEMAC-3D computational domains often include a vertical scale which is much lower than
the horizontal scale. The vertical scale can be distorted by a multiplicative factor so that RUBENS
or any other post-processing tool, can output a clearer display. That can be done using the
keyword DISTORTION BETWEEN VERTICAL AND HORIZONTAL.
WARNING: If one uses a distortion factor, the vertical velocities will themselves be multiplied by
that factor. That is necessary to properly represent the velocity vector directions.
Lastly, note that the starting point is located to the left of the cross section and the ending point to,
the right of it i.e. the length is defined from the bottom to the free surface.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
RELEASE 7.1
Page 97
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 98
APPENDIX N° 6. REFERENCES
[1]
HERVOUET, J.-M. " Hydrodynamics of Free Surface Flows. Modelling with the finite element
method". Wiley, 2007.
[2]
QUETIN, B. Modèles mathématiques de calcul des écoulements induits par le vent", 17e
congrès de l'A.I.R.H, 15-19 août 1977, Baden-Baden.
[3]
TSANIS, I. "Simulation of wind-induced water currents", Journal of hydraulic Engineering,
Vol. 115, n°8, pp1113-1134, 1989.
[4]
SMAGORINSKY, J. "General simulation experiments with the primitive equations". Mon.
Wea. Rev. 91, 99-164, 1963.
[5]
HERVOUET, J.-M. “Guide to programming in the Telemac system version 6.0”. EDF report
H-P74-2009-00801-EN.
[6]
PHAM, C.-T., BOURBAN, S., DURAND, N., TURNBULL, M. "Méthodologie pour la
simulation de la marée avec la version 6.2 de TELEMAC-2D et TELEMAC-3D". EDF report H-P742012-02534-FR.
[7]
JOLY, A., GOEURY, C., HERVOUET, J.-M. "Adding a particle transport module to
TELEMAC-2D with applications to algae blooms and oil spills". EDF report H-P74-2013-02317-EN.
RELEASE 7.1
JULY 2016
TELEMAC modelling system
TELEMAC-3D / Operating manual
Page 99
BINARY DATABASE 2 FOR TIDE .......................................... 46
&
&ETA................................................................................ 20
&FIN ................................................................................ 20
&IND................................................................................ 20
&LIS ................................................................................. 20
&STO ............................................................................... 20
BINARY RESULTS FILE.................................................. 26, 76
BLUE KENUE ...................................................... 8, 15, 18, 21
BORD3D ............................................................... 44, 57, 73
BOTTOM TOPOGRAPHY FILE ............................................. 27
BOUNDARY CONDITION ON THE BOTTOM......................... 40
BOUNDARY CONDITIONS FILE ................................ 21, 41, 73
BTBOR ............................................................................. 37
2
BYPASS VOID VOLUMES.................................................... 69
C
2D CONTINUATION..................................................... 22, 33
2D RESULT FILE................................................................. 25
2D RESULT FILE FORMAT ............................................ 18, 25
CALCOT ...................................................................... 28, 90
CGSTAB ............................................................................ 66
3
CHECKING THE MESH ....................................................... 87
COEFFICIENT FOR HORIZONTAL DIFFUSION OF TRACERS.... 75
3D RESULT FILE........................................................... 24, 81
3D RESULT FILE FORMAT ............................................ 18, 24
A
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES50,
53
COEFFICIENT FOR VERTICAL DIFFUSION OF TRACERS ......... 75
COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES 50, 53
ABSCISSAE OF SOURCES.................................................... 55
COEFFICIENT OF WIND INFLUENCE.................................... 56
ACCURACY FOR DIFFUSION OF K-EPSILON ......................... 67
COEFFICIENT TO CALIBRATE SEA LEVEL.............................. 48
ACCURACY FOR DIFFUSION OF TRACERS............................ 67
COEFFICIENT TO CALIBRATE THE ATMOSPHERE-WATER
ACCURACY FOR DIFFUSION OF VELOCITIES ........................ 67
EXCHANGE MODEL....................................................... 58
ACCURACY FOR PPE .................................................... 65, 67
COEFFICIENT TO CALIBRATE TIDAL RANGE ........................ 48
ACCURACY FOR PROPAGATION ......................................... 67
COEFFICIENT TO CALIBRATE TIDAL VELOCITIES .................. 48
ACCURACY FOR VERTICAL VELOCITY............................ 65, 67
COMPUTATION CONTINUED ....................................... 22, 32
ADD_PARTICLE ....................................................... 76, 81
CONDIM ........................................ 26, 28, 29, 30, 31, 32, 72
Adding new variables ....................................................... 84
CONSTANT DEPTH ............................................................ 31
advection schemes ........................................................... 63
CONSTANT ELEVATION ..................................................... 31
ADVECTION STEP.............................................................. 61
CONSTANT VISCOSITY....................................................... 50
AIR PRESSURE................................................................... 57
CORIOLIS .......................................................................... 59
AIR TEMPERATURE ........................................................... 57
CORIOLIS COEFFICENT ...................................................... 60
ARTEMIS .......................................................................... 59
CORRXY .......................................................................... 83
ASCII DATABASE FOR TIDE ................................................ 46
CORSTR ............................................................................ 85
ATBOR ............................................................................. 37
COUPLING PERIOD............................................................ 86
ATMOSPHERE-WATER EXCHANGE MODEL ........................ 58
COUPLING WITH............................................................... 86
AUBOR ............................................................................. 37
Crout ............................................................................... 68
AVERAGE WATER DENSITY................................................ 74
B
BETA EXPANSION COEFFICIENT FOR TRACERS.................... 74
D
DAMOCLES................................................................. 19, 89
DAMPING FUNCTION........................................................ 51
BINARY DATA FILE 1.................................................... 26, 59
DEBUGGER ....................................................................... 18
BINARY DATA FILE 1 FORMAT ........................................... 18
DECLARATION_TELEMAC3D ........................................... 9
BINARY DATA FILE 2.......................................................... 26
DECLARATIONS_TELEMAC3D ...................................... 26
BINARY DATABASE 1 FOR TIDE .......................................... 46
DELWAQ .................................................................... 87, 93
RELEASE 7.1
JULY 2016
TELEMAC modelling system
100
TELEMAC-3D / Operating manual
DENSITY LAW ................................................................... 74
Page
HYDROSTATIC INCONSISTENCY FILTER............................... 69
DIFFUSION STEP ............................................................... 61
DRIUTI .......................................................................... 51
DROGUES FILE ................................................. 76, 77, 79, 81
DURATION ....................................................................... 28
DYNAMIC BOUNDARY CONDITION .................................... 45
DYNAMIC PRESSURE IN WAVE EQUATION ......................... 62
E
I
I_ORIG ............................................................................. 83
IDP ................................................................................. 16
IELEM2 .......................................................................... 16
IELM2 ............................................................................. 16
IKLE ............................................................................... 21
IKLE2 ............................................................................. 16
ELEMENT ......................................................................... 30
IKLE3 ............................................................................. 16
ELEVATIONS OF SOURCES ................................................. 55
IMPLICITATION FOR DEPTH............................................... 65
EXCHANGE_WITH_ATMOSPHERE...................................... 58
IMPLICITATION FOR DIFFUSION ........................................ 65
IMPLICITATION FOR VELOCITIES........................................ 65
F
FES................................................................................... 46
FILE FOR 2D CONTINUATION.................................. 18, 22, 33
FILE FOR 2D CONTINUATION FORMAT ........................ 18, 33
FILE FOR SCOPE ................................................................ 26
FILPOL.............................................................................. 85
FLOT3D .................................................................... 76, 90
FORMATTED DATA FILE 1.................................................. 26
FORMATTED DATA FILE 2.................................................. 26
FORMATTED RESULTS FILE................................................ 26
INFORMATION ABOUT MASS-BALANCE FOR EACH LISTING
PRINTOUT .................................................................... 25
INITIAL CONDITIONS................................................... 31, 32
INITIAL DEPTH .................................................................. 31
INITIAL ELEVATION ........................................................... 31
INITIAL GUESS FOR DEPTH ................................................ 70
INITIAL TIME SET TO ZERO ................................................ 33
INITIAL VALUES OF TRACERS ............................................. 72
J
Fortran 90 .......................................................................... 9
J_ORIG ............................................................................. 83
FORTRAN FILE ............................................................ 21, 26
JANET......................................................................... 15, 21
FREE SURFACE COMPATIBILITY GRADIENT......................... 62
JMJ .................................................................................. 46
FRICTION COEFFICIENT FOR LATERAL SOLID BOUNDARIES . 54
FRICTION COEFFICIENT FOR THE BOTTOM......................... 54
FUDAA-PREPRO............................................8, 18, 19, 21, 22
G
GEOGRAPHIC SYSTEM....................................................... 47
GEOMETRY FILE..................................................... 16, 21, 85
GEOMETRY FILE FORMAT ........................................... 18, 21
GLOBAL NUMBER OF THE POINT TO CALIBRATE HIGH WATER
.................................................................................... 47
GMRES ............................................................................. 66
GRAPHIC PRINTOUT PERIOD ....................................... 23, 24
GRAVITY ACCELERATION................................................... 60
Guide for programming in the Telemac system ................... 8
H
K
k-ε ............................................................................. 14, 52
KARMAN CONSTANT ........................................................ 49
L
LATITU .......................................................................... 83
LATITUDE OF ORIGIN POINT.............................................. 84
Launching the computation .............................................. 89
LAW OF BOTTOM FRICTION .............................................. 54
LAW OF FRICTION ON LATERAL BOUNDARIES .................... 54
LEO .................................................................................. 85
LIHBOR .................................................................... 36, 37
LIMI3D ......................................................... 35, 38, 39, 41
LINEARIZED PROPAGATION............................................... 62
LIQUID BOUNDARIES FILE ...................................... 22, 42, 73
HARMONIC CONSTANTS FILE ............................................ 47
liquid boundary ................................................................ 36
HBOR ............................................................................... 37
LISTING PRINTOUT PERIOD ............................................... 25
HORIZONTAL CROSS SECTIONS ......................................... 94
LITBOR ............................................................... 36, 37, 72
HORIZONTAL TURBULENCE MODEL............................. 49, 52
LIUBOR ............................................................... 36, 37, 40
RELEASE 7.1
JULY 2016
TELEMAC modelling system
101
TELEMAC-3D / Operating manual
Page
LIVBOR ....................................................36, 37, 38, 39, 40
NEA.................................................................................. 46
LOCAL NUMBER OF THE POINT TO CALIBRATE HIGH WATER
NELEM2 .................................................................... 16, 21
.................................................................................... 47
NELEM3 ................................................................. 9, 16, 91
LONGITUDE OF ORIGIN POINT .......................................... 59
NOMVAR_2D_IN_3D ........................................................ 84
NOMVAR_TELEMAC3D ............................................... 84, 85
M
MARDAT .......................................................................... 28
MARTIM........................................................................... 28
MASS-BALANCE ................................................................ 25
MASS-LUMPING FOR DEPTH ............................................. 70
MASS-LUMPING FOR DIFFUSION....................................... 70
MASS-LUMPING FOR VELOCITIES ...................................... 70
MASS-LUMPING FOR WEAK CHARACTERISTICS............ 64, 70
MATISSE.......................................................... 15, 21, 27, 38
MATRIX STORAGE............................................................. 70
MAXIMUM NUMBER OF BOUNDARIES .............................. 34
MAXIMUM NUMBER OF BOUNDARIES ON THE BED .......... 41
MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION
SCHEMES ..................................................................... 67
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF KEPSILON ....................................................................... 67
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF
TRACERS ...................................................................... 67
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF
VELOCITIES................................................................... 67
NON-HYDROSTATIC VERSION...................................... 49, 64
NPLAN .................................................................. 16, 29, 30
NPOIN2 ............................................................... 16, 21, 85
NPOIN3 ................................................................. 9, 16, 91
NPTFR2 ...................................................................... 9, 44
NSOM .............................................................................. 85
NUMBER OF 2D PRIVATE ARRAYS ............................... 25, 85
NUMBER OF BOTTOM SMOOTHINGS .......................... 27, 83
NUMBER OF DROGUES ............................................... 76, 78
NUMBER OF FIRST TIME STEP FOR GRAPHIC PRINTOUTS .. 23,
24
NUMBER OF FIRST TIME STEP FOR LISTING PRINTOUTS ..... 25
NUMBER OF GAUSS POINTS FOR WEAK CHARACTERISTICS 64
NUMBER OF HORIZONTAL LEVELS ..................................... 28
NUMBER OF PRIVATE ARRAYS..................................... 24, 85
NUMBER OF SUB ITERATIONS FOR NON LINEARITIES ......... 64
NUMBER OF TIME STEPS................................................... 28
NUMBER OF TRACERS....................................................... 72
O
MAXIMUM NUMBER OF ITERATIONS FOR PPE ............ 65, 67
OIL SPILL MODEL .............................................................. 78
MAXIMUM NUMBER OF ITERATIONS FOR PROPAGATION . 67
OIL SPILL STEERING FILE.............................................. 78, 79
MAXIMUM NUMBER OF ITERATIONS FOR VERTICAL
OIL_FLOT ................................................................ 78, 80
VELOCITY ............................................................... 65, 67
OPEN BOUNDARY CONDITIONS ON THE BED ..................... 40
MAXIMUM NUMBER OF ITERATIONS IS REACHED ............. 67
OPTION FOR CHARACTERISTICS......................................... 64
MAXIMUM NUMBER OF SOURCES .................................... 56
OPTION FOR LIQUID BOUNDARIES .................................... 45
MAXIMUM NUMBER OF TRACERS..................................... 72
OPTION FOR THE BOUNDARY CONDITIONS OF K-EPSILON . 53
MEAN DEPTH FOR LINEARIZATION .................................... 62
OPTION FOR THE DIFFUSION............................................. 62
MESH ......................................................................... 15, 85
OPTION FOR THE TREATMENT OF TIDAL FLATS.................. 69
MESH TRANSFORMATION...................................... 16, 28, 30
OPTION FOR TIDAL BOUNDARY CONDITIONS .............. 46, 47
METEO ............................................................................. 57
OPTION OF SOLVER FOR DIFFUSION OF K-EPSILON ............ 66
MINIMUM VALUE FOR DEPTH........................................... 69
OPTION OF SOLVER FOR DIFFUSION OF TRACERS .............. 66
MINOR CONSTITUENTS INFERENCE ................................... 46
OPTION OF SOLVER FOR DIFFUSION OF VELOCITIES........... 66
MIXING LENGTH ............................................................... 50
OPTION OF SOLVER FOR PPE ....................................... 65, 66
MIXING LENGTH MODEL................................................... 50
OPTION OF SOLVER FOR PROPAGATION ............................ 66
MODIFICATION OF BOTTOM TOPOGRAPHY....................... 83
ORDINATES OF SOURCES .................................................. 55
MODIFYING COORDINATES............................................... 83
ORIGIN COORDINATES ...................................................... 83
ORIGINAL DATE OF TIME ................................. 28, 47, 59, 91
N
NAMES OF 2D PRIVATE VARIABLES ................................... 85
NAMES OF TRACERS ......................................................... 72
NDP ...................................................................... 16, 21, 91
RELEASE 7.1
ORIGINAL HOUR OF TIME ................................ 28, 47, 59, 91
P
parallel ...................................................................... 88
JULY 2016
TELEMAC modelling system
102
TELEMAC-3D / Operating manual
Page
PARALLEL PROCESSORS .................................................... 88
SCHEME FOR ADVECTION OF TRACERS........................ 61, 75
PARALLELISM ................................................................... 88
SCHEME FOR ADVECTION OF VELOCITIES .................... 61, 64
PARTICULAR ..................................................................... 31
SCHEME FOR DIFFUSION OF K-EPSILON............................. 62
PARTITIONING TOOL......................................................... 88
SCHEME FOR DIFFUSION OF TRACERS ......................... 61, 75
POSTEL-3D ....................................................................... 94
SCHEME FOR DIFFUSION OF VELOCITIES............................ 61
PRANDTL NUMBER ........................................................... 49
SCOPE ............................................................................. 26
PRECONDITIONING FOR DIFFUSION OF K-EPSILON ............ 68
Secchi .............................................................................. 58
PRECONDITIONING FOR DIFFUSION OF TRACERS............... 68
SERAFIN ........................................................................... 18
PRECONDITIONING FOR DIFFUSION OF VELOCITIES ........... 68
SERAFIND......................................................................... 18
PRECONDITIONING FOR PPE ....................................... 65, 68
SINUSX............................................................................. 27
PRECONDITIONING FOR PROPAGATION ............................ 68
SISYPHE............................................................................ 86
PRECONDITIONING FOR VERTICAL VELOCITY ............... 65, 68
SISYPHE STEERING FILE ..................................................... 87
PRERES_TELEMAC3D .................................................. 84, 85
SL3 ................................................................................. 43
PRESCRIBED ELEVATIONS............................................ 40, 73
Smagorinsky..................................................................... 52
PRESCRIBED FLOWRATES............................................ 40, 73
SOLVER FOR DIFFUSION OF K-EPSILON .............................. 66
PRESCRIBED FLOWRATES ON THE BED .............................. 41
SOLVER FOR DIFFUSION OF TRACERS ................................ 66
PRESCRIBED TRACERS VALUES .......................................... 73
SOLVER FOR DIFFUSION OF VELOCITIES............................. 66
PRESCRIBED VELOCITIES ................................................... 40
SOLVER FOR PPE......................................................... 65, 66
PREVIMER ........................................................................ 46
SOLVER FOR PROPAGATION.............................................. 66
PREVIOUS COMPUTATION FILE ........................ 22, 31, 32, 33
SOLVER FOR VERTICAL VELOCITY ................................ 65, 66
PREVIOUS COMPUTATION FILE FORMAT ..................... 19, 33
SORG2D ........................................................................... 85
PRINTOUT PERIOD FOR DROGUES ............................... 76, 78
SORG3D ........................................................................... 85
PRISM .............................................................................. 30
SOURCE .......................................................................... 56
PRIVE ............................................................................... 84
SOURCE_TRAC .................................................... 58, 75, 90
PROPAGATION STEP ......................................................... 62
SOURCES FILE ............................................................. 23, 55
SPATIAL PROJECTION TYPE ............................................... 84
Q
Q3 ................................................................................... 43
SPHERICAL COORDINATES........................................... 60, 83
STAGE-DISCHARGE CURVES .............................................. 43
STAGE-DISCHARGE CURVES FILE ................................. 22, 43
R
STANDARD VALUES FOR TRACERS ..................................... 74
STBTEL ............................................................ 15, 21, 27, 36
RAIN OR EVAPORATION .................................................... 57
SUPG OPTION................................................................... 64
RAIN OR EVAPORATION IN MM PER DAY........................... 57
RECORD NUMBER FOR RESTART ................................. 32, 33
T
RECORD NUMBER IN WAVE FILE ....................................... 59
REFERENCE FILE.......................................................... 22, 86
REFERENCE FILE FORMAT ...................................... 19, 22, 86
References ....................................................................... 98
RESTART FILE.................................................................... 33
RESTART FILE FORMAT ............................................... 19, 33
RESTART MODE ................................................................ 33
Rouse profile.................................................................... 73
RUBENS...................................................................... 18, 21
T3D_CORFON ............................................................. 27, 83
T3D_DEBSCE .................................................................... 55
T3D_TRSCE................................................................. 55, 75
T3DBI1 ............................................................................. 26
T3DBI2 ............................................................................. 26
T3DFO1............................................................................ 26
T3DFO2............................................................................ 26
T3DRBI ............................................................................. 26
T3DRBO ........................................................................... 26
S
TAGFLO .................................................................... 76, 81
TBOR ......................................................................... 37, 73
SCHEME FOR ADVECTION OF DEPTH ................................. 61
TETRAHEDRON ........................................................... 15, 30
SCHEME FOR ADVECTION OF K-EPSILON ........................... 61
TEXTPR............................................................................. 85
RELEASE 7.1
JULY 2016
TELEMAC modelling system
103
TELEMAC-3D / Operating manual
Page
THE STEERING FILE ........................................................... 19
VBOR ......................................................................... 37, 44
THRESHOLD DEPTH FOR WIND.......................................... 57
VEL_PROF_Z ................................................................. 45
THRESHOLD FOR VISCOSITY CORRECTION ON TIDAL FLATS 69
VELOCITIES OF THE SOURCES ALONG X ............................. 56
TIDAL DATA BASE ............................................................. 46
VELOCITIES OF THE SOURCES ALONG Y ............................. 56
TIDAL FLATS ..................................................................... 68
VELOCITIES OF THE SOURCES ALONG Z.............................. 56
TIDAL MODEL FILE ............................................................ 46
VELOCITY PROFILES .......................................................... 44
TIDE GENERATING FORCE ................................................. 59
VELOCITY PROJECTED ON BOTTOM ................................... 70
TIME STEP ............................................................. 23, 25, 28
VELOCITY PROJECTED ON SOLID LATERAL BOUNDARIES .... 70
TITLE ................................................................................ 28
VELOCITY VERTICAL PROFILES ........................................... 45
TOMAWAC................................................................. 59, 86
VERTICAL CROSS SECTIONS............................................... 95
TOMAWAC STEERING FILE ................................................ 87
VERTICAL TURBULENCE MODEL .................................. 49, 52
TPXO................................................................................ 46
VERTICAL VELOCITY DERIVATIVES ..................................... 52
TPXO SATELLITE ALTIMETRY.............................................. 31
VISCLIP............................................................................. 69
TRACERS VERTICAL PROFILES ............................................ 73
VISCLM .......................................................................... 52
TRANSF_PLANE ..................................................... 28, 29, 30
VISCOS .......................................................................... 50
TREATMENT OF FLUXES AT THE BOUNDARIES ................... 73
VIT3 ............................................................................... 43
TREATMENT OF NEGATIVE DEPTHS................................... 69
VSCE ............................................................................... 56
TREATMENT ON TIDAL FLATS FOR K-EPSILON .................... 69
TREATMENT ON TIDAL FLATS FOR TRACERS ...................... 69
TREATMENT ON TIDAL FLATS FOR VELOCITIES................... 69
TRISOU............................................................................. 59
TURBULENCE REGIME FOR LATERAL SOLID BOUNDARIES... 54
TURBULENCE REGIME FOR THE BOTTOM .......................... 54
TYPE OF SOURCES............................................................. 56
U
UBOR .................................................................... 37, 44, 45
USCE ............................................................................... 56
USE BIEF ........................................................................ 9
W
WATER DISCHARGE OF SOURCES ...................................... 55
WAVE DRIVEN CURRENTS ................................................. 59
WIND ............................................................................... 56
WIND VELOCITY ALONG X ................................................. 56
WIND VELOCITY ALONG Y ................................................. 56
X
XFLOT ....................................................................... 76, 81
Y
USE DECLARATIONS_TELEMAC3D ............................... 9
YFLOT ....................................................................... 76, 81
V
Z
VALIDA............................................................................. 86
VALIDATING A COMPUTATION.......................................... 86
ZERO DEPTH..................................................................... 31
VALIDATION ..................................................................... 86
ZERO ELEVATION.............................................................. 31
VALUE OF THE TRACERS AT THE SOURCES ......................... 74
ZFLOT ............................................................................. 76
VALUES OF TRACERS IN THE RAIN ..................................... 57
ZONE NUMBER IN GEOGRAPHIC SYSTEM .......................... 48
VARIABLES FOR 2D GRAPHIC PRINTOUTS .................... 24, 84
ZPLANE ............................................................................ 30
VARIABLES FOR 3D GRAPHIC PRINTOUTS .......................... 23
ZSTAR ........................................................................ 29, 30
VARIABLES FOR LISTING PRINTOUTS ................................. 84
RELEASE 7.1
JULY 2016
Source Exif Data:
File Type : PDF File Type Extension : pdf MIME Type : application/pdf PDF Version : 1.6 Linearized : No Author : J97554 Create Date : 2016:07:04 16:16:02+02:00 Modify Date : 2018:01:15 17:55:53Z Subject : Has XFA : No XMP Toolkit : Adobe XMP Core 5.6-c015 84.159810, 2016/09/10-02:41:30 Producer : GPL Ghostscript 9.05 Keywords : Creator Tool : PDFCreator Version 1.6.2 Metadata Date : 2018:01:15 17:55:53Z Document ID : uuid:5a6b26c1-444d-11e6-0000-2ffba6aedf52 Instance ID : uuid:964c0f18-9c89-334c-a2d3-29251ed1e795 Format : application/pdf Title : telemac3d_user_manual_v7p1_ap_verif_JF Creator : J97554 Description : Description (x-repair) : Page Count : 108EXIF Metadata provided by EXIF.tools