BC inlet for SST solver based on local conditions at the Inlet

Common/src/CConfig.cpp

@@ -1222,7 +1222,7 @@ void CConfig::SetConfig_Options() {
   addDoubleListOption("MU_T_REF", nMu_Temperature_Ref, Mu_Temperature_Ref);
   /* DESCRIPTION: Sutherland constant, default value for AIR SI */
   addDoubleListOption("SUTHERLAND_CONSTANT", nMu_S, Mu_S);
-  
+
   /*--- Options related to Viscosity Model ---*/
   /*!\brief MIXINGVISCOSITY_MODEL \n DESCRIPTION: Mixing model of the viscosity \n OPTIONS: See \link ViscosityModel_Map \endlink \n DEFAULT: DAVIDSON \ingroup Config*/
   addEnumOption("MIXING_VISCOSITY_MODEL", Kind_MixingViscosityModel, MixingViscosityModel_Map, MIXINGVISCOSITYMODEL::DAVIDSON);
@@ -3798,7 +3798,7 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i
 
     if (Kind_FluidModel == FLUID_MIXTURE) {
       /*--- Check whether the number of entries of each specified fluid property equals the number of transported scalar
-       equations solved + 1. nMolecular_Weight and nSpecific_Heat_Cp are used because it is required for the fluid mixing models. 
+       equations solved + 1. nMolecular_Weight and nSpecific_Heat_Cp are used because it is required for the fluid mixing models.
        * Cp is required in case of MIXTURE_FLUID_MODEL because the energy equation needs to be active.--- */
       if (nMolecular_Weight != nSpecies_Init + 1 || nSpecific_Heat_Cp != nSpecies_Init + 1) {
         SU2_MPI::Error(
@@ -6054,10 +6054,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::LM:    cout << "Transition model: Langtry and Menter's 4 equation model (2009)" << endl; break;        
+          case TURB_TRANS_MODEL::LM:    cout << "Transition model: Langtry and Menter's 4 equation model (2009)" << endl; break;
         }
         cout << "Hybrid RANS/LES: ";
         switch (Kind_HybridRANSLES) {

SU2_CFD/include/fluid/CFluidModel.hpp

@@ -83,7 +83,7 @@ class CFluidModel {
    * \brief Instantiate the right type of conductivity model based on config.
    */
   static unique_ptr<CConductivityModel> MakeThermalConductivityModel(const CConfig* config, unsigned short iSpecies);
-  
+
   /*!
    * \brief Instantiate the right type of mass diffusivity model based on config.
    */
@@ -167,7 +167,7 @@ class CFluidModel {
     mass_diffusivity = MassDiffusivity->GetDiffusivity();
     return mass_diffusivity;
   }
-  
+
   /*!
    * \brief Get fluid pressure partial derivative.
    */
@@ -242,7 +242,7 @@ class CFluidModel {
    * \brief Set thermal conductivity model.
    */
   virtual void SetThermalConductivityModel(const CConfig* config);
-  
+
   /*!
    * \brief Set mass diffusivity model.
    */

SU2_CFD/include/variables/CEulerVariable.hpp

@@ -39,7 +39,7 @@
  * \author F. Palacios, T. Economon
  */
 class CEulerVariable : public CFlowVariable {
-public:
+ public:
   static constexpr size_t MAXNVAR = 12;
 
   template <class IndexType>

SU2_CFD/include/variables/CPrimitiveIndices.hpp

@@ -0,0 +1,92 @@
+/*!
+ * \file CPrimitiveIndices.hpp
+ * \brief Abstract representation of flow primitive variable indices that tries to be efficient.
+ * \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 <cstdint>
+#include <cstring>
+#include "CEulerVariable.hpp"
+#include "CIncEulerVariable.hpp"
+#include "CNEMOEulerVariable.hpp"
+
+/*!
+ * \class CPrimitiveIndices
+ * \brief Abstract wrapper for the CIndices classes of CEulerVariable, CIncEulerVariable, etc.
+ * This should only be used when the concrete type of solver variable is not known.
+ * \note We are using some tricks to try to make this more efficient than virtual functions.
+ */
+template <class IndexType>
+struct CPrimitiveIndices {
+  CPrimitiveIndices(bool incompressible, bool nemo, IndexType nDim, IndexType nSpecies) {
+    /*--- Instantiate the right type of CIndices, without heap allocations. ---*/
+    if (incompressible) {
+      type_ = 0;
+      Construct<CIncEulerVariable::template CIndices<IndexType>>(nDim, nSpecies);
+    } else if (nemo) {
+      type_ = 1;
+      Construct<CNEMOEulerVariable::template CIndices<IndexType>>(nDim, nSpecies);
+    } else {
+      type_ = 2;
+      Construct<CEulerVariable::template CIndices<IndexType>>(nDim, nSpecies);
+    }
+  }
+
+  /*--- To use the CIndices object we stored on the buffer, we cast to the
+   * concrete type using the same logic used during construction. ---*/
+#define GET_INDEX_IMPL(NAME)                                                                            \
+  inline IndexType NAME() const {                                                                       \
+    switch (type_) {                                                                                    \
+     case 0:                                                                                            \
+      return reinterpret_cast<const CIncEulerVariable::template CIndices<IndexType>*>(data_)->NAME();   \
+     case 1:                                                                                            \
+      return reinterpret_cast<const CNEMOEulerVariable::template CIndices<IndexType>*>(data_)->NAME();  \
+     default:                                                                                           \
+      return reinterpret_cast<const CEulerVariable::template CIndices<IndexType>*>(data_)->NAME();      \
+    }                                                                                                   \
+  }
+  GET_INDEX_IMPL(Temperature)
+  GET_INDEX_IMPL(Velocity)
+  GET_INDEX_IMPL(Pressure)
+  GET_INDEX_IMPL(Density)
+  GET_INDEX_IMPL(Enthalpy)
+  GET_INDEX_IMPL(SoundSpeed)
+  GET_INDEX_IMPL(LaminarViscosity)
+  GET_INDEX_IMPL(EddyViscosity)
+  GET_INDEX_IMPL(ThermalConductivity)
+  GET_INDEX_IMPL(CpTotal)
+  GET_INDEX_IMPL(Temperature_ve)
+
+ private:
+  template <class ConcreteIndices>
+  void Construct(IndexType nDim, IndexType nSpecies) {
+    /*--- Build the indices object in the static buffer owned by this class. ---*/
+    static_assert(sizeof(ConcreteIndices) <= 2 * sizeof(IndexType), "");
+    new(data_) ConcreteIndices(nDim, nSpecies);
+  }
+
+  alignas(sizeof(IndexType)) char data_[2 * sizeof(IndexType)]{};  /*!< \brief Space for the largest CIndices class. */
+  uint8_t type_;
+};

SU2_CFD/src/solvers/CTurbSSTSolver.cpp

@@ -28,6 +28,7 @@
 #include "../../include/solvers/CTurbSSTSolver.hpp"
 #include "../../include/variables/CTurbSSTVariable.hpp"
 #include "../../include/variables/CFlowVariable.hpp"
+#include "../../include/variables/CPrimitiveIndices.hpp"
 #include "../../../Common/include/parallelization/omp_structure.hpp"
 #include "../../../Common/include/toolboxes/geometry_toolbox.hpp"
 
@@ -322,8 +323,8 @@ void CTurbSSTSolver::Source_Residual(CGeometry *geometry, CSolver **solver_conta
     numerics->SetCrossDiff(nodes->GetCrossDiff(iPoint));
 
     /*--- Effective Intermittency ---*/
-    
-    if (TURB_TRANS_MODEL::LM == config->GetKind_Trans_Model()) {      
+
+    if (TURB_TRANS_MODEL::LM == config->GetKind_Trans_Model()) {
       numerics->SetIntermittencyEff(solver_container[TRANS_SOL]->GetNodes()->GetIntermittencyEff(iPoint));
     }
 
@@ -578,6 +579,9 @@ void CTurbSSTSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container, C
 
   const bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
 
+  const CPrimitiveIndices<unsigned short> prim_idx(config->GetKind_Regime() == ENUM_REGIME::INCOMPRESSIBLE,
+                                                   config->GetNEMOProblem(), nDim, config->GetnSpecies());
+
   /*--- Loop over all the vertices on this boundary marker ---*/
 
   SU2_OMP_FOR_STAT(OMP_MIN_SIZE)
@@ -608,13 +612,38 @@ void CTurbSSTSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container, C
 
       conv_numerics->SetPrimitive(V_domain, V_inlet);
 
-      /*--- Non-dimensionalize Inlet_TurbVars if Inlet-Files are used. ---*/
       su2double Inlet_Vars[MAXNVAR];
-      Inlet_Vars[0] = Inlet_TurbVars[val_marker][iVertex][0];
-      Inlet_Vars[1] = Inlet_TurbVars[val_marker][iVertex][1];
       if (config->GetInlet_Profile_From_File()) {
-        Inlet_Vars[0] /= pow(config->GetVelocity_Ref(), 2);
-        Inlet_Vars[1] *= config->GetViscosity_Ref() / (config->GetDensity_Ref() * pow(config->GetVelocity_Ref(), 2));
+        /*--- Non-dimensionalize Inlet_TurbVars if Inlet-Files are used. ---*/
+        Inlet_Vars[0] = Inlet_TurbVars[val_marker][iVertex][0] / pow(config->GetVelocity_Ref(), 2);
+        Inlet_Vars[1] = Inlet_TurbVars[val_marker][iVertex][1] * config->GetViscosity_Ref() /
+                        (config->GetDensity_Ref() * pow(config->GetVelocity_Ref(), 2));
+      } else {
+        /*--- Obtain fluid model for computing the  kine and omega to impose at the inlet boundary. ---*/
+        CFluidModel* FluidModel = solver_container[FLOW_SOL]->GetFluidModel();
+
+        /*--- Obtain flow velocity vector at inlet boundary node ---*/
+
+        const su2double* Velocity_Inlet = &V_inlet[prim_idx.Velocity()];
+        su2double Density_Inlet;
+        if (config->GetKind_Regime() == ENUM_REGIME::COMPRESSIBLE) {
+          Density_Inlet = V_inlet[prim_idx.Density()];
+          FluidModel->SetTDState_Prho(V_inlet[prim_idx.Pressure()], Density_Inlet);
+        } else {
+          const su2double* Scalar_Inlet = nullptr;
+          if (config->GetKind_Species_Model() != SPECIES_MODEL::NONE) {
+            Scalar_Inlet = config->GetInlet_SpeciesVal(config->GetMarker_All_TagBound(val_marker));
+          }
+          FluidModel->SetTDState_T(V_inlet[prim_idx.Temperature()], Scalar_Inlet);
+          Density_Inlet = FluidModel->GetDensity();
+        }
+        const su2double Laminar_Viscosity_Inlet = FluidModel->GetLaminarViscosity();
+        const su2double Intensity = config->GetTurbulenceIntensity_FreeStream();
+        const su2double viscRatio = config->GetTurb2LamViscRatio_FreeStream();
+        const su2double VelMag2 = GeometryToolbox::SquaredNorm(nDim, Velocity_Inlet);
+
+        Inlet_Vars[0] = 3.0 / 2.0 * (VelMag2 * pow(Intensity, 2));
+        Inlet_Vars[1] = Density_Inlet * Inlet_Vars[0] / (Laminar_Viscosity_Inlet * viscRatio);
       }
 
       /*--- Set the turbulent variable states. Use free-stream SST

TestCases/hybrid_regression.py

@@ -254,7 +254,7 @@ def main():
     axi_rans_air_nozzle_restart.cfg_dir   = "axisymmetric_rans/air_nozzle"
     axi_rans_air_nozzle_restart.cfg_file  = "air_nozzle_restart.cfg"
     axi_rans_air_nozzle_restart.test_iter = 10
-    axi_rans_air_nozzle_restart.test_vals         = [-12.094405, -6.649212, -8.899548, -2.491859, -1924.900000]
+    axi_rans_air_nozzle_restart.test_vals         = [-12.085749, -7.439410, -8.772091, -4.045091, -1924.800000]
     axi_rans_air_nozzle_restart.test_vals_aarch64 = [-12.093539, -6.630357, -8.798732, -2.399130, -1938.200000]
     test_list.append(axi_rans_air_nozzle_restart)
 
@@ -633,7 +633,7 @@ def main():
     bars_SST_2D.cfg_dir   = "sliding_interface/bars_SST_2D"
     bars_SST_2D.cfg_file  = "bars.cfg"
     bars_SST_2D.test_iter = 13
-    bars_SST_2D.test_vals = [13.000000, -0.609170, -1.523885]
+    bars_SST_2D.test_vals = [13.000000, -0.604409, -1.523885]
     bars_SST_2D.multizone = True
     test_list.append(bars_SST_2D)
 

TestCases/parallel_regression.py

@@ -394,7 +394,7 @@ def main():
     axi_rans_air_nozzle_restart.cfg_dir   = "axisymmetric_rans/air_nozzle"
     axi_rans_air_nozzle_restart.cfg_file  = "air_nozzle_restart.cfg"
     axi_rans_air_nozzle_restart.test_iter = 10
-    axi_rans_air_nozzle_restart.test_vals = [-12.098072, -6.657157, -8.872565, -2.491820, -1924.900000]
+    axi_rans_air_nozzle_restart.test_vals = [-12.089268, -7.493381, -8.716391, -4.021218, -1924.800000]
     axi_rans_air_nozzle_restart.tol       = 0.0001
     test_list.append(axi_rans_air_nozzle_restart)
 
@@ -919,7 +919,7 @@ def main():
     coolprop_nozzle.cfg_dir   = "nicf/coolprop"
     coolprop_nozzle.cfg_file  = "coolprop_nozzle.cfg"
     coolprop_nozzle.test_iter = 10
-    coolprop_nozzle.test_vals = [-4.692515, -3.779307, 3.419174, 0.000000, 0.000000]
+    coolprop_nozzle.test_vals = [-4.692149, -1.965079, 3.428130, 0.000000, 0.000000]
     test_list.append(coolprop_nozzle)
 
     ######################################
@@ -1025,7 +1025,7 @@ def main():
     bars_SST_2D.cfg_dir   = "sliding_interface/bars_SST_2D"
     bars_SST_2D.cfg_file  = "bars.cfg"
     bars_SST_2D.test_iter = 13
-    bars_SST_2D.test_vals = [13.000000, -0.609170, -1.523885]
+    bars_SST_2D.test_vals = [13.000000, -0.604409, -1.523885]
     bars_SST_2D.command   = TestCase.Command(exec = "SU2_CFD")
     bars_SST_2D.multizone = True
     test_list.append(bars_SST_2D)
@@ -1244,7 +1244,7 @@ def main():
     pywrapper_unsteadyCHT.cfg_dir       = "py_wrapper/flatPlate_unsteady_CHT"
     pywrapper_unsteadyCHT.cfg_file      = "unsteady_CHT_FlatPlate_Conf.cfg"
     pywrapper_unsteadyCHT.test_iter     = 5
-    pywrapper_unsteadyCHT.test_vals     = [-1.614167, 2.245726, -0.001241, 0.175715]
+    pywrapper_unsteadyCHT.test_vals     = [-1.614167, 2.245725, -0.001241, 0.175713]
     pywrapper_unsteadyCHT.command       = TestCase.Command("mpirun -np 2", "python", "launch_unsteady_CHT_FlatPlate.py --parallel -f")
     pywrapper_unsteadyCHT.unsteady      = True
     pywrapper_unsteadyCHT.new_output    = True
@@ -1330,7 +1330,7 @@ def main():
     species2_primitiveVenturi_mixingmodel.cfg_dir   = "species_transport/venturi_primitive_3species"
     species2_primitiveVenturi_mixingmodel.cfg_file  = "species2_primitiveVenturi_mixingmodel.cfg"
     species2_primitiveVenturi_mixingmodel.test_iter = 50
-    species2_primitiveVenturi_mixingmodel.test_vals = [-5.476757, -4.585053, -4.579453, -5.693611, -0.069507, -5.630568, 5.000000, -1.667485, 5.000000, -4.865481, 5.000000, -1.168603, 0.000425, 0.000405, 0.000020, 0.000000]
+    species2_primitiveVenturi_mixingmodel.test_vals = [-5.477173, -4.589710, -4.582462, -5.690334, -0.069005, -5.631182, 5.000000, -1.668297, 5.000000, -4.866834, 5.000000, -1.168385, 0.000425, 0.000404, 0.000021, 0.000000]
     species2_primitiveVenturi_mixingmodel.new_output = True
     test_list.append(species2_primitiveVenturi_mixingmodel)
 
@@ -1339,7 +1339,7 @@ def main():
     species2_primitiveVenturi_mixingmodel_viscosity.cfg_dir   = "species_transport/venturi_primitive_3species"
     species2_primitiveVenturi_mixingmodel_viscosity.cfg_file  = "species2_primitiveVenturi_mixingmodel_viscosity.cfg"
     species2_primitiveVenturi_mixingmodel_viscosity.test_iter = 50
-    species2_primitiveVenturi_mixingmodel_viscosity.test_vals = [-4.421969, -4.027438, -4.346703, -5.465691, 0.436661, -4.685235, 5.000000, -1.816533, 5.000000, -5.298965, 5.000000, -1.724682, 2.287735, 0.971641, 0.607256, 0.708838]
+    species2_primitiveVenturi_mixingmodel_viscosity.test_vals = [-4.415103, -4.019010, -4.341952, -5.471507, 0.437029, -4.685683, 5.000000, -1.814481, 5.000000, -5.282862, 5.000000, -1.717346, 2.280789, 0.971585, 0.607256, 0.701947]
     species2_primitiveVenturi_mixingmodel_viscosity.new_output = True
     test_list.append(species2_primitiveVenturi_mixingmodel_viscosity)
 
@@ -1348,7 +1348,7 @@ def main():
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg_dir   = "species_transport/venturi_primitive_3species"
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg_file  = "species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg"
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_iter = 50
-    species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_vals = [-6.123576, -4.998143, -4.888374, -7.379916, 2.465109, -5.639096, 30.000000, -5.719962, 12.000000, -8.161358, 10.000000, -8.602778, 2.084606, 1.000000, 0.600000, 0.484606]
+    species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_vals = [-6.119300, -4.997334, -4.886952, -7.382047, 2.439615, -5.627780, 30.000000, -5.723010, 12.000000, -8.145296, 9.000000, -8.075817, 2.084563, 1.000000, 0.600000, 0.484563]
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2.su2_exec  = "mpirun -n 2 SU2_CFD"
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2.timeout   = 1600
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2.new_output = True
@@ -1360,7 +1360,7 @@ def main():
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg_dir   = "species_transport/venturi_primitive_3species"
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg_file  = "species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg"
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_iter = 50
-    species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_vals = [-5.729433, -5.302965, -5.193202, -8.383664, 2.160286, -5.244935, 30.000000, -5.719915, 12.000000, -8.161218, 10.000000, -8.602842, 2.084608, 1.000000, 0.600000, 0.484608]
+    species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_vals = [-5.725157, -5.302154, -5.191782, -8.385791, 2.134792, -5.233618, 30.000000, -5.722966, 12.000000, -8.145150, 9.000000, -8.075867, 2.084565, 1.000000, 0.600000, 0.484565]
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.su2_exec  = "mpirun -n 2 SU2_CFD"
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.timeout   = 1600
     species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.new_output = True

TestCases/parallel_regression_AD.py

@@ -150,7 +150,7 @@ def main():
     discadj_axisymmetric_rans_nozzle.cfg_dir    = "axisymmetric_rans/air_nozzle"
     discadj_axisymmetric_rans_nozzle.cfg_file   = "air_nozzle_restart.cfg"
     discadj_axisymmetric_rans_nozzle.test_iter  = 10
-    discadj_axisymmetric_rans_nozzle.test_vals  = [9.524990, 5.023431, 9.417118, 2.545450, 0.000000, -246129999999999991087104.000000]
+    discadj_axisymmetric_rans_nozzle.test_vals  = [9.523445, 5.023948, 7.509000, 2.841876, 0.000000, -246139999999999992659968.000000]
     discadj_axisymmetric_rans_nozzle.no_restart = True
     test_list.append(discadj_axisymmetric_rans_nozzle)
 

TestCases/serial_regression.py

@@ -264,7 +264,7 @@ def main():
     turb_wallfunction_flatplate_sst.cfg_dir   = "wallfunctions/flatplate/compressible_SST"
     turb_wallfunction_flatplate_sst.cfg_file  = "turb_SST_flatplate.cfg"
     turb_wallfunction_flatplate_sst.test_iter = 10
-    turb_wallfunction_flatplate_sst.test_vals = [-4.229952, -1.930555, -1.998476, 1.250387, -1.635646, 1.462394, 10.000000, -2.151907, 0.072873, 0.002514] #last 10 columns
+    turb_wallfunction_flatplate_sst.test_vals = [-4.229999, -1.930714, -1.998514, 1.250334, -1.635256, 1.474343, 10.000000, -2.152260, 0.072873, 0.002514] #last 10 columns
     test_list.append(turb_wallfunction_flatplate_sst)
 
     # FLAT PLATE, WALL FUNCTIONS, COMPRESSIBLE SA
@@ -353,7 +353,7 @@ def main():
     axi_rans_air_nozzle_restart.cfg_dir   = "axisymmetric_rans/air_nozzle"
     axi_rans_air_nozzle_restart.cfg_file  = "air_nozzle_restart.cfg"
     axi_rans_air_nozzle_restart.test_iter = 10
-    axi_rans_air_nozzle_restart.test_vals = [-12.098379, -6.656873, -8.869081, -2.491818, -1924.900000]
+    axi_rans_air_nozzle_restart.test_vals = [-12.088204, -7.497812, -8.714688, -4.019578, -1924.800000]
     axi_rans_air_nozzle_restart.tol       = 0.0001
     test_list.append(axi_rans_air_nozzle_restart)
 
@@ -1023,7 +1023,7 @@ def main():
     bars_SST_2D.cfg_dir   = "sliding_interface/bars_SST_2D"
     bars_SST_2D.cfg_file  = "bars.cfg"
     bars_SST_2D.test_iter = 13
-    bars_SST_2D.test_vals = [13.000000, -0.609170, -1.523885]
+    bars_SST_2D.test_vals = [13.000000, -0.604409, -1.523885]
     bars_SST_2D.multizone = True
     test_list.append(bars_SST_2D)