[WIP] Implementation of Simplified LM Transition model

Common/include/option_structure.hpp

@@ -1174,6 +1174,7 @@ static const MapType<std::string, TURB_TRANS_MODEL> Trans_Model_Map = {
 enum class LM_OPTIONS {
   NONE,         /*!< \brief No option / default. */
   LM2015,       /*!< \brief Cross-flow corrections. */
+  SLM,          /*!< \brief Simplified version. */
   MALAN,        /*!< \brief Kind of transition correlation model (Malan). */
   SULUKSNA,     /*!< \brief Kind of transition correlation model (Suluksna). */
   KRAUSE,       /*!< \brief Kind of transition correlation model (Krause). */
@@ -1187,6 +1188,7 @@ enum class LM_OPTIONS {
 static const MapType<std::string, LM_OPTIONS> LM_Options_Map = {
   MakePair("NONE", LM_OPTIONS::NONE)
   MakePair("LM2015", LM_OPTIONS::LM2015)
+  MakePair("SLM", LM_OPTIONS::SLM)
   MakePair("MALAN", LM_OPTIONS::MALAN)
   MakePair("SULUKSNA", LM_OPTIONS::SULUKSNA)
   MakePair("KRAUSE", LM_OPTIONS::KRAUSE)
@@ -1216,6 +1218,7 @@ enum class TURB_TRANS_CORRELATION {
 struct LM_ParsedOptions {
   LM_OPTIONS version = LM_OPTIONS::NONE;  /*!< \brief LM base model. */
   bool LM2015 = false;                    /*!< \brief Use cross-flow corrections. */
+  bool SLM = false;                    /*!< \brief Use simplified version. */
   TURB_TRANS_CORRELATION Correlation = TURB_TRANS_CORRELATION::DEFAULT;
 };
 
@@ -1235,6 +1238,7 @@ inline LM_ParsedOptions ParseLMOptions(const LM_OPTIONS *LM_Options, unsigned sh
   };
 
   LMParsedOptions.LM2015 = IsPresent(LM_OPTIONS::LM2015);
+  LMParsedOptions.SLM = IsPresent(LM_OPTIONS::SLM);
 
   int NFoundCorrelations = 0;
   if (IsPresent(LM_OPTIONS::MALAN)) {

SU2_CFD/src/solvers/CSolverFactory.cpp

@@ -381,7 +381,7 @@ CSolver* CSolverFactory::CreateTransSolver(TURB_TRANS_MODEL kindTransModel, CSol
         solver[FLOW_SOL]->Preprocessing(geometry, solver, config, iMGLevel, NO_RK_ITER, RUNTIME_FLOW_SYS, false);
         transSolver->Postprocessing(geometry, solver, config, iMGLevel);
         solver[FLOW_SOL]->Preprocessing(geometry, solver, config, iMGLevel, NO_RK_ITER, RUNTIME_FLOW_SYS, false);
-        break;      
+        break;
       case TURB_TRANS_MODEL::NONE:        
         break;
     }

SU2_CFD/src/solvers/CTransLMSolver.cpp

@@ -60,11 +60,17 @@ CTransLMSolver::CTransLMSolver(CGeometry *geometry, CConfig *config, unsigned sh
 
   nDim = geometry->GetnDim();
 
-  /*--- Define variables needed for transition from config file */
+  /*--- Define variables needed for transition from config file ---*/
   options = config->GetLMParsedOptions();
   TransCorrelations.SetOptions(options);
   TurbFamily = TurbModelFamily(config->GetKind_Turb_Model());
 
+  /*--- Check if Simplified version is used ---*/
+  if (options.SLM) {
+    nVar = 1;
+    nPrimVar = 1;
+  }
+
   /*--- Single grid simulation ---*/
 
   if (iMesh == MESH_0) {
@@ -101,33 +107,41 @@ CTransLMSolver::CTransLMSolver(CGeometry *geometry, CConfig *config, unsigned sh
   lowerlimit[0] = 1.0e-4;
   upperlimit[0] = 5.0;
 
-  lowerlimit[1] = 1.0e-4;
-  upperlimit[1] = 1.0e15;
+  if (!options.SLM) {
+    lowerlimit[1] = 1.0e-4;
+    upperlimit[1] = 1.0e15;
+  }
 
   /*--- Far-field flow state quantities and initialization. ---*/
   const su2double Intensity = config->GetTurbulenceIntensity_FreeStream()*100.0;
 
   const su2double Intermittency_Inf  = 1.0;
-  su2double ReThetaT_Inf = 100.0;
+  Solution_Inf[0] = Intermittency_Inf;
 
-  /*--- Momentum thickness Reynolds number, initialized from freestream turbulent intensity*/
-  if (Intensity <= 1.3) {
-    if(Intensity >=0.027) {
-      ReThetaT_Inf = (1173.51-589.428*Intensity+0.2196/(Intensity*Intensity));
+  su2double ReThetaT_Inf = 0.0;
+
+  if (!options.SLM) {
+    su2double ReThetaT_Inf = 100.0;
+
+    /*--- Momentum thickness Reynolds number, initialized from freestream turbulent intensity*/
+    if (Intensity <= 1.3) {
+      if(Intensity >=0.027) {
+        ReThetaT_Inf = (1173.51-589.428*Intensity+0.2196/(Intensity*Intensity));
+      }
+      else {
+        ReThetaT_Inf = (1173.51-589.428*Intensity+0.2196/(0.27*0.27));
+      }
     }
-    else {
-      ReThetaT_Inf = (1173.51-589.428*Intensity+0.2196/(0.27*0.27));
+    else if(Intensity>1.3) {
+      ReThetaT_Inf = 331.5*pow(Intensity-0.5658,-0.671);
     }
-  }
-  else if(Intensity>1.3) {
-    ReThetaT_Inf = 331.5*pow(Intensity-0.5658,-0.671);
-  }
 
-  Solution_Inf[0] = Intermittency_Inf;
-  Solution_Inf[1] = ReThetaT_Inf;
+    Solution_Inf[1] = ReThetaT_Inf;
+  }
 
   /*--- Initialize the solution to the far-field state everywhere. ---*/
   nodes = new CTransLMVariable(Intermittency_Inf, ReThetaT_Inf, 1.0, 1.0, nPoint, nDim, nVar, config);
+
   SetBaseClassPointerToNodes();
 
   /*--- MPI solution ---*/

SU2_CFD/src/variables/CTransLMVariable.cpp

@@ -31,10 +31,13 @@
 CTransLMVariable::CTransLMVariable(su2double Intermittency, su2double ReThetaT, su2double gammaSep, su2double gammaEff, unsigned long npoint, unsigned long ndim, unsigned long nvar, CConfig *config) 
   : CTurbVariable(npoint, ndim, nvar, config) {
 
+  LM_ParsedOptions options = config->GetLMParsedOptions();
+
   for(unsigned long iPoint=0; iPoint<nPoint; ++iPoint)
   {
     Solution(iPoint,0) = Intermittency;
-    Solution(iPoint,1) = ReThetaT;
+    if(!options.SLM)
+      Solution(iPoint,1) = ReThetaT;
   }
 
   Solution_Old = Solution;