diff --git a/src/thermophysicalModels/thermophysicalFunctions/NSRDSfunctions/NSRDSfunc14/NSRDSfunc14.H b/src/thermophysicalModels/thermophysicalFunctions/NSRDSfunctions/NSRDSfunc14/NSRDSfunc14.H
index 2b4a7cea96d5a72ae9fbc8727f7c00f1a108bf5f..20876115ac69c7961f119210ebf02abb22a2ed80 100644
--- a/src/thermophysicalModels/thermophysicalFunctions/NSRDSfunctions/NSRDSfunc14/NSRDSfunc14.H
+++ b/src/thermophysicalModels/thermophysicalFunctions/NSRDSfunctions/NSRDSfunc14/NSRDSfunc14.H
@@ -105,9 +105,11 @@ public:
         //- Evaluate the function and return the result
         scalar f(scalar, scalar T) const
         {
-            scalar t = 1.0 - T/Tc_;
+            scalar Tdash = min(T, Tc_ - ROOTVSMALL);
+
+            scalar t = 1.0 - Tdash/Tc_;
             return
-                a_*a_/t + b_ - t
+                a_*a_/(t + ROOTVSMALL) + b_ - t
                *(
                    2.0*a_*c_
                  + t*(a_*d_ + t*(c_*c_/3.0 + t*(0.5*c_*d_ + 0.2*d_*d_*t)))