github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-04-23 03:00:17 -04:00
Code Issues Packages Projects Releases Wiki Activity

Update sympy code for cubic EOS solver rho(T,p)

Browse Source
This commit is contained in:
Ian Bell
2016-06-13 19:29:39 -06:00
parent 63c6c184e5
commit b44d02576b
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/Backends/Cubics/CubicBackend.h
View File

@@ -133,7 +133,7 @@ public:
* \f[ B = \frac{bp}{RT} \f]
*
* Sympy code:
* R,T,v,b,a,Delta_1,Delta_2,p,Z = symbols('R,T,v,b,a,Delta_1,Delta_2,p,Z')
* R,T,v,b,a,Delta_1,Delta_2,p,Z,A,B = symbols('R,T,v,b,a,Delta_1,Delta_2,p,Z,A,B')
* eqn = (R*T/(v-b)-a/(v+Delta_1*b)/(v+Delta_2*b)-p).subs(v,Z*R*T/p)
* eqn2 = eqn*(R*T*Z/p-b)*(R*T*Z/p+Delta_2*b)*(R*T*Z/p+Delta_1*b)/(R**3*T**3/p**2)
* collect(expand(factor(-eqn2)),Z).subs(b*p/(R*T),B).subs(a*p/(R**2*T**2),A)