open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone, I'm trying to make a thermal plumes simulation, but I can't find a manual or example that allows me to start. I ask if anyone can help me with some source of information or any example that allows me make my simulation.

Thanks and regards.

Jaime

Hello,

To start with, you can look at cases "stratification" and "violet", which are cases with an active temperature.

With best regards,

Jean-Michel Hervouet

Hello sir. Hervouet thanks for your quick response, but I still have a pair of questions that I hope you can answer me:

1) How could I place an initial temperature to the water body that I'm going to study?

2) How could I tell to the program than by some border from the body of water that I'll study, will be discharged a quantity of water with a temperature higher?

I hope you can help me, I really would appreciate it.

Thanks and regards.

Jaime.

Hello,

If you want to initialize the temperature to a constant value, this is done in the cas file. Assuming you have one tracer:
INITIAL VALUES OF TRACERS: 20.

If you need to do a more complex initialization, this is done in the subroutine condim.f of telemac3d. You can replace the lines:
IF(NTRAC.GT.0) THEN
DO I=1,NTRAC
CALL OS( 'X=C ', X=TA%ADR(I)%P, C=TRAC0(I))
ENDDO
ENDIF
with your own initialisation. For example:
IF(NTRAC.GT.0) THEN
PI=ACOS(-1.D0)
L_BOX = 200.D0
DO IPOIN2=1,NPOIN2
DO IPLAN=1,NPLAN
IPOIN3 = (IPLAN-1)*NPOIN2 + IPOIN2
TA%ADR(IND_T)%P%R(IPOIN3) = (1.D0 - 16.D0*PI**2/L_BOX**2)
& * SIN(2.D0*PI/L_BOX*MESH%X%R(IPOIN3))
& * SIN(2.D0*PI/L_BOX*MESH%Y%R(IPOIN3))
ENDDO
ENDDO
ENDIF

To answer you second question, you can set the temperature at each open boundary in the cas file. Assuming you have two open boundaries this gives:
PRESCRIBED TRACERS VALUES : 20.;20.

Note that at outflow boundaries the tracer value is interpolated from inside the domain so this value is overwritten.

Best regards,
Agnès

Best regards,
Agnès