open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everybody,

I am currently working on a comparison of models representing a current state and a project state. The initial model is validated and stable. The modification made is the exclusion of 4 buildings in the mesh by solid boundaries.
These modifications lead to strong variations in the results even within non-modified and non-hydraulically connected areas (torrential transitions).

The whole model seems to shake (attachments).

I've been going around the parameters for several iterations now, without finding any way to improve it. Your ideas are welcome.

Thank you very much in advance.
Jean

Here is the steering file :

/

---

/ ENVIRONNEMENT INFORMATIQUE
/

---

/
FICHIER DES CONDITIONS AUX LIMITES : './Cli.cli'
FICHIER DE GEOMETRIE : './V2.2_STB.slf'
FICHIER DES RESULTATS : './V2.2_STB.res'
FICHIER DES FRONTIERES LIQUIDES : './1800T_Q10_long.txt'
FICHIER DU CALCUL PRECEDENT : './ini.rep'
FICHIER DES COURBES DE TARAGE : './Barrage_ok.txt'
FICHIER DES SECTIONS DE CONTROLE : './controleG3b.txt'
FICHIER DE SORTIE DES SECTIONS DE CONTROLE : './resu.2.txt'
IMPRESSION DU CUMUL DES FLUX : OUI
NOMBRE MAXIMUM DE FRONTIERES : 9000
FICHIER DE DONNEES DES BRECHES : './breachG3b.txt'
/
/

---

/ OPTIONS GENERALES
/

---

/
/
VARIABLES POUR LES SORTIES GRAPHIQUES :'U,V,S,B,H,Q,
M,L,W,F,MAXV,MAXZ'
VARIABLES A IMPRIMER : ''
PAS DE TEMPS : 1.5
DUREE DU CALCUL : 125000
PERIODE POUR LES SORTIES GRAPHIQUES : 600
PERIODE DE SORTIE LISTING : 600
PROCESSEURS PARALLELES : 6
/
/

---

/ CONDITIONS INITIALES
/

---

/
/
CONDITIONS INITIALES : 'HAUTEUR NULLE'
SUITE DE CALCUL : OUI
REMISE A ZERO DU TEMPS : NON
BRECHE =true
/
/

---

/ CONDITIONS AUX LIMITES
/

---

/
LOI DE FROTTEMENT SUR LE FOND : 3
COEFFICIENT DE FROTTEMENT : 30
DEBITS IMPOSES : 10.0;10;10;0;0;0;10
COTES IMPOSEES : 0;0;0;196.5;195;198.1;0
COEFFICIENT DE DIFFUSION DES VITESSES : 0.01
PROFONDEUR LIMITE POUR PROCEDURE DE RESSUYAGE : 0.05
MODELE DE TURBULENCE : 1
COURBES DE TARAGE : 0;0;0;1;1;0;0
/
/

---

/ OPTIONS NUMERIQUES
/

---

/
BANCS DECOUVRANTS : OUI
BILAN DE MASSE : OUI
PRECISION DU SOLVEUR : 1.D-8
FORME DE LA CONVECTION : 1;5
SOLVEUR : 1
MASS-LUMPING SUR H : 1
MASS-LUMPING SUR LA VITESSE : 1
TRAITEMENT DU SYSTEME LINEAIRE : 2
PROFILS DE VITESSE : 4;4;4;1;1;1;4
COMPATIBILITE DU GRADIENT DE SURFACE LIBRE : 0.8
TRAITEMENT DES HAUTEURS NEGATIVES : 2
CORRECTION DE CONTINUITE : OUI
OPTION DE SUPG : 0;0;0;0
MAXIMUM D'ITERATIONS POUR LE SOLVEUR : 500
IMPLICITATION POUR LA VITESSE : 1
IMPLICITATION POUR LA HAUTEUR : 1
&ETA

Hello Jean,

Would it be possible to have figures of the 2 meshes please? In particular zooms located where you changed the mesh with your 4 excluded buildings (before and after). Please indicate the size of the elements (edge lengths)

Chi-Tuan