open TELEMAC-MASCARET The mathematically superior suite of solvers

Can someone tell me wether the effect of wave rollers has ever been implemented in TOMAWAC? A roller is the turbulent water mass (white water) travelling with the wave after breaking. It causes a delay in the wave energy dissipation.
It is described by an advection-diffusion equation just like the action balance equation (left-hand side advection, right-hand side source terms).

If anyone attempted the implementation of this effect, please notify me.
Also, does anyone know whether there is a detailed quide to the TOMAWAC Source Code?

Kind regards,

Simon Zenner

Hi Simon,

As far as I know, there is no implementation of the wave roller in Tomawac.
I let the developments team answer you about any developments that I am not aware of.

For the guide of the source code, you can have a look at the doxygen page of the wac subroutine (main routine of Tomawac) and see the call graph for this function:
docs.opentelemac.org/wac_8f.html

I hope it helps.

Regards,
Laurent

Hello

You're right Laurent there are no implementation of the wave roller.

As for the guide of the source code i have nothing to propose more than laurent except to help you when you will have question. Though be aware of the fact that this graph is for version 7P0, it gives you a general idea but some things changed since then.

regards

Thierry

Hello Laurent and Thierry,

Thank you for your quick reply.

I would like to know specifically:
1. in which subroutine in TOMAWAC the action balance equation is solved
2. Which general method is used to slove this balance equation?
3. Where the timestep for the advection terms is defined
4. Where advection terms are computed
5. Where source terms are computed

This would be of tremendous help as I am trying to determine how the roller could be implemented in the program. The roller is defined whith an advection-source/sink balance equation, just like the action balance equation.

Kind regards,

Simon

Hello Simon,

The transport equation is solved by the method of characteristics. So you calculate first the foot of characteristics (general subroutine prepro, to prepare the propagation) and at each time step you call propa to propagate.
Those 2 subroutines are general subroutine.
The time step is simply a keyword ot dictionnary TIME STEP which will be stored in variable DT
The advection field is calculated in conwac if you don't have current and conw4d if you have current.
The source terms are all calculated in routine semimp

Hope it helps

regards

Thierry