open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi
I did a simulation by Telemac 2D for a rectangular channel (simple simulation). I use Strickler as a roughness law.
I have two problems when I selected a time step number equal to 0.1 and a number of the time steps equal to 10000.
The calculation stops and I find this notification:
Telemac System 5.6 to 6.1 - Perl scripts version 6.1
=========================================================
starting...
HOSTTYPE : win
PROJECT : C:\TELEMAC\V6P1
BASE DIRECTORY : C:\Users\KOUKA\Desktop\nazu
LAUNCH DIRECTORY : C:\Users\KOUKA\Desktop\nazu
WORK DIRECTORY : C:\Users\KOUKA\Desktop\nazu\cas2.txt6740_tmp
PARAMETER FILE : cas2.txt
*** Fichier de configuration SPECIFIQUE :
C:\TELEMAC\V6P1\config\systel.ini ***
*** Exploitation de la version specifique v6p1 ***
*** TELEMAC2D sur STATION ***
*** Lancement en interactif ***
*** VERSION V6P1 ***
________________________________________________________
Fichier des parametres : cas2.txt
________________________________________________________
Starting execution: telemac2d.bat
_________________________________________________________________________________________________________________________________
*** EXECUTABLE PAR DEFAUT ***
C:\TELEMAC\V6P1\telemac2d\tel2d_v6p1\win\telemac2dv6p1.exe
______________________________________________________________________________
*** ACQUISITION DES FICHIERS ***

- FICHIER DES PARAMETRES : cas2.txt
- DICTIONNAIRE : telemac2dv6p1.dico
- FICHIER DE GEOMETRIE : geo.slf
- FICHIER DES CONDITIONS AUX LIMITES : conlim.cli
______________________________________________________________________________
*** EXECUTION ***
LISTING DE TELEMAC-2D

---

---

TTTTT EEEEE L EEEEE M M AAAAA CCCCC
T E L E MM MM A A C
T EEE L EEE M M M AAAAA C
T E L E M M A A C
T EEEEE LLLLL EEEEE M M A A CCCCC

2D VERSION 6.1 FORTRAN 90
WITH SEVERAL TRACERS
COUPLED WITH SISYPHE AND TOMAWAC
LE MOT CLE : REMISE A ZERO DU TEMPS
EST CITE AU MOINS 2 FOIS, SEULE LA DERNIERE VALEUR EST CONSERVEE...
FIN DU FICHIER POUR DAMOCLES
********************************************
* LECDON: *
* APRES APPEL DE DAMOCLES *
* VERIFICATION DES DONNEES LUES *
* SUR LE FICHIER DES PARAMETRES *
********************************************

