TELEMAC-MASCARET Forum: Adding sediment to transport due to splash erosion. (1/1)

Welcome, Guest

TOPIC: Adding sediment to transport due to splash erosion.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25415

chelobarros

Hi, I am a new user to sishype module on the OpenTELEMAC suite and i have some experience using the telemac2d module.

I need to incorporate erosion due to splash erosion (rain drops removing soils particles).

my formulation for the splash erosion kg/h . m2 requires water depth in each point to account for the damping effect of water in the erosion yield.
So I need to call water depth variable to every node to introduce it in a formula.

Is there anywhere of a sishype subroutines to incorporate mass(in any form) to and individual node ? since splash erosion is proportional to the assigned area to a particular node and water depth, i can not assign a general value of soil mass, but a individual one to each node.

thanks.

Chelo.

The administrator has disabled public write access.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25421

Pablo
OFFLINE
Moderator
Posts: 316
Thank you received: 121

Dear Chelo,
Your problem looks quite similar to one of the topics being studied in our group. Please take a look the the TUC 2015' page 18 [Numerical modeling of sediment transfers at the catchment scale with TELEMAC by Taccone and Antoine, www.opentelemac.org/downloads/Papers%20a...2015-proceedings.pdf].

I think that also there is an article presented in the TUC 2016, but unfortunately the proceedings are not (yet) available for download.

Best wishes,

Pablo

The administrator has disabled public write access.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25451

chelobarros	Thank you. I will revise it and report back. Is it possible to contact them later on? Greetings, Chelo.
	The administrator has disabled public write access.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25703

chelobarros

Hi again Pablo, after some research I have been able to locate the variable FLUER (erosion flux). This is the variable I have added the splash erosion to.

I added g/DT due to splash erosion to this section in the suspension_computation.f subroutine.

!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
! 3.  EROSION FLUX   : FLUER
!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!
!     SKIN FRICTION TAUP  --> T4
!
      CALL OS('X=CYZ   ', X= T4, Y= TOB, Z= MU, C=1.D0)
      CALL OS('X=+(Y,C)', X=T4, Y=T4, C=ZERO)
!
!     SAND ONLY
!
      IF(.NOT.MIXTE) THEN
        IF(.NOT.SEDCO) THEN
          IF (DEBUG > 0) WRITE(LU,*) 'SUSPENSION_EROSION'
          CALL SUSPENSION_EROSION(T4,HN,FDM,FD90,AVA,NPOIN,CHARR,XMVE,
     &                            XMVS,VCE,GRAV,HMIN,XWC,ZERO,
     &                            ZREF,AC,FLUER,CSTAEQ,QS_C,ICQ,U2D,V2D,
     &                            CSRATIO,T14,DEBUG)
          IF (DEBUG > 0) WRITE(LU,*) 'END_SUSPENSION_EROSION'
!
!         NOTE JMH : THIS SHOULD BE INCLUDED IN SUSPENSION_EROSION
!
          DO I=1,NPOIN
            FLUER%R(I)=MIN(FLUER%R(I),ELAY%R(I)*AVA(I)/DT*CSF_SABLE) +
     &      splash_ero
          ENDDO
!

but doing some test I realized that NPOIN (NUMBER OF POINTS IN THE MESH) value is less that the total number of points in my mesh. I was able to associate a number point of the mesh whit an area corresponding to that point in the mesh. but the DO I=1,NPOIN goes to a lesser number that my actual number of points.
This led me to believe that I adding this term to the wrong part of the code or subroutine.
And does anyone know why NPOIN is less than the points in my mesh?.
What units should is use to add the splash erosion to FLUER?

Grettings,
Chelo.

The administrator has disabled public write access.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25704

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	Hi$ You probably run in parallel so NPOIN is the number of point inside the subdomain... Regards
	Christophe The administrator has disabled public write access.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25705

Pablo OFFLINE Moderator Posts: 316 Thank you received: 121	Hello, As in Sisyphe the depth-averaged concentration is expressed in volume, to be consistent with fluer the units of splash erosion should be in m/s. Cheers Pablo
	The administrator has disabled public write access.

Adding sediment to transport due to splash erosion. 7 years 8 months ago #25707

chelobarros

Hi Pablo, thank you for the response.

I was wondering about the units of the FLUER , in some part of the posterior processing this value get multiply by the 'assigned' area of the node? Where can I see this?

Because if not, this would mean that the value of the volume that enters that point is the one that comes from an area of 1 m2, which is significantly less than the real 'assigned' area of each node.

I have already calculated the 'assigned' area of each node whit the intention to multiply them by the value of the splash erosion value(m/s) to account for the entire volume of the soil detachment due to rain. But this would not be concordant whit the required units (m/s != m3/s)

I am sorry if this is a silly confusion on my part.

Greetings,
Chelo.

The administrator has disabled public write access.