Merge pull request #24 from su2code/feature_nlctrl

Tutorial Updates for v7

Author: talbring

Inc_Inviscid_Hydrofoil/inv_naca0012.cfg

@@ -155,11 +155,10 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 % Convergence criteria (CAUCHY, RESIDUAL)
 %
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_PRESSURE
 %
 % Min value of the residual (log10 of the residual)
-CONV_RESIDUAL_MINVAL= -18
+CONV_RESIDUAL_MINVAL= -15
 %
 % Start Cauchy criteria at iteration number
 CONV_STARTITER= 10
@@ -204,4 +203,4 @@ WRT_SOL_FREQ= 50
 WRT_CON_FREQ= 1
 % 
 % Screen output
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, WALL_TIME, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)

Inc_Laminar_Cavity/lam_buoyancy_cavity.cfg

@@ -15,7 +15,6 @@
 %                               POISSON_EQUATION)
 SOLVER= INC_NAVIER_STOKES
 %
-%
 % Specify turbulent model (NONE, SA, SA_NEG, SST)
 KIND_TURB_MODEL= NONE
 %
@@ -123,11 +122,11 @@ NUM_METHOD_GRAD= GREEN_GAUSS
 CFL_NUMBER= 50.0
 %
 % Adaptive CFL number (NO, YES)
-CFL_ADAPT= NO
+CFL_ADAPT= YES
 %
 % Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
 %                                        CFL max value )
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 15.0, 1e10)
+CFL_ADAPT_PARAM= ( 0.1, 2.0, 50.0, 1e10)
 %
 % Maximum Delta Time in local time stepping simulations
 MAX_DELTA_TIME= 1E6
@@ -153,18 +152,18 @@ LINEAR_SOLVER_ILU_FILL_IN= 0
 LINEAR_SOLVER_ERROR= 1E-15
 %
 % Max number of iterations of the linear solver for the implicit formulation
-LINEAR_SOLVER_ITER= 10
+LINEAR_SOLVER_ITER= 30
 
 % -------------------------- MULTIGRID PARAMETERS -----------------------------%
 %
 % Multi-Grid Levels (0 = no multi-grid)
-MGLEVEL= 2
+MGLEVEL= 0
 %
 % Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE)
 MGCYCLE= V_CYCLE
 %
 % Multi-grid pre-smoothing level
-MG_PRE_SMOOTH= ( 1, 1, 1, 1 )
+MG_PRE_SMOOTH= ( 1, 2, 3, 3 )
 %
 % Multi-grid post-smoothing level
 MG_POST_SMOOTH= ( 0, 0, 0, 0 )
@@ -173,10 +172,10 @@ MG_POST_SMOOTH= ( 0, 0, 0, 0 )
 MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 )
 %
 % Damping factor for the residual restriction
-MG_DAMP_RESTRICTION= 0.8
+MG_DAMP_RESTRICTION= 0.95
 %
 % Damping factor for the correction prolongation
-MG_DAMP_PROLONGATION= 0.8
+MG_DAMP_PROLONGATION= 0.95
 
 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
 %
@@ -199,11 +198,10 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Convergence criteria (CAUCHY, RESIDUAL)
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_TEMPERATURE
 %
 % Min value of the residual (log10 of the residual)
-CONV_RESIDUAL_MINVAL= -12
+CONV_RESIDUAL_MINVAL= -8
 %
 % Start convergence criteria at iteration number
 CONV_STARTITER= 10
@@ -213,7 +211,6 @@ CONV_CAUCHY_ELEMS= 100
 %
 % Epsilon to control the series convergence
 CONV_CAUCHY_EPS= 1E-6
-%
 
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
@@ -266,12 +263,4 @@ WRT_SOL_FREQ= 100
 WRT_CON_FREQ= 1
 %
 % Screen output
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)
-% ----------------------- GEOMETRY EVALUATION PARAMETERS ----------------------%
-%
-% Extract a slice using GEO_BOUNDS (works with 2D cartesian grids only)
-WRT_SLICE= YES
-%
-% Geometrical bounds (Y coordinate) for the wing geometry analysis or
-% fuselage evaluation (X coordinate).
-GEO_BOUNDS= ( 0.499, 0.501)
+SCREEN_OUTPUT= (INNER_ITER, WALL_TIME, RMS_PRESSURE, RMS_TEMPERATURE, LIFT, DRAG)

Inc_Laminar_Flat_Plate/lam_flatplate.cfg

@@ -135,11 +135,11 @@ NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
 CFL_NUMBER= 15.0
 %
 % Adaptive CFL number (NO, YES)
