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

negative displacement thickness

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

Like Tree2Likes
  • 2 Post By 3dfernando

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 18, 2015, 11:40
Default negative displacement thickness
  #1
Member
 
robo
Join Date: May 2013
Posts: 47
Rep Power: 13
robo is on a distinguished road
Hi all,
I'm working on an incompressible axisymmetric boundary layer code for some bodies of revolution. For some of my results in turbulent flow, I'm getting negative values for the displacement and momentum thicknesses. The rest of the flow field looks correct, and the code appears to be functioning properly in other cases and validated well for flat plate flow and laminar axisymmetric flows. I'm not sure if these negative values are physical, and if they are, how to interpret them. I've done a bit of googling and it seems like they might be indicative of strong streamline convergence leading to an excess of mass and momentum flow in some region or somesuch but I'm not 100% on that. Can these values be negative? What is the best way to interpret that result physically if they can? Finally, are there any good references on this? I'm using Schlichting's and White's books, but haven't been able to find any references in there to negative values for these parameters, at least not so far.
Thanks!
robo is offline   Reply With Quote

Old   December 19, 2015, 05:24
Default
  #2
Senior Member
 
Troy Snyder
Join Date: Jul 2009
Location: Akron, OH
Posts: 219
Rep Power: 18
tas38 is on a distinguished road
If the flow is incompressible, then I believe that the only way for these values to become 'negative' is if the bulk flow direction is reversed.
tas38 is offline   Reply With Quote

Old   December 19, 2015, 08:22
Default
  #3
Member
 
robo
Join Date: May 2013
Posts: 47
Rep Power: 13
robo is on a distinguished road
Thanks for the reply. I've since realised that the displacement and momentum thickness formula's need to be modified for axisymmetric systems and that has resolved my problem.
robo is offline   Reply With Quote

Old   February 22, 2021, 13:46
Default
  #4
New Member
 
Fernando
Join Date: Apr 2013
Posts: 1
Rep Power: 0
3dfernando is on a distinguished road
I know this is an old topic, but I think people could use some actual answer here. As an experimentalist, I've seen momentum thickness (as well as displacement thickness) go negative many times.

Just look at the formulation of displacement thickness. You're integrating:

\theta(x)=\int_0^\infty \frac{u(x,y)}{u_\infty(x)} \bigg(1-\frac{u(x,y)}{u_\infty(x)}\bigg) dy

The term in the parenthesis can go negative if u(x,y)>u_\infty(x). This happens quite often in boundary layers where there is a highly favorable pressure gradient near the surface (but not in the free stream). For example, cases where I encountered negative \theta usually occurred near a sharp edge of a bluff body. Say for example, the Ahmed body or a 3D backwards step with low aspect ratio. The suction under the separation bubble after the step "propagates" upstream to the boundary layer prior separation, accelerating it. Since the free stream does not see this effect, you get a greater velocity close to the surface (in this strange boundary layer), than in the free stream. Then the integral becomes negative.

I see people in the internet say that the momentum thickness becomes negative when the flow separates. I think that is not correct, and I have never personally encountered this particular case in my experiments. If the flow separates, sure, u(x,y) becomes negative in the region of reverse flow, but usually velocity magnitudes at this region are rather small compared to the free stream, making the ratio u(x,y)/u_\infty too small to overcome the momentum defect in the remainder of the separation bubble. In all separation bubbles I've analyzed I've never encountered a negative displacement thickness. It could be a limitation of my experience, though (mostly bluff body separation).

With regards to interpretation, I think the negative momentum thickness has a fairly reasonable one: If a boundary layer has negative momentum thickness because u(x,y)>u_\infty(x), it means it's more energized than the free stream. Thus, it is "highly energized" and will be quite resistant to separation. Obviously, a sharp edge will still separate it.

Thanks,
-FZ
aero_head and Yadira like this.
3dfernando is offline   Reply With Quote

Old   May 23, 2022, 06:29
Default
  #5
Axd
New Member
 
Pauper
Join Date: May 2022
Posts: 1
Rep Power: 0
Axd is on a distinguished road
Recirculation is possible
1) if u changes its direction then the integral of momentum defect will be negative.
2) or if u> U, the momentum thickness will be negative too.
But the chances of the 2) option are quite low (0.00001%)
Reason: 'u' can only energize if the tendency of boundary separation exists.
And we should not forget the pressure gradient along the length of the body.
Axd 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
[Other] mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 13 March 7, 2022 17:22
[snappyHexMesh] determining displacement for added points CFDnewbie147 OpenFOAM Meshing & Mesh Conversion 1 October 22, 2013 09:53
[blockMesh] error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Meshing & Mesh Conversion 2 March 14, 2012 09:56
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 14:11
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 14:00


All times are GMT -4. The time now is 02:55.