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 × 10 8 . The second option is the Sutherland law, which abides the following expression:(2)
μ = μ 0 ( T T 0 ) 3 2 ( T 0 + S ) T + S
Where,
μ 0
Is the reference viscosity.
T 0
Is the reference temperature.
S MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaam4uaaaa@36E1@
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 ( 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.