Merge pull request #1751 from sun5k/develop

Langtry and Menter transition model

Author: sun5k

AUTHORS.md

@@ -119,6 +119,7 @@ Samet Cakmakcioglu
 Scott Imlay
 Steffen Schotthöfer
 Steven Endres
+Sunoh Kang
 Teus van der Stelt
 Thomas D. Economon
 Tim Albring

Common/include/CConfig.hpp

@@ -848,6 +848,7 @@ class CConfig {
   Viscosity_FreeStream,            /*!< \brief Free-stream viscosity of the fluid.  */
   Tke_FreeStream,                  /*!< \brief Total turbulent kinetic energy of the fluid.  */
   Intermittency_FreeStream,        /*!< \brief Freestream intermittency (for sagt transition model) of the fluid.  */
+  ReThetaT_FreeStream,             /*!< \brief Freestream Transition Momentum Thickness Reynolds Number (for LM transition model) of the fluid.  */
   TurbulenceIntensity_FreeStream,  /*!< \brief Freestream turbulent intensity (for sagt transition model) of the fluid.  */
   Turb2LamViscRatio_FreeStream,    /*!< \brief Ratio of turbulent to laminar viscosity. */
   NuFactor_FreeStream,             /*!< \brief Ratio of turbulent to laminar viscosity. */
@@ -1926,6 +1927,12 @@ class CConfig {
    */
   su2double GetIntermittency_FreeStream(void) const { return Intermittency_FreeStream; }
 
+  /*!
+   * \brief Get the value of the freestream momentum thickness Reynolds number.
+   * \return Freestream momentum thickness Reynolds number.
+   */
+  su2double GetReThetaT_FreeStream() const { return ReThetaT_FreeStream; }
+
   /*!
    * \brief Get the value of the non-dimensionalized freestream turbulence intensity.
    * \return Non-dimensionalized freestream intensity.
@@ -2564,6 +2571,12 @@ class CConfig {
    */
   void SetOmega_FreeStream(su2double val_omega_freestream) { Omega_FreeStream = val_omega_freestream; }
 
+  /*!
+   * \brief Set the freestream momentum thickness Reynolds number.
+   * \param[in] val_ReThetaT_freestream - Value of the freestream momentum thickness Reynolds number.
+   */
+  void SetReThetaT_FreeStream(su2double val_ReThetaT_freestream) { ReThetaT_FreeStream = val_ReThetaT_freestream; }
+
   /*!
    * \brief Set the non-dimensional freestream energy.
    * \param[in] val_energy_freestreamnd - Value of the non-dimensional freestream energy.

Common/include/option_structure.hpp

@@ -478,7 +478,7 @@ enum RUNTIME_TYPE {
   RUNTIME_ADJFEA_SYS = 30,    /*!< \brief One-physics case, the code is solving the adjoint FEA equation. */
   RUNTIME_HEAT_SYS = 21,      /*!< \brief One-physics case, the code is solving the heat equation. */
   RUNTIME_ADJHEAT_SYS = 31,   /*!< \brief One-physics case, the code is solving the adjoint heat equation. */
-  RUNTIME_TRANS_SYS = 22,     /*!< \brief One-physics case, the code is solving the turbulence model. */
+  RUNTIME_TRANS_SYS = 22,     /*!< \brief One-physics case, the code is solving the transition model. */
   RUNTIME_RADIATION_SYS = 23, /*!< \brief One-physics case, the code is solving the radiation model. */
   RUNTIME_ADJRAD_SYS = 24,    /*!< \brief One-physics case, the code is solving the adjoint radiation model. */
   RUNTIME_SPECIES_SYS = 25,   /*!< \brief One-physics case, the code is solving the species model. */

Common/src/CConfig.cpp

@@ -4729,8 +4729,8 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i
     for (int i=0; i<7; ++i) eng_cyl[i] /= 12.0;
   }
 
-  if (Kind_Trans_Model == TURB_TRANS_MODEL::LM) {
-    SU2_MPI::Error("The LM transition model is under maintenance.", CURRENT_FUNCTION);
+  if ((Kind_Turb_Model != TURB_MODEL::SST) && Kind_Trans_Model == TURB_TRANS_MODEL::LM) {
+    SU2_MPI::Error("LM transition model currently only available in combination with SST turbulence model!", CURRENT_FUNCTION);
   }
 
   if(Turb_Fixed_Values && !OptionIsSet("TURB_FIXED_VALUES_DOMAIN")){
@@ -6047,6 +6047,10 @@ void CConfig::SetOutput(SU2_COMPONENT val_software, unsigned short val_izone) {
             }
             cout << "." << endl;
             break;
+        }        
+        switch (Kind_Trans_Model) {
+          case TURB_TRANS_MODEL::NONE:  break;
+          case TURB_TRANS_MODEL::LM:    cout << "Transition model: Langtry and Menter's 4 equation model (2009)" << endl; break;        
         }
         cout << "Hybrid RANS/LES: ";
         switch (Kind_HybridRANSLES) {

SU2_CFD/include/drivers/CDriver.hpp

@@ -221,6 +221,12 @@ class CDriver {
   void InstantiateTurbulentNumerics(unsigned short nVar_Turb, int offset, const CConfig *config,
                                     const CSolver* turb_solver, CNumerics ****&numerics) const;
 
+  /*!
+   * \brief Helper to instantiate transition numerics specialized for different flow solvers.
+   */
+  template <class FlowIndices>
+  void InstantiateTransitionNumerics(unsigned short nVar_Trans, int offset, const CConfig *config,
+                                    const CSolver* trans_solver, CNumerics ****&numerics) const;
   /*!
    * \brief Helper to instantiate species transport numerics specialized for different flow solvers.
    */

SU2_CFD/include/numerics/CNumerics.hpp

@@ -82,6 +82,8 @@ class CNumerics {
   turb_ke_i,  /*!< \brief Turbulent kinetic energy at point i. */
   turb_ke_j;  /*!< \brief Turbulent kinetic energy at point j. */
   su2double
+  intermittency_eff_i;  /*!< \brief effective intermittency at point i. */
+  su2double
   Pressure_i,  /*!< \brief Pressure at point i. */
   Pressure_j;  /*!< \brief Pressure at point j. */
   su2double
@@ -698,6 +700,14 @@ class CNumerics {
    */
   virtual void SetCrossDiff(su2double val_CDkw_i) {/* empty */};
 
+  /*!
+   * \brief Set the value of the effective intermittency for the LM model.
+   * \param[in] intermittency_eff_i - Value of the effective intermittency at point i.
+   */
+  void SetIntermittencyEff(su2double val_intermittency_eff_i) {
+    intermittency_eff_i = val_intermittency_eff_i;
+  }
+
   /*!
    * \brief Set the gradient of the auxiliary variables.
    * \param[in] val_auxvar_grad_i - Gradient of the auxiliary variable at point i.