open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear all

I'm trying to simulate thermal plume using TELEMAC-3D.

My conditons are belows.
Intake and Outfall Q = 30 CMS
Excess temperature = 10 deg.C

So, There is information of my steering file.
NUMBER OF TRACERS = 1
NAMES OF TRACERS : 'TEMPERATURE'
INITIAL VALUES OF TRACERS = 5.0

ABSCISSAE OF SOURCES : 80.0; 420.0
ORDINATES OF SOURCES : 40.0; 40.0
ELEVATIONS OF SOURCES : 0.0; 0.0
WATER DISCHARGE OF SOURCES : 30.0; -30.0
VALUE OF THE TRACERS AT THE SOURCES : 10.0; 10.0

However, I have no idea how to calculate excess temerature and re-circulation.
Could you explain how to simulate thermal discharge, pleae?

Thank you.
Cho.

Hi,

You have to put in our fortran file (inside the time loop of telemac3d.f for example) the lines below:

ITRAC = 1 / the temperature

TA_SCE2%ADR(ITRAC)%P%R(1) = TA%ADR(ITRAC)%P%R(ISCE(2)+(KSCE(2)-1)*NPOIN2) + 10.D0 / For DT = 10°C, the first source the outfull & the 2nd the intake

Regards,
Olivier

Thank you, Olivier

It was very useful to me your answer.
But I need more deail about your description.
(source code of telemac3d.f replace? or modify? or add where?)

The only way is source code modifying?
Is there any other way to simulate thermal discharge?

Regards,
Cho.

Hi Cho,

You have to put in your steering file the key word :
FORTRAN FILE = fortran.f
Copy telemac3d.f in your working directory and rename it in fortran.f

You can simulate thermal discharge without code modifying but in this case without recirculation!

Regards,
Olivier

Hi Everyone,

So I am using telemac3d for a thermal plume simulation and I have a small problem.

For my project I have 3 intake nodes. In the fortran code in the BORD3D Subroutine I implement the Temperature at the outfall as the average of the three temperatures at the intake and add the deltaT (which is 5°C) : Tout = ((T1+T2+T3)/3)+5°C (you can find the code below).

But at the end of the model, I notice that the temperature at the outfall is not 5°C but less (see figure below)

Does anyone know why I have this strange results?

Thank you for the help

Mehdi Dadou

CODE :


IF(IFRLIQ.EQ.2) THEN
TABORL%ADR(ITRAC)%P%R(IBORD) =
& ((TA%ADR(1)%P%R(18651)+TA%ADR(1)%P%R(18679)
& +TA%ADR(1)%P%R(18707))/3
& +5.0D0)
ENDIF

Hi Mehdi,

If you use the sources, the intake is not a boundary node and you haven't to modify the bord3d surbroutine.
If you want to use a boundary node, you have also the key word 'TREATMENT OF FLUXES AT THE BOUNDARIES'.
The outfall temperature can be also mixed with the surrounding temperature inside the mesh.

Regards,
Olivier

Thank you Olivier,

One last thing : I want to know how does the code calculate this? I mean is there a subroutine/function that explain how it calculates the temperature according to the mixing inside the mesh?

Thank you again

B. regards,

Mehdi Dadou

No subroutine!
The temperature inside de mesh is calculated by the thermal equation.

Regards,
Olivier

Ah Great!

Thank you

Kind Regards,

Mehdi Dadou