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

Momentum and displacement thickness: du/dy calculation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 29, 2010, 20:02
Default Momentum and displacement thickness: du/dy calculation
  #1
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
Hi all -

I am calculating the momentum and displacement thicknesses along an airfoil polyline I created. I require the displacement thickness to calculate the momentum thickness. I am using the integral von Karman momentum equation, found here on page 2 (I will be solving for theta, the momentum thickness):
http://www.owlnet.rice.edu/~mech371/...t_integ_BL.pdf

I know how to export the wall shear, density, and velocity data along the polyline to an Excel file. However, I am unsure of how to properly integrate the displacement thickness at each point along the curve. I could create lines at various points along the curve, integrate along those lines to calculate the displacement thicknesses, and then use the integral von Karman momentum equation to find the momentum thicknesses corresponding to those points. However, that would be very time consuming, and would only give me select points rather than a smooth curve.

Alternatively, if there is an easier method of finding the momentum thickness in CFX, I'd love to hear it.

Any help is appreciated.
Josh is offline   Reply With Quote

Old   April 30, 2010, 07:36
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,726
Rep Power: 77
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
This type of analysis is not easily done in CFD-Post. I think your best bet is to do as you describe, but you should be able to script it so that will involve some setup time but once done it should be easy to process data sets.
ghorrocks is offline   Reply With Quote

Old   April 30, 2010, 13:48
Default
  #3
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
As usual, thanks so much for the reply, Glenn.

One problem I am having with creating the lines (to integrate the displacement thickness on) is that I'd like them to be orthogonal to the surface at each point. Since I don't have an equation for the airfoil curvature, I am unsure of how to go about this.

Any ideas, or should I just guess-and-check?
Josh is offline   Reply With Quote

Old   May 2, 2010, 08:17
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,726
Rep Power: 77
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
You might be able to get the surface normal from the normal variable which I think is defined on wall boundaries. Have you tried that?
ghorrocks is offline   Reply With Quote

Old   May 4, 2010, 13:08
Default
  #5
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
I haven't tried that recently, Glenn, but I recall using it at another time with little success. I may try it again.

However, after contacting the authors of the paper, we have discovered that they assumed a constant pressure within the boundary layer, meaning that the momentum thickness calculation is only based on the wall shear, velocity, and density. So much for integrating the displacement thickness.

Thanks for your help, Glenn. I really ought to put you as an author on our research paper.
Josh is offline   Reply With Quote

Old   June 30, 2010, 14:22
Default How to calcualate Momentum thicknees in CFX-post
  #6
New Member
 
Join Date: Jun 2010
Posts: 1
Rep Power: 0
Raj007 is on a distinguished road
Dear Joshua, if you are able to find the momentum thickness in CFX-post then please let me know the processor. waiting for your reply. Thanks

Last edited by Raj007; June 30, 2010 at 14:25. Reason: spelling mistake
Raj007 is offline   Reply With Quote

Old   June 30, 2010, 14:51
Default
  #7
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
Hi Raj -

Unfortunately, I never successfully found a way to calculate momentum/displacement thickness without it being a huge headache. The above methods are the only ideas we've had.

Good luck.
Josh is offline   Reply With Quote

Old   June 30, 2010, 20:20
Default
  #8
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
For the record, here is what I sent to Raj:

I'm not sure what the form factor is. Could you tell me what equation you're using for total pressure loss coefficient is?

As for the momentum thickness...

In Post, I created a polyline to wrap around the airfoil. Then, I exported the wall shear data from this polyline to a .csv (comma separated values) file. The file comes out in this format:

X Y Z Wall Shear
. . . .
. . . .
. . . .

I converted the coordinates (x, y, z) into a curve length, s. I found the difference between each curve value, delta_s. I multiplied each shear value by each delta_s value and summed up the totals to give me the accumulated drag force at each point. I defined momentum thickness using the Blasius similarity solution:

mom. thickness = D(x) / (density * b * U^2)

where D(x) is the drag force in the x direction, b is the span, and U is the freestream velocity. b = 1 for me since my simulations were 2D. Now, I had a format like this:

