adding eddy diffusivity in enthalpy diffusion

Author: Cristopher-Morales

SU2_CFD/include/fluid/CFluidScalar.hpp

@@ -45,6 +45,7 @@ class CFluidScalar final : public CFluidModel {
   const su2double Ref_Temperature;        /*!< \brief Reference temperature. */
   const su2double GasConstant_Ref;        /*!< \brief Gas constant reference needed for Nondimensional problems. */
   const su2double Prandtl_Number;         /*!< \brief Prandlt number.*/
+  const su2double Schmidt_Turb_Number;         /*!< \brief Prandlt number.*/
 
   const bool wilke;
   const bool davidson;

SU2_CFD/src/fluid/CFluidScalar.cpp

@@ -48,6 +48,7 @@ CFluidScalar::CFluidScalar(su2double value_pressure_operating, const CConfig* co
       Ref_Temperature(config->GetStandard_RefTemperatureND()),
       GasConstant_Ref(config->GetGas_Constant_Ref()),
       Prandtl_Number(config->GetPrandtl_Turb()),
+      Schmidt_Turb_Number(config->GetSchmidt_Number_Turbulent()),
       wilke(config->GetKind_MixingViscosityModel() == MIXINGVISCOSITYMODEL::WILKE),
       davidson(config->GetKind_MixingViscosityModel() == MIXINGVISCOSITYMODEL::DAVIDSON) {
   if (n_species_mixture > ARRAYSIZE) {
@@ -228,6 +229,7 @@ void CFluidScalar::GetEnthalpyDiffusivity(su2double* enthalpy_diffusions) {
     const su2double enthalpy_species_i = specificHeat[iVar] * (Temperature - Ref_Temperature);
     enthalpy_diffusions[iVar] = Density * (enthalpy_species_i * massDiffusivity[iVar] -
                                            enthalpy_species_N * massDiffusivity[n_species_mixture - 1]);
+    enthalpy_diffusions[iVar] += Mu_Turb * (enthalpy_species_i - enthalpy_species_N) / Schmidt_Turb_Number;
   }
 }
 
@@ -236,6 +238,8 @@ void CFluidScalar::GetGradEnthalpyDiffusivity(su2double* grad_enthalpy_diffusion
     grad_enthalpy_diffusions[iVar] = Density *
                                (specificHeat[iVar] * massDiffusivity[iVar] -
                                 specificHeat[n_species_mixture - 1] * massDiffusivity[n_species_mixture - 1]);
+    grad_enthalpy_diffusions[iVar] +=
+        Mu_Turb * (specificHeat[iVar] - specificHeat[n_species_mixture - 1]) / Schmidt_Turb_Number;
   }
 }
 

SU2_CFD/src/solvers/CIncNSSolver.cpp

@@ -329,6 +329,7 @@ void CIncNSSolver::Compute_Enthalpy_Diffusion(unsigned long iEdge, CGeometry* ge
   su2double EnthalpyDiffusion_i[MAXNVAR_SPECIES]{0.0};
   su2double GradEnthalpyDiffusion_i[MAXNVAR_SPECIES]{0.0};
   FluidModel->SetTDState_T(nodes->GetPrimitive(iPoint)[prim_idx.Temperature()], Species_i);
+  FluidModel->SetEddyViscosity(nodes->GetPrimitive(iPoint)[prim_idx.EddyViscosity()]);
   FluidModel->GetEnthalpyDiffusivity(EnthalpyDiffusion_i);
   if (implicit) FluidModel->GetGradEnthalpyDiffusivity(GradEnthalpyDiffusion_i);
 
@@ -337,6 +338,7 @@ void CIncNSSolver::Compute_Enthalpy_Diffusion(unsigned long iEdge, CGeometry* ge
   su2double EnthalpyDiffusion_j[MAXNVAR_SPECIES]{0.0};
   su2double GradEnthalpyDiffusion_j[MAXNVAR_SPECIES]{0.0};
   FluidModel->SetTDState_T(nodes->GetPrimitive(jPoint)[prim_idx.Temperature()], Species_j);
+  FluidModel->SetEddyViscosity(nodes->GetPrimitive(jPoint)[prim_idx.EddyViscosity()]);
   FluidModel->GetEnthalpyDiffusivity(EnthalpyDiffusion_j);
   if (implicit) FluidModel->GetGradEnthalpyDiffusivity(GradEnthalpyDiffusion_j);