open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

Could anyone help me calculate energy extraction in T2D?
I have attempted it by adapting CORSTR.f, and updating the Mannings friction coefficient, based on power uptake proportional to velocity cubed.
I want to simulate multiple turbines from one node.
I'm getting confused when I think about area upon which the turbine(s) act.

Any example code would be much appreciated.

Thanks,
Peter

Hi
You should use (and adapt) the subroutine dragfo.f (example commented in the subroutine) which is best suited for such case.

regards

Thanks for that,

I've now tried adapting DRAGFO.f, but can't seem to generate any changes to the velocities - altering FUDRAG%R(I) and FVDRAG%R(I) do not change anything (I have also modified the quadrilateral zone of interest).

Do I need to modify any other code in the subroutine in order to affect velocities locally?

Thanks,
Peter

Hi
You should activate the function with a keyword (vertical structure or something similar).
You could have a look into the dictionnary

Hope this helps

Thanks!

I am having a similar problem, I am trying to calculate the energy extraction in a defined area in a simple open channel flow. The area contains a defined number of tidal turbines, I have edited the dragfo subroutine to include a calculation for power, in the same format that is used for calculating the drag force, however when I run teh code I get an error that 'A dummy argument name is req'd in this context. [POWER] type(BIEF_OBJ), intent(INOUT) :: POWER'

Does this mean I need to define POWER somewhere else?

I have attached my dragfo subroutine, any help would be much appreciated.

Cheers

Danny