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_{1} + C_{2} T + C_{3} T^{2} + C_{4} T^{3} + C_{5} T^{4}

With $T$ 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.

The 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 (Viscosity-Temperature Coupling).

Note: If the viscTempCoupling is turned on, all of the fluid phase viscosities in the case have to be defined through the Viscosity-temperature coupling parameters section.