[tareqkh]

Hello Foamers,

I am learning OpenFOAM for my PhD. I will be doing wake effect analysis of the shrouded wind turbine by using actuator disc model. I just found a good paper that explains the implementations of an actuator disk in OpenFoam (See the links below).

http://www.tfd.chalmers.se/~hani/kur...ningReport.pdf

Furthermore, I found the actuator disk example (attached below) by the same author Erik Svenning posted in upper link, and it was made in version 1.5. I am currently using version 2.3.2. Moreover, I have been able to run blockMesh successfully, however, I can't view the mesh in paraFoam. And I get errors when trying to run simpleFoam of the example.
I am asking if anybody of you could find some time to explain what is exactly wrong with version 2.3.2.

Here is "the error"what I got when I ran blockMesh and simpleFoam:

"--> FOAM FATAL IO ERROR: 
Cannot find patchField entry for cyclic fan_half0
Is your field uptodate with split cyclics?
Run foamUpgradeCyclics to convert mesh and fields to split cyclics.

file: /home/tariq/windturbine/project3/erikSvenningFiles/fan/0/p.boundaryField from line 34 to line 59.

From function GeometricField<Type, PatchField, GeoMesh>::GeometricBoundaryField::readField(const DimensionedField<Type, GeoMesh>&, const dictionary&)
in file /home/tariq/OpenFOAM/OpenFOAM-2.3.1/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 192.

FOAM exiting
"

Thank in advance
fan.zip

[eruwaedhiel]

This might not be much help, but I was trying to get that code working a few weeks ago for my own project, and came across this: http://www.cfd-online.com/Forums/ope...f2-1-help.html, which basically says not to try using that code.

What I wound up doing in the end was making a pressure jump baffle for the actuator disk following these tutorial: https://www.youtube.com/watch?v=7Ex1RHDQqIc and https://www.youtube.com/watch?v=qTYAr8jgViM

I'm currently trying to get it to converge with a helicopter fuselage mesh in simpleFoam, but the actuator disk portion at least seems to be doing what I expect it to.

[tareqkh]

Hello eruwaedhiel.

Thank you for your response. Is there any way to start the same example from scratch? I feel that it would be perfect.

Thanks in advance.

[eruwaedhiel]

Hi tareqkh,

Sorry for the delay, my computer died last week and I just got it back running yesterday. The two tutorials I posted actually provide all input files associated with them in the descriptions, so you can probably start from that (that's what I did).

[AnasCFD]

In your case you are dealing with "cyclic" boundary conditions. This means that you have, after running the blockMesh, two faces named as "fan_half0" and "fan_half1". So you should define them in the 0 directory. If the two faces have the same boundary condition you can set, for example, in the 0/p file "fan.+" instead of "fan".

Enjoy

[tareqkh]

Hello Anas and Eruwaedhiel,

Thank you for your time clarifying my confusion in the first part of the problem. I have tried to follow the steps about creating the solver in order to solve both thrust forces and swirl in the same project. I have successfully implemented the solver, however, for some reason I am not getting any velocities in Ux, Uy, and Uz. How can I make sure that the implementations are made correctly. I am really stuck at this solver.

Regards,