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

Fix tau derivative of n2Aijk

Browse Source
This commit is contained in:
Ian Bell
2015-08-12 23:42:46 -06:00
parent ad643aff7e
commit 14b1bb829e
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/Backends/Helmholtz/MixtureDerivatives.h
View File

@@ -194,8 +194,8 @@ class MixtureDerivatives{
}
static CoolPropDbl d_n2Aijk_dX(HelmholtzEOSMixtureBackend &HEOS, std::size_t i, std::size_t j, std::size_t k, x_N_dependency_flag xN_flag, parameters WRT){
CoolPropDbl RT = HEOS.gas_constant()*HEOS.T();
CoolPropDbl summer = -d_nAij_dX(HEOS, i, j, xN_flag, WRT);
if (WRT == iTau){
double summer = 0;
summer += d2_ndln_fugacity_i_dnj_dtau2__constdelta_x(HEOS, i, j, xN_flag)*ndtaudni__constT_V_nj(HEOS, k, xN_flag);
summer += d_ndln_fugacity_i_dnj_dtau__constdelta_x(HEOS, i, j, xN_flag)*d_ndtaudni_dTau(HEOS, k, xN_flag);
summer += d2_ndln_fugacity_i_dnj_ddelta_dtau__constx(HEOS, i, j, xN_flag)*nddeltadni__constT_V_nj(HEOS, k, xN_flag);
@@ -205,11 +205,11 @@ class MixtureDerivatives{
for (std::size_t m = 0; m < mmax; ++m){
summer -= HEOS.mole_fractions[m]*d2_ndln_fugacity_i_dnj_dxk_dTau__constdelta(HEOS, i, j, m, xN_flag);
}
return 1/RT*summer + 1/(HEOS.gas_constant()*HEOS.T_reducing())*nd_ndln_fugacity_i_dnj_dnk__constT_V_xi(HEOS,i,j,k,xN_flag) - d_nAij_dX(HEOS, i, j, xN_flag, WRT);
}
else if (WRT== iDelta){
return -1000;
}
return summer;
}
static Eigen::MatrixXd Lstar(HelmholtzEOSMixtureBackend &HEOS, x_N_dependency_flag xN_flag){
std::size_t N = HEOS.mole_fractions.size();