[WIP] Implementation of eN transition model

Common/include/CConfig.hpp

@@ -850,6 +850,7 @@ class CConfig {
   Tke_FreeStream,                  /*!< \brief Total turbulent kinetic energy of the fluid.  */
   Intermittency_FreeStream,        /*!< \brief Freestream intermittency (for sagt transition model) of the fluid.  */
   TurbulenceIntensity_FreeStream,  /*!< \brief Freestream turbulent intensity (for sagt transition model) of the fluid.  */
+  AmplificationFactor_FreeStream,  /*!< \brief Freestream amplifictation factor for the eN 1 equation transition model.  */
   Turb2LamViscRatio_FreeStream,    /*!< \brief Ratio of turbulent to laminar viscosity. */
   NuFactor_FreeStream,             /*!< \brief Ratio of turbulent to laminar viscosity. */
   NuFactor_Engine,                 /*!< \brief Ratio of turbulent to laminar viscosity at the engine. */
@@ -1934,6 +1935,12 @@ class CConfig {
    * \return Non-dimensionalized freestream intermittency.
    */
   su2double GetIntermittency_FreeStream(void) const { return Intermittency_FreeStream; }
+
+  /*!
+   * \brief Get the value of the non-dimensionalized freestream amplification factor.
+   * \return Non-dimensionalized freestream amplification factor.
+   */
+  su2double GetAmplificationFactor_FreeStream(void) const { return AmplificationFactor_FreeStream; }
 
   /*!
    * \brief Get the value of the non-dimensionalized freestream turbulence intensity.

Common/include/option_structure.hpp

@@ -1166,10 +1166,12 @@ inline SA_ParsedOptions ParseSAOptions(const SA_OPTIONS *SA_Options, unsigned sh
 enum class TURB_TRANS_MODEL {
   NONE,  /*!< \brief No transition model. */
   LM,    /*!< \brief Kind of transition model (Langtry-Menter (LM) for SST and Spalart-Allmaras). */
+  EN,	 /*!< \brief Kind of transition model using Amplification factor purely for Spalart-Allmaras*/
 };
 static const MapType<std::string, TURB_TRANS_MODEL> Trans_Model_Map = {
   MakePair("NONE", TURB_TRANS_MODEL::NONE)
   MakePair("LM", TURB_TRANS_MODEL::LM)
+  MakePair("EN", TURB_TRANS_MODEL::EN)
 };
 
 /*!

Common/src/CConfig.cpp

@@ -4722,6 +4722,10 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i
   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::SA) && Kind_Trans_Model == TURB_TRANS_MODEL::EN) {
+      SU2_MPI::Error("eN transition model currently only available in combination with SA (basic/FT2) turbulence model!", CURRENT_FUNCTION);
+  }
 
   if(Turb_Fixed_Values && !OptionIsSet("TURB_FIXED_VALUES_DOMAIN")){
     SU2_MPI::Error("TURB_FIXED_VALUES activated, but no domain set with TURB_FIXED_VALUES_DOMAIN.", CURRENT_FUNCTION);
@@ -6032,6 +6036,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::EN:    cout << "Low-turbulence Transition model: eN 1 equation model (2014)" << endl; break;
+		}
         cout << "Hybrid RANS/LES: ";
         switch (Kind_HybridRANSLES) {
           case NO_HYBRIDRANSLES: cout << "No Hybrid RANS/LES" << endl; break;

SU2_CFD/include/numerics/turbulent/transition/trans_convection.hpp

@@ -0,0 +1,83 @@
+/*!
+ * \file trans_convection.hpp
+ * \brief Delarations of numerics classes for discretization of
+ *        convective fluxes in transition problems.
+ * \author R. Roos
+ * \version 7.4.0 "Blackbird"
+ *
+ * SU2 Project Website: https://su2code.github.io
+ *
+ * The SU2 Project is maintained by the SU2 Foundation
+ * (http://su2foundation.org)
+ *
+ * Copyright 2012-2022, SU2 Contributors (cf. AUTHORS.md)
+ *
+ * SU2 is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2.1 of the License, or (at your option) any later version.
+ *
+ * SU2 is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with SU2. If not, see <http://www.gnu.org/licenses/>.
+ */
+
+#pragma once
+
+#include "../../scalar/scalar_convection.hpp"
+
+/*!
+ * \class CUpwSca_TransEN
+ * \brief Class for doing a scalar upwind solver for the e^N transition model equations.
+ * \ingroup ConvDiscr
+ * \author R. Roos.
+ */
+template <class FlowIndices>
+class CUpwSca_TransEN final : public CUpwScalar<FlowIndices> {
+private:
+  using Base = CUpwScalar<FlowIndices>;
+  using Base::nDim;
+  using Base::V_i;
+  using Base::V_j;
+  using Base::a0;
+  using Base::a1;
+  using Base::Flux;
+  using Base::Jacobian_i;
+  using Base::Jacobian_j;
+  using Base::ScalarVar_i;
+  using Base::ScalarVar_j;
+  using Base::implicit;
+  using Base::idx;
+
+  /*!
+   * \brief Adds any extra variables to AD.
+   */
+  void ExtraADPreaccIn() override {}
+
+  /*!
+   * \brief SA specific steps in the ComputeResidual method
+   * \param[in] config - Definition of the particular problem.
+   */
+  void FinishResidualCalc(const CConfig* config) override {
+	  Flux[0] = a0*V_i[idx.Density()]*ScalarVar_i[0] + a1*V_j[idx.Density()]*ScalarVar_j[0];
+
+    if (implicit) {
+      Jacobian_i[0][0] = a0;
+      Jacobian_j[0][0] = a1;
+    }
+  }
+
+public:
+  /*!
+   * \brief Constructor of the class.
+   * \param[in] val_nDim - Number of dimensions of the problem.
+   * \param[in] val_nVar - Number of variables of the problem.
+   * \param[in] config - Definition of the particular problem.
+   */
+  CUpwSca_TransEN(unsigned short val_nDim, unsigned short val_nVar, const CConfig* config)
+    : CUpwScalar<FlowIndices>(val_nDim, val_nVar, config) {}
+};

SU2_CFD/include/numerics/turbulent/transition/trans_diffusion.hpp

@@ -0,0 +1,103 @@
+/*!
+ * \file trans_diffusion.hpp
+ * \brief Declarations of numerics classes for discretization of
+ *        viscous fluxes in transition problems.
+ * \author R. Roos
+ * \version 7.4.0 "Blackbird"
+ *
+ * SU2 Project Website: https://su2code.github.io
+ *
+ * The SU2 Project is maintained by the SU2 Foundation
+ * (http://su2foundation.org)
+ *
+ * Copyright 2012-2022, SU2 Contributors (cf. AUTHORS.md)
+ *
+ * SU2 is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2.1 of the License, or (at your option) any later version.
+ *
+ * SU2 is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with SU2. If not, see <http://www.gnu.org/licenses/>.
+ */
+#pragma once
+
+
+#include "../../scalar/scalar_diffusion.hpp"
+
+/*!
+ * \class CAvgGrad_TransEN
+ * \brief Class for computing viscous term using average of gradient with correction (e^N transition model).
+ * \ingroup ViscDiscr
+ * \author R. Roos
+ */
+template <class FlowIndices>
+class CAvgGrad_TransEN final : public CAvgGrad_Scalar<FlowIndices> {
+private:
+  using Base = CAvgGrad_Scalar<FlowIndices>;
+  using Base::Laminar_Viscosity_i;
+  using Base::Laminar_Viscosity_j;
+  using Base::Eddy_Viscosity_i;
+  using Base::Eddy_Viscosity_j;
+  using Base::Density_i;
+  using Base::Density_j;
+  using Base::ScalarVar_i;
+  using Base::ScalarVar_j;
+  using Base::Proj_Mean_GradScalarVar;
+  using Base::proj_vector_ij;
+  using Base::implicit;
+  using Base::Flux;
+  using Base::Jacobian_i;
+  using Base::Jacobian_j;
+
+  const su2double sigma_n = 1.0;
+
+  /*!
+   * \brief Adds any extra variables to AD
+   */
+  void ExtraADPreaccIn() override {}
+
+  /*!
+   * \brief SA specific steps in the ComputeResidual method
+   * \param[in] config - Definition of the particular problem.
+   */
+  void FinishResidualCalc(const CConfig* config) override {
+
+	/*--- Compute mean effective dynamic viscosity ---*/
+	const su2double diff_i_amplification = (Laminar_Viscosity_i + Eddy_Viscosity_i)/sigma_n;
+	const su2double diff_j_amplification = (Laminar_Viscosity_j + Eddy_Viscosity_j)/sigma_n;
+
+	const su2double diff_amplification = 0.5*(diff_i_amplification + diff_j_amplification);
+
+    Flux[0] = diff_amplification*Proj_Mean_GradScalarVar[0];
+
+    /*--- For Jacobians -> Use of TSL approx. to compute derivatives of the gradients ---*/
+
+    if (implicit) {
+      Jacobian_i[0][0] = (0.5*Proj_Mean_GradScalarVar[0]-diff_amplification*proj_vector_ij);
+      //const su2double proj_on_rho_i = proj_vector_ij/Density_i;
+      //Jacobian_i[0][0] = -diff_amplification*proj_on_rho_i;
+
+      Jacobian_j[0][0] = (0.5*Proj_Mean_GradScalarVar[0]+diff_amplification*proj_vector_ij);
+      //const su2double proj_on_rho_j = proj_vector_ij/Density_j;
+      //Jacobian_j[0][0] = diff_amplification*proj_on_rho_j;
+    }
+  }
+
+public:
+  /*!
+   * \brief Constructor of the class.
+   * \param[in] val_nDim - Number of dimensions of the problem.
+   * \param[in] val_nVar - Number of variables of the problem.
+   * \param[in] correct_grad - Whether to correct gradient for skewness.
+   * \param[in] config - Definition of the particular problem.
+   */
+  CAvgGrad_TransEN(unsigned short val_nDim, unsigned short val_nVar,
+                  bool correct_grad, const CConfig* config)
+    : CAvgGrad_Scalar<FlowIndices>(val_nDim, val_nVar, correct_grad, config) {}
+};

SU2_CFD/include/numerics/turbulent/transition/trans_sources.hpp

@@ -0,0 +1,93 @@
+/*!
+ * \file trans_sources.hpp
+ * \brief Numerics classes for integration of source terms in transition problems.
+ * \author R. Roos
+ * \version 7.4.0 "Blackbird"
+ *
+ * SU2 Project Website: https://su2code.github.io
+ *
+ * The SU2 Project is maintained by the SU2 Foundation
+ * (http://su2foundation.org)
+ *
+ * Copyright 2012-2022, SU2 Contributors (cf. AUTHORS.md)
+ *
+ * SU2 is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2.1 of the License, or (at your option) any later version.
+ *
+ * SU2 is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with SU2. If not, see <http://www.gnu.org/licenses/>.
+ */
+
+#pragma once
+#include "../../../../../Common/include/toolboxes/geometry_toolbox.hpp"
+#include "../../scalar/scalar_sources.hpp"
+
+
+/*!
+ * \class CSourcePieceWise_TranEN
+ * \brief Class for integrating the source terms of the e^N transition model equations.
+ * \ingroup SourceDiscr
+ * \author R. Roos
+ */
+template <class FlowIndices>
+class CSourcePieceWise_TransEN final : public CNumerics {
+ private:
+  const FlowIndices idx; /*!< \brief Object to manage the access to the flow primitives. */
+
+  su2double g_eff_i,
+  g_eff_j,
+  g_sep_i,
+  g_sep_j;
+
+  su2double Vorticity;
+  su2double Residual, *Jacobian_i;
+  su2double Jacobian_Buffer; /*!< \brief Static storage for the Jacobian (which needs to be pointer for return type). */
+
+ public:
+  /*!
+   * \brief Constructor of the class.
+   * \param[in] val_nDim - Number of dimensions of the problem.
+   * \param[in] val_nVar - Number of variables of the problem.
+   * \param[in] config - Definition of the particular problem.
+   */
+  CSourcePieceWise_TransEN(unsigned short val_nDim, unsigned short val_nVar, const CConfig* config)
+      : CNumerics(val_nDim, 1, config),
+        idx(val_nDim, config->GetnSpecies()) {
+
+    /*--- "Allocate" the Jacobian using the static buffer. ---*/
+	Jacobian_i = &Jacobian_Buffer;
+  }
+
+  /*!
+   * \brief Residual for source term integration.
+   * \param[in] config - Definition of the particular problem.
+   * \return A lightweight const-view (read-only) of the residual/flux and Jacobians.
+   */
+  ResidualType<> ComputeResidual(const CConfig* config) override {
+
+    Residual = 0.0;
+    Jacobian_i[0] = 0.0;
+
+    if (dist_i > 1e-10) {
+
+      /*--- Source ---*/
+      Residual = Volume;
+
+      /*--- Implicit part ---*/
+      Jacobian_i[0] *= Volume;
+    }
+
+    AD::SetPreaccOut(Residual);
+    AD::EndPreacc();
+
+    return ResidualType<>(&Residual, &Jacobian_i, nullptr);
+  }
+};
+

SU2_CFD/include/solvers/CSolver.hpp

@@ -3145,6 +3145,12 @@ class CSolver {
    * \return Value of the turbulent frequency.
    */
   inline virtual su2double GetOmega_Inf(void) const { return 0; }
+
+  /*!
+   * \brief Get value of the Amplification Factor.
+   * \return Value of the Amplification Factor.
+   */
+  inline virtual su2double GetAmplificationFactor_Inf() const { return 0; }
 
   /*!
    * \brief A virtual member.