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