CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

How to model an internal fan

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 8, 2008, 09:27
Default Dear OpenFOAMers: I am tryi
  #1
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 317
Rep Power: 18
hsieh is on a distinguished road
Dear OpenFOAMers:

I am trying to model a heat transfer problem where inside the fluid domain (air), there is a fan with specified speed.

The internal fan has an inlet side (into the fan) and an outlet side (back into the fluid domain). Both the inlet side and the outlet side will have velocity specified. How should pressure and temperature BCs at the inlet and outlet be applied?

For temperature, I was thinking
at inlet: zeroGradient
at outlet: compute average temperature at the inlet, and, assign this to the outlet

Is this reasonable? How about pressure? zeroGradient for both?

By the way, this is a transient problem.

Can someone shed some light on this? Thanks!

Pei
hsieh is offline   Reply With Quote

Old   July 8, 2008, 10:12
Default Hi Pei, The pressure at a f
  #2
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 77
Rep Power: 17
mathieu is on a distinguished road
Hi Pei,

The pressure at a fan BC is related to the fan performance curve which specifies the pressure difference between both sides of the fan as a function of the flow (or the velocity...) that crosses the fan. For your problem, you should determine the pressure-velocity relation for your fan first. Then you have two choices :

1. You compute the pressure difference with the fan curve and you use it as BC.

2. You implement a BC that compute the velocity depending on the pressure difference (I don't think that a fan BC already exists in OpenFOAM... correct me if I'm wrong).

You should use the first option only if you know that the fan produces a given velocity with a given pressure difference in a given domain (that means that you know the fan behavior in a particular geometry from experimental results for example). You can also use this option if you are sure that the pressure is not an important issue of the problem.

If you really want to know the effect of a specific fan in a specific geometry, then use the option 2.

Hope this helps.

Mathieu

PS : I assume you work with an incompressible flow.

PS2 : For the temperature, I would guess the same as you.
mathieu is offline   Reply With Quote

Old   July 9, 2008, 09:24
Default Thanks Mathieu! I have thin
  #3
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 317
Rep Power: 18
hsieh is on a distinguished road
Thanks Mathieu!

I have think about this more as I am still not quite clear how to implement this in the OF code.

Yes, I am considering incompressible flow.

You refer pressure difference, is it the pressure difference between the fan inlet and fan outlet or?

Pei
hsieh is offline   Reply With Quote

Old   July 9, 2008, 13:02
Default Hi ! Since I am not an exp
  #4
Member
 
Mathieu Olivier
Join Date: Mar 2009
Location: Quebec City, Canada
Posts: 77
Rep Power: 17
mathieu is on a distinguished road
Hi !

Since I am not an expert with C++ and the OF code, I can't help you very much with the programming aspects. However, you can take a look in your "~/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/fields/fvPatchFields/derived" folder for some ideas.

And yes, the pressure difference means the pressure difference between the fan inlet and fan outlet.

Mathieu
mathieu is offline   Reply With Quote

Old   August 24, 2008, 09:02
Default Hello All, I am working on
  #5
Senior Member
 
Mark Couwenberg
Join Date: Mar 2009
Location: Netherlands
Posts: 130
Rep Power: 17
markc is on a distinguished road
Hello All,

I am working on something similar. I am trying to model a ships propeller using the actuator disc theory in order to investigate interaction with nozzle and ships hull.
I modelled the propeller by two discs which are very close to each other.
The picture shows a cut of my model:



The vertical lines represent the propeller disc (left: suction, right: pressure side.
The wing profile represents a cut of a nozzle in which the propeller is positioned.

We more or less know the water velocity through the disc so I apply this velocity. Doing so I am able to run icoFoam. Solutions are however not very accurate yet on summed pressure results. This is due to a lot of reasons (mesh, turbulence, etc).

However I am also eager to run it by applying pressure on the discs (on each disc half the calculated pressure increase, suction minus, pressure side positive) In my undersanding this is physically meaningful. However so far I did not succeed. Velocities continue to grow until solution explodes. I am still trying to find out what I am doing wrong here.
BC's:
p:
>>>
Prop_pres
{
type fixedValue;
value uniform 13.4;
}

Prop_suct
{
type fixedValue;
value uniform -13.4;
}
<<<
U:
>>>
Prop_pres
{
type pressureInletVelocity;
value uniform (0 0 0);
}

Prop_suct
{
type zeroGradient;
}
<<<
To my understanding pressureInletVelocity is specially designed for pressure Inlet BC's. The value (0 0 0) are not intended to be physical values but are only used for calculations. Please anyone correct me if I am wrong.

Concluding:
-Modeling an internal fan (or propeller) can be done by using actuator disc theory and applying velocity BC's
-applying p BC's did not work yet in my case. I am very eager for reactions on what I did wrong here.

Brgds,

Mark
markc is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Commercial meshers] TGridFluent mesh with internal by prism layer and internal face for diagnostic sponiar OpenFOAM Meshing & Mesh Conversion 2 March 30, 2009 16:02
thermal analysis - how to model internal fan? Pei-Ying Hsieh Main CFD Forum 6 March 20, 2008 11:35
how to get fluid internal internal energy vbk FLUENT 0 December 17, 2007 02:54
turbulence model in internal combustion engine George Main CFD Forum 0 December 4, 2006 19:01
How to model internal Wall in ICEM5? Fredrick CFX 4 January 29, 2005 18:01


All times are GMT -4. The time now is 10:13.