open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everyone, I'm a new Telemac 2D user who working for my thesis about Effect of Climate Change on River Morphology. The river been simulated is assumed to be steady state river. I have created the boundary conditions file and geometry file using Blue Kenue. For the Boundary Conditions File, Prescribed Q is used at the inlet while Prescribed H is indicated at the outlet. The cas file is created by using notepad.

But there are errors happened when simulating in Telemac 2D. I'm not sure what the problems are, either steering file or the meshing files. I wish to ask advises from everyone. Please help me find out what I'm missing. Your responses will be really appreciated. Thanks.

I have attached the slf, cli and cas files.
And the log is shown below:


Loading Options and Configurations
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

_____ __ __
|___ | /_ | /_ |
__ __ _/ / _ __ | | _ __ | |
\ \ / / |_ _| | '_ \ | || '__| | |
\ V / / / | |_) | | || | | |
\_/ /_/ | .__/ |_||_| |_|
| |
|_|
_ _ __ ___ _ _ __ _____
_| || |_ /_ | / _ \ | || | /_ ||___ |
_ __ ___ __ __ |_ __ _| | || | | || || |_ | | _/ /
| '__| / _ \\ \ / / _| || |_ | || | | ||__ _| | | |_ _|
| | | __/ \ V / _ |_ __ _| | || |_| | | | | | / /
|_| \___| \_/ (_) |_||_| |_| \___/ |_| |_| /_/


... parsing configuration file: C:\opentelemac-mascaret\v7p1\configs\systel_v7p1.cfg


Running your CAS file for:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

+> configuration: win7gfors
+> root: C:\opentelemac-mascaret\v7p1


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


... reading the main module dictionary

... processing the main CAS file(s)
+> running in English

... handling temporary directories

... checking coupling between codes

... checking parallelisation

... first pass at copying all input files
copying: Tutorial_BK_geo.slf C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\T2DGEO
copying: FOND_BC.cli C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\T2DCLI
re-copying: C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\T2DCAS
copying: telemac2d.dico C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\T2DDICO

... checking the executable
re-copying: telemac2d.exe C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\out_telemac2d.exe

... handling sortie file(s)


Running your simulation(s) :
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\out_telemac2d.exe


