open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I'm modifying source terms in continuity and momentum equations of fluid motion. I'm currently modifying PROSOU subroutine. FU and FV are the for the force terms and SMH is for continuity equation. While the addition of momentum source terms is easy, continuity equation influences the tracer transport equation. If I add something to SMH (eg. 0.0001HN%R in every point of the mesh), that can do damage to the tracer transport equation when I add sources with tracer inflow. Just to make things clear, my addition will be present no matter if there are tracer inflow sources or not.

Is this the proper way to alter source terms for continuity and momentum equations or is there something better?

Regards,
Budo

Hello Budo,

Yes, SMH is the good way to add sources of water, like rain (you can see how rain is treated). SMH is taken into account in the tracer equation, so this should not alter the tracer treatment (see also logical value YASMH). You can also look at punctual sources, which are given through their coordinates, discharge, values of tracers, etc. You can see in DIFSOU how they are taken into account before a call to CVDFTR for advection-difusion. You will see also that depending on the fact that you bring or take water, the corresponding value of tracers change.

I hope this helps,

With best regards,

Jean-Michel

Hello Jean-Michel,

Thank you for pointing out the YASMH variable . I was aware that adding the source term in continuity equation influences the tracer equation. Is the logical variable YASMH used to control only point sources of tracer or is it used in general whenever I try to add a source term in continuity equation?

Regards,
Budo

Hello Budo,

Well the YASMH variable just tells that there is something into SMH, wherever it comes from. Later we added a component in structure BIEF_OBJ : TYPR, that may also tell what is inside an array, but it is not used yet for SMH.

With best regards,

Jean-Michel

Thank you Jean-Michel on your prompt response.

Best regards,
Budo Zindovic