[Nek5000-users] Nek5000 and Non-Newtonian Fluids

[nek5000-users at lists.mcs.anl.gov]

Paul,

I switched this to the nek5000-users channel  in case we discuss some thing relevant to others.

For the Bingham model, it seems that all I need is to calculate the second invariant of the strain rate tensor, which is then used to calculate the time-varying apparent viscosity.   I suspect the routine building-blocks for generating this tensor are already in the nek distribution.  If so, could you point me in that direction?

Once those routines are in order, all the calculation will be directed out of uservp, correct?

Thanks.
--Mike


On Jan 14, 2011, at 7:27 AM, Paul Fischer wrote:

> Hi Mike,
> 
> [...]
> The "peris" case uses the stress formulation, as do the "fs"
> cases.  These have the correct SiZE files (resp. for 3D and 2D)
> and flags in the .rea file, save that you would set "IFMVBD"
> to false, since you don't want moving boundaries.
> 
> [...]
> The code was used
> for a lot of non-Newtonian stuff by Lee Ho and Einar, however,
> so it has been exercised significantly in the past.  
> [...]