[WIP] NEMO - New Symmetry BC

Common/include/option_structure.hpp

@@ -64,12 +64,12 @@ using MapType = CEmptyMap<T,U>;
  * \brief Different software components of SU2
  */
 enum SU2_COMPONENT {
-  SU2_CFD = 1,  /*!< \brief Running the SU2_CFD software. */
-  SU2_DEF = 2,  /*!< \brief Running the SU2_DEF software. */
-  SU2_DOT = 3,  /*!< \brief Running the SU2_DOT software. */
-  SU2_MSH = 4,  /*!< \brief Running the SU2_MSH software. */
-  SU2_GEO = 5,  /*!< \brief Running the SU2_GEO software. */
-  SU2_SOL = 6   /*!< \brief Running the SU2_SOL software. */
+  SU2_CFD = 1,   /*!< \brief Running the SU2_CFD software. */
+  SU2_DEF = 2,   /*!< \brief Running the SU2_DEF software. */
+  SU2_DOT = 3,   /*!< \brief Running the SU2_DOT software. */
+  SU2_MSH = 4,   /*!< \brief Running the SU2_MSH software. */
+  SU2_GEO = 5,   /*!< \brief Running the SU2_GEO software. */
+  SU2_SOL = 6,   /*!< \brief Running the SU2_SOL software. */
 };
 
 const unsigned int EXIT_DIVERGENCE = 2;   /*!< \brief Exit code (divergence). */

Common/src/CConfig.cpp

