Date: 24th May 2001

The principal features of this version v5p1 of TELEMAC-3D are:

• A number of new options and modifications programmed by Sogreah engineers (turbulence, source terms, etc.)
• A re-writing of the input/outputs part, to use the common subroutines from BIEF.
• Source terms moved to the diffusion step.
• Diffusion step in tetrahedrons, discretized in the physical space.
• New options of Telemac-2D available (edge-based storage, wave equation).

The keyword TURBULENCE MODEL has been split into two: HORIZONTAL TURBULENCE MODEL and VERTICAL TURBULENCE MODEL. There are two main innovations for turbulence:

- The introduction of two new mixing length model for the vertical turbulence model: Quetin and Tsanis models which are very useful for wind forcing.

- A new turbulence model for the horizontal turbulence :Smagorinsky.

Integer value. The default value is 1

Six choices are currently possible:

1 Constant viscosity

2 Classical mixing length

3 k-epsilon model

4 Mixing length with jet

5 Mixing length “Quetin model”

6 Mixing length “Tsanis model”

Integer value. The default is 1

This is the previous NPRIV in the PRINCI program.

Three choices are currently possible:

1 Constant viscosity

3 k-epsilon model

7 Smagorinsky model

Some improvement in the management of sources. We can now introduce sources in the domain like in TELEMAC-2D.

Real value. The default value is None

Real numbers giving the abscissae of possible water inflow (also called sources) in the domain (in meters).

Real value. The default value is None

Real numbers giving the ordinate of possible water inflow (also called sources) in the domain (in meters).

Real value. The default value is None

Real numbers giving the elevation of possible water inflow (also called sources) in the domain (in meters).

Real value. The default value is None

Real numbers giving the flowrates at the sources in m3/s.

Real value. The default value is None

Real numbers giving the value of the active tracers at the source.

Real value. The default value is None

Real numbers giving the value of the passive tracers at the source.

Real value. The default value is None

Real numbers specified the component of the velocity ejection at each source (m/s) along X.

Real value. The default value is None

Real numbers specified the component of the velocity ejection at each source (m/s) along Y.

There is now only tracers, not passive and active tracers. Consequently all the duplicated key-words have been simplified and do not mention 'active' or 'passive'. All the tracers are potentially active, depending on the density law used.

Example :

ADVECTION OF PASSIVE TRACERS and

ADVECTION OF ACTIVE TRACERS

Are now only :

ADVECTION OF TRACERS

For density laws, see the key-word DENSITY LAW below

Character strings with ”;” as a separator. The default value is None

Names of tracers as in NOMVAR_TELEMAC3D subroutine are given in 32 characters (16 for the name, 16 for the unit).

When not specified by the user, the names are:

“TRACER ” and ”?? ” for unit (active tracers) (english)

“TRACEUR ” and ”?? ” (french)

The MNEMO strings used in key-word VARIABLES FOR GRAPHIC PRINTOUTS are: TA1, TA2, … for tracers

The names of viscosities associated with each tracer are called: NUX NAME_OF_TRACER (unit M2/S), NUY NAME_OF_TRACER (unit M2/S), and NUZ NAME_OF_TRACER (unit M2/S)

And the MNEMO strings are, for example, NAX1, NAY1, NAZ1 for tracer number 1.

Integer value. The default value is 0

Replaces the key-word DENSITY ACCORDING TO T AND S (DENSITE FONCTION DE T ET S)

if 0 : no density variations taken into account

if 1 : density function of temperature in which case, a tracer must be called TEMPERATURE (use key-word NAMES OF TRACERS)

if 2 : density function of salinity in which case, a tracer must be called SALINITY

if 3 : density function of both temperature and salinity in which case, a tracer must be called TEMPERATURE and a tracer must be called SALINITY

if 4 : density function of all tracers, in which case, the density will be computed with all tracers and the “BETA EXPANSION COEFFICIENT FOR ACTIVE TRACERS” will be used (it is now an array and the SUBROUTINE INBETA has been suppressed). If it is different between tracers, use user SUBROUTINE INBETA. The beta coefficients are used only in case 4 (built in coefficients are used in cases 1, 2 and 3, in SUBROUTINE DRSURR)

