[vitorspadetoventurin]

I am simulating a non-newtonian flow, I know that openfoam can simulate with Power law model. But how can I define powerlaw coeffs on the transport properties file? Just have coeffs to crosspowerlaw and Carreau model..

transportModel  Newtonian;

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

CrossPowerLawCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
    m               m [ 0 0 1 0 0 0 0 ] 1;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}

BirdCarreauCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
    k               k [ 0 0 1 0 0 0 0 ] 0;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}

thanks

[vitorspadetoventurin]

Please, someone? It is important for my project. Thanks!

[kmooney]

Originally Posted by vitorspadetoventurin

I am simulating a non-newtonian flow, I know that openfoam can simulate with Power law model. But how can I define powerlaw coeffs on the transport properties file? Just have coeffs to crosspowerlaw and Carreau model..

transportModel  Newtonian;

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

CrossPowerLawCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
    m               m [ 0 0 1 0 0 0 0 ] 1;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}

BirdCarreauCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
    k               k [ 0 0 1 0 0 0 0 ] 0;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}

thanks

Hi Vitor,

You can set a power law fluid in a manner similar to this:

transportModel  powerLaw;


powerLawCoeffs
{
    nuMax             nuMax [ 0 2 -1 0 0 0 0 ] 1e-04;
    nuMin           nuMin [ 0 2 -1 0 0 0 0 ] 1e-06;
    k               k [ 0 0 1 0 0 0 0 ] 1;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}

I might be wrong about the dimensions on this.

Cheers,
Kyle

[vitorspadetoventurin]

thanks kmooney, but I have a question.... I paper show me the power law definition:
http://postimg.org/image/54gbo0nf3/

But what would be the nuMin and nuMax? I did not undertand.....


Thank you, If I get success, I will post my results (validation)

Best Regards,
Vitor Spadeto

[shikamaru]

nuMin and nuMax are values that limit the viscosity in case the value for viscosity calculated by the power law equation gives something that may cause an instability.

if nu<nu_min, nu = nu_min and if nu>nu_max, nu = nu_max

This is common to all power-law implementation I have seen in commercial CFD softwares. If you know the limit of your shear rate, then you can calculate them. Otherwise, I guess you take an educated guess (I am sure there is a better way).

Could you please share the reference you posted picture of?

Regards.

[vitorspadetoventurin]

Dear shikamaru, my reference:

http://www.comsol.se/papers/6540/download/Fodor.pdf

thanks!