s Wall Shear delta_s (Wall Shear * delta_s) D(x) mom. thick.
. . . . . .
. . . . . .
. . . . . .

Here's a sample calculation:

s Wall Shear delta_s (Wall Shear * delta_s) D(x)
s1 w1 d1=s2-s1 wd1=w1*d1 D1=wd1
s2 w2 d2=s3-s2 wd2=w2*d2 D2=wd1+wd2.
s3 w3 d3=s4-s3 wd3=w3*d3 D3=wd1+wd2+wd3

momentum thickness
m1=D1/(density * U^2)
m2=D2/(density * U^2)
etc.

Using this method, I predicted very poor values of momentum thickness, so I suspect I did something incorrectly.

Good luck!
Josh is offline   Reply With Quote

Old   June 30, 2010, 20:23
Default
  #9
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
Guess the html format treats multiple spaces as one unit, so the formatting is poor. Sorry about that.
Josh is offline   Reply With Quote

Old   April 10, 2011, 16:09
Default
  #10
New Member
 
Irfan
Join Date: Jan 2011
Location: Netherlands
Posts: 16
Rep Power: 5
iffy is on a distinguished road
Hi,
Does anybody knows how to create lines normal to the airfoil surface in CFX-Post? I have to compute displacement and momentum thickness for an airfoil for which I need lines normal to the airfoil surface?
If anyone here can help?

Regards,

Irfan
iffy is offline   Reply With Quote

Old   April 11, 2011, 03:11
Default
  #11
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 8
Josh is on a distinguished road
Hi iffy -

Unfortunately, it's not easy. For my airfoil, the NACA 0012, I did the following:

Obtain the equation describing your airfoil. If the airfoil does not have one specific equation describing it, subdivide the airfoil into sections of similar slope and create individual equations for each section. y1 = (m)(x1) + b1

Use the equation of the airfoil curve to find the negative, inverse slope (-m^-1) at every point of interest (I had 10 evenly spaced points from x/c = 0 to 1).

Determine how long you want each line to be (I used a length of L = 2 m).

Find the angle between your new, normal lines and the horizontal.

Use L, the angles, and this equation (y2 = (-m^-1)(x2) + b2) to find your new points.

Create your lines.

Hope this helps.
Josh is offline   Reply With Quote

Old   April 11, 2011, 04:01
Default
  #12
New Member
 
Irfan
Join Date: Jan 2011
Location: Netherlands
Posts: 16
Rep Power: 5
iffy is on a distinguished road
Hi Josh,
Thank you for your reply. I am also simulating naca 0012 airfoil. I will try the method which you speicified.

Irfan
iffy is offline   Reply With Quote

Old   October 8, 2013, 06:38
Default
  #13
Member
 
Ye Zhang
Join Date: Dec 2009
Location: Delft,Netherland
Posts: 86
Rep Power: 6
kiddmax is on a distinguished road
Hello, guys.

Did you find good methods to calculate the displacement thickness? How is the method working as Josh mentioned? is it working well?

Best regards,
Ye
kiddmax is offline   Reply With Quote

Old   October 27, 2013, 16:38
Default
  #14
Senior Member
 
Syavash
Join Date: Apr 2010
Location: Tehran, Karaj
Posts: 177
Rep Power: 6
syavash is on a distinguished road
Hello,

I think there might be a way! In Fluent, we create a line normal to the airfoil surface. Then we export the velocity magnitude along that line to a file. Now we should export cell wall distance to another file. By having cell wall distance, we can calculate each cell "dy" via writing a simple script. Now we are able to compute thickness and momentum thickness by integrating the following equations:





Goodluck
syavash 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
general momentum source/hydraulic pressure jump/positive displacement machine ckleanth CFX 10 August 25, 2012 07:23
Displacement and momentum thickness ganesh FLUENT 1 April 15, 2011 11:47
UDF to compute displacement thickness zitzitoun FLUENT 0 October 6, 2009 11:48


All times are GMT -4. The time now is 20:12.