-CFL_ADAPT= NO
+CFL_ADAPT= YES
 %
 % Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
 %                                        CFL max value )
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 1.0, 100.0 )
+CFL_ADAPT_PARAM= ( 0.1, 2.0, 1.0, 1e10 )
 %
 % Number of total iterations
 ITER= 99999
@@ -164,7 +164,7 @@ LINEAR_SOLVER_ITER= 25
 % -------------------------- MULTIGRID PARAMETERS -----------------------------%
 %
 % Multi-Grid Levels (0 = no multi-grid)
-MGLEVEL= 3
+MGLEVEL= 0
 %
 % Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE)
 MGCYCLE= W_CYCLE
@@ -209,12 +209,10 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Convergence criteria (CAUCHY, RESIDUAL)
-%
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_PRESSURE
 %
 % Min value of the residual (log10 of the residual)
-CONV_RESIDUAL_MINVAL= -14
+CONV_RESIDUAL_MINVAL= -16
 %
 % Start convergence criteria at iteration number
 CONV_STARTITER= 10
@@ -224,7 +222,6 @@ CONV_CAUCHY_ELEMS= 100
 %
 % Epsilon to control the series convergence
 CONV_CAUCHY_EPS= 1E-6
-%
 
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
@@ -277,12 +274,4 @@ WRT_SOL_FREQ= 250
 WRT_CON_FREQ= 1
 %
 % Screen output
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)
-
-% ----------------------- GEOMETRY EVALUATION PARAMETERS ----------------------%
-%
-WRT_SLICE= YES
-%
-% Geometrical bounds (Y coordinate) for the wing geometry analysis or
-% fuselage evaluation (X coordinate).
-GEO_BOUNDS= ( 0.304, 0.305)
+SCREEN_OUTPUT= (INNER_ITER, WALL_TIME,  RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)

Inc_Laminar_Step/lam_backwardstep.cfg

@@ -194,9 +194,7 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Convergence criteria (CAUCHY, RESIDUAL)
-%
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_PRESSURE
 %
 % Min value of the residual (log10 of the residual)
 CONV_RESIDUAL_MINVAL= -16
@@ -209,7 +207,6 @@ CONV_CAUCHY_ELEMS= 100
 %
 % Epsilon to control the series convergence
 CONV_CAUCHY_EPS= 1E-6
-%
 
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
@@ -262,14 +259,4 @@ WRT_SOL_FREQ= 250
 WRT_CON_FREQ= 1
 %
 % Screen output
-SCREEN_OUTPUT=(INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)
-
-% ----------------------- GEOMETRY EVALUATION PARAMETERS ----------------------%
-%
-% Extract a slice using GEO_BOUNDS (works with 2D cartesian grids only)
-WRT_SLICE= YES
-%
-% Geometrical bounds (Y coordinate) for the wing geometry analysis or
-% fuselage evaluation (X coordinate).
-%GEO_BOUNDS= ( 6.99, 7.01)
-GEO_BOUNDS= ( 14.99, 15.01)
+SCREEN_OUTPUT=(INNER_ITER, WALL_TIME, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)

Inc_Turbulent_Flat_Plate/turb_flatplate.cfg

@@ -13,7 +13,7 @@
 % Physical governing equations (EULER, NAVIER_STOKES,
 %                               WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
 %                               POISSON_EQUATION)
-SOLVER= INC_NAVIER_STOKES
+SOLVER= INC_RANS
 %
 % If Navier-Stokes, kind of turbulent model (NONE, SA)
 KIND_TURB_MODEL= SA
@@ -100,11 +100,11 @@ NUM_METHOD_GRAD= GREEN_GAUSS
 CFL_NUMBER= 100.0
 %
 % Adaptive CFL number (NO, YES)
-CFL_ADAPT= NO
+CFL_ADAPT= YES
 %
 % Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
 %                                        CFL max value )
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 1.0, 100.0 )
+CFL_ADAPT_PARAM= ( 0.1, 2.0, 100.0, 1e3 )
 %
 % Runge-Kutta alpha coefficients
 RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 )
@@ -201,8 +201,7 @@ TIME_DISCRE_TURB= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Convergence criteria (CAUCHY, RESIDUAL)
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_PRESSURE
 %
 % Min value of the residual (log10 of the residual)
 CONV_RESIDUAL_MINVAL= -14
@@ -268,13 +267,5 @@ WRT_SOL_FREQ= 250
 WRT_CON_FREQ= 1
 % 
 % Screen output
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)
-
-% ----------------------- GEOMETRY EVALUATION PARAMETERS ----------------------%
-%
-WRT_SLICE= YES
-%
-% Geometrical bounds (Y coordinate) for the wing geometry analysis or
-% fuselage evaluation (X coordinate).
-GEO_BOUNDS= ( 0.970, 0.9701)
+SCREEN_OUTPUT= (INNER_ITER, WALL_TIME, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)
 

