
[Sponsors] 
March 14, 2021, 19:02 
expression force_z()@

#1 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
Hello,
I´ve got problems with the calculating of the force of my fan (metro). F= massFlow * velocity areaAve(Velocity w)@Outlet 40,79 m s^1 massFlow()@Outlet 60,07 kg s^1 massFlow()@Outlet * areaAve(Velocity w)@Outlet 2450,37 N force_z()@Outlet  2254,37 N I´m wondering about the difference of 200N. It would be great to know how the expression force_z()@ is defined. There is no hint in de Ansys Post manual, or something else in the internet. Can me pls someone help Thanks eska 

March 15, 2021, 03:26 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,700
Rep Power: 143 
Why are you calculating the force on the outlet? Don't you want the force on the fan blades? So shouldn't you be calculating force_z()@WallBlades?
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

March 15, 2021, 04:46 

#3 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
Hi,
thanks for the answer. I want to calculate the force of the whole fan. So I think I need it at the Outlet? I ve got a value of a metro fan (experimentel) which I have to validate. eska 

March 15, 2021, 04:59 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,700
Rep Power: 143 
The fluid in the outlet does not touch the fan. So how can it apply a force to the fan? It is the fluid next to the fan which applies the force to the fan. So what is the point of measuring the force at the outlet? And what does force at the outlet mean anyway? If you want the force on the whole fan then do force_z()@TheWholeFan. As simple as that. (where "TheWholeFan" is a wall boundary defined as the entire outer surface of the fan.)
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

March 15, 2021, 06:41 

#5 
Senior Member
M
Join Date: Dec 2017
Posts: 642
Rep Power: 12 
Or are you looking for thrust?


March 15, 2021, 06:45 

#6 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
yes, I m looking for the thrust.
In den experimental test they installed a meassuring outside of the fan, to calculate the thrust(/the force?). On some pages I saw the word force as a definition for thrust. 

March 15, 2021, 06:47 

#7  
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
Quote:
My Outlet is at the end of the fan, but still inside. So the expression force_z() is in my case wrong? 

March 15, 2021, 07:14 

#8 
Senior Member
Join Date: Jun 2009
Posts: 1,803
Rep Power: 32 
Please review the equations you are using, keeping in mind the following
integral of ( a * b) =! integral of (a) * integral of (b) and recall massFlow(a)@Boundary is an "integral of", similar to areaInt(b)@Boundary
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

March 15, 2021, 07:15 

#9 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
Im trying to explain it better:
In a company they calculated the thrust of a metro fan, with measuring equipment on the outside of the fan. The fan is standing on a frame with boxes which can measure forces. At least they got the thrust of the hole fan (during the test). Now I have to simulate the fan, and the target is to get the same thrust. In the Internet I saw different ways for calculating the thrust, but I have some works where they work it out with this formula: F(thrust) = massflow * velocity(outlet) = area*density*velocity^2 The massflow is calculated by the same velocity, also at the outlet. This formula works great, I have a thrust with 2450 N, and I need 2221 N. But now I found the expression force_z(). With this one I got 2250 N. Even better then with my formula. But which is now the right one? After ur answer (thank you!!) I think the expression force_z() is wrong. Maybe u can give me a hint where the difference could be. It would be great eska 

March 15, 2021, 07:22 

#10 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
thank you, I will check this now, I m a bit confused, but I found some things in this forum.


March 15, 2021, 09:54 

#11  
Senior Member
Join Date: Jun 2009
Posts: 1,803
Rep Power: 32 
Quote:
The total force on the outlet is by definition an integral, and it cannot be decomposed in the product of the integral of the independent arguments. Can you in the lab measure the independent quantities? Sure you can but it is a lot more work (and money) to do so. You need to measure a distribution of the quantities of interest in the crosssection, and add them up appropriately. Which is better for your final goals? Depends on what you want to rely on. Recall the measurements have an error (this is just on of them), and the simulation has another (truncation error and model error). You must plot them on the same graph with some estimates of their errors (error bars), and you can then start making conclusions. Hope the above helps,
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

