From 38449f47ce5a20525b62a1f85f8cca5149ddd07a Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Wed, 27 Nov 2019 14:35:32 +0100
Subject: [PATCH] [MatL] Adsorption; Remove unused getDReactionRate.
---
MaterialLib/Adsorption/Adsorption.cpp | 23 -----------------------
MaterialLib/Adsorption/Adsorption.h | 11 -----------
2 files changed, 34 deletions(-)
diff --git a/MaterialLib/Adsorption/Adsorption.cpp b/MaterialLib/Adsorption/Adsorption.cpp
index 7a2f2624d1e..a5f257acd0c 100644
--- a/MaterialLib/Adsorption/Adsorption.cpp
+++ b/MaterialLib/Adsorption/Adsorption.cpp
@@ -118,29 +118,6 @@ double AdsorptionReaction::getReactionRate(const double p_Ads, const double T_Ad
// this the rate in terms of loading!
}
-void AdsorptionReaction::getDReactionRate(const double p_Ads, const double T_Ads,
- const double M_Ads, const double /*loading*/,
- std::array<double, 3> &dqdr) const
-{
- const double A = getPotential(p_Ads, T_Ads, M_Ads);
- const double p_S = getEquilibriumVapourPressure(T_Ads);
- const double dAdT = MaterialLib::PhysicalConstant::IdealGasConstant * log(p_S/p_Ads) / (M_Ads*1.e3);
- const double dAdp = - MaterialLib::PhysicalConstant::IdealGasConstant * T_Ads / M_Ads / p_Ads;
-
- const double W = characteristicCurve(A);
- const double dWdA = dCharacteristicCurve(A);
-
- const double rho_Ads = getAdsorbateDensity(T_Ads);
- const double drhodT = - rho_Ads * getAlphaT(T_Ads);
-
- dqdr = std::array<double, 3>{{
- rho_Ads*dWdA*dAdp,
- drhodT*W + rho_Ads*dWdA*dAdT,
- -k_rate
- }};
-}
-
-
// Evaluate adsorbtion potential A
double AdsorptionReaction::getPotential(const double p_Ads, double T_Ads, const double M_Ads) const
{
diff --git a/MaterialLib/Adsorption/Adsorption.h b/MaterialLib/Adsorption/Adsorption.h
index 7f2c51739f1..8f139642f51 100644
--- a/MaterialLib/Adsorption/Adsorption.h
+++ b/MaterialLib/Adsorption/Adsorption.h
@@ -39,17 +39,6 @@ public:
double getReactionRate(const double p_Ads, const double T_Ads,
const double M_Ads,
const double loading) const override;
- /**
- * @brief get_d_reaction_rate
- * @param p_Ads
- * @param T_Ads
- * @param M_Ads
- * @param loading
- * @param dqdr array containing the differentials wrt: p, T, C
- */
- virtual void getDReactionRate(const double p_Ads, const double T_Ads,
- const double M_Ads, const double loading,
- std::array<double, 3>& dqdr) const;
protected:
virtual double getAdsorbateDensity(const double T_Ads) const = 0;
--
GitLab