# Viscosity Models

nanoFluidX now supports both temperature-viscosity and non-Newtonian modeling behavior. For temperature-viscosity, three models were implemented: polynomial, Sutherland, and power law. For non-Newtonian, the Cross model is available, which can be used to approximate power law behavior without risking instabilities due to viscosity unboundness of the power model.

## Non-Newtonian Viscosity Model (Cross model)

Newtonian fluids, where shear and strain rate have a linear relationship through a single constant called viscosity, are one of the simplest models used to describe fluid behavior. However, some common fluids exhibit more complex behavior where the shearing resistance changes depending on the duration of the shear or the rate of strain applied to the fluid. Some lubricants, many colloidal suspensions, yogurt, ketchup and blood are a few examples of such fluids, collectively called non-Newtonian fluids.

- $\tau $
- Stress tensor
- $\dot{\epsilon}$
- Strain (shear) rate tensor
- ${\mu}_{e}\left(\dot{\epsilon}\right)$
- Effective viscosity as a function of strain rate

In Cross fluids, the effective viscosity ${\mu}_{e}$ varies smoothly from a predefined zero shear rate viscosity ${\mu}_{0}$ in Pa.s to an infinite shear rate viscosity ${\mu}_{\infty}$ in Pa.s. To control the transition between these two ends, a coefficient $k$ [s] and a dimensionless power of $n$ are used in the following form to formalize the relationship with shear rate $\dot{\epsilon}$ in 1/s.

In Cross fluids, shear thickening or shear thinning behavior is achieved by changing ${\mu}_{0}$ and ${\mu}_{\infty}$ while $k$ and $n$ control the curve connecting these two values.

In the absence of an explicit definition for reference viscosity, the code will use the maximum value specified for ${\mu}_{0}$ and ${\mu}_{\infty}$ . This may not be ideal for shear thickening flows as shear rate levels producing the maximum viscosity may not appear in the simulation. It is not recommended to set ${\mu}_{0}$ and/or ${\mu}_{\infty}$ to zero as this may result in inviscid fluid at certain shear rates. In addition, while it may be possible to use negative ${\mu}_{\infty}$ , $k$ or $n$ values without jeopardizing the simulation, these scenarios were not tested and are not recommend as they may result in negative viscosity or viscosity values beyond the interval defined by ${\mu}_{0}$ and ${\mu}_{\infty}$ . Setting $k$ to zero recovers Newtonian behavior with a viscosity value equal to ${\mu}_{0}$ .

## Temperature-viscosity Models

`varvisc_temperature`and

`energy_transport`switches in the Simulation Parameters section must be set to true.

- PolynomialSpecified by five coefficients (4th order polynomial), such that the dynamic viscosity is given by:
(3) Where,$$\mu ={C}_{1}+{C}_{2}T+{C}_{3}{T}^{2}+{C}_{4}{T}^{3}+{C}_{5}{T}^{4}$$- $T$
- Temperature of the particle
- ${C}_{n}$
- Coefficients

For air, the viscosity can be approximated by a linear function with ${C}_{1}=5\times {10}^{-8}$ .

- Sutherland lawAbides by the following expression:
(4) Where,$$\mu =\frac{{\mu}_{0}{\left(\frac{T}{{T}_{0}}\right)}^{\frac{3}{2}}\left({T}_{0}+S\right)}{T+S}$$- ${\mu}_{0}$
- Reference viscosity
- ${T}_{0}$
- Reference temperature
- $S$
- Sutherland temperature (coefficient)

For air, these values are: ${\mu}_{0}$ = 1.72 x 10

^{-5}Pas, ${T}_{0}$ = 273.15 K and S = 110.4. - Power lawDefined by:
(5) $$\mu ={\mu}_{0}{\left(\frac{T}{{T}_{0}}\right)}^{n}$$Where $n$ is the exponent. For air, the power law values are: ${\mu}_{0}$ = 1.72 x 10

^{-5}Pas, ${T}_{0}$ = 273.15 K and $n$ = 0.66.

- Since the viscosity field is updated after establishing the time step,
reference viscosity must be specified in the Domain Parameter section as
`ref_visc`. The reference viscosity should be the highest expected viscosity during the simulation. - After the reference viscosity has been set, the Viscosity Models section must be defined.
- If the
`varvisc_temperature`or`varvisc_non_Newtonian`options are turned on, it is recommended that all of the fluid phase viscosities in the case be defined in the viscosity models section of the configuration file.