Welcome, Guest
Username: Password: Remember me

TOPIC: Suspended sediment point sources

Suspended sediment point sources 13 years 8 months ago #1194

  • murphyenda
  • murphyenda's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 26
  • Thank you received: 4
Hello,

Can anyone advise the best way to introduce time-varying point sources of suspended sediment in Sisyphe? Can this be implemented through CONDIM_SUSP, or is this subroutine only called at initialization?

Thanks in advance.

Enda
The administrator has disabled public write access.

Re:Suspended sediment point sources 13 years 8 months ago #1209

  • jmhervouet
  • jmhervouet's Avatar
Hello Enda,

There is nothing programmed yet. We are preparing a link to software DredgeSim of BAW and University of Munich, to add and remove sediment according to a given scenario, but this will be for bedload. So far you can only increase the concentration by a delta(C), for example in subroutine suspension_computation, knowing that it is a depth-averaged concentration (CST%R(I) for a point I), so that the added mass during the time step is delta(C)*depth*(integral of test function which is here V2DPAR%R(I)). This can be done after the call to CVDFTR.
This is however a rough approximation, the correct way to do it would be adding also a quantity of water containing the extra sediment, and activating arrays TEXP (explicit source term, here T11), SMH (source in the continuity equation) and YASMH (then set to .TRUE.) in the call to CVDFTR, but it implies that also the hydrodynamics sees the source. This latter solution is a rather heavy development that will imply transmitting all the data of sources from Telemac-2D to Sisyphe.

With best regards,

Jean-Michel Hervouet
The administrator has disabled public write access.
Moderators: Pablo, pavans

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.