%------------------------------------------------------------------------------%
% NACA 0021 at 60 deg, Re = 2.7e5: hybrid RANS/LES (V&V page NACA0021_DDES).     %
% Settings of the SST-EDDES run. For the other runs change only:                %
%   SA-DDES  : KIND_TURB_MODEL= SA,  HYBRID_RANSLES= SA_DDES  (remove SST_OPTIONS) %
%   SA-EDDES : KIND_TURB_MODEL= SA,  HYBRID_RANSLES= SA_EDDES (remove SST_OPTIONS) %
%   SST-DDES : KIND_TURB_MODEL= SST, HYBRID_RANSLES= SST_DDES                    %
% The SST variants and the current SA_EDDES length scale need SU2 PR #2150.     %
%------------------------------------------------------------------------------%

% ------------- Problem --------------------------------------------------------%
SOLVER= INC_RANS
KIND_TURB_MODEL= SST
SST_OPTIONS= ( V2003m )
HYBRID_RANSLES= SST_EDDES
MATH_PROBLEM= DIRECT

% Start from a converged steady RANS solution of the same model
% (restart_flow.dat). RESTART_ITER= 2 for the first unsteady run.
RESTART_SOL= YES
RESTART_ITER= 2

% ------------- Free stream: U = 4.1625 m/s at 60 deg in the x-z plane ---------%
% Chord c = 1 m, Re = rho U c / mu = 2.7e5. Span along y.
INC_DENSITY_INIT= 1.2
INC_VELOCITY_INIT= ( 2.08125, 0.0, 3.6048 )
VISCOSITY_MODEL= CONSTANT_VISCOSITY
MU_CONSTANT= 1.85e-5
% Inflow eddy viscosity equal to the laminar viscosity, as in Guseva et al. (2017)
FREESTREAM_TURBULENCEINTENSITY= 0.018
FREESTREAM_TURB2LAMVISCRATIO= 1.0

REF_ORIGIN_MOMENT_X= 0.25
REF_ORIGIN_MOMENT_Y= 0.00
REF_ORIGIN_MOMENT_Z= 0.00
REF_LENGTH= 1.0
REF_AREA= 1.0

% ------------- Time integration: dt = 0.025 c/U --------------------------------%
TIME_DOMAIN= YES
TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
TIME_STEP= 0.006
TIME_ITER= 80000
INNER_ITER= 50

% ------------- Boundaries -------------------------------------------------------%
MARKER_HEATFLUX= ( Wall, 0.0 )
MARKER_PERIODIC= ( Periodic0Y, PeriodicPlusY, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0 )
MARKER_FAR= ( Farfield )
MARKER_PLOTTING= ( Wall )
MARKER_MONITORING= ( Wall )

% ------------- Numerics ---------------------------------------------------------%
NUM_METHOD_GRAD= GREEN_GAUSS
CFL_NUMBER= 10.0
CFL_ADAPT= NO
LINEAR_SOLVER= RESTARTED_FGMRES
LINEAR_SOLVER_ERROR= 1E-6
LINEAR_SOLVER_ITER= 60

CONV_NUM_METHOD_FLOW= JST
MUSCL_FLOW= NO
TIME_DISCRE_FLOW= EULER_IMPLICIT
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
MUSCL_TURB= NO
TIME_DISCRE_TURB= EULER_IMPLICIT

CONV_RESIDUAL_MINVAL= -8
CONV_FIELD= DENSITY
CONV_STARTITER= 50
CONV_CAUCHY_ELEMS= 10
CONV_CAUCHY_EPS= 1E-5

% ------------- Input / output ---------------------------------------------------%
MESH_FILENAME= Mesh.cgns
MESH_FORMAT= CGNS
SOLUTION_FILENAME= restart_flow.dat
RESTART_FILENAME= restart_flow.dat
TABULAR_FORMAT= CSV
CONV_FILENAME= history
VOLUME_FILENAME= flow
SURFACE_FILENAME= surface_flow
% Surface files every 5 steps: they give the sectional forces and the mean Cp
OUTPUT_WRT_FREQ= ( 5, 500, 500 )
OUTPUT_FILES= ( SURFACE_CGNS, RESTART, CGNS )
SCREEN_OUTPUT= ( TIME_ITER, INNER_ITER, WALL_TIME, RMS_RES, LIFT, DRAG, LINSOL_ITER )
HISTORY_OUTPUT= ( ITER, LINSOL, RMS_RES, AERO_COEFF, AERO_COEFF_SURF )
VOLUME_OUTPUT= ( COORDINATES, SOLUTION, PRIMITIVE, VORTEX_IDENTIFICATION, TIME_AVERAGE )