Inc_Turbulent_NACA0012/turb_naca0012.cfg

@@ -14,7 +14,7 @@
 % Physical governing equations (EULER, NAVIER_STOKES,
 %                               WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
 %                               POISSON_EQUATION)
-SOLVER= INC_NAVIER_STOKES
+SOLVER= INC_RANS
 %
 % Specify turbulent model (NONE, SA, SA_NEG, SST)
 KIND_TURB_MODEL= SA
@@ -101,7 +101,7 @@ CFL_ADAPT= NO
 %
 % Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
 %                                        CFL max value )
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 25.0, 10000.0 )
+CFL_ADAPT_PARAM= ( 0.1, 1.2, 25.0, 1e3 )
 %
 % Maximum Delta Time in local time stepping simulations
 MAX_DELTA_TIME= 1E6
@@ -203,9 +203,7 @@ TIME_DISCRE_TURB= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Convergence criteria (CAUCHY, RESIDUAL)
-%
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_PRESSURE
 %
 % Min value of the residual (log10 of the residual)
 CONV_RESIDUAL_MINVAL= -14
@@ -218,7 +216,6 @@ CONV_CAUCHY_ELEMS= 100
 %
 % Epsilon to control the series convergence
 CONV_CAUCHY_EPS= 1E-6
-%
 
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
@@ -271,5 +268,5 @@ WRT_SOL_FREQ= 250
 WRT_CON_FREQ= 1
 %
 % Screen output
-SCREEN_OUTPUT = (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT = (INNER_ITER, WALL_TIME, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)
 

Inviscid_Bump/inv_channel.cfg

@@ -89,11 +89,11 @@ NUM_METHOD_GRAD= GREEN_GAUSS
 CFL_NUMBER= 50.0
 %
 % Adaptive CFL number (NO, YES)
-CFL_ADAPT= NO
+CFL_ADAPT= YES
 %
 % Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
 %                                        CFL max value )
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 1.0, 100.0 )
+CFL_ADAPT_PARAM= ( 0.1, 2.0, 50.0, 1e10 )
 %
 % Runge-Kutta alpha coefficients
 RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 )
@@ -110,18 +110,18 @@ LINEAR_SOLVER= FGMRES
 LINEAR_SOLVER_PREC= ILU
 %
 % Minimum error of the linear solver for implicit formulations
-LINEAR_SOLVER_ERROR= 1E-4
+LINEAR_SOLVER_ERROR= 1E-10
 %
 % Max number of iterations of the linear solver for the implicit formulation
-LINEAR_SOLVER_ITER= 10
+LINEAR_SOLVER_ITER= 20
 
 % -------------------------- MULTIGRID PARAMETERS -----------------------------%
 %
 % Multi-Grid Levels (0 = no multi-grid)
 MGLEVEL= 3
 %
 % Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE)
-MGCYCLE= V_CYCLE
+MGCYCLE= W_CYCLE
 %
 % Multi-grid pre-smoothing level
 MG_PRE_SMOOTH= ( 1, 2, 3, 3 )
@@ -133,27 +133,17 @@ MG_POST_SMOOTH= ( 0, 0, 0, 0 )
 MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 )
 %
 % Damping factor for the residual restriction
-MG_DAMP_RESTRICTION= 0.9
+MG_DAMP_RESTRICTION= 1.0
 %
 % Damping factor for the correction prolongation
-MG_DAMP_PROLONGATION= 0.9
+MG_DAMP_PROLONGATION= 1.0
 
 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
 %
 % Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
-%           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
-%
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
-% Coefficient for the limiter (smooth regions)
-VENKAT_LIMITER_COEFF= 0.03
-%
 % 2nd and 4th order artificial dissipation coefficients
 JST_SENSOR_COEFF= ( 0.5, 0.02 )
 %
@@ -163,12 +153,10 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Convergence criteria (CAUCHY, RESIDUAL)
-%
-CONV_CRITERIA= RESIDUAL
-%
+CONV_FIELD= RMS_DENSITY
 %
 % Min value of the residual (log10 of the residual)
-CONV_RESIDUAL_MINVAL= -12
+CONV_RESIDUAL_MINVAL= -10
 %
 % Start convergence criteria at iteration number
 CONV_STARTITER= 10
@@ -178,7 +166,6 @@ CONV_CAUCHY_ELEMS= 100
 %
 % Epsilon to control the series convergence
 CONV_CAUCHY_EPS= 1E-10
-%
 
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
@@ -231,4 +218,4 @@ WRT_SOL_FREQ= 250
 WRT_CON_FREQ= 1
 % 
 % Screen output
-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, WALL_TIME, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)