SORTIE DE LECDON. TITRE DE L'ETUDE :
Couplage avec sisyphe
OUVERTURE DES FICHIERS POUR TELEMAC2D
*****************************
* ALLOCATION DE LA MEMOIRE *
*****************************
READGEO1 : TITRE= newSelafin
NOMBRE D'ELEMENTS: 8509
NOMBRE REEL DE POINTS: 4529
MXPTEL (BIEF) : NOMBRE MAXIMUM D'ELEMENTS VOISINS D'UN POINT : 8
NOMBRE MAXIMUM DE POINTS VOISINS D'UNPOINT : 8
CORRXY (BIEF) : PAS DE MODIFICATION DES COORDONNEES
MAILLAGE : MESH ALLOUE
****************************************
* FIN DE L'ALLOCATION DE LA MEMOIRE : *
****************************************
INBIEF (BIEF) : MACHINE NON VECTORIELLE (SELON VOS DONNEES)
STRCHE (BIEF) : PAS DE MODIFICATION DU FROTTEMENT
IL Y A 2 FRONTIERE(S) LIQUIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 250 , DE NUMERO GLOBAL 4506
ET DE COORDONNEES : 4.000000 0.8000000E-02
FIN AU POINT DE BORD 273 , DE NUMERO GLOBAL 4529
ET DE COORDONNEES : 4.000000 0.3840000
FRONTIERE 2 :
DEBUT AU POINT DE BORD 524 , DE NUMERO GLOBAL 25
ET DE COORDONNEES : 0.000000 0.3880000
FIN AU POINT DE BORD 1 , DE NUMERO GLOBAL 1
ET DE COORDONNEES : 0.000000 0.000000
IL Y A 2 FRONTIERE(S) SOLIDE(S) :
FRONTIERE 1 :
DEBUT AU POINT DE BORD 1 , DE NUMERO GLOBAL 1
ET DE COORDONNEES : 0.000000 0.000000
FIN AU POINT DE BORD 250 , DE NUMERO GLOBAL 4506
ET DE COORDONNEES : 4.000000 0.8000000E-02
FRONTIERE 2 :
DEBUT AU POINT DE BORD 273 , DE NUMERO GLOBAL 4529
ET DE COORDONNEES : 4.000000 0.3840000
FIN AU POINT DE BORD 524 , DE NUMERO GLOBAL 25
ET DE COORDONNEES : 0.000000 0.3880000
CORFON (TELEMAC2D) : PAS DE MODIFICATION DU FOND
================================================================================
ITERATION 0 TEMPS : 0.0000 S
PRERES : NOMBRE DE COURANT MAXIMUM : 0.000000
STREAMLINE: USING PARALLEL VERSION OF CHARACTERISTICS 6.1
forrtl: severe (165): Program Exception - integer overflow
Image PC Routine Line Source
out6740_win.exe 004BABCB CHAR11 124 char11.f
Duree du calcul : 7 secondes ( 0:0:7 ) (systeme=0.047 sec)
______________________________________________________________________________
*** RESTITUTION DES FICHIERS ***
- FICHIER DES RESULTATS : resultatstrickler
- resultatstrickler
______________________________________________________
Execution terminee: telemac2d.bat
________________________________________________________
No compilation/linking/file errors detected.
No execution errors detected.
'.' n'est pas reconnu en tant que commande interne
ou externe, un programme exécutable ou un fichier de commandes.
## Error : System command failed for ./delete_cas2.txt6740.bat :256
then I change the time step of 0.01 with the same number of time steps 10000. Here there is a calculation but I don’t know why doesn’t work with a step of 0.1?? and I also I find a strange results of the water depth and velocity. I used a water depth equal to 0.15 m as a boundary condition in the downstream but the result of the water depth is between 0-0.01 m!!!
And the mean measured velocity was 0.2 m/s, but the calculated is between: 0.02-0.15
thank you in advince

best regrads

Hello,

withou some details about your case it is difficult to trace the error.
Can you send the case?

Best regards,
Clemens

Hi
thank you for your reply
you find in the attached document my case file

also the boundary conditions file and the geometry?

and these are the boundary conditions and the geometry

Ok,
according to your boundary condition file, at the east side of your domain you specify an imposed water level and at the west side the inflow. Since the outlet ist the FIRST liquid boundary and the inlet the SECOND liquid boundary, you should invert the values in your steering files in order to represent the right order of your liquid boundaries. With your definition e.g. you impose a water level of 0.0 m which can cause troubles, depending on your bottom heights etc..

Instead of:
COTES IMPOSEES : 0.0;0.15
DEBITS IMPOSES : 0.012;0.0

->
COTES IMPOSEES : 0.15;0.00
DEBITS IMPOSES : 0.00;0.012

And for river applications usually you don't need the Thompson boundary condition, see the Telemac manual.. So I would set: OPTION POUR LES FRONTIERES LIQUIDES :1

The case works but an advice for your mesh quality:
you should modify the edge elements at the liquid / solid boundaries in order to avoid that 3 points of the same element lie on a boundary. This can cause instabilities at the boundaries which propagate into your domain and they can be really annoying since they can prevent the simulation to get steady. You can do this in BlueKenue by -> Swap edge. After that you must create a new selafin object because you have a new connectivity table. Or in FUDAA (STBTEL) which has a function which does this automatically (Erase overstressed mesh).

Best regards,
Clemens

Hi Clemens

Thank you so much for your help and for your explanation, but just a question, I try to rectify the mesh by blue kenue, the swap edge is disabled, I didn’t use this option in T3 mesh toolbar, I don’t know how to activate this option
by fudaa, it’s simple, but, I find an extension of the new geo was .ser?? should I change the extension

best regards
MORRI

Hi,
in BlueKenue you should work directly on the mesh, in your case it has the name "FOND". E.g. you can save the FOND as t3s mesh file, open it, do what you want with it and after you are sure that all is right, create a new selafin object and drag the mesh file into it.

The extension doesn't matter but you can rename it in order to get conform with Telemac naming.

Best regards,
Clemens