-
Notifications
You must be signed in to change notification settings - Fork 968
Description
What should the default value of RefElemLength in the code below be? The comment says it is 0.1 but the code says 1.0, which are inconsistent.
/*--- Length based parameter for slope limiters uses a default value of
0.1m ---*/
RefElemLength = 1.0;
if (SystemMeasurements == US) RefElemLength /= 0.3048;Line 4696 in 9f082de
| RefElemLength = 1.0; |
For shape optimization, when I set ADJ_SHARP_LIMITER_COEFF to default (3.0) and set SENS_REMOVE_SHARP to YES, surface sensitivity becomes zero within fairly large distance from the sharp edges. The cutoff distance is defined by the product of ADJ_SHARP_LIMITER_COEFF (AdjSharp_LimiterCoeff) and RefElemLength, so RefElemLength seems to be too large.
SU2/SU2_CFD/src/solvers/CDiscAdjSolver.cpp
Line 488 in 9f082de
| const su2double eps = config->GetAdjSharp_LimiterCoeff()*config->GetRefElemLength(); |
Because of this, I'm assuming that the default value of RefElemLength should have been 0.1 instead of 1.0. However, I realized that RefElemLength is used by the slope limiters as well, and I'm currently not sure about the impact of changing this variable.
I am investigating myself but if someone knows the intent of this variable or the implementation of the slope limiters, I would like to ask for your help before I open PR.