Fix biparabolic regression

TestCases/euler/biparabolic/BIPARABOLIC.cfg

@@ -30,47 +30,24 @@ REF_ORIGIN_MOMENT_Z = 0.00
 REF_LENGTH= 1.0
 REF_AREA= 1.0
 
-% ------------------------ SUPERSONIC SIMULATION ------------------------------%
-%
-EQUIV_AREA= YES
-EA_INT_LIMIT= ( -100.0, 100.0, 2.0 )
-
 % ----------------------- BOUNDARY CONDITION DEFINITION -----------------------%
 %
 MARKER_EULER= ( airfoil )
-MARKER_FAR= ( farfield )
-MARKER_NEARFIELD= ( upper_nearfield, lower_nearfield )
+MARKER_FAR= ( farfield, upper_nearfield, lower_nearfield )
 MARKER_PLOTTING= ( airfoil )
-MARKER_MONITORING= ( airfoil, upper_nearfield )
+MARKER_MONITORING= ( airfoil )
 
 % ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
 %
 NUM_METHOD_GRAD= GREEN_GAUSS
-CFL_NUMBER= 25.0
-CFL_ADAPT= NO
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 1.0, 100.0 )
-RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 )
-EXT_ITER= 99999999
-
-% -------------------------- MULTIGRID PARAMETERS -----------------------------%
-%
-MGLEVEL= 3
-MGCYCLE= V_CYCLE
-MG_PRE_SMOOTH= ( 4, 4, 4, 4 )
-MG_POST_SMOOTH= ( 4, 4, 4, 4 )
-MG_CORRECTION_SMOOTH= ( 1, 1, 1, 1 )
-MG_DAMP_RESTRICTION= 0.5
-MG_DAMP_PROLONGATION= 0.5
-
-% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
-%
-CONV_NUM_METHOD_FLOW= JST
-VENKAT_LIMITER_COEFF= 0.03
-JST_SENSOR_COEFF= ( 0.5, 0.02 )
-TIME_DISCRE_FLOW= EULER_IMPLICIT
+CFL_NUMBER= 2.0
+CFL_ADAPT= YES
+CFL_ADAPT_PARAM= ( 0.5, 1.5, 1.0, 50.0 )
+ITER= 25
 
 % ---------------- ADJOINT-FLOW NUMERICAL METHOD DEFINITION -------------------%
 %
+
 OBJECTIVE_FUNCTION= EQUIVALENT_AREA
 CONV_NUM_METHOD_ADJFLOW= JST
 ADJ_SHARP_LIMITER_COEFF= 3.0
@@ -80,19 +57,44 @@ TIME_DISCRE_ADJFLOW= EULER_IMPLICIT
 
 % ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
 %
+
 DV_KIND= PARABOLIC
 DV_MARKER= ( airfoil )
 DV_PARAM= ( 0.5, 5 )
 DV_VALUE= 0.01
 HOLD_GRID_FIXED= YES
 HOLD_GRID_FIXED_COORD= ( -0.2, -0.4, 0.0, 2.2, 0.4, 0.0 )
 
+% -------------------------- MULTIGRID PARAMETERS -----------------------------%
+%
+
+MGLEVEL= 0
+
+% ------------------------ SUPERSONIC SIMULATION ------------------------------%
+%
+
+% Equivalent-area / nearfield setup for sonic-boom studies. Disabled for the
+% direct solve: enabling these produces 760 non-physical points and rms[Rho]
+% stalls at +3.4 instead of -3.9. Re-enable together with MARKER_NEARFIELD
+% below (and remove those two markers from MARKER_FAR) for adjoint EA work.
+%
+% EQUIV_AREA= YES
+% EA_INT_LIMIT= ( -100.0, 100.0, 2.0 )
+% MARKER_NEARFIELD= ( upper_nearfield, lower_nearfield )
+
+
+% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
+%
+CONV_NUM_METHOD_FLOW= ROE
+MUSCL_FLOW= YES
+SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
+VENKAT_LIMITER_COEFF= 0.05
+TIME_DISCRE_FLOW= EULER_IMPLICIT
+
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 CONV_RESIDUAL_MINVAL= -6
 CONV_STARTITER= 10
-CONV_CAUCHY_ELEMS= 100
-CONV_CAUCHY_EPS= 1E-10
 
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
@@ -110,4 +112,4 @@ VOLUME_ADJ_FILENAME= volume_adjoint
 GRAD_OBJFUNC_FILENAME= of_grad
 SURFACE_FILENAME= surface_flow
 SURFACE_ADJ_FILENAME= surface_adjoint
-OUTPUT_WRT_FREQ= 100
+OUTPUT_WRT_FREQ= 500

TestCases/serial_regression.py

@@ -112,6 +112,14 @@ def main():
     wedge.test_vals = [-3.701409, 2.023011, -0.249531, 0.043953]
     test_list.append(wedge)
 
+    # Supersonic biparabolic airfoil
+    biparabolic           = TestCase('biparabolic')
+    biparabolic.cfg_dir   = "euler/biparabolic"
+    biparabolic.cfg_file  = "BIPARABOLIC.cfg"
+    biparabolic.test_iter = 25
+    biparabolic.test_vals = [-2.348749, 0.298540, 0.112150, 3.213351]
+    test_list.append(biparabolic)
+
     # ONERA M6 Wing
     oneram6           = TestCase('oneram6')
     oneram6.cfg_dir   = "euler/oneram6"