Viscosity-temperature Dependence Models

New viscosity-temperature coupling (viscTempCoupling) has been introduced into the nanoFluidX code as an option. Three models were implemented: polynomial, Sutherland and power law.

The polynomial dependency can be specified by five coefficients (4th order polynomial), such that the dynamic viscosity is given by:(1)

μ = C_{0} + C_{1} T + C_{2} T^{2} + C_{3} T^{3} + C_{4} T^{4}

With

T

of course being the temperature of the particle and

C_{n}

being the coefficients. For air the viscosity can be approximated by a linear function, with

C_{1} = 5 \times 10^{- 8}

. The second option is the Sutherland law, which abides the following expression:(2)

μ = \frac{μ_{0} {(\frac{T}{T_{0}})}^{\frac{3}{2}} (T_{0} + S)}{T + S}

Where,

$μ_{0}$: Is the reference viscosity.
$T_{0}$: Is the reference temperature.
$S$: Is the Sutherland coefficient.

For air, these values are: $μ_{0}$ = 1.72 x 10^-5Pas, $T_{0}$ = 273.15 K and S = 110.4.

The last option is the power law, which is defined by:(3)

μ = μ_{0} {(\frac{T}{T_{0}})}^{n}

With $n$ being the exponent. For air, the power law values are: $μ_{0}$ = 1.72 x 10^-5 Pas, $T_{0}$ = 273.15 K and $n$ = 0.66.

First thing that needs to be specified is the viscTempCoupling switch in the Simulation parameters section.

Since the viscosity field is updated after establishing the time step, reference viscosity has to 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-temperature coupling parameter section needs to be defined: see Viscosity-temperature Coupling Parameters.

Important to note is that if the viscTempCoupling is turned on, then all the fluid phase viscosities in the case have to be defined through the Viscosity-temperature coupling parameters section.