CFD Online Discussion Forums

CFD Online Discussion Forums (
-   SU2 (
-   -   How to run a wind turbine blade simulation (

gerritgroot May 13, 2013 09:16

How to run a wind turbine blade simulation

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?



taxalian May 13, 2013 16:15

Hi gerritgroot,
First of all related to the PBC boundary condition, i took the extracts from the thread

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)
% 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)
% 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:

hope this helps.

gerritgroot May 14, 2013 04:25

Thanks a lot!!

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

economon May 21, 2013 03:17

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.


All times are GMT -4. The time now is 16:22.