From 5e86a1171f9b5787899e347f0d5857dba79d97f7 Mon Sep 17 00:00:00 2001
From: Ian Bell <ian.bell@nist.gov>
Date: Thu, 14 Oct 2021 15:53:50 -0400
Subject: [PATCH] Also time Peng-Robinson
---
src/time_Ar0n.cpp | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/src/time_Ar0n.cpp b/src/time_Ar0n.cpp
index 82bc69e..df8b6b7 100644
--- a/src/time_Ar0n.cpp
+++ b/src/time_Ar0n.cpp
@@ -20,6 +20,7 @@
#include "teqp/models/multifluid.hpp"
#include "teqp/models/eos.hpp"
+#include "teqp/models/cubics.hpp"
#include "teqp/models/pcsaft.hpp"
#include "teqp/derivs.hpp"
@@ -250,11 +251,19 @@ int main()
};
auto SAFT = build_PCSAFT(Ncomp);
+ std::valarray<double> Tc_K = { 369.89 }, pc_Pa = { 4251200.0 }, acentric = { 0.1521 };
+ auto PR = canonical_PR(Tc_K, pc_Pa, acentric);
+
outputs.push_back(one_deriv<0, 0>(thing, Ncomp, taus, deltas, vdW, "vdW", Ts, rhos)); append_Ncomp();
outputs.push_back(one_deriv<0, 1>(thing, Ncomp, taus, deltas, vdW, "vdW", Ts, rhos)); append_Ncomp();
outputs.push_back(one_deriv<0, 2>(thing, Ncomp, taus, deltas, vdW, "vdW", Ts, rhos)); append_Ncomp();
outputs.push_back(one_deriv<0, 3>(thing, Ncomp, taus, deltas, vdW, "vdW", Ts, rhos)); append_Ncomp();
+ outputs.push_back(one_deriv<0, 0>(thing, Ncomp, taus, deltas, PR, "PR", Ts, rhos)); append_Ncomp();
+ outputs.push_back(one_deriv<0, 1>(thing, Ncomp, taus, deltas, PR, "PR", Ts, rhos)); append_Ncomp();
+ outputs.push_back(one_deriv<0, 2>(thing, Ncomp, taus, deltas, PR, "PR", Ts, rhos)); append_Ncomp();
+ outputs.push_back(one_deriv<0, 3>(thing, Ncomp, taus, deltas, PR, "PR", Ts, rhos)); append_Ncomp();
+
outputs.push_back(one_deriv<0, 0>(thing, Ncomp, taus, deltas, SAFT, "PCSAFT", Ts, rhos)); append_Ncomp();
outputs.push_back(one_deriv<0, 1>(thing, Ncomp, taus, deltas, SAFT, "PCSAFT", Ts, rhos)); append_Ncomp();
outputs.push_back(one_deriv<0, 2>(thing, Ncomp, taus, deltas, SAFT, "PCSAFT", Ts, rhos)); append_Ncomp();
--
GitLab