add switch

Author: PENGYAN777

SU2_CFD/src/fluid/CCoolProp.cpp

@@ -59,8 +59,12 @@ void CCoolProp::SetTDState_rhoe(su2double rho, su2double e) {
     // assume it is pure gas
     fluid_entity->specify_phase(CoolProp::iphase_gas);
     if (abs(Pressure / Pressure_Critical - 1) < epsilon) {
-      // if P is very cloase to Pc, the CoolProp has issue
-      Pressure = Pressure_Critical * (1-epsilon);
+      if (Pressure >= Pressure_Critical){
+        Pressure = Pressure_Critical * (1 + epsilon);
+      }
+      else {
+        Pressure = Pressure_Critical * (1 - epsilon);
+      }
     }
     fluid_entity->update(CoolProp::PT_INPUTS, Pressure, Temperature);
     if (abs(fluid_entity->rhomass() / Density - 1) < epsilon) {
@@ -80,7 +84,12 @@ void CCoolProp::SetTDState_rhoe(su2double rho, su2double e) {
 
 void CCoolProp::SetTDState_PT(su2double P, su2double T) {
   if (abs(P / Pressure_Critical - 1) < epsilon) {
-    P = Pressure_Critical * (1-epsilon);
+    if (P >= Pressure_Critical){
+      P = Pressure_Critical * (1 + epsilon);
+    }
+    else {
+      P = Pressure_Critical * (1 - epsilon);
+    }
   }
   fluid_entity->update(CoolProp::PT_INPUTS, P, T);
   su2double rho = fluid_entity->rhomass();
@@ -90,7 +99,12 @@ void CCoolProp::SetTDState_PT(su2double P, su2double T) {
 
 void CCoolProp::SetTDState_Prho(su2double P, su2double rho) {
   if (abs(P / Pressure_Critical - 1) < epsilon) {
-    P = Pressure_Critical * (1-epsilon);
+    if (P >= Pressure_Critical){
+      P = Pressure_Critical * (1 + epsilon);
+    }
+    else {
+      P = Pressure_Critical * (1 - epsilon);
+    }
   }
   fluid_entity->update(CoolProp::DmassP_INPUTS, rho, P);
   su2double e = fluid_entity->umass();
@@ -99,7 +113,12 @@ void CCoolProp::SetTDState_Prho(su2double P, su2double rho) {
 
 void CCoolProp::SetEnergy_Prho(su2double P, su2double rho) {
   if (abs(P / Pressure_Critical - 1) < epsilon) {
-    P = Pressure_Critical * (1-epsilon);
+    if (P >= Pressure_Critical){
+      P = Pressure_Critical * (1 + epsilon);
+    }
+    else {
+      P = Pressure_Critical * (1 - epsilon);
+    }
   }
   fluid_entity->update(CoolProp::DmassP_INPUTS, rho, P);
   StaticEnergy = fluid_entity->umass();
@@ -114,7 +133,12 @@ void CCoolProp::SetTDState_hs(su2double h, su2double s) {
 
 void CCoolProp::SetTDState_Ps(su2double P, su2double s) {
   if (abs(P / Pressure_Critical - 1) < epsilon) {
-    P = Pressure_Critical * (1-epsilon);
+    if (P >= Pressure_Critical){
+      P = Pressure_Critical * (1 + epsilon);
+    }
+    else {
+      P = Pressure_Critical * (1 - epsilon);
+    }
   }
   fluid_entity->update(CoolProp::PSmass_INPUTS, P, s);
   su2double rho = fluid_entity->rhomass();