open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello

I want to model one breach with sediment transport.
But without Gaia it's ok but with gaia i have topographic anomalies (see join image).

I try to modify lot of variables: solid discharge, classes sediment diameters, SUSPENSION TRANSPORT FORMULA FOR ALL SANDS, and others but all tests are no good.

Do you have an idea ? may be the problem could be that i have the same friction coeffcient on all mesh?

Thank you

Hi,

without discussing on my difficulty , do you think it is possible to model a breach taking into account sediment transport (with gaia)?

Hello Gregory,

SISYPHE/GAIA has been recently used to reproduce breaching processes in a PhD thesis at EDF (title: "Physical and numerical modeling of overtopping induced levee breaching", by Lydia Kheloui).

It's difficult to answer just looking to this figure and the GAIA's steering file (where I can't notice anything strange).
In these cases I usually suggest to simplify your problem, e.g. testing just suspension, then adding bedload or viceversa; looking at the listing if there is something strange, looking at the time step influence, looking at the values of velocities, shear stress in your domain...

By the way, I don't think that the problem is related to the same friction coefficient.

Best regards,
Sara

Thanks you Sara.

In fact, i test MPM with differents shield paramaters.(same results with formules 10 and 30 for "transport for alls sands") .

I continue my tests, simplifying the model even more as you suggest. If I find a solution I will post it here.

good day

Hi,

After lot of tests, i can't find a solution for get a normal breach with gaia. I would like to test with differents friction coefficient. But, the E annexe or "friction_explanation.pdf" no gives example to assign an ID friction at each nodes in friction_user.f. (why??? )

I already create a list with all nodes+id friction but i don't know integrate it in friction_user.f .

Please, does anyone have an example friction_user.f file where they have assigned a friction id for each node?

In fact, my friction_user.f isn't explicit for me (see join file).

Thanking you

Hi Gregory,

it seems that the friction_explanation.pdf is not up to date.
If you have a ZONES FILE with friction ID and node number you don't need to change the friction_user.f. You can directly change the friction values in the FRICTION DATA FILE. See the donau and vegetation validation example of telemac2d. The donau case works without vegetation friction and uses a ZONES FILE. The vegetation example uses also vegetation friction and the friction IDs are written in the geometry file so that a ZONES FILE is not needed anymore.

Best regards,
Rebekka

Hi Rebeka,

A big thank you for your reply. I don't forget to explain a solution for use Gaia wwith Breach if i find a solution.

Good afternoon

Hi,

In conclusion, with v8p2 , i think (beginner user) that it isn't possible to use the Breach and Gaia modules at the same time during a river flood.
Maybe it works with Sisyphe module, or user_confon.f, or some ingenious subroutine?

Cordially

Hi (-:

I test to modifiy the bottom with user_conf.f :

_____________________________________________________________________
USE BIEF
USE DECLARATIONS_TELEMAC2D
!
USE DECLARATIONS_SPECIAL
IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
INTEGER I
! CCV MODIFICATION DU FOND FRONTIERE 2
! ++++++++++++++++++++++++++++++++++++
I = GLOBAL_TO_LOCAL_POINT(313669, MESH)
IF(I.NE.0) ZF%R(I) = 10.50D0
I = GLOBAL_TO_LOCAL_POINT(313671, MESH)
IF(I.NE.0) ZF%R(I) = 10.50D0
.....etc

RETURN
END
______________________________________________________________________

And i add in cas file the repository where:

FICHIER FORTRAN = 'user_fortran' (.../user_fortran/user_corfon.f)

I still have at the beginning of the calculation:

"CORDON (TELEMAC 2D): NO MODIFICATION OF BOTTOM"

How should I proceed to integrate my user_corfon.f in the computation?

Thanks

Hi Gregory,

the message from corfon.f referred to BOTTOM SMOOTHING which you don't have. Nevertheless your user_corfon.f is read. You can do a print message in user_corfon.f in order to be sure that your subroutine is read.
Best regards,
Rebekka