March 16, 2021, 04:53 

#12 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
thank you for your great answer!
but there is one more question: what does this expression force_z() mean? can I take it as the thrust or not. where is the difference? 

March 16, 2021, 06:12 

#13 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,700
Rep Power: 143 
If you use force_z() on a wall boundary it gives you the total force acting on the surface, both pressure force and shear force, in the z direction. I don't know exactly what it does if you ask for force_z() on a fluid region (eg an outlet boundary) and it is unlikely to be meaningful.
This is why I suggested you work out the force on the blade wall boundaries, not the fluid.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

March 16, 2021, 06:25 

#14 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
ok, thank you, this helps.
I would like to calculate it on the blade wall, but I have only this one value for the whole fan 

March 16, 2021, 07:51 

#15 
Senior Member
GertJan
Join Date: Oct 2012
Location: Europe
Posts: 1,827
Rep Power: 27 
The fastest route would be to define a Surface Group in CFDPost, put all surfaces of the blades in that group and determine the force on this Surface Group. How difficult can it be?


March 16, 2021, 08:06 

#16  
Senior Member
Join Date: Jun 2009
Posts: 1,803
Rep Power: 32 
Quote:
From the ANSYS CFX Reference Guide: Quote:
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

March 16, 2021, 17:07 

#17 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 

March 17, 2021, 08:45 

#18 
New Member
Sebastian K
Join Date: Jan 2020
Posts: 29
Rep Power: 6 
Thank you for all the helpful answers. Now I want to present you how I calculated the whole force on the fan.
Because of the experimental value(2221N  measured in a company), which is measured with a frame outside of the fan, with force measurement boxes, the formula F(thrust)= massflow * velocity in my eyes is not right. So I tried it with the expression force_z() (z=flow direction). 10*(force_z()@R1 Blade )+ force_z()@R1 Hub+ force_z()@R1 Shroud+ force_z()@M Blades+ force_z()@M Hub + force_z()@M Shroud + force_z()@S1 Blades+ force_z()@S1 Hub+ force_z()@S1 Shroud+ force_z()@S2 Hub+ force_z()@S2 Shroud + force_z()@K1_wall + force_z()@K2_wall I added my 10 blades and all wall boundaries to my expression. Im getting 2237N (2221 needed). For me its a great value, but is this calculation in your eyes right? Would be great to see some oppinions. Thanks 

March 17, 2021, 16:43 

#19 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,700
Rep Power: 143 
I don't know how you defined your wall boundaries and named them, so I cannot confirm your equation is correct. But if I guess what they represent it looks correct. In short the force measured by your load cell holding the entire fan is the sum of the forces acting on each of the fan elements.
So it looks good to me, and if the result is accurate to your expected value then that is even better.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. 

April 8, 2021, 10:06 

#20 
Senior Member
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,167
Rep Power: 23 
"Force_z()@" at non wall surface is as the literature states, pressure and momentum force:
The closest equivalent expression to make this would be Pressure Force: areaAve(Pressure)@UserSurface Momentum Force: areaInt(Density * Velocity w^2)@UserSurface. (slightly different than Mass flow * velocity, which is only for a uniform velocity profile) Opague Quoted: CFDPost calculates the approximate force as follows: If the locator is a wall boundary, the force is equal to the pressure force. I believe it is actually the pressure force plus the wall shear force. So actually force_z is calculated the same everywhere: Pressure + Wall shear + Momentum. but shear is zero on a nonwall, and momentum is zero at a wall. 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Outlet boundary condition in interFoam  Andrea_85  OpenFOAM Running, Solving & CFD  51  July 20, 2017 13:31 
writing execFlowFunctionObjects  immortality  OpenFOAM PostProcessing  30  September 15, 2013 06:16 
How to install CGNS under windows xp?  lzgwhy  Main CFD Forum  1  January 11, 2011 18:44 
CGNS lib and Fortran compiler  manaliac  Main CFD Forum  2  November 29, 2010 06:25 
Lift, Drag Vs time chart,calculations  Jamesd69climber  CFX  8  February 17, 2005 17:23 