more outputs for debugging

Author: Bot-Enigma-0

SU2_CFD/include/numerics/CNumerics.hpp

@@ -190,7 +190,7 @@ class CNumerics {
 
   bool bounded_scalar = false;    /*!< \brief Flag for bounded scalar problem */
 
-  su2double DiffCoeff_kw[2];      /*!< \brief Storage for diffusion coefficient*/
+  su2double DiffCoeff_kw[6];      /*!< \brief Storage for diffusion coefficient*/
 
 public:
   /*!

SU2_CFD/include/variables/CTurbVariable.hpp

@@ -45,6 +45,27 @@ class CTurbVariable : public CScalarVariable {
   VectorType intermittency;         /*!< \brief Value of the intermittency for the trans. model. */
   VectorType DC_TKE;    /*!< \brief TKE Diffusion Coefficient. */ 
   VectorType DC_Omega;  /*!< \brief Omega Diffusion Coefficient. */
+  VectorType DC_OmegaT1;
+  VectorType DC_OmegaT2;
+  VectorType DC_OmegaT3;
+  VectorType DC_Omega_kc;
+  VectorType DC_Omega_so2i;
+  VectorType DC_Omega_Coi;
+  VectorType DC_Omega_Coj;
+  VectorType DC_Omega_F1i;
+  VectorType DC_Omega_F1j;
+  VectorType DC_Omega_so2j;
+  VectorType DC_Omega_mui;
+  VectorType DC_Omega_muj;
+  VectorType DC_Omega_muti;
+  VectorType DC_Omega_mutj;
+  VectorType DC_Omega_rhoi;
+  VectorType DC_Omega_rhoj;
+  VectorType DC_Omega_ki;
+  VectorType DC_Omega_kj;
+  VectorType DC_Omega_wi;
+  VectorType DC_Omega_wj;
+
   /*!
    * \brief Constructor of the class.
    * \param[in] npoint - Number of points/nodes/vertices in the domain.
@@ -109,6 +130,10 @@ class CTurbVariable : public CScalarVariable {
   inline void SetDiffCoeff_kw(unsigned long iPoint, su2double* val_DC_kw) { 
     DC_TKE(iPoint) = val_DC_kw[0];
     DC_Omega(iPoint) = val_DC_kw[1];
+    DC_OmegaT1(iPoint) = val_DC_kw[2];
+    DC_OmegaT2(iPoint) = val_DC_kw[3];
+    DC_OmegaT3(iPoint) = val_DC_kw[4];
+    DC_Omega_kc(iPoint) = val_DC_kw[5];
   }
 
   /*!
@@ -118,13 +143,46 @@ class CTurbVariable : public CScalarVariable {
   inline su2double GetTKE_DC(unsigned long iPoint) const final { 
     return DC_TKE(iPoint); 
   }  
+
   /*!
    * \brief Get the diffusion coefficient value of omega.
    * \param[in] iPoint - Point index.
    */
   inline su2double GetOmega_DC(unsigned long iPoint) const final { 
     return DC_Omega(iPoint); 
   }
+  
+  /*!
+   * \brief Get omega diffusion coefficient Term 1.
+   * \param[in] iPoint - Point index.
+   */
+  inline su2double GetOmega_DCT1(unsigned long iPoint) const final { 
+    return DC_OmegaT1(iPoint); 
+  }  
+  
+  /*!
+   * \brief Get omega diffusion coefficient Term 2.
+   * \param[in] iPoint - Point index.
+   */
+  inline su2double GetOmega_DCT2(unsigned long iPoint) const final { 
+    return DC_OmegaT2(iPoint); 
+  }
+
+  /*!
+   * \brief Get omega diffusion coefficient Term 3.
+   * \param[in] iPoint - Point index.
+   */
+  inline su2double GetOmega_DCT3(unsigned long iPoint) const final { 
+    return DC_OmegaT3(iPoint); 
+  }  
+  
+  /*!
+   * \brief Get the clip value kc used in omega diffusion coeffient.
+   * \param[in] iPoint - Point index.
+   */
+  inline su2double GetOmega_DC_kc(unsigned long iPoint) const final { 
+    return DC_Omega_kc(iPoint); 
+  }
 
   /*!
    * \brief Register eddy viscosity (muT) as Input or Output of an AD recording.

SU2_CFD/include/variables/CVariable.hpp

@@ -1672,6 +1672,10 @@ class CVariable {
   inline virtual void SetDiffCoeff_kw(unsigned long iPoint, su2double* val_DiffCoeff) const { }
   inline virtual su2double GetTKE_DC(unsigned long iPoint) const {return 0.0; }
   inline virtual su2double GetOmega_DC(unsigned long iPoint) const {return 0.0; }
+  inline virtual su2double GetOmega_DCT1(unsigned long iPoint) const {return 0.0; }
+  inline virtual su2double GetOmega_DCT2(unsigned long iPoint) const {return 0.0; }
+  inline virtual su2double GetOmega_DCT3(unsigned long iPoint) const {return 0.0; }
+  inline virtual su2double GetOmega_DC_kc(unsigned long iPoint) const {return 0.0; }
 
   /*!
    * \brief Get the value of the eddy viscosity.

SU2_CFD/src/output/CFlowOutput.cpp

@@ -1257,7 +1257,13 @@ void CFlowOutput::SetVolumeOutputFieldsScalarSolution(const CConfig* config){
       AddVolumeOutput("TKE", "Turb_Kin_Energy", "SOLUTION", "Turbulent kinetic energy");
       AddVolumeOutput("DISSIPATION", "Omega", "SOLUTION", "Rate of dissipation");
       AddVolumeOutput("DIFF_COEFF_TKE", "TKE_Diff_Coeff", "SOLUTION", "Diffusion Coefficient for TKE equation");
-      AddVolumeOutput("DIFF_COEFF_Omega", "Omega_Diff_Coeff", "SOLUTION", "Diffusion Coefficient for Omega equation");
+      AddVolumeOutput("DIFF_COEFF_Omega", "Omega_Diff_Coeff", "SOLUTION", "Toal Diffusion Coefficient for Omega equation");
+      AddVolumeOutput("DIFF_COEFF_OmegaT1", "w_DiffCoeff_T1", "SOLUTION", "Term 1 in Diffusion Coefficient for Omega equation");
+      AddVolumeOutput("DIFF_COEFF_OmegaT2", "w_DiffCoeff_T2", "SOLUTION", "Term 2 in Diffusion Coefficient for Omega equation");
+      AddVolumeOutput("DIFF_COEFF_OmegaT3", "w_DiffCoeff_T3", "SOLUTION", "Term 3 in Diffusion Coefficient for Omega equation");
+      AddVolumeOutput("DIFF_COEFF_Omegakc", "w_DiffCoeff_kc", "SOLUTION", "kc value in Diffusion Coefficient for Omega equation");
+      AddVolumeOutput("CDkw", "CDkw", "SOLUTION", "Cross-Diffusion term");
+      AddVolumeOutput("F1", "F1", "SOLUTION", "F1 blending function");
       break;
 
     case TURB_FAMILY::NONE:
@@ -1540,6 +1546,12 @@ void CFlowOutput::LoadVolumeDataScalar(const CConfig* config, const CSolver* con
       SetVolumeOutputValue("DISSIPATION", iPoint, Node_Turb->GetSolution(iPoint, 1));
       SetVolumeOutputValue("DIFF_COEFF_TKE", iPoint, Node_Turb->GetTKE_DC(iPoint));
       SetVolumeOutputValue("DIFF_COEFF_Omega", iPoint, Node_Turb->GetOmega_DC(iPoint));
+      SetVolumeOutputValue("CDkw", iPoint, Node_Turb->GetCrossDiff(iPoint));
+      SetVolumeOutputValue("F1", iPoint, Node_Turb->GetF1blending(iPoint));
+      SetVolumeOutputValue("DIFF_COEFF_OmegaT1", iPoint, Node_Turb->GetOmega_DCT1(iPoint));
+      SetVolumeOutputValue("DIFF_COEFF_OmegaT2", iPoint, Node_Turb->GetOmega_DCT2(iPoint));
+      SetVolumeOutputValue("DIFF_COEFF_OmegaT3", iPoint, Node_Turb->GetOmega_DCT3(iPoint));
+      SetVolumeOutputValue("DIFF_COEFF_Omegakc", iPoint, Node_Turb->GetOmega_DC_kc(iPoint));
       SetVolumeOutputValue("RES_TKE", iPoint, turb_solver->LinSysRes(iPoint, 0));
       SetVolumeOutputValue("RES_DISSIPATION", iPoint, turb_solver->LinSysRes(iPoint, 1));
       if (limiter) {

SU2_CFD/src/variables/CTurbVariable.cpp

@@ -38,6 +38,10 @@ CTurbVariable::CTurbVariable(unsigned long npoint, unsigned long ndim, unsigned
     if (TurbModelFamily(config->GetKind_Turb_Model()) == TURB_FAMILY::KW) {
       DC_TKE.resize(nPoint) = su2double(1.0);
       DC_Omega.resize(nPoint) = su2double(1.0);
+      DC_OmegaT1.resize(nPoint) = su2double(1.0);
+      DC_OmegaT2.resize(nPoint) = su2double(1.0);
+      DC_OmegaT3.resize(nPoint) = su2double(1.0);
+      DC_Omega_kc.resize(nPoint) = su2double(1.0);
     }
    }