Enable the use of conductivity models different from CONSTANT_PRANDTL for compressible flows.

Common/src/CConfig.cpp

@@ -4163,10 +4163,6 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i
           SU2_MPI::Error("Only SUTHERLAND viscosity model can be used with US Measurement", CURRENT_FUNCTION);
         }
       }
-      if (Kind_ConductivityModel != CONDUCTIVITYMODEL::CONSTANT_PRANDTL) {
-        SU2_MPI::Error("Only CONSTANT_PRANDTL thermal conductivity model can be used with STANDARD_AIR and IDEAL_GAS",
-                       CURRENT_FUNCTION);
-      }
     }
     /*--- Check for Boundary condition option agreement ---*/
   if (Kind_InitOption == REYNOLDS){

SU2_CFD/include/numerics/flow/flow_diffusion.hpp

@@ -355,15 +355,10 @@ class CGeneralAvgGrad_Flow final : public CAvgGrad_Base {
    * \brief Compute the heat flux due to molecular and turbulent diffusivity
    * \param[in] val_gradprimvar - Gradient of the primitive variables.
    * \param[in] val_laminar_viscosity - Laminar viscosity.
-   * \param[in] val_eddy_viscosity - Eddy viscosity.
    * \param[in] val_thermal_conductivity - Thermal Conductivity.
-   * \param[in] val_heat_capacity_cp - Heat Capacity at constant pressure.
    */
-  void SetHeatFluxVector(const su2double* const *val_gradprimvar,
-                         su2double val_laminar_viscosity,
-                         su2double val_eddy_viscosity,
-                         su2double val_thermal_conductivity,
-                         su2double val_heat_capacity_cp);
+  void SetHeatFluxVector(const su2double* const* val_gradprimvar, su2double val_laminar_viscosity,
+                         su2double val_thermal_conductivity);
 
   /*!
    * \brief Compute the Jacobian of the heat flux vector
@@ -373,17 +368,11 @@ class CGeneralAvgGrad_Flow final : public CAvgGrad_Base {
    *
    * \param[in] val_Mean_PrimVar - Mean value of the primitive variables.
    * \param[in] val_Mean_SecVar - Mean value of the secondary variables.
-   * \param[in] val_eddy_viscosity - Value of the eddy viscosity.
    * \param[in] val_thermal_conductivity - Value of the thermal conductivity.
-   * \param[in] val_heat_capacity_cp - Value of the specific heat at constant pressure.
    * \param[in] val_dist_ij - Distance between the points.
    */
-  void SetHeatFluxJacobian(const su2double *val_Mean_PrimVar,
-                           const su2double *val_Mean_SecVar,
-                           su2double val_eddy_viscosity,
-                           su2double val_thermal_conductivity,
-                           su2double val_heat_capacity_cp,
-                           su2double val_dist_ij);
+  void SetHeatFluxJacobian(const su2double* val_Mean_PrimVar, const su2double* val_Mean_SecVar,
+                           su2double val_thermal_conductivity, su2double val_dist_ij);
 
 public:
 

SU2_CFD/src/numerics/flow/flow_diffusion.cpp

@@ -797,11 +797,9 @@ CGeneralAvgGrad_Flow::CGeneralAvgGrad_Flow(unsigned short val_nDim,
 
 void CGeneralAvgGrad_Flow::SetHeatFluxVector(const su2double* const *val_gradprimvar,
                                              const su2double val_laminar_viscosity,
-                                             const su2double val_eddy_viscosity,
-                                             const su2double val_thermal_conductivity,
-                                             const su2double val_heat_capacity_cp) {
+                                             const su2double val_thermal_conductivity) {
 
-  const su2double heat_flux_factor = val_thermal_conductivity + val_heat_capacity_cp*val_eddy_viscosity/Prandtl_Turb;
+  const su2double heat_flux_factor = val_thermal_conductivity;
 
   /*--- Gradient of primitive variables -> [Temp vel_x vel_y vel_z Pressure] ---*/
   for (unsigned short iDim = 0; iDim < nDim; iDim++) {
@@ -811,9 +809,7 @@ void CGeneralAvgGrad_Flow::SetHeatFluxVector(const su2double* const *val_gradpri
 
 void CGeneralAvgGrad_Flow::SetHeatFluxJacobian(const su2double *val_Mean_PrimVar,
                                                const su2double *val_Mean_SecVar,
-                                               const su2double val_eddy_viscosity,
                                                const su2double val_thermal_conductivity,
-                                               const su2double val_heat_capacity_cp,
                                                const su2double val_dist_ij) {
   /* Viscous flux Jacobians for arbitrary equations of state */
 
@@ -835,8 +831,7 @@ void CGeneralAvgGrad_Flow::SetHeatFluxJacobian(const su2double *val_Mean_PrimVar
   su2double dTdu1= dTde_rho*(-val_Mean_PrimVar[1])*(1/rho);
   su2double dTdu2= dTde_rho*(-val_Mean_PrimVar[2])*(1/rho);
 
-  su2double total_conductivity = val_thermal_conductivity + val_heat_capacity_cp*val_eddy_viscosity/Prandtl_Turb;
-  su2double factor2 = total_conductivity/val_dist_ij;
+  su2double factor2 = val_thermal_conductivity/val_dist_ij;
 
   heat_flux_jac_i[0] = factor2*dTdu0;
   heat_flux_jac_i[1] = factor2*dTdu1;
@@ -903,7 +898,6 @@ CNumerics::ResidualType<> CGeneralAvgGrad_Flow::ComputeResidual(const CConfig* c
   Laminar_Viscosity_i = V_i[nDim+5];    Laminar_Viscosity_j = V_j[nDim+5];
   Eddy_Viscosity_i = V_i[nDim+6];       Eddy_Viscosity_j = V_j[nDim+6];
   Thermal_Conductivity_i = V_i[nDim+7]; Thermal_Conductivity_j = V_j[nDim+7];
-  Cp_i = V_i[nDim+8]; Cp_j = V_j[nDim+8];
 
   /*--- Mean secondary variables ---*/
 
@@ -917,7 +911,6 @@ CNumerics::ResidualType<> CGeneralAvgGrad_Flow::ComputeResidual(const CConfig* c
   Mean_Eddy_Viscosity       = 0.5*(Eddy_Viscosity_i + Eddy_Viscosity_j);
   Mean_turb_ke              = 0.5*(turb_ke_i + turb_ke_j);
   Mean_Thermal_Conductivity = 0.5*(Thermal_Conductivity_i + Thermal_Conductivity_j);
-  Mean_Cp                   = 0.5*(Cp_i + Cp_j);
 
   /*--- Mean gradient approximation ---*/
 
@@ -954,8 +947,7 @@ CNumerics::ResidualType<> CGeneralAvgGrad_Flow::ComputeResidual(const CConfig* c
   if (config->GetSAParsedOptions().qcr2000) AddQCR(nDim, &Mean_GradPrimVar[1], tau);
   if (Mean_TauWall > 0) AddTauWall(UnitNormal, Mean_TauWall);
 
-  SetHeatFluxVector(Mean_GradPrimVar, Mean_Laminar_Viscosity,
-                    Mean_Eddy_Viscosity, Mean_Thermal_Conductivity, Mean_Cp);
+  SetHeatFluxVector(Mean_GradPrimVar, Mean_Laminar_Viscosity, Mean_Thermal_Conductivity);
 
   GetViscousProjFlux(Mean_PrimVar, Normal);
 
@@ -975,8 +967,7 @@ CNumerics::ResidualType<> CGeneralAvgGrad_Flow::ComputeResidual(const CConfig* c
 
       SetTauJacobian(Mean_PrimVar, Mean_Laminar_Viscosity, Mean_Eddy_Viscosity, dist_ij, UnitNormal);
 
-      SetHeatFluxJacobian(Mean_PrimVar, Mean_SecVar, Mean_Eddy_Viscosity,
-                          Mean_Thermal_Conductivity, Mean_Cp, dist_ij);
+      SetHeatFluxJacobian(Mean_PrimVar, Mean_SecVar, Mean_Thermal_Conductivity, dist_ij);
 
       GetViscousProjJacs(Mean_PrimVar, Area, Proj_Flux_Tensor, Jacobian_i, Jacobian_j);
     }

SU2_CFD/src/numerics/flow/flow_sources.cpp

@@ -140,10 +140,8 @@ void CSourceAxisymmetric_Flow::ResidualDiffusion(){
   su2double laminar_viscosity_i    = V_i[nDim+5];
   su2double eddy_viscosity_i       = V_i[nDim+6];
   su2double thermal_conductivity_i = V_i[nDim+7];
-  su2double heat_capacity_cp_i     = V_i[nDim+8];
 
   su2double total_viscosity_i = laminar_viscosity_i + eddy_viscosity_i;
-  su2double total_conductivity_i = thermal_conductivity_i + heat_capacity_cp_i*eddy_viscosity_i/Prandtl_Turb;
 
   su2double u = U_i[1]/U_i[0];
   su2double v = U_i[2]/U_i[0];
@@ -156,7 +154,7 @@ void CSourceAxisymmetric_Flow::ResidualDiffusion(){
   residual[3] -= Volume*(yinv*(total_viscosity_i*(u*(PrimVar_Grad_i[2][0]+PrimVar_Grad_i[1][1])
                                                  +v*TWO3*(2*PrimVar_Grad_i[2][1]-PrimVar_Grad_i[1][0]
                                                  -v*yinv+U_i[0]*turb_ke_i))
-                                                 +total_conductivity_i*PrimVar_Grad_i[0][1])
+                                                 +thermal_conductivity_i*PrimVar_Grad_i[0][1])
                                                  -TWO3*(AuxVar_Grad_i[1][1]+AuxVar_Grad_i[2][0]));
 }
 

SU2_CFD/src/solvers/CAdjNSSolver.cpp

@@ -1641,10 +1641,8 @@ void CAdjNSSolver::BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_cont
       }
 
       /*--- Get transport coefficient information ---*/
-      Laminar_Viscosity    = solver_container[FLOW_SOL]->GetNodes()->GetLaminarViscosity(iPoint);
       Eddy_Viscosity       = solver_container[FLOW_SOL]->GetNodes()->GetEddyViscosity(iPoint);
-      Thermal_Conductivity = Cp * ( Laminar_Viscosity/Prandtl_Lam
-                                   +Eddy_Viscosity/Prandtl_Turb);
+      Thermal_Conductivity = solver_container[FLOW_SOL]-> GetNodes()->GetThermalConductivity(iPoint);
 
 //      GradV = solver_container[FLOW_SOL]->GetNodes()->GetGradient_Primitive(iPoint);
 
@@ -1660,8 +1658,7 @@ void CAdjNSSolver::BC_Isothermal_Wall(CGeometry *geometry, CSolver **solver_cont
         GradT = solver_container[FLOW_SOL]->GetNodes()->GetGradient_Primitive(iPoint)[0];
         kGTdotn = 0.0;
         for (iDim = 0; iDim < nDim; iDim++)
-          kGTdotn += Cp * Laminar_Viscosity/Prandtl_Lam*GradT[iDim]*Normal[iDim]/Area;
-        // Cp * Viscosity/Prandtl_Lam matches term used in solver_direct_mean
+          kGTdotn += (Thermal_Conductivity - Cp * Eddy_Viscosity/Prandtl_Turb) * GradT[iDim]*Normal[iDim]/Area;
         /*--- constant term to multiply max heat flux objective ---*/
         Xi = solver_container[FLOW_SOL]->GetTotal_HeatFlux(); // versions for max heat flux
         Xi = pow(Xi, 1.0/pnorm-1.0)/pnorm;

SU2_CFD/src/solvers/CNSSolver.cpp

@@ -619,10 +619,7 @@ void CNSSolver::BC_Isothermal_Wall_Generic(CGeometry *geometry, CSolver **solver
 
   const bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
   const su2double Temperature_Ref = config->GetTemperature_Ref();
-  const su2double Prandtl_Lam = config->GetPrandtl_Lam();
-  const su2double Prandtl_Turb = config->GetPrandtl_Turb();
   const su2double Gas_Constant = config->GetGas_ConstantND();
-  const su2double Cp = (Gamma / Gamma_Minus_One) * Gas_Constant;
 
   /*--- Identify the boundary and retrieve the specified wall temperature from
    the config (for non-CHT problems) as well as the wall function treatment. ---*/
@@ -689,9 +686,7 @@ void CNSSolver::BC_Isothermal_Wall_Generic(CGeometry *geometry, CSolver **solver
 
     /*--- Get transport coefficients ---*/
 
-    su2double laminar_viscosity    = nodes->GetLaminarViscosity(iPoint);
-    su2double eddy_viscosity       = nodes->GetEddyViscosity(iPoint);
-    su2double thermal_conductivity = Cp * (laminar_viscosity/Prandtl_Lam + eddy_viscosity/Prandtl_Turb);
+    su2double thermal_conductivity = nodes->GetThermalConductivity(iPoint);
 
     // work in progress on real-gases...
     //thermal_conductivity = nodes->GetThermalConductivity(iPoint);

SU2_CFD/src/variables/CNSVariable.cpp

@@ -186,9 +186,10 @@ bool CNSVariable::SetPrimVar(unsigned long iPoint, su2double eddy_visc, su2doubl
 
   SetLaminarViscosity(iPoint, FluidModel->GetLaminarViscosity());
 
-  /*--- Set eddy viscosity ---*/
+  /*--- Set eddy viscosity locally and in the fluid model.---*/
 
   SetEddyViscosity(iPoint, eddy_visc);
+  FluidModel->SetEddyViscosity(eddy_visc);
 
   /*--- Set thermal conductivity ---*/
 

TestCases/hybrid_regression.py

@@ -250,7 +250,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.070954, -7.407644, -8.698118, -4.008751, 0]
+    axi_rans_air_nozzle_restart.test_vals = [-12.070958, -7.407694, -8.697876, -4.008653, 0]
     test_list.append(axi_rans_air_nozzle_restart)
 
     #################################
@@ -573,8 +573,8 @@ def main():
     transonic_stator_restart.cfg_dir   = "turbomachinery/transonic_stator_2D"
     transonic_stator_restart.cfg_file  = "transonic_stator_restart.cfg"
     transonic_stator_restart.test_iter = 20
-    transonic_stator_restart.test_vals         = [-5.007735, -3.099310, -2.751696, 1.091966, -3.542819, 2.163237, -471630.000000, 94.866000, -0.035738]
-    transonic_stator_restart.test_vals_aarch64 = [-5.007735, -3.099310, -2.751696, 1.091966, -3.542819, 2.163237, -471630, 94.866, -0.035738]
+    transonic_stator_restart.test_vals         = [-5.007734, -3.099310, -2.751691, 1.091967, -3.542818, 2.163238, -471630.000000, 94.866000, -0.036113]
+    transonic_stator_restart.test_vals_aarch64 = [-5.007734, -3.099310, -2.751691, 1.091967, -3.542818, 2.163238, -471630, 94.866, -0.036113]
     test_list.append(transonic_stator_restart)
 
     # Multiple turbomachinery interface restart

TestCases/parallel_regression.py

@@ -392,6 +392,14 @@ def main():
     turb_flatplate_species.test_vals = [-4.243064, -0.634797, -1.706652, 1.231264, -3.266203, 9.000000, -6.632972, 5.000000, -6.985977, 10.000000, -6.007208, 0.996237, 0.996237]
     test_list.append(turb_flatplate_species)
 
+    # Flat plate (compressible) species transport using constant conductivity model
+    turb_flatplate_species_ConstConductivity = TestCase('turb_flatplate_species_ConstConductivity')
+    turb_flatplate_species_ConstConductivity.cfg_dir   = "rans/flatplate"
+    turb_flatplate_species_ConstConductivity.cfg_file  = "turb_SA_flatplate_species_ConstConductivity.cfg"
+    turb_flatplate_species_ConstConductivity.test_iter = 20
+    turb_flatplate_species_ConstConductivity.test_vals = [-4.385761, -0.753729, -2.143649, 1.114937, -4.088701, 5.000000, -2.508344, 4.000000, -5.119461, 5.000000, -2.723197, 0.999990, 0.999990]
+    test_list.append(turb_flatplate_species_ConstConductivity)
+
     # Flat plate SST compressibility correction Wilcox
     turb_flatplate_CC_Wilcox = TestCase('turb_flatplate_CC_Wilcox')
     turb_flatplate_CC_Wilcox.cfg_dir   = "rans/flatplate"
@@ -1102,8 +1110,8 @@ def main():
     transonic_stator_restart.cfg_dir   = "turbomachinery/transonic_stator_2D"
     transonic_stator_restart.cfg_file  = "transonic_stator_restart.cfg"
     transonic_stator_restart.test_iter = 20
-    transonic_stator_restart.test_vals = [-5.011834, -3.091110, -2.757795, 1.087934, -3.544707, 2.166101, -471630, 94.868, -0.035888]
-    transonic_stator_restart.test_vals_aarch64 = [-5.011834, -3.091110, -2.757795, 1.087934, -3.544707, 2.166101, -471630, 94.868, -0.035888]
+    transonic_stator_restart.test_vals = [-5.011833, -3.091100, -2.757789, 1.087935, -3.544705, 2.166101, -471630, 94.868, -0.036253]
+    transonic_stator_restart.test_vals_aarch64 = [-5.011833, -3.091100, -2.757789, 1.087935, -3.544705, 2.166101, -471630, 94.868, -0.036253]
     test_list.append(transonic_stator_restart)
 
     # Multiple turbomachinery interface restart

TestCases/parallel_regression_AD.py

@@ -230,8 +230,8 @@ def main():
     discadj_trans_stator.cfg_dir   = "disc_adj_turbomachinery/transonic_stator_2D"
     discadj_trans_stator.cfg_file  = "transonic_stator.cfg"
     discadj_trans_stator.test_iter = 79
-    discadj_trans_stator.test_vals = [79, 0.769967, 0.374299, 0.474436, -0.996528, 2.154053, -4.446085]
-    discadj_trans_stator.test_vals_aarch64 = [79, 0.769967, 0.374299, 0.474436, -0.996528, 2.154053, -4.446085]
+    discadj_trans_stator.test_vals = [79, 0.769976, 0.374299, 0.474435, -0.9965283, 2.154068, -4.446072]
+    discadj_trans_stator.test_vals_aarch64 = [79, 0.769976, 0.374299, 0.474435, -0.996523, 2.154068, -4.446072]
     test_list.append(discadj_trans_stator)
 
     ###################################

TestCases/rans/flatplate/turb_SA_flatplate_species_ConstConductivity.cfg

@@ -0,0 +1,140 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: Turbulent flow over flat plate with zero pressure gradient %
+%                   using conductivity model CONSTANT_CONDUCTIVITY             %
+% Author: C. Morales Ubal                                                      %
+% Institution: Eindhoven University of Technology                              %
+% Date: 2025.02.22                                                             %
+% File Version 8.1.0 "Harrier"                                                 %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= RANS
+KIND_TURB_MODEL= SA
+MATH_PROBLEM= DIRECT
+RESTART_SOL= NO
+
+% ----------- COMPRESSIBLE AND INCOMPRESSIBLE FREE-STREAM DEFINITION ----------%
+%
+MACH_NUMBER= 0.2
+AOA= 0.0
+SIDESLIP_ANGLE= 0.0
+FREESTREAM_TEMPERATURE= 300.0
+REYNOLDS_NUMBER= 5000000.0
+REYNOLDS_LENGTH= 1.0
+
+% ---------------------- REFERENCE VALUE DEFINITION ---------------------------%
+%
+REF_ORIGIN_MOMENT_X = 0.25
+REF_ORIGIN_MOMENT_Y = 0.00
+REF_ORIGIN_MOMENT_Z = 0.00
+REF_LENGTH= 1.0
+REF_AREA= 2.0
+
+% -------------------- FLUID PROPERTIES ------------------------------------- %
+%
+FLUID_MODEL= IDEAL_GAS
+%
+MOLECULAR_WEIGHT= 28.960
+%
+SPECIFIC_HEAT_CP = 1009.39
+%
+CONDUCTIVITY_MODEL= CONSTANT_CONDUCTIVITY
+THERMAL_CONDUCTIVITY_CONSTANT= 0.0258 
+%
+PRANDTL_LAM= 0.72
+TURBULENT_CONDUCTIVITY_MODEL= CONSTANT_PRANDTL_TURB
+PRANDTL_TURB= 0.90
+%
+% --------------------------- VISCOSITY MODEL ---------------------------------%
+%
+VISCOSITY_MODEL= CONSTANT_VISCOSITY
+%
+MU_CONSTANT= 1.8551E-05 
+%
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_HEATFLUX= ( wall, 0.0 )
+SPECIFIED_INLET_PROFILE= NO
+INLET_FILENAME= inlet.dat
+INC_INLET_TYPE=  VELOCITY_INLET
+MARKER_INLET= ( inlet, 302.4, 118309.784, 1.0, 0.0, 0.0 )
+MARKER_INLET_SPECIES= ( inlet, 1.0 )
+MARKER_OUTLET= ( outlet, 115056.0, farfield, 115056.0 )
+MARKER_SYM= ( symmetry )
+MARKER_PLOTTING= ( wall )
+MARKER_MONITORING= ( wall )
+
+% -------------------- SCALAR TRANSPORT ---------------------------------------%
+%
+KIND_SCALAR_MODEL= SPECIES_TRANSPORT
+DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
+DIFFUSIVITY_CONSTANT= 0.001
+CONV_NUM_METHOD_SPECIES= SCALAR_UPWIND
+MUSCL_SPECIES= NO
+SLOPE_LIMITER_SPECIES = NONE
+TIME_DISCRE_SPECIES= EULER_IMPLICIT
+SPECIES_INIT= 1.0
+SPECIES_CLIPPING= YES
+SPECIES_CLIPPING_MIN= 0.0
+SPECIES_CLIPPING_MAX= 1.0
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+%
+NUM_METHOD_GRAD= GREEN_GAUSS
+CFL_NUMBER= 10.0
+CFL_ADAPT= YES
+CFL_ADAPT_PARAM= ( 0.5, 2.0, 1.0, 1000.0 )
+ITER= 10000
+% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
+%
+LINEAR_SOLVER= FGMRES
+LINEAR_SOLVER_PREC= ILU
+LINEAR_SOLVER_ERROR= 1E-5 
+LINEAR_SOLVER_ITER= 5
+%
+% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
+%
+CONV_NUM_METHOD_FLOW= ROE
+MUSCL_FLOW= NO
+SLOPE_LIMITER_FLOW= NONE
+JST_SENSOR_COEFF= ( 0.5, 0.02 )
+TIME_DISCRE_FLOW= EULER_IMPLICIT
+
+% -------------------- TURBULENT NUMERICAL METHOD DEFINITION ------------------%
+%
+CONV_NUM_METHOD_TURB= SCALAR_UPWIND
+MUSCL_TURB= NO
+SLOPE_LIMITER_TURB= VENKATAKRISHNAN
+TIME_DISCRE_TURB= EULER_IMPLICIT
+
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+CONV_RESIDUAL_MINVAL= -15
+CONV_STARTITER= 10
+CONV_CAUCHY_ELEMS= 100
+CONV_CAUCHY_EPS= 1E-6
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+MESH_FILENAME= mesh_flatplate_turb_137x97.su2
+SCREEN_OUTPUT= INNER_ITER WALL_TIME \
+               RMS_DENSITY RMS_MOMENTUM-X RMS_MOMENTUM-Y RMS_ENERGY RMS_TKE RMS_DISSIPATION RMS_SPECIES_0 \
+               LINSOL_ITER LINSOL_RESIDUAL \
+               LINSOL_ITER_TURB LINSOL_RESIDUAL_TURB \
+               LINSOL_ITER_SPECIES LINSOL_RESIDUAL_SPECIES \
+               SURFACE_SPECIES_0
+SCREEN_WRT_FREQ_INNER= 10
+HISTORY_OUTPUT= ITER RMS_RES LINSOL SPECIES_COEFF SPECIES_COEFF_SURF
+CONV_FILENAME= history
+MARKER_ANALYZE= outlet
+MARKER_ANALYZE_AVERAGE= AREA
+OUTPUT_FILES= RESTART_ASCII, PARAVIEW_MULTIBLOCK
+VOLUME_OUTPUT= RESIDUAL, PRIMITIVE
+OUTPUT_WRT_FREQ= 100
+READ_BINARY_RESTART= NO
+RESTART_FILENAME= restart
+SOLUTION_FILENAME= solution
+