open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

With reference to post www.opentelemac.org/index.php/kunena/17-...in-sisyphe-v6p2#6945 , currently sediment source is not implemented in SISYPHE.

However has anyone tried to run a simulation with sediment source which is not at the boundary nodes, (For example dumping of dredged material in water). Is there any other method to accomplish the task, such as coupling with Delft3D modules.

Will the functionality of sediment source be implemented in forthcoming version of TELEMAC.

Thanks,
Gourish

Hello,

So far nothing is done for this at EDF in the next version 6.3, to my knowledge (but HRW also work on Sisyphe). Another approach on the way is the use of DredgeSim. In the while if it is one class of sediment, would it be just a matter of changing the bottom elevation ?

Regards,

JMH

Hi JMH,

Sorry for my very late reply.

Yes currently I working with just one class of sediments, and apparently it is just as you mentioned the matter of changing of bottom elevation. I there any method to achieve it ?

I think this post is also relevant to a very recent question posted at
www.opentelemac.org/index.php/kunena/17-...as-point-source#9323

With Best Regards,
Gourish

Hello,

In sisyphe.f after the lines:

IF(SUSP)...
IF(SLIDE)...
IF(TASS)...

You could put the following lines :

Assuming that you want to add 2 cubic meters of sediment on point 3521

CALL OS('X=0 ',X=T3)
T3%R(3521)=2.D0*UNSV2D%R(3521)
CALL OS('X=X+Y ',X=ZF,Y=T3)
CALL OS('X=X+Y ',X=E,Y=T3)

or more simply if it is only on one point:

ZF%R(3521)=ZF%R(3521)*UNSV2D%R(3521)
E%R(3521)=E%R(3521)*UNSV2D%R(3521)

UNSV2D is 1/V2DPAR, V2DPAR is the area associated to nodes (finite volumes vision), or the integral of test functions (finite element vision).

Not tested ! And maybe my colleagues will have better ideas...

Regards,

Jean-Michel Hervouet

Hello,gourish.

Have you figured out solution for the point source of sediment?

Actually, I have the same problem with you.

I read the subroutine suspension_computation.f and thought the source can be introduced by giving the variable "MASSOU" some value.

I wonder if my thoughs is right or not.

Hope we can have a discussion.

Best wishes

Hi Hong,

I have currently started working on the problem based on Jean-Michel Hervouet's suggestion. I shall post the program if I get the solution.

I shall check the subroutine suspension_computation.f

To be continued...

With Best Regards,
Gourish