CFD Online Discussion Forums - Openfoam 10/11

Hi everyone,

I am struggling to setup a non-newtonian fluid within multiphaseInterFoam. It was straightforward before Openfoam 10 was released, were it was setup within transportProperties.

This new method where the non-newtonian model is setup within momentum transport is too confusing.

Can someone please help me figure it out. I found a tutorial implementing a non-newtonian fluid within pimpleFoam. But that also makes no sense as the viscosity model is set to "constant"

Here is what I have, it managed to run but I am not sure if the crossPowerLaw is applied.

MomentumTransport:

Code:

simulationType laminar;



laminar

{

    model               generalisedNewtonian;

    viscosityModel      CrossPowerLaw;

    CrossPowerLawCoeffs

    {

        nu0         [0 2 -1 0 0 0 0]  1-05;

        nuInf       [0 2 -1 0 0 0 0]  1E-12;

        m           [0 0 1 0 0 0 0]   1.00E-03;

        n           [0 0 0 0 0 0 0]   0.7;

    }

}

physicalProperties (Non-newtonian Fluid):

Code:

viscosityModel  constant;



rho           1000;

m             1.00E-03;

n             0.7;

nu            1E-05;

nuInf         1E-12;

I manage to run the simulation just fine with the following output, but I do not get any output from openfoam to show me the non-newtonian viscosity so I am not sure if that is even implemented:

Code:

PIMPLE: Operating solver in transient mode with 6 outer correctors





Reading field p_rgh



Reading field U



Reading/calculating face flux field phi



Selecting viscosity model Newtonian

Selecting viscosity model Newtonian

Selecting viscosity model constant

Selecting turbulence model type laminar

Selecting laminar stress model generalisedNewtonian

Selecting generalised Newtonian model CrossPowerLaw



Reading g



Reading hRef

Calculating field g.h



No MRF models present



No fvModels present

No fvConstraints present

I would appreciate any help, thank you!