open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everyone,

I tried to simulate an overtopping embankment flow by using 2D Telemac FV scheme. Whew I am running without fortran file is fine, but when coupling with my fortran file which i copies from Malpasset 2D case and changes the coordinate. However, when I use the 2D malpasset fortran file in my 3D telemac case everything going well. Another problem what I have with 2D case is not consider for my prescribed inflow rate while my boundary file is fine in 3D and working well.

Please any suggestion?

I will attach the log and other files

Hello,

What happens is that this Malpasset Fortran file contains an old subroutine condin which is correct but does not consider keywords in English. So either you change your keyword INITIAL CONDITION and put :

INITIAL CONDITION : 'COTE CONSTANTE'

or you change subroutine condin and replace the line :

ELSEIF(CDTINI(1:14).EQ.'COTE CONSTANTE') THEN

by:

ELSEIF(CDTINI(1:18).EQ.'CONSTANT ELEVATION'.OR.
& CDTINI(1:14).EQ.'COTE CONSTANTE' ) THEN

So this does not seem to be linked with Finite Volumes, unless there is another problem behind, but let's see.

With best regards,

Jean-Michel Hervouet

Hello dear Jean-Michel Hervouet,

Thank you so much for your help. When I changed the initial condition into COTE CONSTANTE it works well, but still I have a problem of my boundary inflow. I have two open boundaries, one with prescribed discharge U/S and the other head .D/S. I set the values in the steering value as follow: for depth: 0.001,0 and for flow: 0,0.001
but when I run my model, doesnt consider for these value. I tried the same geometry mesh and boundary file in Telemac3D, it works very well and considered for inflow conditions.

Any idea please.
I have all files in the previous post.

Best Regards,
AMANJ

/

---

/ TELEMAC2D Version v6p3 6-Dec-2014
/ amanj
/

---

/

---

/ EQUATIONS
/

---

FRICTION COEFFICIENT =92

LAW OF BOTTOM FRICTION =3


/

---

/ EQUATIONS, BOUNDARY CONDITIONS
/

---

PRESCRIBED FLOWRATES =0.0;0.001

PRESCRIBED ELEVATIONS =0.01;0.0


/

---

/ EQUATIONS, INITIAL CONDITIONS
/

---

INITIAL ELEVATION =0.01

INITIAL CONDITIONS ='COTE CONSTANTE'


/

---

/ INPUT-OUTPUT, FILES
/

---

FORTRAN FILE ='caseone_amanj0.f'

GEOMETRY FILE ='Geometry.slf'

STEERING FILE ='TT2D.cas'

BOUNDARY CONDITIONS FILE ='BOTTOM_BC2.cli'

RESULTS FILE ='TT2dresult.slf'


/

---

/ INPUT-OUTPUT, GRAPHICS AND LISTING
/

---

LISTING PRINTOUT PERIOD =20

GRAPHIC PRINTOUT PERIOD =10

VARIABLES FOR GRAPHIC PRINTOUTS =U,V,B,H,S,F

MASS-BALANCE =true


/

---

/ INPUT-OUTPUT, INFORMATION
/

---

TITLE ='amanj'


/

---

/ NUMERICAL PARAMETERS
/

---

NUMBER OF TIME STEPS =100

NUMBER OF PRIVATE ARRAYS =2

DISCRETIZATIONS IN SPACE =11;11

FINITE VOLUME SCHEME =2

SUPG OPTION =2;2

TIME STEP =0.01

VARIABLE TIME-STEP =true

DURATION =5

TYPE OF ADVECTION =7;5

EQUATIONS ='SAINT-VENANT FV'


/

---

/ NUMERICAL PARAMETERS, SOLVER
/

---

SOLVER OPTION =3

SOLVER =7


/

---

/ NUMERICAL PARAMETERS, VELOCITY-CELERITY-HIGHT
/

---

IMPLICITATION FOR DEPTH =1

IMPLICITATION FOR VELOCITY =0.6

Hello,

Before going further, could you test in finite elements with 3 modified keywords and a new one:

VARIABLE TIME-STEP =.FALSE.
TYPE OF ADVECTION =1;5
EQUATIONS ='SAINT-VENANT EF'
TREATMENT OF THE LINEAR SYSTEM : 2

Your elevation seems very small, do you have dry zones ?

With best regards,

Jean-Michel Hervouet

Hello,

I tried with what you mentioned before, but still i have problem with the boundary inflow at inlet (u/s). Yes, I have dry zone at the crest and downstream slope of the embankment.
/

---

/ TELEMAC2D Version v6p3 10-Dec-2014
/ amanj
/

---

/

---

/ EQUATIONS
/

---

FRICTION COEFFICIENT =92

LAW OF BOTTOM FRICTION =3


/

---

/ EQUATIONS, BOUNDARY CONDITIONS
/

---

PRESCRIBED FLOWRATES =0.0;0.01

PRESCRIBED ELEVATIONS =0.01;0.0


/

---

/ EQUATIONS, INITIAL CONDITIONS
/

---

INITIAL ELEVATION =0.01

INITIAL CONDITIONS ='COTE CONSTANTE'


/

---

/ INPUT-OUTPUT, FILES
/

---

FORTRAN FILE ='caseone_amanj0.f'

GEOMETRY FILE ='Geometry.slf'

STEERING FILE ='TT2D.cas'

BOUNDARY CONDITIONS FILE ='BOTTOM_BC2.cli'

RESULTS FILE ='TT2dresult.slf'


/

---

/ INPUT-OUTPUT, GRAPHICS AND LISTING
/

---

LISTING PRINTOUT PERIOD =20

GRAPHIC PRINTOUT PERIOD =10

VARIABLES FOR GRAPHIC PRINTOUTS =U,V,B,H,S,F

MASS-BALANCE =true


/

---

/ INPUT-OUTPUT, INFORMATION
/

---

TITLE ='amanj'


/

---

/ NUMERICAL PARAMETERS
/

---

NUMBER OF TIME STEPS =100

NUMBER OF PRIVATE ARRAYS =2

DISCRETIZATIONS IN SPACE =11;11

FINITE VOLUME SCHEME =2

TREATMENT OF THE LINEAR SYSTEM =2

SUPG OPTION =2;2

TIME STEP =0.01

DURATION =5

TYPE OF ADVECTION =1;5


/

---

/ NUMERICAL PARAMETERS, SOLVER
/

---

SOLVER OPTION =3

SOLVER =7


/

---

/ NUMERICAL PARAMETERS, VELOCITY-CELERITY-HIGHT
/

---

IMPLICITATION FOR DEPTH =1

IMPLICITATION FOR VELOCITY =0.6

Hello,

If you still use the same coundary conditions file, your boundary 1 is actually the one with prescribed discharge (the first one totally crossed if you follow the file line by line), and your boundary 2 is the one with prescribed elevation. So the values given for PRESCRIBED FLOWRATES and PRESCRIBED ELEVATIONS must be swapped.

Regards,

JMH

Hello dear Jean-Michel Hervouet,

Thank you so much. Now it works.

Kind Regards,
AMANJ

But I have another question. When I have dry zone, does it has solution when using FV scheme?

Sincerely,
AMANJ

Hello,

Yes, the FV options should work also on dry zones.

Rgards,

JMH