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

nonNewtonianIcoFoam, wired nu

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 10, 2013, 23:19
Default nonNewtonianIcoFoam, wired nu
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Hi FOAMERS,
I want to testify a simple case to try nonnowtonian fluid.
transport dic is this:
Code:
transportModel  CrossPowerLaw;

CrossPowerLawCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 5e-5;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 5e-2;
    m               m [ 0 0 1 0 0 0 0 ] 20;
    n               n [ 0 0 0 0 0 0 0 ] 0.2;
}
Its a kind of shear thinning fluid.

But I got the result is like this:

I think the nu in the bottom cornor should be big,but its not. Whats the problem?

Thanks in advance.

just run nonNewtonianIcoFoam



This problem has been handled. nu0 and nuInf is contrary.

Another issue:

What is m? is it consistency index? But the dimension is s which is different with pa.s^n.
Attached Images
File Type: jpg nu.jpg (30.4 KB, 13 views)
Attached Files
File Type: gz cavity-nonnewtonianicofoam.tar.gz (20.0 KB, 10 views)

Last edited by sharonyue; March 11, 2013 at 04:40.
sharonyue is offline   Reply With Quote

Old   March 12, 2013, 09:06
Default
  #2
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Does anyone know whats the mean of "m"?
sharonyue is offline   Reply With Quote

Old   March 12, 2013, 11:06
Default
  #3
New Member
 
Daniel Duque
Join Date: Jan 2011
Location: ETSIN, Madrid
Posts: 28
Rep Power: 15
dduque is on a distinguished road
I found this in file src/transportModels/incompressible/viscosityModels/CrossPowerLaw:

Code:
Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::CrossPowerLaw::calcNu() const
{
    return (nu0_ - nuInf_)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_;
}
m seems to be some relaxation time in this model. Other than that, I have no idea ...
dduque is offline   Reply With Quote

Old   March 12, 2013, 11:11
Default
  #4
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Quote:
Originally Posted by dduque View Post
I found this in file src/transportModels/incompressible/viscosityModels/CrossPowerLaw:

Code:
Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::CrossPowerLaw::calcNu() const
{
    return (nu0_ - nuInf_)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_;
}
m seems to be some relaxation time in this model. Other than that, I have no idea ...
Hi Daniel;



Thanks all the same. in the Cross viscosity equation, it should be like this. and m is the consistency index which dimension is pa*s^n.
Code:
Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::CrossPowerLaw::calcNu() const
{
    return (nu0_ - nuInf_)/(scalar(1) + m_*pow(strainRate(), n_)) + nuInf_;
}
but in Foam, its just like you post the code. So its wired....
Attached Images
File Type: jpg 1.jpg (5.1 KB, 5 views)

Last edited by sharonyue; March 12, 2013 at 23:02.
sharonyue is offline   Reply With Quote

Old   March 12, 2013, 23:16
Default
  #5
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Yesterday I just run a case.Actually I dont know what does m mean,But I set it like this:I think its like a kind of gum.

Code:
phase1
{
    transportModel  HerschelBulkley;
rho             rho [ 1 -3 0 0 0 0 0 ] 1000;
HerschelBulkleyCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 0.05;
    tau0           tau0 [ 0 2 -2 0 0 0 0 ] 0.01;
    k              k [ 0 2 -1 0 0 0 0 ] 0.2;
    n               n [ 0 0 0 0 0 0 0 ] 0.2;
}
}

phase2
{
    transportModel  Newtonian;
    nu              nu [ 0 2 -1 0 0 0 0 ] 1.48e-05;
    rho             rho [ 1 -3 0 0 0 0 0 ] 1;
   
}

sigma           sigma [ 1 0 -2 0 0 0 0 ] 0.07;

Result of alpha looks okay. in 5 sec, there are some bubbles trapped in the gum.

But nu is not good especially in the final time. When its 5 sec, air is in the top and the gum is in the bottom. but nu is the same all the field.Whats going wrong?

Thanks in advance.
Attached Images
File Type: jpg g.jpg (86.8 KB, 14 views)
File Type: jpg h.jpg (97.5 KB, 15 views)
sharonyue is offline   Reply With Quote

Reply

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
dambreak tutorial's weird velocity field. sharonyue OpenFOAM 8 October 18, 2016 15:57
nonNewtonianIcoFoam: errors near boundaries kpax OpenFOAM Running, Solving & CFD 0 January 10, 2013 10:26
wired error in openfoam conceptone OpenFOAM 4 August 7, 2012 16:45
nonNewtonianIcoFoam wmake T.D. OpenFOAM Running, Solving & CFD 1 November 4, 2010 12:36
nonNewtonianIcoFoam - problem with parallel run chris_sev OpenFOAM Bugs 4 April 1, 2009 10:13


All times are GMT -4. The time now is 05:16.