TELEMAC-MASCARET Forum: Problem in defining rigid bed (1/1)

Welcome, Guest

TOPIC: Problem in defining rigid bed

Problem in defining rigid bed 9 years 10 months ago #15503

amanj2013
OFFLINE
Expert Boarder
Posts: 211
Thank you received: 24

Hello everyone,

I have an issue and I would to share with you in order to figure-out and get your opinion on it. I am working on simulating an overtopping non cohesive sand embankment by coupling Telemac 2d FE and Sisyphe. When I am introducing the zero level rigid bed to my sisyphe, my outflow peak discharge reduced by 80% while in experimental physical modelling, the value is much closer to erodable bed but in fact, the model is non erodable bed in the lab model. Hope some one help us to see where is the point that makes this big difference.

Looking forward to hear from you

Sincerely,
AMANJ

Attachments:

27.png (83KB)
28.png (45KB)
29.png (37KB)
30.png (27KB)
31.png (26KB)
32.png (328KB)
33.png (44KB)

The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15505

amanj2013 OFFLINE Expert Boarder Posts: 211 Thank you received: 24	here is the fortran file Thanks AMANJ Attachments: noerod.f (3KB)
	The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15522

cyamin OFFLINE openTELEMAC Guru Posts: 997 Thank you received: 234	Hello Amanj, Uncomment the control at the end of the fortran file and run again. It is this part: ! DO I=1,NPOIN ! IF(ZR(I).GT.ZF(I)) THEN ! WRITE(LU,*) 'POINT ',I,' NON ERODABLE BED HIGHER THAN BED' ! CALL PLANTE(1) ! STOP ! ENDIF ! ENDDO Costas
	The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15527

amanj2013 OFFLINE Expert Boarder Posts: 211 Thank you received: 24	Hello dear Costas, I tried what you mentioned above but I got an error message at the beginning of the run. Pls take a look at the enclosed images. Thanks AMANJ Attachments: 64.png (13KB) 65.png (71KB)
	The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15528

Serrer OFFLINE Moderator Posts: 720 Thank you received: 304	Hi Amanj, Fortran statements MUST start in column 7 Shift the "DO" statement to the right by 1 space and the "ENDDO" statement to the right by 2 spaces. Cheers... Martin
	The administrator has disabled public write access. The following user(s) said Thank You: amanj2013

Problem in defining rigid bed 9 years 10 months ago #15531

amanj2013 OFFLINE Expert Boarder Posts: 211 Thank you received: 24	Thank you dear Martin, Now is running :) AMANJ
	The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15533

amanj2013
OFFLINE
Expert Boarder
Posts: 211
Thank you received: 24

Hello,

Dear Costas,I did what you said and it gives me the same result. I believe something wrong in calculating discharge through my control section because in my control discharge the peak discharge is about 20 l/sec while in the screen printout shows at some certain time that FLUX Boundary 2: about 68 l/sec which more realistic. The question is that how to change or correct the discharge through my control section? or how to print into a file the discharge through my second liquid boundary which is downstream liquid boundary?

Thanks

Attachments:

33_2015-01-20.png (44KB)
control_sections.txt (0KB)

The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15534

amanj2013 OFFLINE Expert Boarder Posts: 211 Thank you received: 24	here is the screen result at Q peak: Attachments: 001.png (53KB) 002_2015-01-20.png (54KB) 003_2015-01-20.png (61KB) 004.png (56KB) 005.png (62KB) 006.png (61KB) 007.png (61KB) 008.png (67KB)
	The administrator has disabled public write access.

Problem in defining rigid bed 9 years 10 months ago #15535

amanj2013 OFFLINE Expert Boarder Posts: 211 Thank you received: 24	more screenshot Attachments: 009.png (63KB) 010.png (60KB) 011.png (61KB)
	The administrator has disabled public write access.