“…The model has been validated against DNS and experimental correlation data, with simulations of Pr = 0.025 fluids for several geometries, involving plane channels, cylindrical ducts and triangle and square bundle lattices [8,19,20]. With the purpose of improving the numerical stability of the model, a change of state variables has been considered in [21], where specific dissipation rates ω and ω θ have been used in place of ε and ε θ , and in [22], where, by following the work of [23], logarithmic values of ω and ω θ , i.e., Ω and Ω θ , have been used. As a consequence of the change of state variables, some model coefficients had to be changed in order to improve the numerical stability, and the resulting models have been validated with simulations of plane channels and cylindrical ducts for Pr = 0.025.…”