open TELEMAC-MASCARET The mathematically superior suite of solvers

hello;
after launching my simulation case (you find the file cas below) the simulation was succefful but it seems without any calculations or itérations in the mass bilan ...and all what is showed is zeros in bilan intial nd final for tracer and water mass
...where the problem is came from ?
DEBUGGER : 1
/
STOCKAGE DES MATRICES : 3
/

FICHIER DES CONDITIONS AUX LIMITES : week.cli
FICHIER DE GEOMETRIE : week.slf
FICHIER DES RESULTATS 3D : r3d_source.slf
FICHIER DES RESULTATS 2D : r2d_source.slf
FICHIER FORTRAN : tide.f
/FICHIER DES SOURCES : src_source.txt

/
NOMBRE DE PLANS HORIZONTAUX : 10
TRANSFORMATION DU MAILLAGE=1

/
/

---

/
/ OPTIONS GENERALES
/

---

/
/TITRE = 'TELEMAC 3D : TEST POUR SOURCE'
/SUITE DE CALCUL : OUI
VARIABLES POUR LES SORTIES GRAPHIQUES 2D : 'U,V,H,S,B,TA1'
VARIABLES POUR LES SORTIES GRAPHIQUES 3D : 'Z,U,V,W,TA1'
PAS DE TEMPS = 100
DUREE DU CALCUL = 4320
PERIODE POUR LES SORTIES GRAPHIQUES : 100
PERIODE POUR LES SORTIES LISTING : 100
VERSION NON-HYDROSTATIQUE = OUI
/TRAITEMENT DES HAUTEURS NEGATIVES = 2
/

---

/
/ PPE
/

---

/
/

---

/
/ OPTIONS GENERALES
/

---

/
SOLVEUR POUR PPE=8
MAXIMUM D'ITERATIONS POUR PPE=1000
PRECISION POUR PPE=0.00000001
/OPTION DU SOLVEUR POUR PPE=
PRECONDITIONNEMENT POUR PPE=34
/

---

/
/ CONDITION INITIALE
/

---

CONDITIONS INITIALES = 'COTE CONSTANTE'
COTE INITIALE =-100
DEBITS IMPOSES =-50000;50000;50000
/

---

/
/ SCHEMAS POUR LA CONVECTION
/

---

/
ETAPE DE CONVECTION = OUI
SCHEMA POUR LA CONVECTION DES VITESSES = 2
SCHEMA POUR LA CONVECTION DES TRACEURS :2
/

---

/
/ SCHEMAS POUR LA DIFFUSION
/

---

/

SCHEMA POUR LA DIFFUSION DES VITESSES = 1
SCHEMA POUR LA DIFFUSION DES VITESSES =1
SCHEMA POUR LA DIFFUSION DES TRACEURS =1
OPTION POUR LA DIFFUSION =1
/MASS-LUMPING POUR LA HAUTEUR =1

/

---

/
/ PROPAGATION
/

---

SOLVEUR POUR LA PROPAGATION : 7
PRECISION POUR LA PROPAGATION = 1.E-6
IMPLICITATION POUR LES VITESSES = 0.6
PRECONDITIONNEMENT POUR LA PROPAGATION : 2

/

---

/
/ TURBULENCE
/

---

LOI DE FROTTEMENT SUR LE FOND : 3
LOI DE FROTTEMENT SUR LES PAROIS LATERALES : 0
COEFFICIENT DE FROTTEMENT POUR LE FOND = 100.
MODELE DE TURBULENCE HORIZONTAL = 1
MODELE DE TURBULENCE VERTICAL = 2
COEFFICIENT DE DIFFUSION HORIZONTAL DES VITESSES : 0.000001
COEFFICIENT DE DIFFUSION VERTICAL DES VITESSES : 0.000001
/

---

/
/ NUMERIQUE
/

---

/
MAXIMUM D'ITERATIONS POUR LA PROPAGATION : 100
IMPLICITATION POUR LA HAUTEUR = 0
IMPLICITATION POUR LES VITESSES = 0
IMPLICITATION POUR LA DIFFUSION =0
ORDRE DU TIR INITIAL POUR LA HAUTEUR : 1
/

---

/
BILAN DE MASSE : OUI
INFORMATION SUR LE BILAN DE MASSE A CHAQUE SORTIE LISTING : OUI
/

---

/
/ TRACEUR
/

---

