CoolProp 8.0.0
An open-source fluid property and humid air property database
MixtureDerivatives.h
Go to the documentation of this file.
1
8// ***************************************************************
9// ***************************************************************
10// ***************** MIXTURE DERIVATIVES ***********************
11// ***************************************************************
12// ***************************************************************
13
14#ifndef MIXTURE_DERIVATIVES_H
15#define MIXTURE_DERIVATIVES_H
16
17#include <Eigen/Core>
18#include "Backends/Helmholtz/HelmholtzEOSMixtureBackend.h"
19#include "ReducingFunctions.h"
20
21namespace CoolProp {
22
23class HelmholtzEOSMixtureBackend;
24
33class MixtureDerivatives
34{
35 public:
44 static CoolPropDbl ndpdV__constT_n(HelmholtzEOSMixtureBackend& HEOS);
45
54 static CoolPropDbl dpdT__constV_n(HelmholtzEOSMixtureBackend& HEOS);
55
56 static CoolPropDbl dpdrho__constT_n(HelmholtzEOSMixtureBackend& HEOS);
57
68 static CoolPropDbl ndpdni__constT_V_nj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
69
80 static CoolPropDbl partial_molar_volume(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
81
89 static CoolPropDbl fugacity_i(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
90
97 static CoolPropDbl ln_fugacity_coefficient(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
98
109 static CoolPropDbl dln_fugacity_i_dT__constrho_n(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
110
121 static CoolPropDbl dln_fugacity_i_drho__constT_n(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
122
134 static CoolPropDbl dln_fugacity_coefficient_dT__constp_n(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
135
136 static CoolPropDbl dln_fugacity_i_dT__constp_n(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
137 static CoolPropDbl dln_fugacity_i_dp__constT_n(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
138 static CoolPropDbl dln_fugacity_dxj__constT_p_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
139 static CoolPropDbl dln_fugacity_i_dtau__constdelta_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
140 static CoolPropDbl dln_fugacity_i_ddelta__consttau_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
141 static CoolPropDbl dln_fugacity_dxj__constT_rho_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
142
143 static CoolPropDbl ndln_fugacity_i_dnj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
144 static CoolPropDbl d_ndln_fugacity_i_dnj_dtau__constdelta_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
145 x_N_dependency_flag xN_flag);
146 static CoolPropDbl d2_ndln_fugacity_i_dnj_dtau2__constdelta_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
147 x_N_dependency_flag xN_flag);
148 static CoolPropDbl d2_ndln_fugacity_i_dnj_ddelta_dtau__constx(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
149 x_N_dependency_flag xN_flag);
150 static CoolPropDbl d2_ndln_fugacity_i_dnj_ddelta2__consttau_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
151 x_N_dependency_flag xN_flag);
152 static CoolPropDbl d_ndln_fugacity_i_dnj_ddelta__consttau_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
153 x_N_dependency_flag xN_flag);
154 static CoolPropDbl d_ndln_fugacity_i_dnj_ddxk__consttau_delta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
155 x_N_dependency_flag xN_flag);
156
157 static CoolPropDbl d2_ndln_fugacity_i_dnj_dxk_dTau__constdelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
158 x_N_dependency_flag xN_flag);
159 static CoolPropDbl d2_ndln_fugacity_i_dnj_dxk_dDelta__consttau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
160 x_N_dependency_flag xN_flag);
161
162 static CoolPropDbl nd_ndln_fugacity_i_dnj_dnk__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
163 x_N_dependency_flag xN_flag);
164
165 static CoolPropDbl nAij(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag) {
166 return ndln_fugacity_i_dnj__constT_V_xi(HEOS, i, j, xN_flag);
167 }
168 static CoolPropDbl n2Aijk(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k, x_N_dependency_flag xN_flag) {
169 return nd_ndln_fugacity_i_dnj_dnk__constT_V_xi(HEOS, i, j, k, xN_flag) - nAij(HEOS, i, j, xN_flag);
170 }
171 static CoolPropDbl d_nAij_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag) {
172 return d_nAij_dX(HEOS, i, j, xN_flag, iTau);
173 }
174 static CoolPropDbl d_nAij_dDelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag) {
175 return d_nAij_dX(HEOS, i, j, xN_flag, iDelta);
176 }
177 static CoolPropDbl d_nAij_dX(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag, parameters WRT) {
178 if (WRT == iTau) {
179 return MixtureDerivatives::d_ndln_fugacity_i_dnj_dtau__constdelta_x(HEOS, i, j, xN_flag);
180 } else if (WRT == iDelta) {
181 return MixtureDerivatives::d_ndln_fugacity_i_dnj_ddelta__consttau_x(HEOS, i, j, xN_flag);
182 } else {
183 throw ValueError(format("wrong WRT"));
184 }
185 }
186 static CoolPropDbl d_n2Aijk_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k, x_N_dependency_flag xN_flag) {
187 return d_n2Aijk_dX(HEOS, i, j, k, xN_flag, iTau);
188 }
189 static CoolPropDbl d_n2Aijk_dDelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k, x_N_dependency_flag xN_flag) {
190 return d_n2Aijk_dX(HEOS, i, j, k, xN_flag, iDelta);
191 }
192 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,
193 parameters WRT) {
194 double summer = 0;
195 if (WRT == iTau) {
196 summer += d2_ndln_fugacity_i_dnj_dtau2__constdelta_x(HEOS, i, j, xN_flag) * ndtaudni__constT_V_nj(HEOS, k, xN_flag);
197 summer += d_ndln_fugacity_i_dnj_dtau__constdelta_x(HEOS, i, j, xN_flag) * d_ndtaudni_dTau(HEOS, k, xN_flag);
198 summer += d2_ndln_fugacity_i_dnj_ddelta_dtau__constx(HEOS, i, j, xN_flag) * nddeltadni__constT_V_nj(HEOS, k, xN_flag);
199 summer += d2_ndln_fugacity_i_dnj_dxk_dTau__constdelta(HEOS, i, j, k, xN_flag);
200 std::size_t mmax = HEOS.mole_fractions.size();
201 if (xN_flag == XN_DEPENDENT) {
202 mmax--;
203 }
204 for (std::size_t m = 0; m < mmax; ++m) {
205 summer -= HEOS.mole_fractions[m] * d2_ndln_fugacity_i_dnj_dxk_dTau__constdelta(HEOS, i, j, m, xN_flag);
206 }
207 } else if (WRT == iDelta) {
208 summer += d2_ndln_fugacity_i_dnj_ddelta_dtau__constx(HEOS, i, j, xN_flag) * ndtaudni__constT_V_nj(HEOS, k, xN_flag);
209 summer += d2_ndln_fugacity_i_dnj_ddelta2__consttau_x(HEOS, i, j, xN_flag) * nddeltadni__constT_V_nj(HEOS, k, xN_flag);
210 summer += d_ndln_fugacity_i_dnj_ddelta__consttau_x(HEOS, i, j, xN_flag) * d_nddeltadni_dDelta(HEOS, k, xN_flag);
211 summer += d2_ndln_fugacity_i_dnj_dxk_dDelta__consttau(HEOS, i, j, k, xN_flag);
212 std::size_t mmax = HEOS.mole_fractions.size();
213 if (xN_flag == XN_DEPENDENT) {
214 mmax--;
215 }
216 for (std::size_t m = 0; m < mmax; ++m) {
217 summer -= HEOS.mole_fractions[m] * d2_ndln_fugacity_i_dnj_dxk_dDelta__consttau(HEOS, i, j, m, xN_flag);
218 }
219 } else {
220 return _HUGE;
221 }
222 return summer - d_nAij_dX(HEOS, i, j, xN_flag, WRT);
223 }
224 static Eigen::MatrixXd Lstar(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag) {
225 std::size_t N = HEOS.mole_fractions.size();
226 Eigen::MatrixXd L;
227 L.resize(N, N);
228 for (std::size_t i = 0; i < N; ++i) {
229 for (std::size_t j = i; j < N; ++j) {
230 L(i, j) = nAij(HEOS, i, j, xN_flag);
231 }
232 }
233 // Fill in the symmetric elements
234 for (std::size_t i = 0; i < N; ++i) {
235 for (std::size_t j = 0; j < i; ++j) {
236 L(i, j) = L(j, i);
237 }
238 }
239 return L;
240 }
241 static Eigen::MatrixXd dLstar_dX(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag, parameters WRT) {
242
243 std::size_t N = HEOS.mole_fractions.size();
244 Eigen::MatrixXd dLstar_dX(N, N);
245 for (std::size_t i = 0; i < N; ++i) {
246 for (std::size_t j = i; j < N; ++j) {
247 dLstar_dX(i, j) = d_nAij_dX(HEOS, i, j, xN_flag, WRT);
248 }
249 }
250 // Fill in the symmetric elements
251 for (std::size_t i = 0; i < N; ++i) {
252 for (std::size_t j = 0; j < i; ++j) {
253 dLstar_dX(i, j) = dLstar_dX(j, i);
254 }
255 }
256 return dLstar_dX;
257 }
258 static Eigen::MatrixXd d2Lstar_dX2(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag, parameters WRT1, parameters WRT2) {
259
260 std::size_t N = HEOS.mole_fractions.size();
261 Eigen::MatrixXd d2Lstar_dX2(N, N);
262 for (std::size_t i = 0; i < N; ++i) {
263 for (std::size_t j = i; j < N; ++j) {
264 if (WRT1 == iTau && WRT2 == iTau) {
265 d2Lstar_dX2(i, j) = MixtureDerivatives::d2_ndln_fugacity_i_dnj_dtau2__constdelta_x(HEOS, i, j, xN_flag);
266 } else {
267 throw ValueError(format("d2Lstar_dX2 invalid WRT"));
268 }
269 }
270 }
271 // Fill in the symmetric elements
272 for (std::size_t i = 0; i < N; ++i) {
273 for (std::size_t j = 0; j < i; ++j) {
274 d2Lstar_dX2(i, j) = d2Lstar_dX2(j, i);
275 }
276 }
277 return d2Lstar_dX2;
278 }
279 static Eigen::MatrixXd Mstar(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag, Eigen::MatrixXd& L) {
280
281 std::size_t N = HEOS.mole_fractions.size();
282 Eigen::MatrixXd M = L, adjL = adjugate(L);
283
284 // Last row
285 for (std::size_t i = 0; i < N; ++i) {
286 Eigen::MatrixXd n2dLdni(N, N);
287 for (std::size_t j = 0; j < N; ++j) {
288 for (std::size_t k = j; k < N; ++k) {
289 n2dLdni(j, k) = n2Aijk(HEOS, j, k, i, xN_flag);
290 // Fill in the symmetric elements
291 n2dLdni(k, j) = n2dLdni(j, k);
292 }
293 }
294 M(N - 1, i) = (adjL * n2dLdni).trace();
295 }
296 return M;
297 }
298 static Eigen::MatrixXd dMstar_dX(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag, parameters WRT, Eigen::MatrixXd& L,
299 Eigen::MatrixXd& dL_dX) {
300
301 std::size_t N = HEOS.mole_fractions.size();
302 Eigen::MatrixXd dMstar = dL_dX, adjL = adjugate(L), d_adjL_dX = adjugate_derivative(L, dL_dX);
303
304 // Last row in the d(Mstar)/d(X) requires derivatives of
305 for (std::size_t i = 0; i < N; ++i) {
306 Eigen::MatrixXd n2dLdni(N, N), d_n2dLdni_dX(N, N);
307 for (std::size_t j = 0; j < N; ++j) {
308 for (std::size_t k = j; k < N; ++k) {
309 n2dLdni(j, k) = n2Aijk(HEOS, j, k, i, xN_flag);
310 d_n2dLdni_dX(j, k) = d_n2Aijk_dX(HEOS, j, k, i, xN_flag, WRT);
311 // Fill in the symmetric elements
312 n2dLdni(k, j) = n2dLdni(j, k);
313 d_n2dLdni_dX(k, j) = d_n2dLdni_dX(j, k);
314 }
315 }
316 dMstar(N - 1, i) = (n2dLdni * d_adjL_dX + adjL * d_n2dLdni_dX).trace();
317 }
318 return dMstar;
319 }
320
337 static CoolPropDbl ndalphar_dni__constT_V_nj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
338
340 static CoolPropDbl dnalphar_dni__constT_V_nj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
341
352 static CoolPropDbl dln_fugacity_coefficient_dp__constT_n(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
353
369 static CoolPropDbl ndln_fugacity_coefficient_dnj__constT_p(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
370 x_N_dependency_flag xN_flag);
371
384 static CoolPropDbl dln_fugacity_coefficient_dxj__constT_p_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
385 x_N_dependency_flag xN_flag);
386
397 static CoolPropDbl dpdxj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t j, x_N_dependency_flag xN_flag);
398
410 static CoolPropDbl d2nalphar_dxj_dni__constT_V(HelmholtzEOSMixtureBackend& HEOS, std::size_t j, std::size_t i, x_N_dependency_flag xN_flag) {
411 return MixtureDerivatives::d_ndalphardni_dxj__constT_V_xi(HEOS, i, j, xN_flag)
412 + MixtureDerivatives::dalphar_dxj__constT_V_xi(HEOS, j, xN_flag);
413 };
414
417 static CoolPropDbl dalphar_dxj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t j, x_N_dependency_flag xN_flag);
418
421 static CoolPropDbl d_ndalphardni_dxj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
422
425 static CoolPropDbl d_dalpharddelta_dxj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t j, x_N_dependency_flag xN_flag);
426
429 static CoolPropDbl ddelta_dxj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t j, x_N_dependency_flag xN_flag);
430
433 static CoolPropDbl dtau_dxj__constT_V_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t j, x_N_dependency_flag xN_flag);
434
448 static CoolPropDbl d2nalphar_dni_dT(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
449
462 static CoolPropDbl d_ndalphardni_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
463
474 static CoolPropDbl d2_ndalphardni_dTau2(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
475
476 static CoolPropDbl d3_ndalphardni_dTau3(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
477
490 static CoolPropDbl d_ndalphardni_dDelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
491
504 static CoolPropDbl d2_ndalphardni_dDelta2(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
505
506 static CoolPropDbl d3_ndalphardni_dDelta3(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
507
521 static CoolPropDbl d2_ndalphardni_dDelta_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
522
523 static CoolPropDbl d3_ndalphardni_dDelta2_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
524 static CoolPropDbl d3_ndalphardni_dDelta_dTau2(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
525
541 static CoolPropDbl d2_ndalphardni_dxj_dDelta__consttau_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
542 x_N_dependency_flag xN_flag);
543
559 static CoolPropDbl d2_ndalphardni_dxj_dTau__constdelta_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
560 x_N_dependency_flag xN_flag);
561 static CoolPropDbl d3_ndalphardni_dxj_dTau2__constdelta_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
562 x_N_dependency_flag xN_flag);
563 static CoolPropDbl d3_ndalphardni_dxj_dDelta2__consttau_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
564 x_N_dependency_flag xN_flag);
565 static CoolPropDbl d3_ndalphardni_dxj_dDelta_dTau__constxi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
566 x_N_dependency_flag xN_flag);
567
590 static CoolPropDbl nd2nalphardnidnj__constT_V(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
591
592 /* \brief GERG 2004 Eqn. 7.47
593 * \f{eqnarray*}{
594 * n\left( \frac{\partial}{\partial n_j}\left(n\left(\frac{\partial \alpha^r}{\partial n_i}\right)_{T,V,n_j} \right) \right)_{T,V,n_i} &=& \left( \frac{\partial}{\partial \delta}\left(n\left(\frac{\partial\alpha^r}{\partial n_i}\right)_{T,V,n_j}\right)\right)_{\tau,\bar x}\cdot n\left(\frac{\partial\delta}{\partial n_j}\right)_{T,V,n_i}\\
595 * &+& \left( \frac{\partial}{\partial \tau}\left(n\left(\frac{\partial\alpha^r}{\partial n_i}\right)_{T,V,n_j}\right)\right)_{\tau,\bar x}\cdot n\left(\frac{\partial\tau}{\partial n_j}\right)_{T,V,n_i}\\
596 * &+& \left( \frac{\partial}{\partial x_j}\left(n\left(\frac{\partial\alpha^r}{\partial n_i}\right)_{T,V,n_j}\right)\right)_{\delta,\tau,x_i}-\sum_{k=1}^{N}x_k \left( \frac{\partial}{\partial x_k}\left(n\left(\frac{\partial\alpha^r}{\partial n_i}\right)_{T,V,n_j}\right)\right)_{\delta,\tau,x_i}\\
597 * \f}
598 * @param HEOS The HelmholtzEOSMixtureBackend to be used
599 * @param i The index of interest
600 * @param j The second index of interest
601 * @param xN_flag A flag specifying whether the all mole fractions are independent or only the first N-1
602 */
603 static CoolPropDbl nd_ndalphardni_dnj__constT_V(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
604
605 /* \brief \f$\tau\f$ derivative of GERG 2004 7.47
606 *
607 */
608 static CoolPropDbl d_nd_ndalphardni_dnj_dTau__constdelta_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
609 x_N_dependency_flag xN_flag);
610 static CoolPropDbl d2_nd_ndalphardni_dnj_dTau2__constdelta_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
611 x_N_dependency_flag xN_flag);
612 static CoolPropDbl d2_nd_ndalphardni_dnj_dDelta2__consttau_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
613 x_N_dependency_flag xN_flag);
614 static CoolPropDbl d2_nd_ndalphardni_dnj_dDelta_dTau__constx(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
615 x_N_dependency_flag xN_flag);
616
617 /* \brief \f$\delta\f$ derivative of GERG 2004 7.47
618 *
619 */
620 static CoolPropDbl d_nd_ndalphardni_dnj_dDelta__consttau_x(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
621 x_N_dependency_flag xN_flag);
622
623 /* \brief \f$x_k\f$ derivative of GERG 2004 7.47
624 *
625 */
626 static CoolPropDbl d_nd_ndalphardni_dnj_dxk__consttau_delta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
627 x_N_dependency_flag xN_flag);
628
629 static CoolPropDbl d2_nd_ndalphardni_dnj_dxk_dTau__constdelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
630 x_N_dependency_flag xN_flag);
631 static CoolPropDbl d2_nd_ndalphardni_dnj_dxk_dDelta__consttau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
632 x_N_dependency_flag xN_flag);
633
644 static CoolPropDbl nddeltadni__constT_V_nj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
645
646 static CoolPropDbl d_nddeltadni_dDelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
647
648 static CoolPropDbl d_nddeltadni_dxj__constdelta_tau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
649
650 static CoolPropDbl d2_nddeltadni_dxj_dDelta__consttau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
651 x_N_dependency_flag xN_flag);
652
663 static CoolPropDbl ndtaudni__constT_V_nj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
664
665 static CoolPropDbl d_ndtaudni_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
666
667 static CoolPropDbl d_ndtaudni_dxj__constdelta_tau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
668
669 static CoolPropDbl d2_ndtaudni_dxj_dTau__constdelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
670
686 static CoolPropDbl d_ndalphardni_dxj__constdelta_tau_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
687 x_N_dependency_flag xN_flag);
688
689 /* \brief \f$x_k\f$ derivative of GERG 2004 7.52
690 *
691 */
692 static CoolPropDbl d2_ndalphardni_dxj_dxk__constdelta_tau_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
693 x_N_dependency_flag xN_flag);
694
695 static CoolPropDbl d3_ndalphardni_dxj_dxk_dTau__constdelta_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
696 x_N_dependency_flag xN_flag);
697 static CoolPropDbl d3_ndalphardni_dxj_dxk_dDelta__consttau_xi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
698 x_N_dependency_flag xN_flag);
699
700 static CoolPropDbl dalpha0_dxi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
701 static CoolPropDbl d2alpha0_dxi_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
702 static CoolPropDbl d2alpha0_dxi_dDelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
703 static CoolPropDbl d2alpha0dxidxj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
704
706 static CoolPropDbl psi(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
707 return (HEOS.alphar() + HEOS.alpha0()) * HEOS.rhomolar() * HEOS.gas_constant() * HEOS.T();
708 }
710 static CoolPropDbl dpsi_dDelta(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
711
713 static CoolPropDbl dpsi_dTau(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
714
716 static CoolPropDbl dpsir_dTau(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
717
719 static CoolPropDbl dpsi_dxi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
720
722 static CoolPropDbl dpsir_dxi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
723
725 static CoolPropDbl d_rhorTr_dxi(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag);
726
728 static CoolPropDbl d2_rhorTr_dxidxj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag);
729
731 static CoolPropDbl d2psi_dDelta2(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
732
734 static CoolPropDbl d2psi_dDelta_dTau(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
735
737 static CoolPropDbl d2psir_dDelta_dTau(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
738
740 static CoolPropDbl d2psi_dTau2(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
741
743 static CoolPropDbl d2psir_dTau2(HelmholtzEOSMixtureBackend& HEOS, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
744
746 static CoolPropDbl d2psi_dxi_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
747
749 static CoolPropDbl d2psir_dxi_dTau(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
750
752 static CoolPropDbl d2psi_dxi_dDelta(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
753
755 static CoolPropDbl d2psi_dxi_dxj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
756
758 static CoolPropDbl d2psir_dxi_dxj(HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, x_N_dependency_flag xN_flag = XN_INDEPENDENT);
759
768
769 static CoolPropDbl PSI_rho(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
770 return HEOS.Reducing->PSI_rho(HEOS.mole_fractions, i, xN_flag);
771 }
772 static CoolPropDbl d_PSI_rho_dxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
773 CoolProp::x_N_dependency_flag xN_flag) {
774 return HEOS.Reducing->d_PSI_rho_dxj(HEOS.mole_fractions, i, j, xN_flag);
775 }
776 static CoolPropDbl d2_PSI_rho_dxj_dxk(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
777 CoolProp::x_N_dependency_flag xN_flag) {
778 return HEOS.Reducing->d2_PSI_rho_dxj_dxk(HEOS.mole_fractions, i, j, k, xN_flag);
779 }
780 static CoolPropDbl PSI_T(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
781 return HEOS.Reducing->PSI_T(HEOS.mole_fractions, i, xN_flag);
782 }
783 static CoolPropDbl d_PSI_T_dxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, CoolProp::x_N_dependency_flag xN_flag) {
784 return HEOS.Reducing->d_PSI_T_dxj(HEOS.mole_fractions, i, j, xN_flag);
785 }
786 static CoolPropDbl d2_PSI_T_dxj_dxk(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
787 CoolProp::x_N_dependency_flag xN_flag) {
788 return HEOS.Reducing->d2_PSI_T_dxj_dxk(HEOS.mole_fractions, i, j, k, xN_flag);
789 }
790
791 static CoolPropDbl alpha0(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
792 return HEOS.calc_alpha0_deriv_nocache(0, 0, HEOS.mole_fractions, HEOS.tau(), HEOS.delta(), HEOS.T_reducing(), HEOS.rhomolar_reducing());
793 }
794
795 static CoolPropDbl alphar(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
796 bool cache_values = false;
797 HelmholtzDerivatives HD = HEOS.residual_helmholtz->all(HEOS, HEOS.mole_fractions, HEOS.tau(), HEOS.delta(), cache_values);
798 return HD.alphar;
799 }
800 static CoolPropDbl dalphar_dxi(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
801 return HEOS.residual_helmholtz->dalphar_dxi(HEOS, i, xN_flag);
802 }
803 static CoolPropDbl d2alphar_dxi_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
804 return HEOS.residual_helmholtz->d2alphar_dxi_dTau(HEOS, i, xN_flag);
805 }
806 static CoolPropDbl d3alphar_dxi_dTau2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
807 return HEOS.residual_helmholtz->d3alphar_dxi_dTau2(HEOS, i, xN_flag);
808 }
809 static CoolPropDbl d4alphar_dxi_dTau3(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
810 return HEOS.residual_helmholtz->d4alphar_dxi_dTau3(HEOS, i, xN_flag);
811 }
812 static CoolPropDbl d2alphar_dxi_dDelta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
813 return HEOS.residual_helmholtz->d2alphar_dxi_dDelta(HEOS, i, xN_flag);
814 }
815 static CoolPropDbl d3alphar_dxi_dDelta2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
816 return HEOS.residual_helmholtz->d3alphar_dxi_dDelta2(HEOS, i, xN_flag);
817 }
818 static CoolPropDbl d4alphar_dxi_dDelta3(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
819 return HEOS.residual_helmholtz->d4alphar_dxi_dDelta3(HEOS, i, xN_flag);
820 }
821 static CoolPropDbl d3alphar_dxi_dDelta_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
822 return HEOS.residual_helmholtz->d3alphar_dxi_dDelta_dTau(HEOS, i, xN_flag);
823 }
824 static CoolPropDbl d4alphar_dxi_dDelta2_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
825 return HEOS.residual_helmholtz->d4alphar_dxi_dDelta2_dTau(HEOS, i, xN_flag);
826 }
827 static CoolPropDbl d4alphar_dxi_dDelta_dTau2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
828 return HEOS.residual_helmholtz->d4alphar_dxi_dDelta_dTau2(HEOS, i, xN_flag);
829 }
830 static CoolPropDbl d2alphardxidxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
831 CoolProp::x_N_dependency_flag xN_flag) {
832 return HEOS.residual_helmholtz->d2alphardxidxj(HEOS, i, j, xN_flag);
833 }
834 static CoolPropDbl d4alphar_dxi_dxj_dTau2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
835 CoolProp::x_N_dependency_flag xN_flag) {
836 return HEOS.residual_helmholtz->d4alphar_dxi_dxj_dTau2(HEOS, i, j, xN_flag);
837 }
838 static CoolPropDbl d4alphar_dxi_dxj_dDelta2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
839 CoolProp::x_N_dependency_flag xN_flag) {
840 return HEOS.residual_helmholtz->d4alphar_dxi_dxj_dDelta2(HEOS, i, j, xN_flag);
841 }
842 static CoolPropDbl d4alphar_dxi_dxj_dDelta_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
843 CoolProp::x_N_dependency_flag xN_flag) {
844 return HEOS.residual_helmholtz->d4alphar_dxi_dxj_dDelta_dTau(HEOS, i, j, xN_flag);
845 }
846 static CoolPropDbl d3alphar_dxi_dxj_dDelta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
847 CoolProp::x_N_dependency_flag xN_flag) {
848 return HEOS.residual_helmholtz->d3alphar_dxi_dxj_dDelta(HEOS, i, j, xN_flag);
849 }
850 static CoolPropDbl d3alphar_dxi_dxj_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
851 CoolProp::x_N_dependency_flag xN_flag) {
852 return HEOS.residual_helmholtz->d3alphar_dxi_dxj_dTau(HEOS, i, j, xN_flag);
853 }
854 static CoolPropDbl d3alphardxidxjdxk(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
855 CoolProp::x_N_dependency_flag xN_flag) {
856 return HEOS.residual_helmholtz->d3alphardxidxjdxk(HEOS, i, j, k, xN_flag);
857 }
858 static CoolPropDbl tau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
859 return HEOS.tau();
860 }
861 static CoolPropDbl Tr(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
862 return HEOS.Reducing->Tr(HEOS.get_mole_fractions());
863 }
864 static CoolPropDbl dTrdxi__constxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
865 return HEOS.Reducing->dTrdxi__constxj(HEOS.get_mole_fractions(), i, xN_flag);
866 }
867 static CoolPropDbl d2Trdxi2__constxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
868 return HEOS.Reducing->d2Trdxi2__constxj(HEOS.get_mole_fractions(), i, xN_flag);
869 }
870 static CoolPropDbl d2Trdxidxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, CoolProp::x_N_dependency_flag xN_flag) {
871 return HEOS.Reducing->d2Trdxidxj(HEOS.get_mole_fractions(), i, j, xN_flag);
872 }
873 static CoolPropDbl d3Trdxidxjdxk(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
874 CoolProp::x_N_dependency_flag xN_flag) {
875 return HEOS.Reducing->d3Trdxidxjdxk(HEOS.get_mole_fractions(), i, j, k, xN_flag);
876 }
877 static CoolPropDbl ndTrdni__constnj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
878 return HEOS.Reducing->ndTrdni__constnj(HEOS.get_mole_fractions(), i, xN_flag);
879 }
880 static CoolPropDbl d_ndTrdni_dxj__constxi(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
881 CoolProp::x_N_dependency_flag xN_flag) {
882 return HEOS.Reducing->d_ndTrdni_dxj__constxi(HEOS.get_mole_fractions(), i, j, xN_flag);
883 }
884 static CoolPropDbl d2_ndTrdni_dxj_dxk__constxi(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
885 CoolProp::x_N_dependency_flag xN_flag) {
886 return HEOS.Reducing->d2_ndTrdni_dxj_dxk__constxi(HEOS.get_mole_fractions(), i, j, k, xN_flag);
887 }
888 static CoolPropDbl dTr_dgammaT(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
889 return HEOS.Reducing->dTr_dgammaT(HEOS.get_mole_fractions());
890 }
891 static CoolPropDbl dTr_dbetaT(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
892 return HEOS.Reducing->dTr_dbetaT(HEOS.get_mole_fractions());
893 }
894 static CoolPropDbl d2Tr_dxidgammaT(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag) {
895 return HEOS.Reducing->d2Tr_dxidgammaT(HEOS.get_mole_fractions(), i, xN_flag);
896 }
897 static CoolPropDbl d2Tr_dxidbetaT(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag) {
898 return HEOS.Reducing->d2Tr_dxidbetaT(HEOS.get_mole_fractions(), i, xN_flag);
899 }
900
901 static CoolPropDbl delta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
902 return HEOS.delta();
903 }
904 static CoolPropDbl rhormolar(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
905 return HEOS.Reducing->rhormolar(HEOS.get_mole_fractions());
906 }
907 static CoolPropDbl drhormolardxi__constxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
908 return HEOS.Reducing->drhormolardxi__constxj(HEOS.get_mole_fractions(), i, xN_flag);
909 }
910 static CoolPropDbl d2rhormolardxidxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
911 CoolProp::x_N_dependency_flag xN_flag) {
912 return HEOS.Reducing->d2rhormolardxidxj(HEOS.get_mole_fractions(), i, j, xN_flag);
913 }
914 static CoolPropDbl d3rhormolardxidxjdxk(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
915 CoolProp::x_N_dependency_flag xN_flag) {
916 return HEOS.Reducing->d3rhormolardxidxjdxk(HEOS.get_mole_fractions(), i, j, k, xN_flag);
917 }
918 static CoolPropDbl drhormolar_dgammaV(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
919 return HEOS.Reducing->drhormolar_dgammaV(HEOS.get_mole_fractions());
920 }
921 static CoolPropDbl drhormolar_dbetaV(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
922 return HEOS.Reducing->drhormolar_dbetaV(HEOS.get_mole_fractions());
923 }
924 static CoolPropDbl d2rhormolar_dxidgammaV(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag) {
925 return HEOS.Reducing->d2rhormolar_dxidgammaV(HEOS.get_mole_fractions(), i, xN_flag);
926 }
927 static CoolPropDbl d2rhormolar_dxidbetaV(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, x_N_dependency_flag xN_flag) {
928 return HEOS.Reducing->d2rhormolar_dxidbetaV(HEOS.get_mole_fractions(), i, xN_flag);
929 }
930 static CoolPropDbl ndrhorbardni__constnj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
931 return HEOS.Reducing->ndrhorbardni__constnj(HEOS.get_mole_fractions(), i, xN_flag);
932 }
933 static CoolPropDbl d_ndrhorbardni_dxj__constxi(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
934 CoolProp::x_N_dependency_flag xN_flag) {
935 return HEOS.Reducing->d_ndrhorbardni_dxj__constxi(HEOS.get_mole_fractions(), i, j, xN_flag);
936 }
937 static CoolPropDbl d2_ndrhorbardni_dxj_dxk__constxi(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j, std::size_t k,
938 CoolProp::x_N_dependency_flag xN_flag) {
939 return HEOS.Reducing->d2_ndrhorbardni_dxj_dxk__constxi(HEOS.get_mole_fractions(), i, j, k, xN_flag);
940 }
941 static CoolPropDbl p(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag) {
942 return HEOS.p();
943 }
944
945 static CoolPropDbl dalpha0_dDelta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
946 return HEOS.dalpha0_dDelta();
947 }
948 static CoolPropDbl d2alpha0_dDelta2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
949 return HEOS.d2alpha0_dDelta2();
950 }
951 static CoolPropDbl dalpha0_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
952 return HEOS.dalpha0_dTau();
953 }
954 static CoolPropDbl d2alpha0_dTau2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
955 return HEOS.d2alpha0_dTau2();
956 }
957
958 static CoolPropDbl dalphar_dDelta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
959 return HEOS.dalphar_dDelta();
960 }
961 static CoolPropDbl d2alphar_dDelta2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
962 return HEOS.d2alphar_dDelta2();
963 }
964 static CoolPropDbl d2alphar_dDelta_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
965 return HEOS.d2alphar_dDelta_dTau();
966 }
967 static CoolPropDbl dalphar_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
968 return HEOS.dalphar_dTau();
969 }
970 static CoolPropDbl d2alphar_dTau2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
971 return HEOS.d2alphar_dTau2();
972 }
973
974 static CoolPropDbl alpha(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
975 return HEOS.alphar() + alpha0(HEOS, xN_flag);
976 }
977 static CoolPropDbl dalpha_dDelta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
978 return HEOS.dalphar_dDelta() + HEOS.dalpha0_dDelta();
979 }
980 static CoolPropDbl dalpha_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
981 return HEOS.dalphar_dTau() + HEOS.dalpha0_dTau();
982 }
983 static CoolPropDbl d2alpha_dDelta2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
984 return HEOS.d2alphar_dDelta2() + HEOS.d2alpha0_dDelta2();
985 }
986 static CoolPropDbl d2alpha_dDelta_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
987 return HEOS.d2alphar_dDelta_dTau() + HEOS.d2alpha0_dDelta_dTau();
988 }
989 static CoolPropDbl d2alpha_dTau2(CoolProp::HelmholtzEOSMixtureBackend& HEOS, CoolProp::x_N_dependency_flag xN_flag = XN_INDEPENDENT) {
990 return HEOS.d2alphar_dTau2() + HEOS.d2alpha0_dTau2();
991 }
992 static CoolPropDbl dalpha_dxi(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
993 return dalphar_dxi(HEOS, i, xN_flag) + dalpha0_dxi(HEOS, i, xN_flag);
994 }
995 static CoolPropDbl d2alpha_dxi_dDelta(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
996 return d2alphar_dxi_dDelta(HEOS, i, xN_flag) + d2alpha0_dxi_dDelta(HEOS, i, xN_flag);
997 }
998 static CoolPropDbl d2alpha_dxi_dTau(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
999 return d2alphar_dxi_dTau(HEOS, i, xN_flag) + d2alpha0_dxi_dTau(HEOS, i, xN_flag);
1000 }
1001 static CoolPropDbl d2alphadxidxj(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, std::size_t j,
1002 CoolProp::x_N_dependency_flag xN_flag) {
1003 return d2alphardxidxj(HEOS, i, j, xN_flag) + d2alpha0dxidxj(HEOS, i, j, xN_flag);
1004 }
1005 static CoolPropDbl ln_fugacity(CoolProp::HelmholtzEOSMixtureBackend& HEOS, std::size_t i, CoolProp::x_N_dependency_flag xN_flag) {
1006 double f_i = MixtureDerivatives::fugacity_i(HEOS, i, xN_flag);
1007 return log(f_i);
1008 }
1009
1010}; /* class MixtureDerivatives */
1011
1012} /* namespace CoolProp*/
1013#endif
Generated by doxygen 1.9.4