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

powerLaw

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 30, 2010, 05:26
Default powerLaw
  #1
New Member
 
Leo Pratama Gani
Join Date: Sep 2010
Posts: 2
Rep Power: 0
lgani is on a distinguished road
Hi everyone.

Thanks in advance for the help.

I am trying to simulate this power law type fluid behaviour in a pipe flow.

However, I am stuck in the transportproperties file where I have to set the properties of the fluid. I found one thread that says

transportModel powerLaw;

nu nu [0 2 -1 0 0 0 0] 1e-05;

powerLawCoeffs
{
k k [0 2 -1 0 0 0 0] 2500;
n n [0 0 0 0 0 0 0] 0.4;
nuMin nuMin [0 2 -1 0 0 0 0] 0;
nuMax nuMax [0 2 -1 0 0 0 0] 100;
}

I am not sure though how to change my unit of k which is in Pa s^1/n to this unit required by openfoam.
Any idea or help with powerLaw fluid flow ?

Cheers
lgani is offline   Reply With Quote

Old   September 30, 2010, 21:04
Default
  #2
New Member
 
Lachlan Graham
Join Date: Dec 2009
Posts: 14
Rep Power: 7
Locky3827 is on a distinguished road
Hi Leo,
OpenFoam uses a kinematic viscosity so you need to divide the k value by the density of your fluid. You will also need to check your pressure results as the pressure is also divided by the density.
You will also need to check the shear rate used by OpenFoam, if you search the forum there is a thread about that as well (don't know how to put the link in from here).
With laminar pipe flow of a power loaw fluid you can check the results against an analytical result.
Let us know how you get on, there only seems to be a few people doing non-Newtonian flow.
Regards,
Locky
Locky3827 is offline   Reply With Quote

Old   October 1, 2010, 02:26
Default
  #3
Member
 
Primoz Ternik
Join Date: Apr 2009
Location: Maribor, Slovenia
Posts: 65
Rep Power: 8
ternik is on a distinguished road
Hi,

your units (normalized with the pressure) for the Power law coefficients are O.K.! It is true that correct unit for the Power law "consistency coefficient" k should be Pa s^(1/n), but I am using the same units and it works!

More important - please note, that the implementation of the shear rate (that is IInd invariant of symmetrical rate of deformation tensor) is not correct (at least in OF-1.5) and it should be:

return sqrt(2.0)*mag(symm(fvc::grad(U_)))

Hope this helps you!

Cheers,
Primoz

Quote:
Originally Posted by lgani View Post
Hi everyone.

Thanks in advance for the help.

I am trying to simulate this power law type fluid behaviour in a pipe flow.

However, I am stuck in the transportproperties file where I have to set the properties of the fluid. I found one thread that says

transportModel powerLaw;

nu nu [0 2 -1 0 0 0 0] 1e-05;

powerLawCoeffs
{
k k [0 2 -1 0 0 0 0] 2500;
n n [0 0 0 0 0 0 0] 0.4;
nuMin nuMin [0 2 -1 0 0 0 0] 0;
nuMax nuMax [0 2 -1 0 0 0 0] 100;
}

I am not sure though how to change my unit of k which is in Pa s^1/n to this unit required by openfoam.
Any idea or help with powerLaw fluid flow ?

Cheers
ternik is offline   Reply With Quote

Old   August 6, 2012, 14:20
Default
  #4
New Member
 
anonymous
Join Date: Jan 2012
Location: Canada
Posts: 24
Rep Power: 5
Mahyar Javidi is on a distinguished road
Hi
Could anybody tell me when we implement power-law coefficients what is the nu in the line before that? What kind of value I have to implement for it while I am using the power-law coefficients

Regards
Mahyar Javidi is offline   Reply With Quote

Old   August 6, 2012, 14:25
Default
  #5
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
hi mahyar
a newbie question :P
when you choose powerlaw, that mentioned nu is trivial, and the program reads the parameters you defined in powerLawCoeffs

P.S: that nu is for Newtonian transport model
nimasam is offline   Reply With Quote

Old   August 6, 2012, 16:46
Default
  #6
New Member
 
anonymous
Join Date: Jan 2012
Location: Canada
Posts: 24
Rep Power: 5
Mahyar Javidi is on a distinguished road
Thanks Nima
Mahyar Javidi 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
Creating a new OpenFOAM case Wesley OpenFOAM 26 May 3, 2011 12:26
Powerlaw for porous medium Shankar FLUENT 0 August 27, 2002 12:54


All times are GMT -4. The time now is 12:45.