addressing PR comments, updating regression

Author: WallyMaier

SU2_CFD/src/fluid/CNEMOGas.cpp

@@ -260,10 +260,6 @@ void CNEMOGas::ComputedTdU(const su2double *V, su2double *val_dTdU){
     val_dTdU[iSpecies]   = (-ef + 0.5*v2 + Cvtrs[iSpecies]*(Ref_Temperature[iSpecies]-T)) / rhoCvtr;
   }
 
-  //if (ionization) {
-  //  SU2_MPI::Error("NEED TO IMPLEMENT dTdU for IONIZED MIX",CURRENT_FUNCTION);
-  //}
-
   /*--- Momentum derivatives ---*/
   for (iDim = 0; iDim < nDim; iDim++)
     val_dTdU[nSpecies+iDim] = -V[VEL_INDEX+iDim] / V[RHOCVTR_INDEX];

TestCases/nonequilibrium/thermalbath/finitechemistry/weakly_ionized.cfg

@@ -4,11 +4,7 @@
 % Case description: 2D 5x5 Square Test Case representing a thermal bath        %
 %                    with thermal nonequilibrium and finite-rate chemistry     %
 %                                                                              %
-%                   This is only serves as a validation case.                  %
-%                                                                              %
-%                                                                              %
-% Author: W. Maier                                                             %
-% File Version 7.3.0 "Blackbird"                                               %
+%                   This only serves as a regression case.                     %
 %                                                                              %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
@@ -23,7 +19,7 @@ GAS_COMPOSITION=(0.0, 0.767, 0.233, 0.0, 0.0, 0.0, 0.0)
 MATH_PROBLEM= DIRECT
 READ_BINARY_RESTART= NO
 
-TIME_MARCHING=TIME_STEPPING
+TIME_MARCHING=DUAL_TIME_STEPPING-2ND_ORDER
 TIME_DOMAIN= YES
 TIME_STEP= 0.000000001 % 1.e-09
 MAX_TIME= 0.001
@@ -47,7 +43,8 @@ MARKER_SYM= (left,bottom,top,right )
 %
 NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
 CFL_NUMBER= 1.0
-TIME_ITER= 11
+TIME_ITER= 3
+INNER_ITER= 10
 
 % ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
 %
@@ -58,7 +55,7 @@ LINEAR_SOLVER_ITER= 10
 
 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
 %
-CONV_NUM_METHOD_FLOW= AUSM
+CONV_NUM_METHOD_FLOW= MSW
 MUSCL_FLOW= NO
 SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
 TIME_DISCRE_FLOW= EULER_EXPLICIT
@@ -90,3 +87,4 @@ SURFACE_FILENAME= surface_flow
 SURFACE_ADJ_FILENAME= surface_adjoint
 SCREEN_WRT_FREQ_INNER= 10
 OUTPUT_WRT_FREQ= 10
+SCREEN_OUTPUT= (ITER, RMS_DENSITY_0, RMS_DENSITY_1, RMS_DENSITY_2, RMS_DENSITY_3, RMS_DENSITY_4, RMS_DENSITY_5, RMS_DENSITY_6, LIFT, DRAG)

TestCases/parallel_regression.py

@@ -59,7 +59,7 @@ def main():
     ionized.cfg_dir   = "nonequilibrium/thermalbath/finitechemistry"
     ionized.cfg_file  = "weakly_ionized.cfg"
     ionized.test_iter = 10
-    ionized.test_vals = [0.444750, 0.788316, 2.290376, -4.057969, -12.041145, -12.168935, -32.000000, 10.324556]
+    ionized.test_vals = [ -8.987250, -1.377981, -0.871645, -1.260237, -1.790041, -0.785945, -4.24928, -3.458887, 5.961210, 0.000000, 0.000000]
     ionized.su2_exec  = "mpirun -n 2 SU2_CFD"
     ionized.timeout   = 1600
     ionized.new_output = True