clean up and regression update

Author: WallyMaier

Common/include/CConfig.hpp

@@ -1657,31 +1657,31 @@ class CConfig {
   su2double GetHeat_Flux_Ref(void) const { return Heat_Flux_Ref; }
 
   /*!
-   * \brief Get the value of the frestream temperature.
+   * \brief Get the value of the freestream temperature.
    * \return Freestream temperature.
    */
   su2double GetTemperature_FreeStream(void) const { return Temperature_FreeStream; }
   /*!
-   * \brief Get the value of the frestream vibrational-electronic temperature.
-   * \return Freestream temperature.
+   * \brief Get the value of the freestream vibrational-electronic temperature.
+   * \return Freestream vibe-el temperature.
    */
   su2double GetTemperature_ve_FreeStream(void) const { return Temperature_ve_FreeStream; }
 
   /*!
-   * \brief Get the value of the frestream temperature.
-   * \return Freestream temperature.
+   * \brief Get the value of the freestream energy.
+   * \return Freestream energy.
    */
   su2double GetEnergy_FreeStream(void) const { return Energy_FreeStream; }
 
   /*!
-   * \brief Get the value of the frestream temperature.
-   * \return Freestream temperature.
+   * \brief Get the value of the freestream viscosity.
+   * \return Freestream viscosity.
    */
   su2double GetViscosity_FreeStream(void) const { return Viscosity_FreeStream; }
 
   /*!
-   * \brief Get the value of the frestream temperature.
-   * \return Freestream temperature.
+   * \brief Get the value of the freestream density.
+   * \return Freestream density.
    */
   su2double GetDensity_FreeStream(void) const { return Density_FreeStream; }
 
@@ -1759,8 +1759,8 @@ class CConfig {
   su2double GetPressure_Ref(void) const { return Pressure_Ref; }
 
   /*!
-   * \brief Get the value of the reference pressure for non-dimensionalization.
-   * \return Reference pressure for non-dimensionalization.
+   * \brief Get the value of the reference energy for non-dimensionalization.
+   * \return Reference energy for non-dimensionalization.
    */
   su2double GetEnergy_Ref(void) const { return Energy_Ref; }
 

SU2_CFD/src/output/CNEMOCompOutput.cpp

@@ -224,12 +224,12 @@ void CNEMOCompOutput::SetVolumeOutputFields(CConfig *config){
   AddVolumeOutput("ENERGY",       "Energy",     "SOLUTION", "Energy");
   AddVolumeOutput("ENERGY_VE",    "Energy_ve",  "SOLUTION", "Energy_ve");
 
+  SetVolumeOutputFields_ScalarSolution(config);
+
   //Auxiliary variables for post-processment
   for(iSpecies = 0; iSpecies < nSpecies; iSpecies++)
     AddVolumeOutput("MASSFRAC_" + std::to_string(iSpecies),  "MassFrac_" + std::to_string(iSpecies),  "AUXILIARY", "MassFrac_" + std::to_string(iSpecies));
 
-  SetVolumeOutputFields_ScalarSolution(config);
-
   // Grid velocity
   if (gridMovement){
     AddVolumeOutput("GRID_VELOCITY-X", "Grid_Velocity_x", "GRID_VELOCITY", "x-component of the grid velocity vector");

SU2_CFD/src/variables/CNEMOEulerVariable.cpp

@@ -46,6 +46,10 @@ CNEMOEulerVariable::CNEMOEulerVariable(su2double val_pressure,
 
   unsigned short iDim, iSpecies;
 
+  const bool dual_time = (config->GetTime_Marching() == TIME_MARCHING::DT_STEPPING_1ST) ||
+                         (config->GetTime_Marching() == TIME_MARCHING::DT_STEPPING_2ND);
+  const bool classical_rk4 = (config->GetKind_TimeIntScheme_Flow() == CLASSICAL_RK4_EXPLICIT);
+
   /*--- Setting variable amounts ---*/
 
   nSpecies        = config->GetnSpecies();
@@ -109,6 +113,16 @@ CNEMOEulerVariable::CNEMOEulerVariable(su2double val_pressure,
   }
 
   Solution_Old = Solution;
+
+  if (classical_rk4) Solution_New = Solution;
+
+  /*--- Allocate and initializate solution for dual time strategy ---*/
+
+  if (dual_time) {
+    Solution_time_n = Solution;
+    Solution_time_n1 = Solution;
+  }
+
 }
 
 void CNEMOEulerVariable::SetVelocity2(unsigned long iPoint) {

TestCases/parallel_regression.py

@@ -54,6 +54,18 @@ def main():
     thermalbath.tol       = 0.00001
     test_list.append(thermalbath)
 
+    # Adiabatic thermal bath
+    ionized           = TestCase('ionized')
+    ionized.cfg_dir   = "nonequilibrium/thermalbath/finitechemistry"
+    ionized.cfg_file  = "weakly_ionized.cfg"
+    ionized.test_iter = 10
+    ionized.test_vals = [0.945997, 0.945997, -12.039262, -12.171767, -32.000000, 10.013239, 0.0, 0.0]
+    ionized.su2_exec  = "mpirun -n 2 SU2_CFD"
+    ionized.timeout   = 1600
+    ionized.new_output = True
+    ionized.tol       = 0.00001
+    test_list.append(ionized)
+
     # Adiabatic frozen thermal bath
     thermalbath_frozen           = TestCase('thermalbath_frozen')
     thermalbath_frozen.cfg_dir   = "nonequilibrium/thermalbath/frozen"