Skip to content

Test cases for turbulence model variants #1418

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
pcarruscag merged 10 commits into from
Nov 6, 2021
Merged

Test cases for turbulence model variants #1418

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
da60eb0
Test cases for Spalart-Allmaras turbulence model variations and corre…
suargi Nov 2, 2021
f01f78b
Better appearance of python scripts
suargi Nov 2, 2021
1a1661f
Update stored residual values for turbulence modeling regression test…
suargi Nov 2, 2021
84c9584
Update serial_regression.py and turbulence models regression tests co…
suargi Nov 3, 2021
ec5e393
Merge branch 'develop' into feature_turbmod_variants_testcases
suargi Nov 3, 2021
faab657
Merge branch 'feature_turbmod_variants_testcases' into develop
suargi Nov 3, 2021
cf12aeb
Merge branch 'feature_turbmod_variants_testcases' of github.com:su2co…
suargi Nov 3, 2021
0e670da
call SU2_CFD directly
pcarruscag Nov 4, 2021
0ea80fc
Merge branch 'feature_turbmod_variants_testcases' of github.com:su2co…
suargi Nov 5, 2021
6bf7b7d
Remove config files header info
suargi Nov 5, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
76 changes: 76 additions & 0 deletions TestCases/parallel_regression.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -620,6 +620,82 @@ def main():
fem_ns_unsteady_cylinder_ader.tol = 0.00001
test_list.append(fem_ns_unsteady_cylinder_ader)

###########################
### Turbulence modeling ###
###########################

# SA Baseline (Identical to RANS SA RAE2822)
turbmod_sa_bsl_rae2822 = TestCase('turbmod_sa_bsl_rae2822')
turbmod_sa_bsl_rae2822.cfg_dir = "turbulence_models/sa/rae2822"
turbmod_sa_bsl_rae2822.cfg_file = "turb_SA_BSL_RAE2822.cfg"
turbmod_sa_bsl_rae2822.test_iter = 20
turbmod_sa_bsl_rae2822.test_vals = [-2.004689, 0.742306, 0.497308, -5.265793, 0.809463, 0.062016]
turbmod_sa_bsl_rae2822.su2_exec = "mpirun -n 2 SU2_CFD"
turbmod_sa_bsl_rae2822.timeout = 1600
turbmod_sa_bsl_rae2822.new_output = True
turbmod_sa_bsl_rae2822.tol = 0.00001
test_list.append(turbmod_sa_bsl_rae2822)

# SA Negative
turbmod_sa_neg_rae2822 = TestCase('turbmod_sa_neg_rae2822')
turbmod_sa_neg_rae2822.cfg_dir = "turbulence_models/sa/rae2822"
turbmod_sa_neg_rae2822.cfg_file = "turb_SA_NEG_RAE2822.cfg"
turbmod_sa_neg_rae2822.test_iter = 20
turbmod_sa_neg_rae2822.test_vals = [-2.004688, 0.742320, 0.497307, -5.265640, 0.809478, 0.061986]
turbmod_sa_neg_rae2822.su2_exec = "mpirun -n 2 SU2_CFD"
turbmod_sa_neg_rae2822.timeout = 1600
turbmod_sa_neg_rae2822.new_output = True
turbmod_sa_neg_rae2822.tol = 0.00001
test_list.append(turbmod_sa_neg_rae2822)

# SA Compressibility Correction
turbmod_sa_comp_rae2822 = TestCase('turbmod_sa_comp_rae2822')
turbmod_sa_comp_rae2822.cfg_dir = "turbulence_models/sa/rae2822"
turbmod_sa_comp_rae2822.cfg_file = "turb_SA_COMP_RAE2822.cfg"
turbmod_sa_comp_rae2822.test_iter = 20
turbmod_sa_comp_rae2822.test_vals = [-2.004687, 0.742304, 0.497309, -5.266081, 0.809467, 0.062029]
turbmod_sa_comp_rae2822.su2_exec = "mpirun -n 2 SU2_CFD"
turbmod_sa_comp_rae2822.timeout = 1600
turbmod_sa_comp_rae2822.new_output = True
turbmod_sa_comp_rae2822.tol = 0.00001
test_list.append(turbmod_sa_comp_rae2822)

