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

Second-Moment Closure in 2D Flow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 18, 2017, 01:17
Default Second-Moment Closure in 2D Flow
  #1
New Member
 
Pawarit
Join Date: May 2017
Posts: 7
Rep Power: 9
pawarit28 is on a distinguished road
Hi,

I was wondering why specifically we require 4 additional transport equations for the Reynolds stresses in 2D flow??

(For the 3D case, it makes perfect sense to take all 6 Reynolds stress components). I understand that turbulence is a 3D phenomenon. Hence, for 2D simulations, we can allow a zero mean W (out of the page) with a fluctuating velocity w'.

But then it raises the question:
Why do we solve for (w'w')?? but not (u'w') or (v'w') then?

Thank you very much everyone!
pawarit28 is offline   Reply With Quote

Old   May 18, 2017, 02:37
Default
  #2
Senior Member
 
Svetlana Tkachenko
Join Date: Oct 2013
Location: Australia, Sydney
Posts: 416
Rep Power: 15
Svetlana is on a distinguished road
Hello, in 3d it is

u'u' u'v' u'w'

v'u' v'v' v'w'

w'u' w'v' w'w'

but w'u' = u'w' and so on so the matrix is symmetric. We are left with 6 unique items which I highlighted in bold. In 2d it is four stresses in total of which only three are unique. The 2d version is underlined.

Last edited by Svetlana; May 18, 2017 at 02:39. Reason: typo in the word 'but'; bold formatting and underline
Svetlana is offline   Reply With Quote

Old   May 18, 2017, 02:51
Default
  #3
New Member
 
Pawarit
Join Date: May 2017
Posts: 7
Rep Power: 9
pawarit28 is on a distinguished road
Hi, thanks for answering!
The Reynolds Stress Models (RSM) in 2D CFD simulations (such as in Fluent) will solve the transport equations for the three unique stresses: u'u', u'v', v'v'
as well as an additional transport equation for: w'w'

This final term in bold is what confuses me. I don't understand why CFD packages solve for w'w' in 2D exactly..Do you have any more insight on this?

