CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

evaluate effective viscosity from U and p already known

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 18, 2012, 04:03
Default evaluate effective viscosity from U and p already known
  #1
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
Hi!

I would like to evaluate an effective viscosity from pressure and velocity field already computed. I mean that I am looking for a field nuEff that obey to:
\nabla \cdot (U U) + \nabla p = \nabla \cdot (\nu_{eff} \nabla U)

In this equation, U and p result from a previous computation.

Does anybody know a method to do such a thing ?

Regards,
Cyp
Cyp is offline   Reply With Quote

Old   January 19, 2012, 04:03
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
If you work out the brackets in the term on the right hand side, can't you just solve it like any other equation in OpenFOAM? So you use fvc instead of fvm on the terms with p and U and solve for nuEff. I can imagine that this will work, but I'll give you now guarantee.
Bernhard is offline   Reply With Quote

Old   January 19, 2012, 04:11
Default
  #3
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
Yes, but in that case, I will have :

T1 = fvc::laplacian(U) which is a vector

T2 = fvc::grad(p)/rho + fvc::div(phi,U) which is a vector as well.

How can I deduce the value of nuEff which should be a scalar (or may be a tensor) ??
Cyp is offline   Reply With Quote

Old   January 19, 2012, 06:53
Default
  #4
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
For the right hand side, don't you get something like

fvm::grad(nuEff) & fvc::grad(U) + nuEff * fvc::laplacian(U)

or something like that, which you can solve if you supply boundary conditions?

But you mean this is an overdetermined system or something? Maybe you can work component-wise?
Bernhard is offline   Reply With Quote

Old   January 19, 2012, 08:31
Default
  #5
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
In fact, what I look for is a method to evaluate an effective viscosity that result from a turbulent calculation (RANS, LES).
In case of RANS simulation, I should recover the value of nut and it will be a test case to validate the approach.

I was thinking that the effective viscosity may satisfy :
\nu_{eff}=(\textbf{U}\textbf{U}+p\textbf{I})\cdot(\nabla U)^{-1}
but I am not sure about that.

I missed the point with your rewritting of the RHS ? if I am right, fvm::grad() has no sense in OF..
Cyp is offline   Reply With Quote

Old   January 19, 2012, 12:22
Default
  #6
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Why don't you just write nut directly then? Seems like a unnecessary exercise to me...
Bernhard is offline   Reply With Quote

Old   January 20, 2012, 04:13
Default
  #7
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
Because the RANS approach are just a particular approach of what I need. It makes sense when LES is considered for example.
Cyp is offline   Reply With Quote

Old   January 20, 2012, 04:41
Default
  #8
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Why? Then you still have an effective viscosity in the way it is implemented in OpenFOAM.
Bernhard is offline   Reply With Quote

Old   January 20, 2012, 04:44
Default
  #9
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
Yes, but what I look for is the effective viscosity that corresponds to time-average velocity and pressure fields. It is not the case of the nuSgs of the LES..
Cyp is offline   Reply With Quote

Old   January 20, 2012, 04:59
Default
  #10
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Why don't you use the time-averaged value of nuSgs/nut? I don't think this will make a very big difference compared to what you want to do. I still don't get why you want to do it this complex...
Bernhard is offline   Reply With Quote

Old   January 20, 2012, 05:22
Default
  #11
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
Because I do not have it and the simulations are too long to be start one more time...
Cyp is offline   Reply With Quote

Old   January 23, 2012, 09:42
Default
  #12
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 227
Rep Power: 9
Cyp is on a distinguished road
I keep working on this problem and I identified some issues.

I need to calculate
Cste*x*\textbf{I}

where I is the identity tensor and x is the position in the x-direction.

How can I access to the value of x ?
Cyp 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
effective viscosity sajad_abasi CFX 3 October 7, 2011 12:40
Effective thermal conductivity and viscosity Rashmi FLUENT 0 October 31, 2007 07:01
effective? turbulent? laminar? viscosity MIssNancy FLUENT 3 December 3, 2002 00:53
effective viscosity in RNG k-e turbulence model Ellen Main CFD Forum 1 August 20, 2002 02:20
effective viscosity in RNG k-e model Ellen FLUENT 0 August 16, 2002 17:46


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