Radial Equilibrium equation fix for Turbo solver (#1008)

* Radial Equilibrium equation fix for Turbo solver

* fixing spelling mistakes

* adding Radius1 condition with EPS

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>

Author: auzbaig

SU2_CFD/src/solvers/CEulerSolver.cpp

@@ -7441,7 +7441,7 @@ void CEulerSolver::BC_Riemann(CGeometry *geometry, CSolver **solver_container,
 
       case STATIC_PRESSURE:
 
-        /*--- Retrieve the staic pressure for this boundary. ---*/
+        /*--- Retrieve the static pressure for this boundary. ---*/
         Pressure_e = config->GetRiemann_Var1(Marker_Tag);
         Pressure_e /= config->GetPressure_Ref();
         Density_e = Density_i;
@@ -7920,7 +7920,7 @@ void CEulerSolver::BC_TurboRiemann(CGeometry *geometry, CSolver **solver_contain
 
           case MIXING_OUT:
 
-            /*--- Retrieve the staic pressure for this boundary. ---*/
+            /*--- Retrieve the static pressure for this boundary. ---*/
             Pressure_e = ExtAveragePressure[val_marker][iSpan];
             Density_e = Density_i;
 
@@ -7936,7 +7936,7 @@ void CEulerSolver::BC_TurboRiemann(CGeometry *geometry, CSolver **solver_contain
 
           case STATIC_PRESSURE:
 
-            /*--- Retrieve the staic pressure for this boundary. ---*/
+            /*--- Retrieve the static pressure for this boundary. ---*/
             Pressure_e = config->GetRiemann_Var1(Marker_Tag);
             Pressure_e /= config->GetPressure_Ref();
             Density_e = Density_i;
@@ -7954,7 +7954,7 @@ void CEulerSolver::BC_TurboRiemann(CGeometry *geometry, CSolver **solver_contain
 
           case RADIAL_EQUILIBRIUM:
 
-            /*--- Retrieve the staic pressure for this boundary. ---*/
+            /*--- Retrieve the static pressure for this boundary. ---*/
             Pressure_e = RadialEquilibriumPressure[val_marker][iSpan];
             Density_e = Density_i;
 
@@ -12446,7 +12446,8 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
                   Radius2    = geometry->GetTurboRadius(iMarker,iSpan);
                   Radius1    = geometry->GetTurboRadius(iMarker,iSpan-1);
                   Vt2        = AverageTurboVelocity[iMarker][iSpan - 1][1]*AverageTurboVelocity[iMarker][iSpan - 1][1];
-                  RadialEquilibriumPressure[iMarker][iSpan -1] =  RadialEquilibriumPressure[iMarker][iSpan] - AverageDensity[iMarker][iSpan -1]*Vt2/Radius2*(Radius2 - Radius1);
+                  Radius1    = (Radius1 > EPS)? Radius1 : Radius2;
+                  RadialEquilibriumPressure[iMarker][iSpan -1] =  RadialEquilibriumPressure[iMarker][iSpan] - AverageDensity[iMarker][iSpan -1]*Vt2/Radius1*(Radius2 - Radius1);
                 }
               }
             }
@@ -12463,7 +12464,8 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
                   Radius2    = geometry->GetTurboRadius(iMarker,iSpan);
                   Radius1    = geometry->GetTurboRadius(iMarker,iSpan-1);
                   Vt2        = AverageTurboVelocity[iMarker][iSpan -1][1]*AverageTurboVelocity[iMarker][iSpan - 1][1];
-                  RadialEquilibriumPressure[iMarker][iSpan -1] =  RadialEquilibriumPressure[iMarker][iSpan] - AverageDensity[iMarker][iSpan -1]*Vt2/Radius2*(Radius2 - Radius1);
+                  Radius1    = (Radius1 > EPS)? Radius1 : Radius2;
+                  RadialEquilibriumPressure[iMarker][iSpan -1] =  RadialEquilibriumPressure[iMarker][iSpan] - AverageDensity[iMarker][iSpan -1]*Vt2/Radius1*(Radius2 - Radius1);
                 }
               }
             }