Merge remote-tracking branch 'upstream/develop' into heat_solver_fixes

Author: pcarruscag

Common/include/CConfig.hpp

@@ -1162,7 +1162,9 @@ class CConfig {
   su2double *Gas_Composition,               /*!< \brief Initial mass fractions of flow [dimensionless] */
   pnorm_heat;                               /*!< \brief pnorm for heat-flux. */
   bool frozen,                              /*!< \brief Flag for determining if mixture is frozen. */
-  ionization;                               /*!< \brief Flag for determining if free electron gas is in the mixture. */
+  ionization,                               /*!< \brief Flag for determining if free electron gas is in the mixture. */
+  vt_transfer_res_limit,                    /*!< \brief Flag for determining if residual limiting for source term VT-transfer is used. */
+  monoatomic;                               /*!< \brief Flag for monoatomic mixture. */
   string GasModel,                          /*!< \brief Gas Model. */
   *Wall_Catalytic;                          /*!< \brief Pointer to catalytic walls. */
 
@@ -5262,6 +5264,16 @@ class CConfig {
    */
   bool GetIonization(void) const { return ionization; }
 
+  /*!
+   * \brief Indicates whether the VT source residual is limited.
+   */
+  bool GetVTTransferResidualLimiting(void) const { return vt_transfer_res_limit; }
+
+  /*!
+   * \brief Indicates if mixture is monoatomic.
+   */
+  bool GetMonoatomic(void) const { return monoatomic; }
+
   /*!
    * \brief Information about computing and plotting the equivalent area distribution.
    * \return <code>TRUE</code> or <code>FALSE</code>  depending if we are computing the equivalent area.

Common/include/option_structure.hpp

@@ -200,7 +200,7 @@ enum ENUM_MAIN_SOLVER {
   FEM_LES = 29,                     /*!< \brief Definition of the finite element Large Eddy Simulation Navier-Stokes' (LES) solver. */
   MULTIPHYSICS = 30,
   NEMO_EULER = 41,                  /*!< \brief Definition of the NEMO Euler solver. */
-  NEMO_NAVIER_STOKES = 42          /*!< \brief Definition of the NEMO NS solver. */
+  NEMO_NAVIER_STOKES = 42           /*!< \brief Definition of the NEMO NS solver. */
 };
 static const MapType<string, ENUM_MAIN_SOLVER> Solver_Map = {
   MakePair("NONE", NO_SOLVER)

Common/src/CConfig.cpp

@@ -140,7 +140,6 @@ CConfig::CConfig(istream &case_buffer, unsigned short val_software, bool verb_hi
 
 }
 
-
 CConfig::CConfig(CConfig* config, char case_filename[MAX_STRING_SIZE], unsigned short val_software, unsigned short val_iZone, unsigned short val_nZone, bool verb_high) {
 
   caseName = config->GetCaseName();
@@ -356,7 +355,6 @@ void CConfig::addEnumOption(const string name, unsigned short & option_field, co
   return;
 }
 
-
 // input_size is the number of options read in from the config file
 template <class Tenum>
 void CConfig::addEnumListOption(const string name, unsigned short & input_size, unsigned short * & option_field, const map<string, Tenum> & enum_map) {
@@ -1165,6 +1163,10 @@ void CConfig::SetConfig_Options() {
   addBoolOption("FROZEN_MIXTURE", frozen, false);
   /* DESCRIPTION: Specify if there is ionization */
   addBoolOption("IONIZATION", ionization, false);
+  /* DESCRIPTION: Specify if there is VT transfer residual limiting */
+  addBoolOption("VT_RESIDUAL_LIMITING", vt_transfer_res_limit, false);
+  /* DESCRIPTION: Specify if the gas is monoatomic */
+  addBoolOption("MONOATOMIC", monoatomic, false);
   /* DESCRIPTION: List of catalytic walls */
   addStringListOption("CATALYTIC_WALL", nWall_Catalytic, Wall_Catalytic);
   /*!\brief MARKER_MONITORING\n DESCRIPTION: Marker(s) of the surface where evaluate the non-dimensional coefficients \ingroup Config*/
@@ -3561,6 +3563,10 @@ void CConfig::SetPostprocessing(unsigned short val_software, unsigned short val_
     SU2_MPI::Error("Only USER_DEFINED_NONEQ fluid model can be used with the NEMO solver. Mutation++ library will soon be available.", CURRENT_FUNCTION);
   }
 
+  if (nemo && Kind_TransCoeffModel != WILKE ) {
+    SU2_MPI::Error("Only WILKE transport model is stable for the NEMO solver.", CURRENT_FUNCTION);
+  }
+
   if (!ideal_gas && !nemo) {
     if (Kind_Upwind_Flow != ROE && Kind_Upwind_Flow != HLLC && Kind_Centered_Flow != JST) {
       SU2_MPI::Error("Only ROE Upwind, HLLC Upwind scheme, and JST scheme can be used for Non-Ideal Compressible Fluids", CURRENT_FUNCTION);
@@ -4990,6 +4996,8 @@ void CConfig::SetPostprocessing(unsigned short val_software, unsigned short val_
                    (Kind_Solver == INC_RANS) ||
                    (Kind_Solver == EULER) ||
                    (Kind_Solver == NAVIER_STOKES) ||
+                   (Kind_Solver == NEMO_EULER) ||
+                   (Kind_Solver == NEMO_NAVIER_STOKES) ||
                    (Kind_Solver == RANS) ||
                    (Kind_Solver == DISC_ADJ_EULER) ||
                    (Kind_Solver == DISC_ADJ_RANS) ||
@@ -5600,15 +5608,15 @@ void CConfig::SetOutput(unsigned short val_software, unsigned short val_izone) {
       case NEMO_EULER:
         if (Kind_Regime == COMPRESSIBLE) cout << "Compressible two-temperature thermochemical non-equilibrium Euler equations." << endl;
         if(Kind_FluidModel == USER_DEFINED_NONEQ){
-          if ((GasModel != "N2") && (GasModel != "AIR-5"))
-          SU2_MPI::Error("The GAS_MODEL given as input is not valid. Choose one of the options: N2, AIR-5.", CURRENT_FUNCTION);
+          if ((GasModel != "N2") && (GasModel != "AIR-5") && (GasModel != "ARGON"))
+          SU2_MPI::Error("The GAS_MODEL given as input is not valid. Choose one of the options: N2, AIR-5, ARGON.", CURRENT_FUNCTION);
         }
         break;
       case NEMO_NAVIER_STOKES:
         if (Kind_Regime == COMPRESSIBLE) cout << "Compressible two-temperature thermochemical non-equilibrium Navier-Stokes equations." << endl;
         if(Kind_FluidModel == USER_DEFINED_NONEQ){
-          if ((GasModel != "N2") && (GasModel != "AIR-5"))
-          SU2_MPI::Error("The GAS_MODEL given as input is not valid. Choose one of the options: N2, AIR-5.", CURRENT_FUNCTION);
+          if ((GasModel != "N2") && (GasModel != "AIR-5") && (GasModel != "ARGON"))
+          SU2_MPI::Error("The GAS_MODEL given as input is not valid. Choose one of the options: N2, AIR-5, ARGON.", CURRENT_FUNCTION);
         }
         break;
       case FEM_LES:
@@ -6126,7 +6134,7 @@ void CConfig::SetOutput(unsigned short val_software, unsigned short val_izone) {
         if (Kind_Upwind_Flow == SLAU2)  cout << "Simple Low-Dissipation AUSM 2 solver for the flow inviscid terms."<< endl;
         if (Kind_Upwind_Flow == FDS)    cout << "Flux difference splitting (FDS) upwind scheme for the flow inviscid terms."<< endl;
         if (Kind_Upwind_Flow == AUSMPLUSUP)  cout << "AUSM+-up solver for the flow inviscid terms."<< endl;
-        if (Kind_Upwind_Flow == AUSMPLUSUP2)  cout << "AUSM+-up2 solver for the flow inviscid terms."<< endl;
+        if (Kind_Upwind_Flow == AUSMPLUSUP2) cout << "AUSM+-up2 solver for the flow inviscid terms."<< endl;
         if (Kind_Upwind_Flow == AUSMPWPLUS)  cout << "AUSMPWPLUS solver for the flow inviscid terms."<< endl;
 
         if (Kind_Solver == EULER         || Kind_Solver == DISC_ADJ_EULER ||

SU2_CFD/include/CMarkerProfileReaderFVM.hpp

@@ -3,7 +3,7 @@
  * \brief Header file for the class CMarkerProfileReaderFVM.
  *        The implementations are in the <i>CMarkerProfileReaderFVM.cpp</i> file.
  * \author T. Economon
- * \version 6.2.0 "Falcon"
+ * \version 7.0.7 "Blackbird"
  *
  * The current SU2 release has been coordinated by the
  * SU2 International Developers Society <www.su2devsociety.org>

SU2_CFD/include/fluid/

@@ -38,14 +38,19 @@ class CMutationTCLib : public CNEMOGas {
 
 private:
 
-  string GasModel;
+  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. */
+
+  su2double Tref;                         /*!< \brief Reference temperature. */
 
 public:
 
   /*!
    * \brief Constructor of the class.
    */
-  CMutationTCLib(const CConfig* config);
+  CMutationTCLib(const CConfig* config, unsigned short val_nDim);
 
   /*!
    * \brief Destructor of the class.
@@ -93,7 +98,7 @@ class CMutationTCLib : public CNEMOGas {
   /*!
    * \brief Get species enthalpies.
    */
-  vector<su2double>& GetSpeciesEnthalpy(su2double val_T, su2double *val_eves) final;
+  vector<su2double>& GetSpeciesEnthalpy(su2double val_T, su2double val_Tve, su2double *val_eves) final;
 
   /*!
    * \brief Get species diffusion coefficients.
@@ -105,7 +110,6 @@ class CMutationTCLib : public CNEMOGas {
    */
   su2double GetViscosity() final;
 
-  
   /*!
    * \brief Get T-R and V-E thermal conductivities vector.
    */
@@ -114,21 +118,21 @@ class CMutationTCLib : public CNEMOGas {
   /*!
    * \brief Get translational and vibrational temperatures vector.
    */
-  vector<su2double>& GetTemperatures(vector<su2double>& rhos, su2double rhoEmix, su2double rhoEve, su2double rhoEvel) final;
+  vector<su2double>& GetTemperatures(vector<su2double>& val_rhos, su2double rhoE, su2double rhoEve, su2double rhoEvel) final;
 
   /*!
-   * \brief Get derivative of pressure w.r.t. conservative variables.
+   * \brief Get species molar mass.
    */
-  void GetdPdU(su2double *V, vector<su2double>& val_eves, su2double *val_dPdU) final;
+  vector<su2double>& GetSpeciesMolarMass() final;
 
   /*!
-   * \brief Get derivative of temperature w.r.t. conservative variables.
+   * \brief Get reference temperature.
    */
-  void GetdTdU(su2double *V, su2double *val_dTdU) final;
+  vector<su2double>& GetRefTemperature() final;
 
   /*!
-   * \brief Get derivative of vibrational temperature w.r.t. conservative variables.
+   * \brief Get species formation enthalpy.
    */
-  void GetdTvedU(su2double *V, vector<su2double>& val_eves, su2double *val_dTvedU) final;
+  vector<su2double>& GetSpeciesFormationEnthalpy() final;
 
-};
+};

SU2_CFD/include/fluid/CMutationTCLib.hpp

@@ -39,8 +39,10 @@ class CNEMOGas : public CFluidModel {
 
 protected:
 
-  bool frozen;                           /*!< \brief Indicates if mixture is frozen. */
-  
+  bool frozen,                           /*!< \brief Indicates if mixture is frozen. */
+  ionization;                            /*!< \brief Presence of charged species in gas mixture. */
+
+  string Kind_GasModel;                  /*!< \brief String gas model. */
 
   unsigned short nSpecies,               /*!< \brief Number of species in the gas mixture. */
   nHeavy,                                /*!< \brief Number of heavy particles in gas */
@@ -77,14 +79,16 @@ class CNEMOGas : public CFluidModel {
   hs,                                    /*!< \brief Species enthalpies */
   MolarFractions,                        /*!< \brief Species molar fractions */
   ws,                                    /*!< \brief Species net production rates */
-  DiffusionCoeff;                        /*!< \brief Species diffusion coefficients*/
+  DiffusionCoeff,                        /*!< \brief Species diffusion coefficients*/
+  Enthalpy_Formation,                    /*!< \brief Enthalpy of formation */
+  Ref_Temperature;                       /*!< \brief Reference temperature for thermodynamic relations */
 
 public:
 
   /*!
    * \brief Constructor of the class.
    */
-  CNEMOGas(const CConfig* config);
+  CNEMOGas(const CConfig* config, unsigned short val_nDim);
 
   /*!
    * \brief Set mixture thermodynamic state.
@@ -97,16 +101,16 @@ class CNEMOGas : public CFluidModel {
   /*!
    * \brief Set mixture thermodynamic state.
    * \param[in] P    - Pressure.
+   * \param[in] Ms   - Mass fraction of the gas.
    * \param[in] T    - Translational/Rotational temperature.
    * \param[in] Tve  - Vibrational/Electronic temperature.
    */
   void SetTDStatePTTv(su2double P, const su2double *val_massfrac, su2double val_temperature, su2double val_temperature_ve);
 
-  
   /*!
    * \brief Get species T-R specific heats at constant volume.
    */
-  virtual vector<su2double>& GetSpeciesCvTraRot() = 0;
+  virtual const vector<su2double>& GetSpeciesCvTraRot() = 0;
 
   /*!
    * \brief Get species V-E specific heats at constant volume.
@@ -136,7 +140,7 @@ class CNEMOGas : public CFluidModel {
   /*!
    * \brief Get species enthalpies.
    */
-  virtual vector<su2double>& GetSpeciesEnthalpy(su2double val_T, su2double *val_eves) = 0;
+  virtual vector<su2double>& GetSpeciesEnthalpy(su2double val_T, su2double val_Tve, su2double *val_eves) = 0;
 
   /*!
    * \brief Get species diffusion coefficients.
@@ -156,32 +160,37 @@ class CNEMOGas : public CFluidModel {
   /*!
    * \brief Get translational and vibrational temperatures vector.
    */
-  virtual vector<su2double>& GetTemperatures(vector<su2double>& rhos, su2double rhoEmix, su2double rhoEve, su2double rhoEvel) = 0;
+  virtual vector<su2double>& GetTemperatures(vector<su2double>& val_rhos, su2double rhoEmix, su2double rhoEve, su2double rhoEvel) = 0;
 
   /*!
    * \brief Get speed of sound.
    */
-  su2double GetSoundSpeed();
+  su2double ComputeSoundSpeed();
 
   /*!
    * \brief Get pressure.
    */
-  su2double GetPressure();
+  su2double ComputePressure();
+
+  /*!
+   * \brief Get gas constant.
+   */
+  su2double ComputeGasConstant();
 
   /*!
    * \brief Get derivative of pressure w.r.t. conservative variables.
    */
-  virtual void GetdPdU(su2double *V, vector<su2double>& val_eves, su2double *val_dPdU){}
+  void ComputedPdU(su2double *V, vector<su2double>& val_eves, su2double *val_dPdU);
 
   /*!
    * \brief Get derivative of temperature w.r.t. conservative variables.
    */
-  virtual void GetdTdU(su2double *V, su2double *val_dTdU){}
+  void ComputedTdU(su2double *V, su2double *val_dTdU);
 
   /*!
    * \brief Get derivative of vibrational temperature w.r.t. conservative variables.
    */
-  virtual void GetdTvedU(su2double *V, vector<su2double>& val_eves, su2double *val_dTvedU){}
+  void ComputedTvedU(su2double *V, vector<su2double>& val_eves, su2double *val_dTvedU);
 
   /*!
    * \brief Set the translational temperature.
@@ -198,11 +207,6 @@ class CNEMOGas : public CFluidModel {
    */
   inline void SetEves(vector<su2double>& val_eves) { eves = val_eves; }
 
-  /*!
-   * \brief Get gas constant.
-   */
-  su2double GetGasConstant();
-
   /*!
    * \brief Get rhoCvtr.
    */
@@ -221,7 +225,15 @@ class CNEMOGas : public CFluidModel {
   /*!
    * \brief Get species molar mass.
    */
-  inline vector<su2double>& GetMolarMass() { return MolarMass; }
+  virtual const vector<su2double>& GetSpeciesMolarMass() = 0;
 
+  /*!
+   * \brief Get reference temperature.
+   */
+  virtual const vector<su2double>& GetRefTemperature() = 0;
 
-};
+  /*!
+   * \brief Get species formation enthalpy.
+   */
+  virtual const vector<su2double>& GetSpeciesFormationEnthalpy() = 0;  
+};