add option specific to turbulence for accurate jacobians

Author: Bot-Enigma-0

Common/include/CConfig.hpp

@@ -590,6 +590,7 @@ class CConfig {
   MUSCL_AdjTurb;           /*!< \brief MUSCL scheme for the adj turbulence equations.*/
   bool MUSCL_Species;      /*!< \brief MUSCL scheme for the species equations.*/
   bool Use_Accurate_Jacobians;  /*!< \brief Use numerically computed Jacobians for AUSM+up(2) and SLAU(2). */
+  bool Use_Accurate_Turb_Jacobians; /*!< \brief Use numerically computed Jacobians for standard SA turbulence model. */
   bool EulerPersson;       /*!< \brief Boolean to determine whether this is an Euler simulation with Persson shock capturing. */
   bool FSI_Problem = false,/*!< \brief Boolean to determine whether the simulation is FSI or not. */
   Multizone_Problem;       /*!< \brief Boolean to determine whether we are solving a multizone problem. */
@@ -4499,6 +4500,12 @@ class CConfig {
    */
   bool GetUse_Accurate_Jacobians(void) const { return Use_Accurate_Jacobians; }
 
+  /*!
+   * \brief Get whether to "Use Accurate Jacobians" for Standard SA turbulence model.
+   * \return yes/no.
+   */
+  bool GetUse_Accurate_Turb_Jacobians(void) const { return Use_Accurate_Turb_Jacobians; }
+
   /*!
    * \brief Get the kind of integration scheme (explicit or implicit)
    *        for the flow equations.

Common/src/CConfig.cpp

@@ -1996,6 +1996,8 @@ void CConfig::SetConfig_Options() {
   /*!\brief CONV_NUM_METHOD_TURB
    *  \n DESCRIPTION: Convective numerical method \ingroup Config*/
   addConvectOption("CONV_NUM_METHOD_TURB", Kind_ConvNumScheme_Turb, Kind_Centered_Turb, Kind_Upwind_Turb);
+  /*!\brief USE_ACCURATE_TURB_JACOBIANS \n DESCRIPTION: Use numerically computed Jacobians for Standard SA model \ingroup Config*/
+  addBoolOption("USE_ACCURATE_TURB_JACOBIANS", Use_Accurate_Turb_Jacobians, false);
 
   /*!\brief MUSCL_ADJTURB \n DESCRIPTION: Check if the MUSCL scheme should be used \ingroup Config*/
   addBoolOption("MUSCL_ADJTURB", MUSCL_AdjTurb, false);

SU2_CFD/include/numerics/turbulent/turb_diffusion.hpp

@@ -93,7 +93,7 @@ class CAvgGrad_TurbSA final : public CAvgGrad_Scalar<FlowIndices> {
 
     /*--- For Jacobians -> Use of TSL approx. to compute derivatives of the gradients ---*/
     su2double diffusion_coefficient = term_1 - term_2;
-    bool UseExactJacobians = config->GetUse_Accurate_Jacobians();
+    bool UseExactJacobians = config->GetUse_Accurate_Turb_Jacobians();
 
     if (implicit) {
       if (UseExactJacobians) {

config_template.cfg

@@ -1849,6 +1849,10 @@ CFL_REDUCTION_TURB= 1.0
 LOWER_LIMIT_K_FACTOR= 1e-15
 LOWER_LIMIT_OMEGA_FACTOR= 1e-5
 
+% Use numerically computed exact Jacobians for standard SA turbulence model
+% Slower per iteration but potentialy more stable and capable of higher CFL
+USE_ACCURATE_FLUX_JACOBIANS= NO
+%
 % --------------------- HEAT NUMERICAL METHOD DEFINITION ----------------------%
 %
 % Value of the thermal diffusivity