[gerritgroot]

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

[taxalian]

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


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:
rotating problem using NS equation

hope this helps.