How to run a wind turbine blade simulation
 

Hi,

First of all, thanks for the great work done on SU2. I'm very happy to have a free solver that runs on windows and doesn't need any complicated compilation with gnu tools (this means I'm running the executables without having compiling capabilities).

I would like to run a Navier-Stokes simulation on a wind turbine blade in a sector of 120° (three bladed turbine).

Supposing I have a valid grid...

1. How do I set the blade to move (rotate) in absolute coordinates?
(to get a moving wall, without moving the whole flow field)

I also read something about SU2_PBC in other posts, it seems to be a preprocessor for grids with periodicity.

2. Is that right?

Thanks,

Gerrit

Hi gerritgroot,
First of all related to the PBC boundary condition, i took the extracts from the thread http://www.cfd-online.com/Forums/su2...-tutorial.html


The procedure to include PBC boundary conditions in your simulation are:

1) Create a regular mesh in .su2 or .cgns format.
2) Use a configuration file which has the option for the periodic marker.

For example,
MARKER_PERIODIC= ( SideWall1, SideWall2, 0,1.0, 0, -60, 0, 0, 0, 0,0, SideWall2, SideWall1, 0,1.0,0, 60, 0,0, 0,0,0 )

where the format of the above quantities is:
% Format: ( periodic marker, donor marker, rotation_center_x, rotation_center_y,
% rotation_center_z, rotation_angle_x-axis, rotation_angle_y-axis,
% rotation_angle_z-axis, translation_x, translation_y, translation_z, ... ).


Use this mesh file and configuration file to run SU2_PBC as:

SU2_PBC configuration.cfg

This will output a new mesh (by default named mesh_out.su2) with a set of halo cells.

Use this mesh file to run a flow simulation. This can be done by just renaming the input mesh file in the configuration file as:
MESH_FILENAME = mesh_out.su2


1. How do I set the blade to move (rotate) in absolute coordinates?
(to get a moving wall, without moving the whole flow field)

% --------------------------- ROTATING FRAME ----------------------------------%
%
% Rotating frame problem (NO, YES)
ROTATING_FRAME= YES
%
% Origin of the axis of rotation
ROTATIONAL_ORIGIN= ( 0.5, -32.0, 0.0 )
%
% Angular velocity vector (rad/s)
ROTATION_RATE= ( 0.0, 0.0, 8.25 )
%
% Reference length (m) for computing force coefficients (e.g. rotor radius)
ROT_RADIUS= 32.0
or
% Reference speed (m/s) for computing force coefficients (e.g. tip speed)
ROT_VEL_REF= 230.0

----------------------------------

For your information, at the moment, only the Euler equations in a rotating frame have been implemented and validated in SU2. For further details, see the following thread:
http://www.cfd-online.com/Forums/su2...-equation.html

hope this helps.

Thanks a lot!!

I'll wait for the NS version and make the mesh in the mean time. :)

Thanks for sharing that information, taxalian.

I just wanted to follow up here and say that we are indeed working on the viscous terms for rotating frame problems (direct and adjoint for design), and this functionality will most likely be available with the next major release of the code, V2.1.

Cheers,
Tom

Hi Gerrit,

Could you please share any information on any success you had with performing wind turbine simulation with SU2? Currently the rotating reference frame formulation has only been applied to the compressible solver and for freestream mach numbers in 0.01, this solver always seems to diverge or not even start in most cases :( Could you please post your config file that is working for you...

Any information is greatly appreciated,
Thanks in advance,
Kind regards,
Kishore

Similar Problem
 

Hello Guys,

I am doing my final year project (Mechanical engineering, be Hons) on ducted wind tubines. I am mainly interested in increasing the power generation (obviously). I have got also some data from an actual turbine built by a company;such as power generated at each wind speed in the venturi.

I managed to run a cfd model on this and got pretty close results to what the company achieved. Using their data (power generated and power coefficient) and my cfd results, I managed to find the pressure difference caused by the turbine.

The study is a 2d study at first which later on i will be doing also a 3d and instead of the turbine i set a porous region.

My question is even if i do get this right, and lets say i change geometry or some other stuff to increase the mass flow rate, how can I know the new pressure difference caused by the turbine? Because, Based on what I think, turbines would obviously give different pressure difference at different rotational speed (that is how it generates different power at different speed).

Pretty much, I am interested to say that the power generated is increased by certain percentage etc.

Any guides ?

Final report on modelling Wind turbine using SU2 and OpenFOAM
 

Hi All,

I've managed to put together the steps and settings for wind turbine simulation using SU2 and OpenFOAM. It also includes some validation for the mexico rotor.

http://publications.tno.nl/publicati...018-R11648.pdf

https://www.slideshare.net/slideshow...n9VMov6E11jeAd

I hope you find it useful,
Kind regards,
Kishore