Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
TELEMAC-3D
storm and tide problem
Flat Threaded

TOPIC: storm and tide problem

storm and tide problem 6 years 4 months ago #30893

  • mcRove
  • mcRove's Avatar
I wanted to simulate a storm by applying a variable wind on my mesh ... I take into account the tide with tpxo ... it seems to me that the effect of the wind is really zero!

This is the .cas file

/---------------------------------------------------------------------
/ TELEMAC3D Version v7p1 Jul 7, 2018
/ VENTO PORTO DI GENOVA
/---------------------------------------------------------------------



/---------------------------------------------------------------------
/ EQUATIONS
/---------------------------------------------------------------------

/---------------------------------------------------------------------
/ EQUATIONS
/---------------------------------------------------------------------
WIND                                =TRUE

OPTION FOR WIND                     = 2

AIR PRESSURE                        =TRUE

COEFFICIENT OF WIND INFLUENCE       =1.625E-6

HORIZONTAL TURBULENCE MODEL         =3

FRICTION COEFFICIENT FOR THE BOTTOM =35.

VERTICAL TURBULENCE MODEL           =3


/---------------------------------------------------------------------
/ EQUATIONS, INITIAL CONDITIONS
/---------------------------------------------------------------------


INITIAL CONDITIONS ='TPXO SATELLITE ALTIMETRY'



/---------------------------------------------------------------------
/ INPUT-OUTPUT, FILES
/---------------------------------------------------------------------

BINARY DATABASE 2 FOR TIDE ='u_tpxo7.2'

STEERING FILE              ='marea_scirocco.cas'

GEOMETRY FILE              ='geometry_60m.slf'

BINARY DATABASE 1 FOR TIDE ='h_tpxo7.2'

3D RESULT FILE             ='r3d_simulazione_scirocco_t.slf'

2D RESULT FILE             ='r2d_simulazione_scirocco_t.slf'

FORTRAN FILE               ='meteo.f'

BOUNDARY CONDITIONS FILE   ='bc_mesh_60m.cli'

ASCII ATMOSPHERIC DATA FILE     = 'evento_di_scirocco_01_11_2014'

/---------------------------------------------------------------------
/ INPUT-OUTPUT, GRAPHICS AND LISTING
/---------------------------------------------------------------------

VARIABLES FOR 3D GRAPHIC PRINTOUTS =U,V,W,Z,K

LISTING PRINTOUT PERIOD            =20

VARIABLES FOR 2D GRAPHIC PRINTOUTS =U,V,H,B,C,S

GRAPHIC PRINTOUT PERIOD            =20


/---------------------------------------------------------------------
/ INPUT-OUTPUT, INFORMATION
/---------------------------------------------------------------------

TITLE ='VENTO PORTO DI GENOVA'


/---------------------------------------------------------------------
/ NUMERIC PARAMETERS
/---------------------------------------------------------------------

NUMBER OF TIME STEPS        =62880

NUMBER OF HORIZONTAL LEVELS =7


/---------------------------------------------------------------------
/ NUMERICAL PARAMETERS
/---------------------------------------------------------------------

ORIGINAL DATE OF TIME                =2014;11;01

OPTION FOR TIDAL BOUNDARY CONDITIONS =1

TIDAL DATA BASE                      =2

GEOGRAPHIC SYSTEM                    =2

PRESCRIBED ELEVATIONS                =0.0

ZONE NUMBER IN GEOGRAPHIC SYSTEM     =32

TIME STEP                            =15


/---------------------------------------------------------------------
/ PHYSICAL CONSTANTS
/---------------------------------------------------------------------

CORIOLIS COEFFICIENT =119.E-6
The tide is simulated well but lacks the storm... what's wrong?
The administrator has disabled public write access.

storm and tide problem 6 years 4 months ago #30905

  • riadh
  • riadh's Avatar
Hello

Sorry for this late reply

You have to add the subroutine meteo.f to you rfortran file and check if the option 2 for wind is working well or not. This option means that the wind is uniform in the space but varying in time.

I hope that this helps.

with my best regards

Riadh
The administrator has disabled public write access.

storm and tide problem 6 years 2 months ago #31447

  • mcRove
  • mcRove's Avatar
the file fortran meteo.f work well but i can see the action of the wind only if i set WIND COEFFICIENT = 1.62*10-4 and not the right value in the order of 10-6!

why??

thanks a lot

best regards
The administrator has disabled public write access.

storm and tide problem 6 years 2 months ago #31473

  • mcRove
  • mcRove's Avatar
I tried to lengthen the simulation to see how it evolv.

These are the results... the simulation explodes in this way.

The COEFFICIENT OF WIND INFLUENCE was 1.65*10^-6.

suggestions?
Attachments:
The administrator has disabled public write access.

storm and tide problem 6 years 2 months ago #31486

  • Leballeur
  • Leballeur's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 434
  • Thank you received: 163
Hello,

To be sure your wind data are correctly read and interpolated on your mesh during the simulation, can you also output the wind velocity (VARIABLES FOR 2D GRAPHIC PRINTOUTS = X,Y) and plot it at the same location than your plots in the previous post ?

Regards,
Laurent
The administrator has disabled public write access.

storm and tide problem 6 years 2 months ago #31498

  • mcRove
  • mcRove's Avatar
