CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

non uniform nu model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 7, 2012, 12:19
Default non uniform nu model
  #1
Member
 
Simone
Join Date: Sep 2012
Posts: 95
Rep Power: 4
batta31 is on a distinguished road
Hi guys, I've got a question:
I would like to set up a simulation with a "laminar" turbulence model, but with non uniform values of nu over the domain. Is that possible in openFoam?

Thanks
Simone
batta31 is offline   Reply With Quote

Old   November 8, 2012, 10:52
Default
  #2
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 280
Rep Power: 9
kmooney is on a distinguished road
Do you mean that you want a non-uniform Newtonian viscosity field?
kmooney is offline   Reply With Quote

Old   November 8, 2012, 11:53
Default
  #3
Member
 
Simone
Join Date: Sep 2012
Posts: 95
Rep Power: 4
batta31 is on a distinguished road
Yes, actually I want to use the distribution of nut of a different simulation and apply it to my new case, which has to run as a laminar case but with non uniform viscosity.
batta31 is offline   Reply With Quote

Old   November 8, 2012, 12:18
Default
  #4
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 280
Rep Power: 9
kmooney is on a distinguished road
Quote:
Originally Posted by batta31 View Post
Yes, actually I want to use the distribution of nut of a different simulation and apply it to my new case, which has to run as a laminar case but with non uniform viscosity.
If you can create a volScalarField out of your old nut distrubution perhaps like this:

Code:
volScalarField myNuField = nut;
You should be able to drop it into the diffusion term in most any solver as such:

Code:
        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(myNuField, U)
        );
Cheers!
kmooney is offline   Reply With Quote

Old   November 8, 2012, 12:30
Default
  #5
Member
 
Simone
Join Date: Sep 2012
Posts: 95
Rep Power: 4
batta31 is on a distinguished road
That's a very clever idea! Thank you very much kmooney, I'll try as soon as I can. If it doesn't work I'll come back.

Cheers!
batta31 is offline   Reply With Quote

Old   January 31, 2013, 10:27
Default
  #6
Member
 
Simone
Join Date: Sep 2012
Posts: 95
Rep Power: 4
batta31 is on a distinguished road
Hi again,
I come back to this thread because I've a doubt..

Since I have a non-uniform field of viscosity, besides the term suggested,

fvm::laplacian(myNuField, U)

shouldn't I insert even a term like:

(fvc::grad(U) & fvc::grad(myNuField)


To better understand my reasoning, think you have to derive the NS equation but with non-uniform spatial viscosity..if you then explode the viscous term,

div(nu*grad(U))

you should obtain both the terms above, right?

Cheers
Simone
batta31 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
LES supersonic free jet martyn88 OpenFOAM 22 April 17, 2015 06:00
How to set uniform heating boundary condition? Sargam05 OpenFOAM 0 September 11, 2012 10:09
buoyantSimpleFoam for HVAC application Mojtaba.a OpenFOAM Running, Solving & CFD 5 July 9, 2012 14:35
singularity? mihaipruna OpenFOAM Running, Solving & CFD 5 April 24, 2012 17:18
Need help with boundary conditions: open to atmosphere Wolle OpenFOAM 2 April 11, 2011 07:32


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