open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Everyone,

I understand that sediment in SISYPHE can only be introduced at a boundary. However, I wish to introduce sediments into the water column as a point source. I saw some forum posts that this is not yet included in SISYPHE.

Does anyone know of a workaround for the time being? Could an existing subroutine be modified from SISYPHE so that it could accept a rate of sediment injected into the water column at a particular source point (x,y,z)?

I read some posts about DREDGESIM, and saw that there is an interface file in the source code. Could DREDGESIM be used to introduce sediment into the water column as a point source? If so, does anyone have ideas how?

Thanks in advance,

Pat

Hi Pat,

If you are only interested in introducing a suspended sediment source, you can assign a concentration to a particular node (or multiple nodes), i.e. via

CS%ADR(I)%P%R(NPOIN)=...

where I is your sediment size class and NPOIN is your global node number.

Hope this helps,

Enda

Hi Enda,

Thank you for your reply. I will definitely try this.

Can you confirm that this change is to be done in conlit.f subroutine?

Thanks,

Pat

Hi Pat,

To avoid confusing with conditions imposed on boundary nodes, I would just create a separate user-defined subroutine similar to conlit and call it from the sisyphe subroutine close to where conlit is called, so you know it is called at every time step.

Enda

Hello Enda,

I have set up a TELEMAC-2D model coupled with SISYPHE where the entire bed of the channel is non-erodable. I did this by adding the following in the NOEROD.F

The sediment boundary conditions do not add any sediment to the problem, either as bed load or as suspended load. I have also turned off bed load, as I am only interested in suspended sediment (at a source).

Then, I created a new subroutine called PP_SUSP that adds suspended sediment at a single point in the model as:


The PP_SUSP subroutine is called it from SISYPHE.f right after the call to CONLIT (as previously suggested).

I can't figure out why I am not able to see the suspended sediment plume starting at node 326 and moving downstream.

Maybe I am doing something contrary to the rules programmed. All model files are attached.

Thanks in advance,

Pat

Hi Pat and others!

I would also like to define sediment as a point source and started to search in the forum archive.

Pat, I see that your mesh has 606 nodes and that you run on 4 processors. Your implementation should work in scalar but not in parallel since you want to introduce a concentration at node 326 (which must correspond to the global 2D numbering) whereas you can expect to have approx 606 / 4 = 152 nodes per subdomains.

Did you get your model to work in the end?

A more general question: does someone knows if new telemac versions offer a better solution for defining sediment as a point source?

Thank you in advance!

Best regards
PL

Hello PL,

I see what you mean. The stuff in the above post was back from when I was trying to set up a dummy case to get sediment injected into the model as a point source. I abandoned the dummy case, and never got it working.

But, I did try something similar about a year ago (back in v6p3). I am attaching the latest version of subroutine conlit.f that I have. It does work, but I never used it in a real world project. I never had time to fully explore the results I got from the simulations. Note that the node I selected to introduce the sediment happened to be part of two sub-domains in parallel computation. The code attached doesn't check for this, so you have to either check manually, or write some some code to check for this condition.)

Let me know if you advance this further,

Pat