open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi I have been using TELEMAC 3D to model the effect of a power plant heated discharge into a lake. The lake has two inflow boundary conditions from 2 rivers, a power plant intake and outfall, and a dam that controls the water level. I was able to run it with the hydrostatic version. However, after checking the results, I am not seeing much stratification in the lake. The temperature difference between the Epilimnion and Hypolimnion is very small.
I am not sure if the heat exchange with the atmosphere is been taken into account correctly. I am using the default equation given in BORD3D subroutine which includes a coefficient A. I am using the METEO subroutine and another subroutine (THERMAL) to input the time variable wind and air temperature conditions. I am attaching the fortran file steering file and output file.

I have also try to run the non-hydrostatic version but it fails after the 2 day of simulation due to a problem at the intake boundary condition. So I still need to see how to fix that.

Any suggestion is appreciated

Thanks

Viviana

Here is the output file

Hi Viviana,

I had previously got problem at the intake boundary, where the intake region goes dry the model run used to fail. In order to avoid the problem the intake quantity has to be reduced as per the available water at the intake.

I had handled this problem by modifying the q function in TELEMAC2D, in TELEMAC3D I think you may have to modify the q3 subroutine.

You may refer to the post
www.opentelemac.org/index.php/kunena/16-...me-dry?start=10#8119
where I have attached the fortran program. The program is meant for telemac2d, in telemac3d you have have to similarly modify the q3 subroutine.

Hope this helps,

With best regards,
Gourish

Hi

Thank you so much for your suggestion, I will try to modify the q3 subroutine based on the program you wrote, and see what happens.

For now I was finally able to fix my problem about the heat exchange with the atmosphere. It seems to be working now.


Thanks again, I really appreciate your help

Viviana

Hi

I have question about the BORD3D subroutine, at the end of the HEAT EXCHANGE WITH THE ATMOSPHERE there are two lines that are commented:

! IMPORTANT:
! STATES THAT ATABOS AND BTABOS ARE NOT ZERO (SEE LIMI3D AND DIFF3D)
! OTHERWISE THEY WILL NOT BE CONSIDERED
! ATABOS%ADR(ITEMP)%P%TYPR='Q'
! BTABOS%ADR(ITEMP)%P%TYPR='Q'

If I am using the default example that you provide for the HEAT EXCHANGE WITH THE ATMOSPHERE, should I uncomment these two lines as well?? Is there other subroutine that I should change in order to get the heat fluxes to the atmosphere? maybe LIMI3D?

I am asking because initially I did a simulation with (commented)
! ATABOS%ADR(ITEMP)%P%TYPR='Q'
! BTABOS%ADR(ITEMP)%P%TYPR='Q'
And the results were giving me much higher values of Temperature than the field measurements. Then I uncommented these two lines and now it looks better but maybe the temperature values are lower compare to the field data.

Thanks so much

Viviana

Hello,

Yes, when the component TYPR='0' we do not consider the terms, this is for optimisation, to avoid adding zeros for nothing.

Regards,

JMH

Thanks so much, I will keep working on this, and check the results for the case with the two lines uncommented in other words with heat exchange taken into account.

Kind Regards

Viviana