This is the actual state of the simulation:
/---------------------------------------------------------------------
/ TELEMAC3D Version v7p1 Jul 7, 2018
/ VENTO PORTO DI GENOVA
/---------------------------------------------------------------------

/---------------------------------------------------------------------
/ EQUATIONS
/---------------------------------------------------------------------
WIND                                =YES

OPTION FOR WIND                     = 2

AIR PRESSURE                        =true

COEFFICIENT OF WIND INFLUENCE       =2.1268e-06

HORIZONTAL TURBULENCE MODEL         =3

FRICTION COEFFICIENT FOR THE BOTTOM =78.8

FRICTION COEFFICIENT FOR LATERAL SOLID BOUNDARIES =78.8

LAW OF FRICTION ON LATERAL BOUNDARIES =5

VERTICAL TURBULENCE MODEL           =3


/---------------------------------------------------------------------
/ EQUATIONS, INITIAL CONDITIONS
/---------------------------------------------------------------------
INITIAL CONDITIONS ='TPXO SATELLITE ALTIMETRY'


/---------------------------------------------------------------------
/ INPUT-OUTPUT, FILES
/---------------------------------------------------------------------

BINARY DATABASE 2 FOR TIDE ='u_tpxo7.2'

STEERING FILE              ='solo_marea.cas'

GEOMETRY FILE              ='geometry_60m.slf'

BINARY DATABASE 1 FOR TIDE ='h_tpxo7.2'

3D RESULT FILE             ='r3d_marea_e_vento_libeccio_tpxo2.slf'

2D RESULT FILE             ='r2d_marea_e_vento_libeccio_tpxo2.slf'

FORTRAN FILE               ='meteo.f'

BOUNDARY CONDITIONS FILE   ='bc_mesh_60m.cli'

ASCII ATMOSPHERIC DATA FILE     = 'libeccio_30_giorni'

LIQUID BOUNDARIES FILE = 'current.cli'

/---------------------------------------------------------------------
/ INPUT-OUTPUT, GRAPHICS AND LISTING
/---------------------------------------------------------------------

VARIABLES FOR 3D GRAPHIC PRINTOUTS =U,V,W,Z,K

LISTING PRINTOUT PERIOD            =20

VARIABLES FOR 2D GRAPHIC PRINTOUTS =U,V,H,B,C,S,X,Y

GRAPHIC PRINTOUT PERIOD            =20


/---------------------------------------------------------------------
/ INPUT-OUTPUT, INFORMATION
/---------------------------------------------------------------------

TITLE ='VENTO PORTO DI GENOVA'


/---------------------------------------------------------------------
/ NUMERIC PARAMETERS
/---------------------------------------------------------------------

NUMBER OF TIME STEPS        =172560

NUMBER OF HORIZONTAL LEVELS =7


/---------------------------------------------------------------------
/ NUMERICAL PARAMETERS
/---------------------------------------------------------------------

ORIGINAL DATE OF TIME                =2017;06;29

OPTION FOR TIDAL BOUNDARY CONDITIONS =1

TIDAL DATA BASE                      =2

GEOGRAPHIC SYSTEM                    =2

PRESCRIBED ELEVATIONS                =0.0

PRESCRIBED VELOCITIES 				 =1.0

ZONE NUMBER IN GEOGRAPHIC SYSTEM     =32

TIME STEP                            =15


/---------------------------------------------------------------------
/ PHYSICAL CONSTANTS
/---------------------------------------------------------------------

CORIOLIS COEFFICIENT =119.E-6



-------------------------------------------------------------------
And i attached the results...

1) the strange results are disappeared.
2) the wind haven't effect on the set up of the free surface. why?
3) why the result of the free surface it's the same everywhere?

thank you a lot!!

Best regards

R.
Attachments:
The administrator has disabled public write access.

storm and tide problem 6 years 2 months ago #31499

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
Hi

What are the water depth? the vertical position of plane?
Is the result without wind exactly the same for free surface and velocity at the surface?

Regards
Christophe
The administrator has disabled public write access.

storm and tide problem 6 years 2 months ago #31500

  • mcRove
  • mcRove's Avatar
Hi
thank you for the answer.

the water depth in this place is 20 m, more or less.
the layer is the superficial layer.
for the last question i want check for this case for the velocity, but the free surface it's exsactly equal than with only the tide.
it's too clean the signal of free surface for a wind disturb in my opinion.

Best regards
R.
The administrator has disabled public write access.

storm and tide problem 6 years 1 month ago #31588

  • Leballeur
  • Leballeur's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 434
  • Thank you received: 163
Hello,

Did you manage to solve your problem?
If no, can you send us a plot of time series of surface elevation and velocity for the surface layer of the 3D result file by superimposing the 2 simulations (with and without wind) for the 20m-depth point you mention and for another shallower point (about 5m) please?

Regards,
Laurent
The administrator has disabled public write access.

storm and tide problem 6 years 1 month ago #31642

  • mcRove
  • mcRove's Avatar
Hi,
thanks for the reply.
i didn't solve my problem with the model.
This is the stuff that you asked me.
What you think about?

Best regards
The administrator has disabled public write access.
Flat Threaded
Moderators: pham
Forum
System Modules
TELEMAC-3D
storm and tide problem

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.