open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everyone,

I'm running a 2D simulation that requires sediment to be recirculated from the outflow back into the inflow (so basically sediment at the outlet needs to be summed and evenly distributed at the inlet). I saw that there was a post on the same topic years ago (www.opentelemac.com/index.php/assistance...diment-recirculation), but the issue didn't seem to be completely addressed.

Hints from the old post suggest modifying conlit.f. In the old post there's also a modified conlit.f provided. I tried modifying that provided conlit.f to suit the model versions I'm using (telemac v8p5r0 with gaia), but there was no luck in making it work.

I wondered if there's an existing conlit.f that takes care of sediment recirculation? If not, I'm trying to modify the code, but another question emerges: is there any resource where I can find the definitions of all the variables? I know there are functions and subroutines explained in developer guide, but I'm having trouble understanding some of the variables in conlit.f (e.g., YADEB, QGL, ...).

Any help will be greatly appreciated!
Jason

Hi Jason,

Have you tried the example case in GAIA ($HOMETEL/examples/gaia) with a recirculating flume and the Wilcock & Crowe equation? The exact name escapes me but it contains "wilcock" in the directory name. It is also documented in the GAIA validation manual.

I have tried adapting this scenario for my own use, but it seems to fail upon reaching 10 classes of sediment.

Hope this helps,
André Renault

Hi Mr. Renault,

Thank you so much for your inputs! I gave it a try but it produced Infinity inflow sediment discharge in my case... (output image attached).

I wondered if you needed to change anything in the conlit_gaia.f file. I changed the number of inflow nodes and classes of sediment, but it didn't work (Infinity outflow). Did you modified anything else for your case?

More about my case (if this helps): there's only 1 class of bedload. In my .cli file for sediment, I set LICBOR=5 at the inflow nodes and LICBOR=4 at the outflow nodes. It ran perfectly fine when not using the user-defined conlit file. All files are attached.

I'll really appreciate any input from anyone!

Many thanks,
Jason

I'd like to reply to my post that the conlit.f file in example of wilcock-crowe-t2d actually works. If anyone wants only one class of sediment, besides changing the lines setting QOUT (near line 49 & 94), you have to change the codetoat line 146.

Also it works pretty well in clusters. It computes the total sediment load in the previous time step and set the sediment input of the current time step to that value.

Jason