small fixes in NEMO

Author: CatarinaGarbacz

SU2_CFD/include/fluid/CMutationTCLib.hpp

@@ -45,8 +45,7 @@ class CMutationTCLib : public CNEMOGas {
 
   vector<su2double> Cv_ks,                /*!< \brief Species specific heats at constant volume. */
   es,                                     /*!< \brief Species energies. */
-  omega_vec,                              /*!< \brief Dummy vector for vibrational energy source term. */
-  h_RT;                                   /*!< \brief Enthalpy divided by R*T. */
+  omega_vec;                              /*!< \brief Dummy vector for vibrational energy source term. */
 
   su2double Tref;                         /*!< \brief Reference temperature. */
 

SU2_CFD/src/fluid/CMutationTCLib.cpp

@@ -36,7 +36,6 @@ CMutationTCLib::CMutationTCLib(const CConfig* config, unsigned short val_nDim):
 
   /* Allocating memory*/
   Cv_ks.resize(nEnergyEq*nSpecies,0.0);
-  h_RT.resize(nSpecies,0.0);
   es.resize(nEnergyEq*nSpecies,0.0);
   omega_vec.resize(1,0.0);
 
@@ -148,12 +147,7 @@ su2double CMutationTCLib::ComputeEveSourceTerm(){
 
 vector<su2double>& CMutationTCLib::ComputeSpeciesEnthalpy(su2double val_T, su2double val_Tve, su2double *val_eves){
 
-  su2double RuSI = UNIVERSAL_GAS_CONSTANT;
-  su2double Ru   = 1000.0*RuSI;
-
-  mix->speciesHOverRT(val_T, val_Tve, val_T, val_Tve, val_Tve, h_RT.data(), NULL, NULL, NULL, NULL, NULL);
-
-  for (iSpecies = 0; iSpecies < nSpecies; iSpecies++) hs[iSpecies] = h_RT[iSpecies]*(RuSI*val_T); 
+  mix->getEnthalpiesMass(hs.data());
 
   return hs;
 }
@@ -213,7 +207,7 @@ vector<su2double>& CMutationTCLib::GetSpeciesFormationEnthalpy() {
 
    mix->speciesHOverRT(Tref, Tref, Tref, Tref, Tref, NULL, NULL, NULL, NULL, NULL, hf_RT.data());
 
-   for (iSpecies = 0; iSpecies < nSpecies; iSpecies++) Enthalpy_Formation[iSpecies] = hf_RT[iSpecies]*(RuSI*Tref);
+   for (iSpecies = 0; iSpecies < nSpecies; iSpecies++) Enthalpy_Formation[iSpecies] = hf_RT[iSpecies]*(RuSI*Tref*1000.0)/MolarMass[iSpecies];
 
    return Enthalpy_Formation;  
 }

SU2_CFD/src/numerics/NEMO/NEMO_sources.cpp

@@ -376,7 +376,7 @@ CNumerics::ResidualType<> CSource_NEMO::ComputeAxisymmetric(const CConfig *confi
                                                      +v*TWO3*(2*PrimVar_Grad_i[nSpecies+2][1]-PrimVar_Grad_i[nSpecies+2][0]
                                                      -v*yinv+rho*turb_ke_i))
                                                      -total_conductivity_i*PrimVar_Grad_i[nSpecies][1])
-                                                     -TWO3*(AuxVar_Grad_i[1][1]+AuxVar_Grad_i[2][1]));
+                                                     -TWO3*(AuxVar_Grad_i[1][1]+AuxVar_Grad_i[2][0]));
     residual[nSpecies+3] -= Volume*(yinv*(sumJeve_y -qy_ve));
   }