Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
TELEMAC-2D
PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY
Flat Threaded

TOPIC: PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY

PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY 3 weeks 2 days ago #44780

  • Laurie
  • Laurie's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 81
  • Thank you received: 3
Hi everybody !

It's not the first time this subject is presented. I've read various posts in this forum, but I didn't find my solution.

I'm using a TELEMAC2D modelling V7P1 (yes... old one :dry: ).
Before I modify loccaly the mesh, adding a new open boundary, the model run correctly.
Now ... I try to debog it.

I have an error message about a boundaty I didn't modify.

================================================================================
 ITERATION        0  TEMPS :   0.0000 S
 LE FICHIER DES FRONTIERES LIQUIDES CONTIENT
         640  LIGNES AVEC :
 Q(2)     Q(3)     
 USING STREAMLINE VERSION 7.2 FOR CHARACTERISTICS
 DEBIMP : PROBLEME SUR LA FRONTIERE     22
          DONNER UN PROFIL DE VITESSES       
          DANS LE :                          
          FICHIER DES CONDITIONS AUX LIMITES 
          OU VERIFIER LES HAUTEURS D'EAU.   
          AUTRE CAUSE POSSIBLE :             
          ENTREE TORRENTIELLE A HAUTEUR LIBRE
          METTRE UNE HAUTEUR NON NULLE       
          DANS LES CONDITIONS INITIALES      
          OU IMPOSER LA HAUTEUR D'EAU EN ENTREE.



 PLANTE : ARRET DU PROGRAMME APRES ERREUR
 RETURNING EXIT CODE:            2

Sorry it's in French.

I firstly verify that this boundary has water. It seems ok (around 1-2m). It's imposed with a surface elevation in the geometry file and in a *.rep (hotstart).
(I used the *.rep ini, I've deleted the 3 layers with old mesh to add my 3 layers with ma new mesh ... maybe that way is wrong).


It seems my boundary n°22 doesn't work correctly, or I forget to verify something.

Dam/bridge structures boundaries are "handled" with a level/discharge law and a fortran file.

I join the model.

File Attachment:

File Name: cas_2024-05-23.txt
File Size: 3 KB


File Attachment:

File Name: tarage_barrage_V2.f
File Size: 27 KB


File Attachment:

File Name: lois_tarage_13barrages.txt
File Size: 4 KB




File Attachment:

File Name: MeshQV_02_STB.cli
File Size: 964 KB





Maybe someone could help me to figure out what i'm doing wrong.




In the fortran, I've updated :
C****************************** Module VAR_UTIL *****************************
      MODULE VAR_UTIL

      INTEGER, PARAMETER :: NB_barrage = 13
      DOUBLE PRECISION, DIMENSION(NB_barrage) :: 
     &ind_front_tar = (/5,7,9,11,12,14,18,19,21,23,26,27,29/),
     &ind_front_Qimp =(/4,6,8,10,13,15,17,16,20,22,25,28,24/)
      DOUBLE PRECISION, DIMENSION(NB_barrage),SAVE ::Q_barrage

     
      END MODULE VAR_UTIL

I've checked that the upstream limits of my structures are associated with their downstream limits (I've matched the number of the limits "&ind_front_tar" and "&ind_front_Qimp").

In the steering file (CAS), I have (I think) correctly associated the ‘water level’ limits with the imposed level and the calibration curves (apparently it's the Fortran file that associates the water level conditions with the dischare conditions), and the ‘discharge’ with the imposed flows.

I have 29 boundaries.
n°1 : free ouptut.
n°2 et 3 : input discharge.
n°4 to 29 are dams are bridge modelled using water level/discharge law (calibration curves).

DEBITS IMPOSES=
0;2;3;4;0;6;0;8;0;10;0;0;13;0;15;16;17;0;0;20;0;22;0;24;25;0;0;28;0
PROFILS DE VITESSE=1;4;4;4;1;4;1;4;1;4;1;1;4;1;4;4;4;1;1;4;1;4;1;4;4;1;1;4;1

COTES IMPOSEES=
1;0;0;0;5;0;7;0;9;0;11;12;0;14;0;0;0;18;19;0;21;0;23;0;0;26;27;0;29
OPTION POUR LES FRONTIERES LIQUIDES=
1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1

COURBES DE TARAGE =1;0;0;0;1;0;1;0;1;0;1;1;0;1;0;0;0;1;1;0;1;0;1;0;0;1;1;0;1

So, have you any advice helping me ?
Thank you in advance,

Laurie
The administrator has disabled public write access.

PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY 3 weeks 2 days ago #44782

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3653
  • Thank you received: 1015
Hi
Read the user manual, particularly the rules of syntax for the steering file
Christophe
The administrator has disabled public write access.
The following user(s) said Thank You: Laurie

PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY 3 weeks 2 days ago #44787

  • Laurie
  • Laurie's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 81
  • Thank you received: 3
Tkanks c.coulet.

I just done it.

"1. The keywords may be of Integer, Real, Logical or Character type," --> ok (I think)

"2. The order of keywords in the steering file is of no importance, --> ok

"3. Each line is limited to 72 characters. However, it is possible to pass from one line to the next as often as required, provided that the name of the keyword is not split between two lines," --> I corrected it.


"4. For keywords of the array type, the separator between two values is the semi-colon. It is not necessary to give a number of values equal to the size of the array. In this case, DAMOCLES returns the number of read values. For example:
TYPE OF ADVECTION = 1;5
(this keyword is declared as an array of 4 values)" --> ok

But I still have the same error.

Corrected cas :

File Attachment:

File Name: cas_2024-05-23-2.txt
File Size: 3 KB


Unfortunately, that doesn't change anything.
I've double-checked that the limits are associated with flow rate or water level. There is water at the beginning of the model.

In any case, there's something I don't understand: the same initial condition, the same CAS file... works just fine before I modify the mesh locally. The errors concern either limit 2 or 22, depending on whether I take into account the result files from the previous calculation or not. It makes no sense to me.


Do you have any other ideas about what I could check?
The administrator has disabled public write access.

PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY 3 weeks 2 days ago #44788

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3653
  • Thank you received: 1015
Some manipulation you made are unclear for me...
As you modified the mesh, I suspect the initial condition could be the problem...

Try to investigate this way. You could first remove the previous computation file and test a constant value (which ensure water on every boundary) to check if the simulation starts
Regards
Christophe
The administrator has disabled public write access.
The following user(s) said Thank You: Laurie

PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY 3 weeks 2 days ago #44796

  • Laurie
  • Laurie's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 81
  • Thank you received: 3
Thank you very much for your help and advice.

I've taken your comments on board. As the ‘simple’ checks didn't work, I've gone back to basics:
- Simple model, no work. Only the necessary boundary conditions (inlet flow, outlet water level), and constants.
- Then I gradually added the result file from the previous calculation and the structures.
After a few hours of trouble, it worked.

The problem was almost certainly due to incorrect numbering of the boundaries. There are 29 of them. That's getting to be a lot!

In short, there's nothing like working step by step to unblock a situation ;)
The administrator has disabled public write access.
Flat Threaded
Moderators: pham
Forum
System Modules
TELEMAC-2D
PROBLEM ON BOUNDARY NUMBER 2. GIVE A VELOCITY PROFILE IN THE BOUNDARY

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.