open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everyone !

Flocculation effects are coded in sources/telemac3d via the wcturb.f and soulsbyfloc3d.f (the first is Van Leussen, 1994 and the second is Soulsby 2013). Some keywords are relative to these source codes(FLOCCULATION, FLOCCULATION FORMULA, FLOCCULATION COEFFICIEN and COEFFICIENT RELATIVE TO FLOC DESTRUCTION. You can declare these keywords in your t3d.cas but if you declare it in your sis.cas it won't recognize the parameters. So I firstly wonder why these source codes are in telemac3d folder ?

Secondly, it seems that even if the simulation works with these keywords, they are not taken into account. I tested many configurations and i didn't seen any modification of my results. In fact, by adding a PRINT to the source code, I saw that the PRINT didn't appeared in my file results and it's only the vitchu_sisyphe.f that was taken into account in the simulation. So, do you have an idea of how taking into account flocculation effect ? I am working with cohesive sediments and the Stokes formulation coded in vitchu_sisyphe.f isn't convenient for me.

Many others keywords relative to SISYPHE are coded in sources/telemac3d and in particular the ones relative to the hinder of settling so my questions are wide and not only concern flocculation.

Thank you !

Hi,
I would suggest you to implement FLOCMOD (Verney 2009, you can fine the code on Github) instead as I have done:

www.researchgate.net/publication/3275721...during_a_flood_event

www.researchgate.net/publication/3274415...o-sedimentary_fluxes

FLOCMOD is based on a Population Balance Equation so the main advantage is that you can have a distribution of flocs in your model (which is not possible with the Soulsby paramezation).

Cheers

Thank you flanagan,
Having discussed with you and having reviewed the modules associated with sediment in telemac 3d as well as the different related forums, I understood that the routines developed in TELEMAC-3D concerning sediment transport in 3D and morphology belong to the SEDI-3D module. No documentation is available to date so you will have to check related keywords in telemac3d.dico and related subroutines.
After considering sediment as a tracer, I was able to work with the flocculation (soulsby 2013) and entrapment (whitehouse 2000) modules and I have fairly realistic results in my area with only one class of fine sediments.

Hope it can helps

Hello,

I tried introducing the flocculation keywords in a first test in the t3d cas and in a second test in the gaia cas file and see that the subroutine soulsbyfloc3d (where I put a PRINT in the code) is not used (see attached cas file). Could you please share an example of the cas file that worked for you to model the flocculation with T3D? Thank you in advance.

Hello,

I tried introducing the flocculation keywords in a first test in the t3d cas and in a second test in the gaia cas file and see that the subroutine soulsbyfloc3d (where I put a PRINT in the code) is not used (see attached cas file). Could you please share an example of the cas file that worked for you to model the flocculation with T3D? Thank you in advance.

In the attached cas, I have FLOCCULATION FORMULA = 1 but I have tested with FLOCCULATION FORMULA = 2 and the subroutine soulsbyfloc3d.f is not used in the simulation. I did not find any other helpful information neither in TELEMAC3D user/validation manuals nor in GAIA or SISYPHE user/validation manuals.

Hello,

for GAIA the keyword concerning flocculation are not read. This has been fixed in the develop version and so it will be ok for the next version.
In the meantime, if you want to use flocculation with GAIA, you need to add these lines in lecdon_gaia.f:And to complete, if you want to use different hindered formula, you have to delete line 89 in gaia_init:At the moment, there are no test cases using flocculation. The relative user doc will be updated in the next version (GAIA).
Hope it helps,
Sara

Remember that when coupling with GAIA , keywords for sediments (including flocculation) must be set in the GAIA's steering file.

Best regards,
Sara

Hello,

Thanks a lot for your help.

Best regards,
Aymane