typos

ClapeyronThermo · Jan 12, 2024 · eca98e2 · eca98e2 · longemen3000 · Jan 12, 2024
1 parent 344083d
commit eca98e2
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/methods/property_solvers/singlecomponent/saturation/AntoineSat.jl b/src/methods/property_solvers/singlecomponent/saturation/AntoineSat.jl
@@ -101,7 +101,7 @@ function saturation_temperature_impl(model,p,method::AntoineSaturation{TT,VV,CC}
     #one (or two) saturation pressure calculations are normally faster than a crit pure calculation
     (p2,vl2,vv2) = saturation_pressure(model,T2,ChemPotVSaturation(crit_retry = false))
     if !isnan(p2) #nice, psat(T2) exists, we can now produce a really good estimate of the saturation temperature
-        dpdT = (VT_entropy(model,vv2,T2) - VT_entropy(model,vl2,T2))/(vvi - vl2)
+        dpdT = (VT_entropy(model,vv2,T2) - VT_entropy(model,vl2,T2))/(vv2 - vl2)
         dTinvdlnp = -p2/(dpdT*T2*T2)
         Δlnp = log(p/p2)
         Tinv0 = 1/T2

diff --git a/src/methods/property_solvers/singlecomponent/singlecomponent.jl b/src/methods/property_solvers/singlecomponent/singlecomponent.jl
@@ -52,7 +52,7 @@ function μp_equality1_T(model1,model2,v1,v2,p,T,ps,μs)
     p1 = - Av1
     p2 = - Av2
     Fp1 = (p1 - p)*ps
-    Fp2 = (p1 - p)*ps
+    Fp2 = (p2 - p)*ps
     Fμ  = (μ1 - μ2)*μs
     return SVector(Fμ,Fp1,Fp2)
 end