CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Modelling a ducted turbine boundaries

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   March 10, 2008, 11:11
Default Hi all, Firstly, please bea
New Member
Thomas Clark
Join Date: Mar 2009
Location: Cambridge, UK
Posts: 10
Rep Power: 8
t_clark is on a distinguished road
Hi all,

Firstly, please bear with me somewhat, as I'm a total newbie to OpenFOAM, fairly new to C++, and fairly new to RANS!

I'm investigating the effect of a duct on an axial turbine. I've built a grid (2D, which I will soon convert to axisymmetric) around my duct shape.

Inside the duct, there is a cylindrical block representing the turbine position. I now need to model the effect of the turbine as either a pressure drop or (preferably) a force and torque distribution.

To make matters more complicated, the forcing varies with the velocity inflow to the turbine position. The function F(Uduct) is known.

My initial strategy is:
1. Get rid of the block. Convert the upstream patch to an outlet, and the downstream patch to an inlet.
2. Constrain the inlet to have the same velocity as the outlet (mass flow conservation).
3. Evaluate U at the outlet, and therefore determine F.
4. Constrain the inlet pressure to be that of the outlet, minus the appropriate amount: (sum(F)/Area).

Firstly, is there anyone who thinks I'm a loony? Is there a better way of doing it (my inexperience with CFD is showing through...)?

Assuming that the above is sensible, my questions are:

1. How do I set up boundary conditions so that their values update (say) every N timesteps?

2. How do I evaluate the velocities and pressures upstream of the turbine, and send them to a piece of my own code in order to evaluate F(Uduct), before I update boundary conditions?

3. Presumably I'll need to compile a new solver - I guess that I'll adapt turbFOAM... Not having a brilliant grasp of how this works yet (I know, RTFM, I'm reading it now) can anyone point me to the right line of turbFOAM where I should begin my tinkering?

Thank you so much in advance for any help!

Tom Clark
t_clark is offline   Reply With Quote

Old   March 11, 2008, 14:57
Default Solved my own problem. Yes,
New Member
Thomas Clark
Join Date: Mar 2009
Location: Cambridge, UK
Posts: 10
Rep Power: 8
t_clark is on a distinguished road
Solved my own problem.

Yes, it turns out that the above is a crazy approach. I'm simply adding a forcing term to the equation being solved, in the block containing the turbine.

I found out how to do this in the thread entitled:

'SIMPLE + force term isn't converging; what have I missed?'

So thanks to Brooks Moses and Michael Prinkey for unwittingly helping out.

Tom C
t_clark is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Flow recirculation around a ducted fan Jane Main CFD Forum 2 January 24, 2012 06:17
Implemeting ducted fan hiraku OpenFOAM Running, Solving & CFD 3 March 19, 2008 15:00
Modelling a turbine blade David Banks Main CFD Forum 2 June 13, 2007 16:27
CFD for Ducted Fans Nash Main CFD Forum 0 August 1, 2003 05:52
Wind Turbine Modelling Neil Campbell Main CFD Forum 1 November 10, 1998 20:03

All times are GMT -4. The time now is 04:28.