The order of tracers has no importance, in cases 1, 2 and 3, only their names is important.

Integer value, The default value is 1

This is the Telemac-2D key-word with:

1 : classical EBE storage

3 : edge-based storage

So far it is used in Telemac-3D only when calling Telemac-2D (as you will be notified by a message if youtry option 3). Option 3 will save a few percents (up to 15% in some cases) in the computation times.

Character*72 value. The default value is ZERO ELEVATION (COTE NULLE)

This is now exactly like in TELEMAC-2D, choices are :

ZERO ELEVATION (COTE NULLE)

CONSTANT ELEVATION (COTE CONSTANTE), in which case the new key-word INITIAL ELEVATION (COTE INITIALE) must be used.

ZERO DEPTH (HAUTEUR NULLE)

CONSTANT DEPTH (HAUTEUR CONSTANTE), in which case the new key-word INITIAL DEPTH (HAUTEUR INITIALE) must be used.

SPECIAL (PARTICULIERES), in which case SUBROUTINE CONDIM must be modified.

All the choices will act only on the depth, velocities are set to 0.

Another possibility for initial conditions is the key-word: 2D CONTINUATION (SUITE 2D). If set to yes, the BINARY FILE 1 (FICHIER BINAIRE 1) will be considered as a 2D result and read to find the depth and the two components of the horizontal velocity.

Integer array. The default value is 1;1

This is a new key-word for driving different options when TELEMAC-2D is called (hydrostatic mode).

So far two integers, separated by semicolons.

First integer:

if 1 : classical technique, TELEMAC-2D only deals with propagation terms.

if 2 : new option, TELEMAC-2D solves a full Saint-Venant problem.

Second integer:

if 1 : classical technique, primitive equations of TELEMAC-2D

if 2 : new option, “wave equation” in TELEMAC-2D. This latter option will save time, but is rather suited for steady-state computations.

Integer array. The default value is 0;0

This is the new name for NUMBER OF PROCESSORS (NOMBRE DE PROCESSEURS), for coherence with TELEMAC-2D.

0 or 1 means that there is no parallelism. With 1, the parallel libraries are linked, with 0 they are not.

This is now exactly as in TELEMAC-2D, for prescribing boundary conditions variable in time. Functions Q3, SL3 and VIT3 have been created and will act as Q, SL and VIT in 2D.

Here is an example of liquid boundaries file:

#
# Variations in time at liquid boundaries
#
T Q(1) SL(2)
s m3/s m
0. 10. 135.
25. 12. 135.2
100. 20. 136.

Lines beginning with # are comments. Line with T Q(1) SL(2) says that the time is the first column, the discharge at liquid boundary number 1 is the second column and the free surface elevation at liquid boundary 2 is the third column.

The following line is mandatory but not used, you may give the units or explicit names for each column. Then the variations in time are given, T being the first figure of the line. The only requirement for the values of T is that they must be increasing. The time of your computation must be in the range described by the file. Interpolations will be done to cope with your time-steps. A linear variation thus needs only two lines of data.

Important notes:

1. The numbering is the numbering of liquid boundaries, as given at the beginning of the listing (by SUBROUTINE FRONT2, or by the boundary conditions file in parallelism)
2. The key-words PRESCRIBED ELEVATIONS, PRESCRIBED DISCHARGES and PRESCRIBED VELOCITIES must be given in the steering file, with constant values, even if the values will be read after in the liquid boundaries file.
3. All this automatic setting of boundary conditions can be erased in the user SUBROUTINE BORD3D, and changed into a simple setting of prescribed values. Q3, SL3 and VIT3 may also be reprogrammed for complex cases.

The SUBROUTINE BORD3D has undergone many changes, see also Q3, SL3 and VIT3.

The SUBROUTINE VISCO has been modified to take into account the two new keywords for turbulence models.

The SUBROUTINE Q3D and SUBROUTINE SOURCES have been modified to take into account the sources.

The SUBROUTINE LONGML and SUBROUTINE VISCLM were modified by the introduction of the new turbulence models.

The SUBROUTINE TRISOU includes the new improvement for the sources.