_ _ ___ _ _____ __
| | | | |__ \ | | |___ | /_ |
| |_ ___ | | ___ _ __ ___ __ _ ___ ) | __| | ______ __ __ _/ / _ __ | |
| __| / _ \| | / _ \| '_ ` _ \ / _` | / __| / / / _` | |______| \ \ / / |_ _| | '_ \ | |
| |_ | __/| || __/| | | | | || (_| || (__ / /_ | (_| | \ V / / / | |_) | | |
\__| \___||_| \___||_| |_| |_| \__,_| \___||____| \__,_| \_/ /_/ | .__/ |_|
| |
|_|

LISTING OF 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 7.2 FORTRAN 90
WITH SEVERAL TRACERS
COUPLED WITH SISYPHE AND TOMAWAC

END OF FILE FOR DAMOCLES

********************************************
* LECDON: *
* AFTER CALLING DAMOCLES *
* CHECKING OF DATA READ *
* IN THE STEERING FILE *
********************************************

EXITING LECDON. NAME OF THE STUDY:
PAHANG RIVER

OPENING FILES FOR TELEMAC2D

*****************************
* MEMORY ORGANIZATION *
*****************************

READ_MESH_INFO: TITLE= newSelafin
NUMBER OF ELEMENTS: 32444
NUMBER OF POINTS: 17453

FORMAT NOT INDICATED IN TITLE

MXPTEL (BIEF) : MAXIMUM NUMBER OF ELEMENTS AROUND A POINT: 9
MAXIMUM NUMBER OF POINTS AROUND A POINT: 9
CORRXY (BIEF):NO MODIFICATION OF COORDINATES

MESH: MESH ALLOCATED

*************************************
* END OF MEMORY ORGANIZATION: *
*************************************

_____________
runcode::main:
:
|runCode: Fail to run
|C:\opentelemac-mascaret\v7p1\examples\sgphg\sgphg.cas_2018-02-12-14h27min46s\out_telemac2d.exe
|~~~~~~~~~~~~~~~~~~
|Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
|
|Backtrace for this error:
|~~~~~~~~~~~~~~~~~~

Hello

Telemac says in the listing that the boundary condition file does not correspond to the geometry file. So please generate it again based on the geometry file.

with my kind regards

Riadh

Hi Riadh,
Thanks for the reply.
I have tried to regenerate the BC files based on the geometry files but the outcome in Telemac 2D still the same.
I have attached my work in Blue Kenue with the latest geometry, BC and cas files.
Thanks again for the help.

Best regards,
Esther

Hello Ester

you have a problem in the mesh in node number 14315 (see attached photo)



1 node can not be in the same time within 2 different boundaries. You have to refine the mesh in these location in order to separate the two islands.

I hope this helps

with my best regards

Riadh

Thanks Riadh for pointing out the problems, it works!
After running the Telemac2D, the log I got was as above:

================================================================================
ITERATION 8500 TIME: 23 H 36 MIN 40.0000 S ( 85000.0000 S)

---

ADVECTION STEP

---

DIFFUSION-PROPAGATION STEP
RESCJG (BIEF) : 54 ITERATIONS, ABSOLUTE PRECISION: 0.9790493E-04

---

BALANCE OF WATER VOLUME
VOLUME IN THE DOMAIN : 0.1383778E+08 M3
FLUX BOUNDARY 1: -0.3227796 M3/S ( >0 : ENTERING <0 : EXITING )
FLUX BOUNDARY 2: -0.000000 M3/S ( >0 : ENTERING <0 : EXITING )
RELATIVE ERROR IN VOLUME AT T = 0.8500E+05 S : 0.4793199E-09

================================================================================
ITERATION 8600 TIME: 23 H 53 MIN 20.0000 S ( 86000.0000 S)

---

ADVECTION STEP

---

DIFFUSION-PROPAGATION STEP
RESCJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS: 100 RELATIVE PRECISION: 0.1114281E-03

---

BALANCE OF WATER VOLUME
VOLUME IN THE DOMAIN : 0.1383663E+08 M3
FLUX BOUNDARY 1: -1.819715 M3/S ( >0 : ENTERING <0 : EXITING )
FLUX BOUNDARY 2: -0.000000 M3/S ( >0 : ENTERING <0 : EXITING )
RELATIVE ERROR IN VOLUME AT T = 0.8600E+05 S : -0.1246394E-09

================================================================================
ITERATION 8640 TIME: 1 D 0 H 0 MN 0.0000 S ( 86400.0000 S)

---

ADVECTION STEP

---

DIFFUSION-PROPAGATION STEP
RESCJG (BIEF) : 40 ITERATIONS, RELATIVE PRECISION: 0.9990622E-04

---

BALANCE OF WATER VOLUME
VOLUME IN THE DOMAIN : 0.1383584E+08 M3
FLUX BOUNDARY 1: -2.103345 M3/S ( >0 : ENTERING <0 : EXITING )
FLUX BOUNDARY 2: -0.000000 M3/S ( >0 : ENTERING <0 : EXITING )
RELATIVE ERROR IN VOLUME AT T = 0.8640E+05 S : 0.8029971E-10

---

FINAL BALANCE OF WATER VOLUME

RELATIVE ERROR CUMULATED ON VOLUME: -0.2468063E-02

INITIAL VOLUME : 0.3209050E+08 M3
FINAL VOLUME : 0.1383584E+08 M3
VOLUME THAT ENTERED THE DOMAIN: -0.1833386E+08 M3 ( IF <0 EXIT )
TOTAL VOLUME LOST : -79201.36 M3

END OF TIME LOOP

CORRECT END OF RUN

ELAPSE TIME :
23 MINUTES
3 SECONDS


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

... handling result files

+> sgphg.cas
moving: trial_SP_results.slf
copying: sgphg.cas_2018-02-13-15h19min02s.sortie


My work is done



Is that meant my simulation success? Thanks a lot for the reply.

Best regards,

Esther

Hello Ester

yes the simulation is done, but not sure that the result is OK.
I see in the listing that :
RESCJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS: 100 RELATIVE PRECISION:

which means that the solver is not converging.
reduce your time step. and run it again.

kind regards

Riadh