optimized crit_pure for ABCubics

ClapeyronThermo · Feb 4, 2022 · 1df75f3 · 1df75f3
1 parent 1c6aba5
commit 1df75f3
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/src/models/cubic/equations.jl b/src/models/cubic/equations.jl
@@ -78,9 +78,11 @@ function p_scale(model::CubicModel,z=SA[1.0])
     return a/ (b^2) # Pc mean
 end
 
-function x0_crit_pure(model::CubicModel)
-    lb_v = lb_volume(model)
-    (1.0, log10(lb_v/0.3))
+function crit_pure(model::ABCubicModel)
+    Tc = model.params.Tc.values[1]
+    Pc = model.params.Pc.values[1]
+    Vc = volume(model,Pc,Tc,SA[1.],phase=:v)
+    return (Tc,Pc,Vc)
 end
 
 function volume(model::ABCubicModel,p,T,z=SA[1.0];phase=:unknown,threaded=false)