# SA Edwards
turbmod_sa_edw_rae2822 = TestCase('turbmod_sa_edw_rae2822')
turbmod_sa_edw_rae2822.cfg_dir = "turbulence_models/sa/rae2822"
turbmod_sa_edw_rae2822.cfg_file = "turb_SA_EDW_RAE2822.cfg"
turbmod_sa_edw_rae2822.test_iter = 20
turbmod_sa_edw_rae2822.test_vals = [-2.004687, 0.742306, 0.497310, -5.290769, 0.809485, 0.062036]
turbmod_sa_edw_rae2822.su2_exec = "mpirun -n 2 SU2_CFD"
turbmod_sa_edw_rae2822.timeout = 1600
turbmod_sa_edw_rae2822.new_output = True
turbmod_sa_edw_rae2822.tol = 0.00001
test_list.append(turbmod_sa_edw_rae2822)

# SA Compressibility and Edwards
turbmod_sa_comp_edw_rae2822 = TestCase('turbmod_sa_comp_edw_rae2822')
turbmod_sa_comp_edw_rae2822.cfg_dir = "turbulence_models/sa/rae2822"
turbmod_sa_comp_edw_rae2822.cfg_file = "turb_SA_COMP_EDW_RAE2822.cfg"
turbmod_sa_comp_edw_rae2822.test_iter = 20
turbmod_sa_comp_edw_rae2822.test_vals = [-2.004687, 0.742306, 0.497310, -5.290769, 0.809485, 0.062036]
turbmod_sa_comp_edw_rae2822.su2_exec = "mpirun -n 2 SU2_CFD"
turbmod_sa_comp_edw_rae2822.timeout = 1600
turbmod_sa_comp_edw_rae2822.new_output = True
turbmod_sa_comp_edw_rae2822.tol = 0.00001
test_list.append(turbmod_sa_comp_edw_rae2822)

# SA QCR
turbmod_sa_qcr_rae2822 = TestCase('turbmod_sa_qcr_rae2822')
turbmod_sa_qcr_rae2822.cfg_dir = "turbulence_models/sa/rae2822"
turbmod_sa_qcr_rae2822.cfg_file = "turb_SA_QCR_RAE2822.cfg"
turbmod_sa_qcr_rae2822.test_iter = 20
turbmod_sa_qcr_rae2822.test_vals = [-2.004793, 0.742353, 0.497315, -5.265974, 0.807841, 0.062027]
turbmod_sa_qcr_rae2822.su2_exec = "mpirun -n 2 SU2_CFD"
turbmod_sa_qcr_rae2822.timeout = 1600
turbmod_sa_qcr_rae2822.new_output = True
turbmod_sa_qcr_rae2822.tol = 0.00001
test_list.append(turbmod_sa_qcr_rae2822)

############################
### Transition ###
############################
Expand Down
168 changes: 168 additions & 0 deletions TestCases/turbulence_models/sa/rae2822/turb_SA_BSL_RAE2822.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,168 @@
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
%
% Physical governing equations
SOLVER= RANS
%
% Specify turbulent model (NONE, SA, SA_NEG, SST)
KIND_TURB_MODEL= SA
%
% Restart solution (NO, YES)
RESTART_SOL= NO
% -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------%
%
% Mach number (non-dimensional, based on the free-stream values)
MACH_NUMBER= 0.729
%
% Angle of attack (degrees, only for compressible flows)
AOA= 2.31
%
% Free-stream temperature (288.15 K by default)
FREESTREAM_TEMPERATURE= 288.15
%
% Reynolds number (non-dimensional, based on the free-stream values)
REYNOLDS_NUMBER= 6.5E6
%
% Reynolds length (1 m by default)
REYNOLDS_LENGTH= 1.0

