diff --git a/include/teqp/models/multifluid.hpp b/include/teqp/models/multifluid.hpp
index d8986c0e226832b8e168e8ee932e75660f96bb84..95f82ec7f8299e0b2dd3c6a77e29649da0be7da1 100644
--- a/include/teqp/models/multifluid.hpp
+++ b/include/teqp/models/multifluid.hpp
@@ -49,7 +49,7 @@ public:
         for (auto i = 0; i < N; ++i) {
             alphar = alphar + molefracs[i] * EOSs[i].alphar(tau, delta);
         }
-        return alphar;
+        return forceeval(alphar);
     }
 
     template<typename TauType, typename DeltaType>
@@ -82,7 +82,7 @@ public:
                 alphar = alphar + molefracs[i] * molefracs[j] * F(i, j) * funcs[i][j].alphar(tau, delta);
             }
         }
-        return alphar;
+        return forceeval(alphar);
     }
 };