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

Vorticity thickness in temporal mixing layer

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 15, 2019, 12:05
Default Vorticity thickness in temporal mixing layer
  #1
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Hi!

I am struggling with calculating the vorticity thickness in temporal mixing layers. Which is defined as

\delta_\omega=\frac{U_1 - U_2}{\frac{\partial U_0}{\partial y}_{max} }



Where U_1 and U:2 are the lower and upper free stream velocities in opposite directions.
A typical nondimensionalized velocity profile is U/U_1=tanh(y).

To demonstrate the problem:
A simple example is if U_1=1 and U_2=-1. The maximum gradient of tanh(y) is 1. Which yileds a vorticity thickness of 2. But I am supposed to get vorticity thickness of order 1e-4, based on a couple of papers, what am I doing wrong?
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 12:07
Default
  #2
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Vorticity thickness =(U_1-U_2) / (du/dy)_max
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 12:31
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by blackpingu View Post
Hi!

I am struggling with calculating the vorticity thickness in temporal mixing layers. Which is defined as

\delta_\omega=\frac{U_1 - U_2}{\frac{\partial U_0}{\partial y}_{max} }



Where U_1 and U:2 are the lower and upper free stream velocities in opposite directions.
A typical nondimensionalized velocity profile is U/U_1=tanh(y).

To demonstrate the problem:
A simple example is if U_1=1 and U_2=-1. The maximum gradient of tanh(y) is 1. Which yileds a vorticity thickness of 2. But I am supposed to get vorticity thickness of order 1e-4, based on a couple of papers, what am I doing wrong?



Are you working with dimensional or non-dimensional variables? The Reynolds number has the vorticity thickness as characteristic lenght, therefore in non-dimensional variables this thickness is O(1).
FMDenaro is offline   Reply With Quote

Old   February 15, 2019, 12:47
Default
  #4
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Are you working with dimensional or non-dimensional variables? The Reynolds number has the vorticity thickness as characteristic lenght, therefore in non-dimensional variables this thickness is O(1).

Dimensional. For example, one paper claims to have vorticity Reynolds number of 100 with 0.2 mach, which for air is about 70 m/s. which gives a vorticity thickness of 140/70 which is 2. In order to fulfill reynolds number of 100 the viscosity has to be 0.35 which doesn't make any sense.

100=((U_1-U_2)*d)/vis


vis=100/(140*2)=0.35

where d is viscosity thickness.
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 12:50
Default
  #5
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by blackpingu View Post
Dimensional. For example, one paper claims to have vorticity Reynolds number of 100 with 0.2 mach, which for air is about 70 m/s. which gives a vorticity thickness of 140/70 which is 2. In order to fulfill reynolds number of 100 the viscosity has to be 0.35 which doesn't make any sense.

100=((U_1-U_2)*d)/vis


vis=100/(140*2)=0.35

where d is viscosity thickness.



vorticity thickness of 2 meters????
FMDenaro is offline   Reply With Quote

Old   February 15, 2019, 12:59
Default
  #6
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
vorticity thickness of 2 meters????
Yes, I know it doesn't make any sense. There should be something I have misunderstood.
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 13:20
Default
  #7
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by blackpingu View Post
Yes, I know it doesn't make any sense. There should be something I have misunderstood.



If you consider that U1=1m/s and U2=-1m/s and the initial velocity is


u(y)= U1 * tan h (2*y/delta1) the derivative is


du/dy = 2*U1 / (delta1 *cos h (2*y/delta1))


Therefore the delta1 value is the initial condition you have to prescribe.


We did that in Sec.5.2.1 https://www.researchgate.net/publica...-uniform_grids
FMDenaro is offline   Reply With Quote

Old   February 15, 2019, 13:36
Default
  #8
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
If you consider that U1=1m/s and U2=-1m/s and the initial velocity is


u(y)= U1 * tan h (2*y/delta1) the derivative is


du/dy = 2*U1 / (delta1 *cos h (2*y/delta1))


Therefore the delta1 value is the initial condition you have to prescribe.


We did that in Sec.5.2.1 https://www.researchgate.net/publica...-uniform_grids

Thank you, but what if the delta is the vorticity thickness? Which looks like the paper I read did?
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 13:41
Default
  #9
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by blackpingu View Post
Thank you, but what if the delta is the vorticity thickness? Which looks like the paper I read did?

But the vorticity thickness evolves starting from some initial condition... If you are interested in the temporal mixing you have to consider the evolution delta(t) but the thickness at t=0 is what you prescribe from your initial velocity profile
FMDenaro is offline   Reply With Quote

Old   February 15, 2019, 13:47
Default
  #10
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
But the vorticity thickness evolves starting from some initial condition... If you are interested in the temporal mixing you have to consider the evolution delta(t) but the thickness at t=0 is what you prescribe from your initial velocity profile
Ok thanks, so that means I can just reverse the calculations in order to find out what initial delta1 they used ?
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 14:51
Default
  #11
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by blackpingu View Post
Ok thanks, so that means I can just reverse the calculations in order to find out what initial delta1 they used ?

from the Re, prescribing velocity and viscosity
FMDenaro is offline   Reply With Quote

Old   February 15, 2019, 18:22
Default
  #12
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
from the Re, prescribing velocity and viscosity
Problem is, they dont list their viscosity
blackpingu is offline   Reply With Quote

Old   February 15, 2019, 18:28
Default
  #13
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by blackpingu View Post
Problem is, they dont list their viscosity



don't worry about the viscosity value, the key is to use the same Reynolds number...

otherwise, you can simply do the simulation using the non-dimensional equations
FMDenaro is offline   Reply With Quote

Old   February 17, 2019, 11:36
Default
  #14
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 8
blackpingu is on a distinguished road
Thank you
blackpingu is offline   Reply With Quote

Reply


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
[snappyHexMesh] SnappyHexMesh - no layer added bejbro OpenFOAM Meshing & Mesh Conversion 5 February 1, 2020 21:05
[snappyHexMesh] snappyHexMesh Boundary Layer at Corner panpanzhong OpenFOAM Meshing & Mesh Conversion 5 July 3, 2018 06:53
[snappyHexMesh] Problem airfoil snappyhexmesh add layer bye bye my blue OpenFOAM Meshing & Mesh Conversion 1 May 5, 2017 06:56
[snappyHexMesh] SnappyHexMesh: Layer on corners Pj. OpenFOAM Meshing & Mesh Conversion 9 July 17, 2014 05:26
boundary layer thickness 98.4F Main CFD Forum 1 December 29, 2000 08:23


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