Hi everybody,

I’m having some trouble when trying to simulate a power law fluid.

For a newtonian fluid we normally rescale the kinematic viscosity, nu_p, as follows: nu_lb = nu_p / (dx^2/dt).

For a power law fluid, the consistency index, K, has units of [Pa.s^n], where n is the power law index. Thus, knowing that

Pa . s^n = N . s^n/m^2 = kg/m^3 . m^2 /s^(2-n), I arrived to the following relation:

nu_lb = (K/rho_p) * [dt^(2-n)/dx^2]

In this case, for n >= 2 (dilatant fluid) I have enormous values of tau (relaxation factor), which get even bigger when I try to refine the lattice spacing. I read lots of papers, but the authors do not give attention to the way they pass from physical units to lbm units.

Is there anyone who has/had this same problem?

Thanks in advance!

Rodrigo