[titio]

Hi Foamers,

I am trying to model non-newtonian flow and in the model I am being using I have a linear function that depends on the trace of stress tensor field, in the form

func_lin=1+K*trace(StressS)

I am using the PISO algorithm to solve the equations, and I believe to calculate that linear function I need to define a volScalarField, and mapped it in the mesh. I am using the following code to do this

volScalarField func_lin
(
IOobject
(
"func_lin",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
1.0 + K*tr(stressS)
);

K is a scalar constant. Am I right? I think this code initializes the scalar field, that is then recalculated in each iteration.

However, when I compare my predictions with analytical results, I do not get good results. I may be due to the definition of the linear function.

Can someone help me out? Thanks,

Titio

[troosta]

Hi,
I am trying to combine two solvers of openfoam (nonNewtonianIcoFoam anddieselFoam). I need an algorithm to do it.
Roosta