/
BANCS DECOUVRANTS : NON
NOMBRE DE TRACEURS = 1
PRECISION POUR LA DIFFUSION DES TRACEURS = 1.E-7
SOLVEUR POUR LA DIFFUSION DES TRACEURS =7
VALEURS INITIALES DES TRACEURS = 10
ABSCISSES DES SOURCES = 350.102
ORDONNEES DES SOURCES = 30.927
COTES DES SOURCES = -26
DEBITS DES SOURCES = 100
VALEURS DES TRACEURS DES SOURCES =828
/

---

/
/ SOLVEURS
/

---

SOLVEUR POUR LA DIFFUSION DES VITESSES =8
SOLVEUR POUR LA PROPAGATION =8
SOLVEUR POUR PPE =8
SOLVEUR POUR LA VITESSE VERTICALE =8
SOLVEUR POUR LA DIFFUSION DES TRACEURS =8
OPTION DU SOLVEUR POUR LA DIFFUSION DES VITESSES=10
/OPTION DU SOLVEUR POUR LA PROPAGATION =3
/OPTION DU SOLVEUR POUR PPE =3
OPTION DU SOLVEUR POUR LA DIFFUSION DES TRACEURS =20
MAXIMUM D'ITERATIONS POUR LA DIFFUSION DES VITESSES =1000
MAXIMUM D'ITERATIONS POUR LA PROPAGATION =1000
MAXIMUM D'ITERATIONS POUR PPE =1000
MAXIMUM D'ITERATIONS POUR LA VITESSE VERTICALE =1000
MAXIMUM D'ITERATIONS POUR LA DIFFUSION DES TRACEURS =1000
MAXIMUM D'ITERATIONS POUR LA DIFFUSION DU K-EPSILON =1000
/ MAXIMUM D''ITERATIONS POUR LES SCHEMAS DE CONVECTION =1000
/

---

/
/ PRECONDITIONNEMENT
/

---

PRECONDITIONNEMENT POUR LA DIFFUSION DES VITESSES =2
PRECONDITIONNEMENT POUR LA PROPAGATION =2
PRECONDITIONNEMENT POUR PPE =2
PRECONDITIONNEMENT POUR LA VITESSE VERTICALE =2
PRECONDITIONNEMENT POUR LA DIFFUSION DES TRACEURS=2




&ETA
&FIN

Hello,

I see in your steering file :

PAS DE TEMPS = 100
DUREE DU CALCUL = 4320
PERIODE POUR LES SORTIES GRAPHIQUES : 100

A duration of 4320 s with a time step of 100. s makes 43 time steps, with a period for graphic printouts of 100 you will thus have only the initial conditions in the results, idem for the listing printout... The misunderstanding must be that the period is in number of time steps, not in seconds, but then why a duration of 4320 if your period allows to see only 4300 ?

JMH

thank you for your suggestion ,but even I rectify it by :
PAS DE TEMPS =1
DUREE DE CALCUL =1000
and I tested also other values but the lising controle give my the same thing ....."my work is done " but with zeros!!
could the problem came from schemes ?

hello, I joind below also a part of the controle listing ...I can't understand why it works such as it don't do any thing !! appel...retour ...appel ...retour
best regards

Hello,

Remove DEBUGGER : 1 to avoid "retour... appel" etc., the listing will be clearer. Then think of the choice of the time step, is it 100. or 1., what corresponding CFL number it makes ?

Then you have:

IMPLICITATION POUR LA HAUTEUR = 0
IMPLICITATION POUR LES VITESSES = 0
IMPLICITATION POUR LA DIFFUSION =0

whereas these numbers should be greater than 0.5. Please inspire from an existing example. "my work is done" is a message from Python, but the listing must say "fin normale du programme".

With best regards,

Jean-Michel Hervouet

think you very mutch dear Micheal,
the problem of "appel...retour " is solved , then I tried with a new case file basing on an example ,and the simulation runs correctly but it stills one small problem that I can't visualize the evolution of my tracer slowly ,the simulation give a very quick tracer evolution which overruns all my domain.
I joind below my case file if there is some thing wrongly choosen or some remarks
and also images of the BlueKenue results are below .
best regards
think you

HELLO ,
my problem is solved now; but it stills one problem ....when I animate result in 3d I observe that all my domain moves and perturbated . I can't understand this phenomenon!!
please help me to correct this p^roblem in the earlier time.
thank you

Hello,

The 3D domain follows the free surface, so it moves in the vertical direction, this is normal. However if you see also the horizontal coordinates moving, this is very weird...

Regards,

JMH