CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > NUMECA

How to best set boundary conditions?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 22, 2012, 11:38
Default How to best set boundary conditions?
  #1
Member
 
Join Date: Jun 2012
Posts: 32
Rep Power: 5
DonQuijote is on a distinguished road
Hello,

I'm working in FINE/turbo with a fan blade geometry, but I'm encountering a lot of problems in finding a converged solution.

I've reviewed all properties but I'm unabled to find a correct convergence, so I hope someone could help me. I think the problem should be related with boundary conditions.

I imposed a total pressure condition at inlet (since a don't have any other value), and a average static pressure outlet. Both of them are constant values since I don't have any kind of profile information (in fact, I don't know the static pressure outlet value). How important is to begin the computation with profiles??

Furthermore, thinking about the physics, if the fan has x rotational speed, it is supposed to obtain an y outlet pressure depending of it, so how can I pre-set this unknow initial condition?(So: known p. inlet, known rotational speed, unk. outlet)

I hope someone could enlighten to me...

Thank you for your help!
DonQuijote is offline   Reply With Quote

Old   October 23, 2012, 00:35
Default
  #2
Senior Member
 
Hamid Zoka
Join Date: Nov 2009
Posts: 164
Rep Power: 8
Hamidzoka is on a distinguished road
Dear Sergio;
it is not neccessary to begin the computation with profiles and it has no negative effect on convergence. Profiles are just applied to find more realistic results specially when a validation with experimental data is required.
I think you should check the boundary conditions. They should be physically meaningful.
When flow passes thnrough the fan its total pressure increases.
So, you should consider this when setting the boundary conditions.
other possibility is the direction of rotation. check if you have correctly set the direction of rotation. It has a big effect on convergence.
what are your boundary conditions values? what is the level of RPM?

regards
Hamidzoka is offline   Reply With Quote

Old   October 23, 2012, 04:29
Default advice
  #3
New Member
 
john
Join Date: Aug 2010
Location: China
Posts: 8
Rep Power: 7
xyz33uu is on a distinguished road
Hello,
Have U ever checked the message box that might show you the problems why U couldn't make a conerged solution. And show me, then maybe I could help.

Best Regards
xyz33uu is offline   Reply With Quote

Old   October 23, 2012, 09:37
Default
  #4
Member
 
Join Date: Jun 2012
Posts: 32
Rep Power: 5
DonQuijote is on a distinguished road
Thank you Hamidzoka and xyz33uu!

I'll explain better my case to show my main purpose and my problem. I'm studyng a specific fan blade (tunnel fans, not aircraft fan), so I have modelled it in autogrid, and now I want to validate the NUMECA results with information I have. After that, I want to change this blade and optimize it for another design point.

So, I have information about massflow, rpm and thrust. Then I set inlet condition as massflow imposed, and the outlet as averaged static pressure. With this configuration, the solution converges, but the axial thrust doesn't reach the real thrust data. The only value I can change is the outlet static pressure (as I unknow).

Thinking that, I was wondering how this parameter is caculated. Because I thought that axial thrust= (difference between inlet-outlet pressures)* A (Ainlet=Aoutlet). With this formula, I need a pressure difference of about 1500 Pa, but I'm unable to reach the mentioned thrust.

So any ideas that could help me?
Thank you very much for your interest!
Regards
DonQuijote is offline   Reply With Quote

Old   October 23, 2012, 22:05
Default
  #5
New Member
 
john
Join Date: Aug 2010
Location: China
Posts: 8
Rep Power: 7
xyz33uu is on a distinguished road
Maybe we have different understanding about thrust. In my opinion, thrust=H(head)*density*gravity*A=((difference between inlet-outlet pressure)+density*(vout^2-vin^2)/2)*A

Last edited by xyz33uu; October 23, 2012 at 23:16.
xyz33uu is offline   Reply With Quote

Old   October 24, 2012, 03:22
Default
  #6
Member
 
Join Date: Jun 2012
Posts: 32
Rep Power: 5
DonQuijote is on a distinguished road
:O I completely missed this part!!! Thank you john!!.

Then , if the thrust is T=massflow/g*V2+A*(P2-P1) (since the fan is static, V1=0), and I know all variables except P2 (P1=Patm), I can solve it, but when I use this parameter in Numeca, the resulting axial thrust is three times lower than the expected... Even when I consider a higher pressure outlet, the axial thrust is half the 'theoretical' one...(but, since the outlet pressure imposed in Numeca is static pressure, it has to be as Pamb.. is that correct?)

What I'm missing???

Last edited by DonQuijote; October 24, 2012 at 04:35.
DonQuijote is offline   Reply With Quote

Old   October 24, 2012, 06:18
Default
  #7
New Member
 
john
Join Date: Aug 2010
Location: China
Posts: 8
Rep Power: 7
xyz33uu is on a distinguished road
Let's think about the unit of (massflow/g*v2),which is (Kg/s)/(m/s^2)*(m/s) that is equal to Kg. And that is different from the unit of thrust. So I think you have got a wrong equation, right?
xyz33uu is offline   Reply With Quote

Old   October 24, 2012, 07:04
Default
  #8
Member
 
Join Date: Jun 2012
Posts: 32
Rep Power: 5
DonQuijote is on a distinguished road
Independently of the ecuation (I use to obtain T in kg, changing for that the pressure units from N/m^2 to kg/m^2), I have a given 'theoretical' thrust (from technical data), and comparing it with the Numeca axial thrust, they are not corresponding (600 N to 1000 N).

But the problem converges quite well. I am sure about the rpm, and I've checked massflow, density, rpm in blade and hub area under the blade...
DonQuijote is offline   Reply With Quote

Old   October 24, 2012, 08:40
Default
  #9
New Member
 
john
Join Date: Aug 2010
Location: China
Posts: 8
Rep Power: 7
xyz33uu is on a distinguished road
I'm so sorry that I havn't notice you want to calculate the axial thrust. And thrust and axial thrust are different. Axial thrust is a component of thrust.
xyz33uu is offline   Reply With Quote

Old   October 24, 2012, 09:10
Default
  #10
Member
 
Join Date: Jun 2012
Posts: 32
Rep Power: 5
DonQuijote is on a distinguished road
But if the fan axis is aligned with the entire fan, it is not suposed to be the same? How would you calculate the component?
DonQuijote is offline   Reply With Quote

Old   October 26, 2012, 06:31
Default
  #11
New Member
 
john
Join Date: Aug 2010
Location: China
Posts: 8
Rep Power: 7
xyz33uu is on a distinguished road
Have you checked your private message box?
xyz33uu is offline   Reply With Quote

Reply

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
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
Boundary conditions for simpleSRFFoam cnsidero OpenFOAM Running, Solving & CFD 4 June 17, 2010 13:37
boundary conditions and mesh exporting vaina74 Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 May 27, 2010 09:38
How to set the changeable boundary conditions. raintung FLUENT 3 December 14, 2004 07:48


All times are GMT -4. The time now is 09:07.