% ---------------------- REFERENCE VALUE DEFINITION ---------------------------%
%
% Reference origin for moment computation
REF_ORIGIN_MOMENT_X = 0.25
REF_ORIGIN_MOMENT_Y = 0.00
REF_ORIGIN_MOMENT_Z = 0.00
%
% Reference length for pitching, rolling, and yawing non-dimensional moment
REF_LENGTH= 1.0
%
% Reference area for force coefficients (0 implies automatic calculation)
REF_AREA= 1.0

% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
%
% Navier-Stokes wall boundary marker(s) (NONE = no marker)
MARKER_HEATFLUX= ( AIRFOIL, 0.0 )
%
% Farfield boundary marker(s) (NONE = no marker)
MARKER_FAR= ( FARFIELD )
%
% Marker(s) of the surface to be plotted or designed
MARKER_PLOTTING= ( AIRFOIL )
%
% Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated
MARKER_MONITORING= ( AIRFOIL )

% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
%
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
%
% Courant-Friedrichs-Lewy condition of the finest grid
CFL_NUMBER= 2.5
%
% Adaptive CFL number (NO, YES)
CFL_ADAPT= NO
%
% Number of total iterations
ITER= 99999
%
% Linear solver for the implicit formulation (BCGSTAB, FGMRES)
LINEAR_SOLVER= BCGSTAB
%
% Min error of the linear solver for the implicit formulation
LINEAR_SOLVER_ERROR= 1E-6
%
% Max number of iterations of the linear solver for the implicit formulation
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= W_CYCLE
%
% Multi-grid pre-smoothing level
MG_PRE_SMOOTH= ( 1, 2, 3, 3 )
%
% Multi-grid post-smoothing level
MG_POST_SMOOTH= ( 0, 0, 0, 0 )
%
% Jacobi implicit smoothing of the correction
MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 )
%
% Damping factor for the residual restriction
MG_DAMP_RESTRICTION= 0.95
%
% Damping factor for the correction prolongation
MG_DAMP_PROLONGATION= 0.95

% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
%
% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
% TURKEL_PREC, MSW)
CONV_NUM_METHOD_FLOW= JST
%
% 2nd and 4th order artificial dissipation coefficients
JST_SENSOR_COEFF= ( 0.5, 0.02 )
%
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
TIME_DISCRE_FLOW= EULER_IMPLICIT

% -------------------- TURBULENT NUMERICAL METHOD DEFINITION ------------------%
%
% Convective numerical method (SCALAR_UPWIND)
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
%
% Time discretization (EULER_IMPLICIT)
TIME_DISCRE_TURB= EULER_IMPLICIT

% --------------------------- CONVERGENCE PARAMETERS --------------------------%
%
% Min value of the residual (log10 of the residual)
CONV_RESIDUAL_MINVAL= -8
%
% Start convergence criteria at iteration number
CONV_STARTITER= 10
%
% Number of elements to apply the criteria
CONV_CAUCHY_ELEMS= 100
%
% Epsilon to control the series convergence
CONV_CAUCHY_EPS= 1E-6

% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
%
% Mesh input file
MESH_FILENAME= mesh_RAE2822_turb.su2
%
% Mesh input file format (SU2, CGNS, NETCDF_ASCII)
MESH_FORMAT= SU2
%
% Mesh output file
MESH_OUT_FILENAME= mesh_out.su2
%
% Restart flow input file
SOLUTION_FILENAME= solution_flow.dat
%
% Output file format (PARAVIEW, TECPLOT, STL)
TABULAR_FORMAT= CSV
%
% Output file convergence history (w/o extension)
CONV_FILENAME= history
%
% Output file restart flow
RESTART_FILENAME= restart_flow.dat
%
% Output file flow (w/o extension) variables
VOLUME_FILENAME= flow
%
% Output file surface flow coefficient (w/o extension)
SURFACE_FILENAME= surface_flow
%
% Writing solution file frequency
OUTPUT_WRT_FREQ= 250
%
% Screen output fields
SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_NU_TILDE, LIFT, DRAG)
Loading