change default values 0.1 -> 0.01 for more guaranteed convergence

Author: bigfooted

Common/src/CConfig.cpp

@@ -1384,9 +1384,9 @@ void CConfig::SetConfig_Options() {
     /*!\brief INC_INLET_USENORMAL \n DESCRIPTION: Use the local boundary normal for the flow direction with the incompressible pressure inlet. \ingroup Config*/
   addBoolOption("INC_INLET_USENORMAL", Inc_Inlet_UseNormal, false);
   /*!\brief INC_INLET_DAMPING \n DESCRIPTION: Damping factor applied to the iterative updates to the velocity at a pressure inlet in incompressible flow (0.1 by default). \ingroup Config*/
-  addDoubleOption("INC_INLET_DAMPING", Inc_Inlet_Damping, 0.1);
+  addDoubleOption("INC_INLET_DAMPING", Inc_Inlet_Damping, 0.01);
   /*!\brief INC_OUTLET_DAMPING \n DESCRIPTION: Damping factor applied to the iterative updates to the pressure at a mass flow outlet in incompressible flow (0.1 by default). \ingroup Config*/
-  addDoubleOption("INC_OUTLET_DAMPING", Inc_Outlet_Damping, 0.1);
+  addDoubleOption("INC_OUTLET_DAMPING", Inc_Outlet_Damping, 0.01);
 
   /*--- Options related to the species solver. ---*/
 

SU2_CFD/src/solvers/CIncEulerSolver.cpp

@@ -2601,18 +2601,7 @@ void CIncEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container,
         //su2double P_domain_max = config->GetRange_Pressure(1);
         //su2double P_domain_range = P_domain_max - P_domain_min;
 
-        /*--- Do not relax when dP is relatively small compared to the pressure range dp = (P_max-P_min). ---*/
-
-        /*
-        We take the entire range of pressures and we check the relative P and the relative dp
-        1. is it sufficient to subtract the min pressure?
-        range [0..1]
-        range [-1..0]
-        range [0..10]
-        range [-10..0]
-        range [10..11]
-        range [-11..10]]
-        */ 
+        /*--- Do not relax when dP is relatively small compared to the pressure range dp = (P-P_min). ---*/
 
         if (abs(dP) < abs(Damping * (P_domain-P_domain_min)))
           Damping = 1.0;

config_template.cfg

@@ -277,7 +277,7 @@ INC_TEMPERATURE_REF = 1.0
 INC_INLET_TYPE= VELOCITY_INLET
 %
 % Damping coefficient for iterative updates at pressure inlets. (0.1 by default)
-INC_INLET_DAMPING= 0.1
+INC_INLET_DAMPING= 0.01
 %
 % Impose inlet velocity magnitude in the direction of the normal of the inlet face
 INC_INLET_USENORMAL= NO
@@ -287,7 +287,7 @@ INC_INLET_USENORMAL= NO
 INC_OUTLET_TYPE= PRESSURE_OUTLET
 %
 % Damping coefficient for iterative updates at mass flow outlets. (0.1 by default)
-INC_OUTLET_DAMPING= 0.1
+INC_OUTLET_DAMPING= 0.01
 %
 % Bulk Modulus for computing the Mach number
 BULK_MODULUS= 1.42E5