open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear All,

I was wondering, as in Telemac-3D is there a facility in Telemac-2D to program your own reaction term for tracer transport equation. I could not find anything similar to source_trac.f in Telemac-2D.

All help is greatly appreciated.

Best regards,
Sumit

Hello,

The name of the subroutine is difsou.f, and there is also a keyword LAW OF TRACERS DEGRADATION for a first built in law, see the feature logs.

REgards,

JMH

Hello JMH,

Is, "LAW OF TRACERS DEGRADATION" a new addition, I could not find such a keyword anywhere in the reference manual. By the way I am using version v6p1. What are the options for LAW OF TRACER DEGRADATION and what do you mean by feature logs.

Thanks a ton for all your help.

Best regards,
Sumit

Hi Sumit
Yes this is a new functionality of V6P2.
I produce this first developement in order to give the opportunity to new users to have a basic reaction function (often use in france for water qaulity studies) and also to give a first example which should be develop in the future.

When JMH indicate to have a look into the feature logs, he advise you to look on the home page (on the what's new link).

Hope this helps

Thanks Coulet,

I thought so. As you explain, it is water quality issues that I am dealing with hence the search for the implementation of the reaction term. As I have already compiled v6p1 and it has been running successfully, I will first try to program the difsou.f subroutine.

best regards,
Sumit

Coulet,

I am trying to solve advection diffusion reaction equation in 2D. All that I want from reaction term is a first-order decay (-KC). The decay rate is 1 -per day. The code that I have changed in difsou.f is as follows.

DO ITRAC=1,NTRAC
CALL OS('X=0 ',X=TIMP%ADR(ITRAC)%P)
ENDDO
DO I=1,HPROP%DIM1
TIMP%ADR(1)%P%R(I)=-0.0000115D0*HPROP%R(I)
ENDDO
YASMI(1)=.TRUE

Could you kindly take a look and let me know if this is right. I was first trying to use TELEMAC-3D but did not completely succeed in replicating the analytical solution. The detail description of the problem is on the thread
www.opentelemac.org/index.php/kunena/21-...ytical?start=10#7950

Also would like to let you know that I am getting very good results with HEC-RAS. In HEC-RAS the resolution of my successive section is just 20 m, 500 sections in all for 10 KM channel. I wonder is this is the problem. Do I need to have very refined mesh to replicate the analytical solution.

Any and all help is greatly appreciated.

Later on I want to implement some water quality modules in T3D but before I go ahead I wanted to solve simple ADR equation.

Kind regards,
Sumit

Hi

Your code is OK.
It could be more simple because you apply it to all the domain so you could do the same operation in the OS call X=CY instead X=0 with Y=H and C=-1/86400.

Regards

Another point
For the préparation of the validation report of Telemac 6.1, We prepared a test of tracer propagation (with and without diffusion) which i swell fitted with the analytical solution.
I'm not able to see if this validation test case and the validation document is available on the website (but I hope so).

Regards

Hello

Thanks a ton for your kind reply. I am wondering if there is any validation test case for advection-diffusion-reaction equation. I am pretty sure that advection and diffusion of tracer works well, its the reaction component that is creating problems for me.

With kind regards,
Sumit