@@ -2782,8 +2782,8 @@ void CConfig::SetConfig_Options() {
 
   /* DESCRIPTION: Size of the edge groups colored for thread parallel edge loops (0 forces the reducer strategy). */
   addUnsignedLongOption("EDGE_COLORING_GROUP_SIZE", edgeColorGroupSize, 512);
-  /* END_CONFIG_OPTIONS */
 
+  /* END_CONFIG_OPTIONS */
 }
 
 void CConfig::SetConfig_Parsing(char case_filename[MAX_STRING_SIZE]) {
@@ -3058,12 +3058,12 @@ void CConfig::SetHeader(unsigned short val_software) const{
     cout << "|   / __| | | |_  )   Release 7.1.0 \"Blackbird\"                         |" << endl;
     cout << "|   \\__ \\ |_| |/ /                                                      |" << endl;
     switch (val_software) {
-    case SU2_CFD: cout << "|   |___/\\___//___|   Suite (Computational Fluid Dynamics Code)         |" << endl; break;
-    case SU2_DEF: cout << "|   |___/\\___//___|   Suite (Mesh Deformation Code)                     |" << endl; break;
-    case SU2_DOT: cout << "|   |___/\\___//___|   Suite (Gradient Projection Code)                  |" << endl; break;
-    case SU2_MSH: cout << "|   |___/\\___//___|   Suite (Mesh Adaptation Code)                      |" << endl; break;
-    case SU2_GEO: cout << "|   |___/\\___//___|   Suite (Geometry Definition Code)                  |" << endl; break;
-    case SU2_SOL: cout << "|   |___/\\___//___|   Suite (Solution Exporting Code)                   |" << endl; break;
+    case SU2_CFD:    cout << "|   |___/\\___//___|   Suite (Computational Fluid Dynamics Code)         |" << endl; break;
+    case SU2_DEF:    cout << "|   |___/\\___//___|   Suite (Mesh Deformation Code)                     |" << endl; break;
+    case SU2_DOT:    cout << "|   |___/\\___//___|   Suite (Gradient Projection Code)                  |" << endl; break;
+    case SU2_MSH:    cout << "|   |___/\\___//___|   Suite (Mesh Adaptation Code)                      |" << endl; break;
+    case SU2_GEO:    cout << "|   |___/\\___//___|   Suite (Geometry Definition Code)                  |" << endl; break;
+    case SU2_SOL:    cout << "|   |___/\\___//___|   Suite (Solution Exporting Code)                   |" << endl; break;
     }
 
     cout << "|                                                                       |" << endl;
@@ -3550,10 +3550,6 @@ void CConfig::SetPostprocessing(unsigned short val_software, unsigned short val_
       SU2_MPI::Error("ARGON is not working with SU2_NONEQ fluid model!", CURRENT_FUNCTION);
   }
 
-  if (GetKind_FluidModel() == MUTATIONPP && GetFrozen() == true){
-      SU2_MPI::Error("The option of FROZEN_MIXTURE is not yet working with Mutation++ support.", CURRENT_FUNCTION);
-  }
-
   if(GetBoolTurbomachinery()){
     nBlades = new su2double[nZone];
     FreeStreamTurboNormal= new su2double[3];

SU2_CFD/include/solvers/CFVMFlowSolverBase.hpp

@@ -324,7 +324,7 @@ class CFVMFlowSolverBase : public CSolver {
    * \param[in] val_marker - Surface marker where the boundary condition is applied.
    */
   inline void BC_Euler_Wall(CGeometry* geometry, CSolver** solver_container, CNumerics* conv_numerics,
-                            CNumerics* visc_numerics, CConfig* config, unsigned short val_marker) final {
+                            CNumerics* visc_numerics, CConfig* config, unsigned short val_marker) override {
     /*--- Call the equivalent symmetry plane boundary condition. ---*/
     BC_Sym_Plane(geometry, solver_container, conv_numerics, visc_numerics, config, val_marker);
   }

SU2_CFD/include/solvers/CNEMOEulerSolver.hpp

@@ -246,6 +246,18 @@ class CNEMOEulerSolver : public CFVMFlowSolverBase<CNEMOEulerVariable, COMPRESSI
   void BC_Sym_Plane(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics,
                     CNumerics *visc_numerics, CConfig *config, unsigned short val_marker) final;
 
+  /*!
+     * \brief Impose the Euler boundary condition using the residual.
+     * \param[in] geometry - Geometrical definition of the problem.
+     * \param[in] solver_container - Container vector with all the solutions.
+     * \param[in] conv_numerics - Description of the numerical method for convective terms.
+     * \param[in] visc_numerics - Description of the numerical method for viscous terms.
+     * \param[in] config - Definition of the particular problem.
+     * \param[in] val_marker - Surface marker where the boundary condition is applied.
+     */
+  void BC_Euler_Wall(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics,
+                    CNumerics *visc_numerics, CConfig *config, unsigned short val_marker) final;
+
   /*!
      * \brief Impose a subsonic inlet boundary condition.
      * \param[in] geometry - Geometrical definition of the problem.

SU2_CFD/include/variables/CNEMOEulerVariable.hpp

@@ -82,6 +82,8 @@ class CNEMOEulerVariable : public CVariable {
 
   su2double Tve_Freestream; /*!< \brief Freestream vib-el temperature. */
 
+  VectorType symmetry;      /*!< \brief Nodes in symmetry planes. */
+
 public:
 
   /*!
@@ -522,6 +524,16 @@ class CNEMOEulerVariable : public CVariable {
    */
   inline su2double GetGamma(unsigned long iPoint) { return Gamma(iPoint); }
 
+  /*!
+   * \brief Retrieves the number of symmetry planes at the specified node.
+   */  
+  inline su2double GetSymmetry(unsigned long iPoint) { return symmetry[iPoint]; }
+
+  /*!
+   * \brief Increases the number of symmetry planes at the specified node by one.
+   */  
+  inline void SetSymmetry(unsigned long iPoint) {symmetry[iPoint] += 1.0;}
+
   /*---------------------------------------*/
   /*---           NEMO indices          ---*/
   /*---------------------------------------*/

SU2_CFD/src/fluid/CMutationTCLib.cpp

@@ -49,6 +49,7 @@ CMutationTCLib::CMutationTCLib(const CConfig* config, unsigned short val_nDim):
     transport_model = "Gupta-Yos";
 
   opt.setStateModel("ChemNonEqTTv");
+  if (frozen) opt.setMechanism("none");
   opt.setViscosityAlgorithm(transport_model);
   opt.setThermalConductivityAlgorithm(transport_model); 
 

SU2_CFD/src/integration/CIntegration.cpp

@@ -138,7 +138,8 @@ void CIntegration::Space_Integration(CGeometry *geometry,
         solver_container[MainSolver]->BC_Far_Field(geometry, solver_container, conv_bound_numerics, visc_bound_numerics, config, iMarker);
         break;
       case SYMMETRY_PLANE:
-        solver_container[MainSolver]->BC_Sym_Plane(geometry, solver_container, conv_bound_numerics, visc_bound_numerics, config, iMarker);
+        if (!config->GetNEMOProblem())
+          solver_container[MainSolver]->BC_Sym_Plane(geometry, solver_container, conv_bound_numerics, visc_bound_numerics, config, iMarker);
         break;
       case ELECTRODE_BOUNDARY:
         solver_container[MainSolver]->BC_Electrode(geometry, solver_container, conv_bound_numerics, config, iMarker);
@@ -184,6 +185,14 @@ void CIntegration::Space_Integration(CGeometry *geometry,
     solver_container[MainSolver]->BC_Periodic(geometry, solver_container, conv_bound_numerics, config);
   }
 
+  /*--- Placing Symmetry Plane BC last, so we only require to double
+  the residuals to use a ghost nodes approach. --- */
+
+  for (iMarker = 0; iMarker < config->GetnMarker_All(); iMarker++) {
+    KindBC = config->GetMarker_All_KindBC(iMarker);
+      if (KindBC == SYMMETRY_PLANE && config->GetNEMOProblem()) solver_container[MainSolver]->BC_Sym_Plane(geometry, solver_container, conv_bound_numerics, visc_bound_numerics, config, iMarker);
+  }
+
 }
 
 void CIntegration::Time_Integration(CGeometry *geometry, CSolver **solver_container, CConfig *config,

SU2_CFD/src/output/CNEMOCompOutput.cpp

@@ -559,7 +559,7 @@ void CNEMOCompOutput::LoadVolumeData(CConfig *config, CGeometry *geometry, CSolv
 
 void CNEMOCompOutput::LoadSurfaceData(CConfig *config, CGeometry *geometry, CSolver **solver, unsigned long iPoint, unsigned short iMarker, unsigned long iVertex){
 
-  if ((config->GetKind_Solver() == NEMO_NAVIER_STOKES)) {
+  if ((config->GetKind_Solver() == NEMO_NAVIER_STOKES && config->GetMarker_All_KindBC(iMarker) != EULER && config->GetMarker_All_KindBC(iMarker) != SYMMETRY_PLANE)) {
     SetVolumeOutputValue("SKIN_FRICTION-X", iPoint, solver[FLOW_SOL]->GetCSkinFriction(iMarker, iVertex, 0));
     SetVolumeOutputValue("SKIN_FRICTION-Y", iPoint, solver[FLOW_SOL]->GetCSkinFriction(iMarker, iVertex, 1));
     if (nDim == 3)