(PS. and since w' is indeed non-zero..then why can we ignore u'w' and v'w' in 2D simulations but not w'w'?)
pawarit28 is offline   Reply With Quote

Old   May 18, 2017, 03:01
Default
  #4
Senior Member
 
Svetlana Tkachenko
Join Date: Oct 2013
Location: Australia, Sydney
Posts: 416
Rep Power: 15
Svetlana is on a distinguished road
Interesting. I apologize that I did not catch this from your original question. Where is this a quote from?
Svetlana is offline   Reply With Quote

Old   May 18, 2017, 03:12
Default
  #5
New Member
 
Pawarit
Join Date: May 2017
Posts: 7
Rep Power: 9
pawarit28 is on a distinguished road
Apologies if I was unclear at first.

Here's a brief overview:
https://www.sharcnet.ca/Software/Ans...ss_vs_rst.html
"The alternative approach, embodied in the RSM, is to solve transport equations for each of the terms in the Reynolds stress tensor. An additional scale-determining equation (normally for epsilon or omega) is also required. This means that five additional transport equations are required in 2D flows"
(i.e. 4 equations for reynolds stresses + 1 scale-dissipation)

(if you google "Reynolds stress additional transport equations", all sources/literature agree that four unique stresses are solved for 2D flow. We also agree that u'v' = v'u'. Therefore, we've only counted the three that you first mentioned)

Then when running any 2D Reynolds Stress Model (RSM) simulations in Fluent, I always observed that w'w' was always being solved for. Leading me to believe that this is the fourth and final stress solved for in 2D

Thanks again!
pawarit28 is offline   Reply With Quote

Old   May 18, 2017, 03:39
Default
  #6
Senior Member
 
Svetlana Tkachenko
Join Date: Oct 2013
Location: Australia, Sydney
Posts: 416
Rep Power: 15
Svetlana is on a distinguished road
Perhaps due to lack of focus (although I tried to focus for about 2 minutes), I do not see your quote in the page linked. Please highlight it and take a screenshot, or change the link. Thank you.
Svetlana is offline   Reply With Quote

Old   May 18, 2017, 03:44
Default
  #7
New Member
 
Pawarit
Join Date: May 2017
Posts: 7
Rep Power: 9
pawarit28 is on a distinguished road
"The alternative approach, embodied in the RSM, is to solve transport equations for each of the terms in the Reynolds stress tensor. An additional scale-determining equation (normally for epsilon or omega) is also required. This means that five additional transport equations are required in 2D flows"

these five equations consist of: 4 reynolds stresses + 1 scale-dissipation
pawarit28 is offline   Reply With Quote

Old   May 18, 2017, 04:03
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I have not experience in such a type of modelling but I can suppose that the 2D case is actually to be intepreted as a statistically averaged of a 3D turbulent case. The third component of the fluctuations has therefore an impact on the 2D averaged flow.
Someone more expert in this model can highlight the issue
FMDenaro is offline   Reply With Quote

Old   May 18, 2017, 13:58
Default
  #9
Senior Member
 
Hamid Zoka
Join Date: Nov 2009
Posts: 293
Rep Power: 19
Hamidzoka is on a distinguished road
Hi;
have you checked the calculated w'w' value or compared its magnitude with v'v' and u'u'?
To my understanding if you consider w'w' term in a 2D simulation it means that a fraction of turbulent kinetic energy (K=(u'u'^2+v'v'^2+w'w'^2)^0.5) is distributed in w direction and therefore the portions of u'u' and v'v' goes down. in other words turbulent kinetic energy leaks from the 2D plane of the domain and the true kinetic energy level is underestimated. since in a pure 2D turbulent flow no normal to plane velocity component exists.
This does not make sense in general unless the magnitude of w'w' term is quite small in comparison to v'v' and u'u'.
one more point to add is that second moment closure models are numerically unstable and solving the w'w' term may be an strategy to add more stability to the solution.

Regards
Hamidzoka is offline   Reply With Quote

Old   May 19, 2017, 02:11
Default
  #10
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
For 2D, you need to also solve for w'w' because it is part of k since k = 1/2 (u'u' + v'v' + w'w'). You use k as part of your closure model for the Reynolds stresses so that you actually do need (u'u', v'v', and u'v').

Btw, k appears in the energy equation.

Quote:
Originally Posted by Hamidzoka View Post
In other words turbulent kinetic energy leaks from the 2D plane of the domain and the true kinetic energy level is underestimated. since in a pure 2D turbulent flow no normal to plane velocity component exists.
This does not make sense in general....
This "leakage" property is from where we say that turbulence is 3D. Any turbulence that lacks this characteristic we call non-physical turbulence. It's very physical, with kinematic origins. If you squeeze a tube of toothpaste, the toothpaste comes out the end. You apply an in-plane force but you get an out-of-plane motion!

Last edited by LuckyTran; May 22, 2017 at 20:30. Reason: Earlier I wrote that k appears in the momentum equation, it actually appears in the energy equation.
LuckyTran is offline   Reply With Quote

Old   May 19, 2017, 04:23
Default
  #11
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I agree...as I wrote above, 2D turbulence actually does not exist but in same cases it is possible to solve it in a statistical meaning. But you cannot disregard the contribution of the part of the third component on the motion on the plane.
FMDenaro is offline   Reply With Quote

Reply

Tags
reynolds stress model, reynolds stress tensor

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
mass flow inlet and pressure outlet with target mass flow rate Zigainer FLUENT 13 October 26, 2018 06:58
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
Flow meter Design CD adapco Group Marketing Siemens 3 June 21, 2011 09:33
transform navier-stokes eq. to euler-eq. pxyz Main CFD Forum 37 July 7, 2006 09:42
Plug Flow Franck Main CFD Forum 3 September 4, 2003 06:57


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