open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi, I'd like to create flow in the computational domain without any commonly used boundary condition (e.g., constant discharge, velocity, etc). Say to model "push"ing liquid caused by a propeller inside my computational domain. I do not intend to input or remove any mass.

Is there a way for this in Telemac2D?

Cheers,
Hossein

Hi
As a first simple idea I would say you could try to use the wind but you probably have to create a very local wind...
Otherwise It's possible by adapting the subroutine prosou.f

regards

Hi, Thanks a lot for your help. The application of local wind worked.

I am now facing an issue. I have a closed model meaning that there is no open boundary segment and rather a re-circulation, forced by the application of the local wind.

The problem is that I have no control on the initial depth, and regardless of what I input as initial depth it is always initial elevation + 2! Yes, it is 2!

Any thoughts why this happens and any way around it?

Cheers,
Hossein

please join the file in order to lets us see

Thanks Christophe,

It is attached.


More investigations revealed that the problem has probably something to do with one or more of the FORTRAN files in the "user_fortran" directory. The fortran files are coming from "wind_txy" example of telemac2d and I am using v8p0r2 version of telemac2d.

I appreciate your help in advance.

Cheers,
Hossein

Hi
You shouldn't use strictly the user frotran file for example...
In wind_txy the corfon subroutine exist to impose bathymetry. This is not mandatory in your case if the bathymetry is in your geometry file...

Hi Hossein and Christophe,

I played around a bit with the model and corfon. It seems it needs the next lines and the value needs to be negative, otherwise the water doesn't start moving. These lines in the wind_txy are executed whether smooothing is selected or not, commenting them results in no water movement.

DO I=1,NPOIN
ZF%R(I) = -2.0D0
ENDDO

Hope this helps with the trouble shooting.

Patrick

Thanks Patrick for testing and checking,

You are right. For a reason that I don't know, removal of corfon.f from user_fortran directory won't allow the model to produce meaningful results.

On the three-line code that you wrote above. Yes, that is needed to be included in the corfon.f code, but the issue is that it changes the actual depth of the domain. I have been trying to play with the INITIAL ELEVATION in the cas file to "trick" it to force the right depth, but weirdly when the summation of absolute values of -2.0D0 and my initial elevation equals the true depth, the model doesn't produce any velocity.

Any thoughts?

Cheers,
Hossein

Hi
Your geometry file has a flat bottom at 0.95
Your steering file indicate an initial water elevation of 1.9
Which means 0.95m for water height.

To ensure this settings, you should remove corfon.f from the user fortran directory

If there is no velocities, you should check the wind take in accound in your model. For this you could print the wind values on the result file...

Thanks Christophe,

I removed the confon.f file as you suggested. The issue I am facing now is that when the initial depth is less than 1 m, telemac2d outputs zero for velocity. Once the initial depth is even slightly over 1, it works fine.

I checked the wind that telemac2d inputs and it is fine.

It is odd. Any thoughts?

Thanks,
Hossein