[Nidalsb]

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:

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):

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:

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!