pressure threshold dp

Author: PENGYAN777

SU2_CFD/include/fluid/CCoolProp.hpp

@@ -47,20 +47,21 @@ class CCoolProp final : public CFluidModel {
   su2double Pressure_Critical{0.0};    /*!< \brief critical pressure */
   su2double Temperature_Critical{0.0}; /*!< \brief critical temperature */
   su2double acentric_factor{0.0};      /*!< \brief acentric factor */
+  su2double dp{0.01};             /*!< threshold for pressure */
 #ifdef USE_COOLPROP
   std::unique_ptr<CoolProp::AbstractState> fluid_entity; /*!< \brief fluid entity */
-  su2double epsilon{0.01};                               /*!< threshold for pressure */
 #endif
   /*!
    * \brief Avoid critical pressure
    * \return .
    */
   void Check_Pressure(su2double Pressure) {
-    if (abs(Pressure / Pressure_Critical - 1) < epsilon) {
+    if (abs(Pressure / Pressure_Critical - 1) < dp) {
       if (Pressure >= Pressure_Critical) {
-        Pressure = Pressure_Critical * (1 + epsilon);
-      } else {
-        Pressure = Pressure_Critical * (1 - epsilon);
+        Pressure = Pressure_Critical * (1 + dp);
+      }
+      else {
+        Pressure = Pressure_Critical * (1 - dp);
       }
     }
   }

SU2_CFD/src/fluid/CCoolProp.cpp

@@ -60,7 +60,7 @@ void CCoolProp::SetTDState_rhoe(su2double rho, su2double e) {
     fluid_entity->specify_phase(CoolProp::iphase_gas);
     Check_Pressure(Pressure);
     fluid_entity->update(CoolProp::PT_INPUTS, Pressure, Temperature);
-    if (abs(fluid_entity->rhomass() / Density - 1) < epsilon) {
+    if (abs(fluid_entity->rhomass() / Density - 1) < dp) {
       // origial phase is near saturation gas, then just compute sound speed
       SoundSpeed2 = pow(fluid_entity->speed_sound(), 2);
     } else {