diff --git a/TestCases/vandv.py b/TestCases/vandv.py index 7b5082f248ad..7eb0b484e06c 100644 --- a/TestCases/vandv.py +++ b/TestCases/vandv.py @@ -64,6 +64,22 @@ def main(): bump_sst1994m.test_vals = [-13.069139, -10.325333, -5.558903, 0.024576, 0.004967] test_list.append(bump_sst1994m) + # SWBLI SA + swbli_sa = TestCase('swbli_sa') + swbli_sa.cfg_dir = "vandv/rans/swbli" + swbli_sa.cfg_file = "config_sa.cfg" + swbli_sa.test_iter = 20 + swbli_sa.test_vals = [-11.029255, -10.511982, -11.400926, -10.128471, -14.536798, 0.002233, -2.608466, 2.786] + test_list.append(swbli_sa) + + # SWBLI SST + swbli_sst = TestCase('swbli_sst') + swbli_sst.cfg_dir = "vandv/rans/swbli" + swbli_sst.cfg_file = "config_sst.cfg" + swbli_sst.test_iter = 20 + swbli_sst.test_vals = [-11.039738, -10.619752, -11.577767, -10.276345, -11.117432, -3.84103, 0.002315, -2.316445, -5.392453, 2.786] + test_list.append(swbli_sst) + ################# ### RUN TESTS ### ################# diff --git a/TestCases/vandv/rans/swbli/config_sa.cfg b/TestCases/vandv/rans/swbli/config_sa.cfg new file mode 100644 index 000000000000..6d4638468dc6 --- /dev/null +++ b/TestCases/vandv/rans/swbli/config_sa.cfg @@ -0,0 +1,98 @@ +% See https://su2code.github.io/vandv/swbli/ +% +% ------------------------------- SOLVER -------------------------------- % +% +SOLVER= RANS +KIND_TURB_MODEL= SA +SA_OPTIONS= WITHFT2, NEGATIVE +REF_DIMENSIONALIZATION= FREESTREAM_VEL_EQ_MACH +% +% ----------------------------- FREESTREAM ------------------------------ % +% +MACH_NUMBER= 5 +AOA= 0 +INIT_OPTION= TD_CONDITIONS +FREESTREAM_OPTION= TEMPERATURE_FS +FREESTREAM_TEMPERATURE= 68.3333 +FREESTREAM_PRESSURE= 4006.88 +REYNOLDS_LENGTH= 0.523 +% +% -------------------------- REFERENCE VALUES --------------------------- % +% +REF_AREA= 0.523 +REF_LENGTH= 0.523 +REF_ORIGIN_MOMENT_X= 0.0 +REF_ORIGIN_MOMENT_Y= 0.0 +REF_ORIGIN_MOMENT_Z= 0.0 +% +% ------------------------- BOUNDARY CONDITIONS ------------------------- % +% +MARKER_HEATFLUX= ( bottom, 0.0, top, 0.0 ) +MARKER_SUPERSONIC_INLET= ( inlet, 68.3333, 4006.88, 828.58, 0, 0 ) +MARKER_OUTLET= ( outlet, 4006.88 ) +MARKER_SYM= ( sym ) +MARKER_PLOTTING= ( bottom ) +MARKER_MONITORING= ( bottom ) +% +% ---------------------------- FLUID MODELS ----------------------------- % +% +FLUID_MODEL= STANDARD_AIR +GAMMA_VALUE= 1.4 +GAS_CONSTANT= 287.058 +% +VISCOSITY_MODEL= SUTHERLAND +MU_REF= 1.716E-5 +MU_T_REF= 273.15 +SUTHERLAND_CONSTANT= 110.4 +% +CONDUCTIVITY_MODEL= CONSTANT_PRANDTL +PRANDTL_LAM= 0.72 +PRANDTL_TURB= 0.90 +% +% ----------------------- SPATIAL DISCRETIZATION ------------------------ % +% +NUM_METHOD_GRAD= GREEN_GAUSS +CONV_NUM_METHOD_FLOW= ROE +ENTROPY_FIX_COEFF= 1e-5 +MUSCL_FLOW= YES +SLOPE_LIMITER_FLOW= VENKATAKRISHNAN +VENKAT_LIMITER_COEFF= 0.1 +% +CONV_NUM_METHOD_TURB= SCALAR_UPWIND +MUSCL_TURB= NO +% +% ---------- PSEUDOTIME INTEGRATION / CONVERGENCE ACCELERATION ---------- % +% +USE_VECTORIZATION= YES +TIME_DISCRE_FLOW= EULER_IMPLICIT +TIME_DISCRE_TURB= EULER_IMPLICIT +% +CFL_NUMBER= 1 +CFL_REDUCTION_TURB= 1 +CFL_ADAPT= YES +CFL_ADAPT_PARAM= ( 1.0, 1.05, 1, 100, 0.8 ) +% +LINEAR_SOLVER= FGMRES +LINEAR_SOLVER_PREC= ILU +LINEAR_SOLVER_ERROR= 0.0001 +LINEAR_SOLVER_ITER= 3 +% +NEWTON_KRYLOV= YES +NEWTON_KRYLOV_IPARAM= ( 0, 3, 1 ) % n0, np, ft +NEWTON_KRYLOV_DPARAM= ( 0.0, 1e-20, -3, 1e-5 ) % r0, tp, rf, e +% +% ------------------------ CONVERGENCE CRITERIA ------------------------- % +% +ITER= 2000 +CONV_RESIDUAL_MINVAL= -11.5 +% +% --------------------------- INPUT / OUTPUT ---------------------------- % +% +MESH_FILENAME= mesh.su2 +MESH_FORMAT= SU2 +RESTART_SOL= YES +SOLUTION_FILENAME= solution_sa.dat +OUTPUT_WRT_FREQ= 500 +SCREEN_WRT_FREQ_INNER= 1 +HISTORY_OUTPUT= ( ITER, RMS_RES, AERO_COEFF ) +SCREEN_OUTPUT= ( INNER_ITER, RMS_RES, FORCE_X, LINSOL_RESIDUAL, AVG_CFL ) diff --git a/TestCases/vandv/rans/swbli/config_sst.cfg b/TestCases/vandv/rans/swbli/config_sst.cfg new file mode 100644 index 000000000000..5ed0a98448ea --- /dev/null +++ b/TestCases/vandv/rans/swbli/config_sst.cfg @@ -0,0 +1,100 @@ +% See https://su2code.github.io/vandv/swbli/ +% +% ------------------------------- SOLVER -------------------------------- % +% +SOLVER= RANS +KIND_TURB_MODEL= SST +SST_OPTIONS= V2003m +REF_DIMENSIONALIZATION= FREESTREAM_VEL_EQ_MACH +% +% ----------------------------- FREESTREAM ------------------------------ % +% +MACH_NUMBER= 5 +AOA= 0 +INIT_OPTION= TD_CONDITIONS +FREESTREAM_OPTION= TEMPERATURE_FS +FREESTREAM_TEMPERATURE= 68.3333 +FREESTREAM_PRESSURE= 4006.88 +REYNOLDS_LENGTH= 0.523 +FREESTREAM_TURBULENCEINTENSITY= 5e-4 +FREESTREAM_TURB2LAMVISCRATIO= 0.01 +% +% -------------------------- REFERENCE VALUES --------------------------- % +% +REF_AREA= 0.523 +REF_LENGTH= 0.523 +REF_ORIGIN_MOMENT_X= 0.0 +REF_ORIGIN_MOMENT_Y= 0.0 +REF_ORIGIN_MOMENT_Z= 0.0 +% +% ------------------------- BOUNDARY CONDITIONS ------------------------- % +% +MARKER_HEATFLUX= ( bottom, 0.0, top, 0.0 ) +MARKER_SUPERSONIC_INLET= ( inlet, 68.3333, 4006.88, 828.58, 0, 0 ) +MARKER_OUTLET= ( outlet, 4006.88 ) +MARKER_SYM= ( sym ) +MARKER_PLOTTING= ( bottom ) +MARKER_MONITORING= ( bottom ) +% +% ---------------------------- FLUID MODELS ----------------------------- % +% +FLUID_MODEL= STANDARD_AIR +GAMMA_VALUE= 1.4 +GAS_CONSTANT= 287.058 +% +VISCOSITY_MODEL= SUTHERLAND +MU_REF= 1.716E-5 +MU_T_REF= 273.15 +SUTHERLAND_CONSTANT= 110.4 +% +CONDUCTIVITY_MODEL= CONSTANT_PRANDTL +PRANDTL_LAM= 0.72 +PRANDTL_TURB= 0.90 +% +% ----------------------- SPATIAL DISCRETIZATION ------------------------ % +% +NUM_METHOD_GRAD= GREEN_GAUSS +CONV_NUM_METHOD_FLOW= ROE +ENTROPY_FIX_COEFF= 1e-5 +MUSCL_FLOW= YES +SLOPE_LIMITER_FLOW= VENKATAKRISHNAN +VENKAT_LIMITER_COEFF= 0.1 +% +CONV_NUM_METHOD_TURB= SCALAR_UPWIND +MUSCL_TURB= NO +% +% ---------- PSEUDOTIME INTEGRATION / CONVERGENCE ACCELERATION ---------- % +% +USE_VECTORIZATION= YES +TIME_DISCRE_FLOW= EULER_IMPLICIT +TIME_DISCRE_TURB= EULER_IMPLICIT +% +CFL_NUMBER= 1 +CFL_REDUCTION_TURB= 1 +CFL_ADAPT= YES +CFL_ADAPT_PARAM= ( 1.0, 1.05, 1, 100, 0.8 ) +% +LINEAR_SOLVER= FGMRES +LINEAR_SOLVER_PREC= ILU +LINEAR_SOLVER_ERROR= 0.0001 +LINEAR_SOLVER_ITER= 3 +% +NEWTON_KRYLOV= YES +NEWTON_KRYLOV_IPARAM= ( 0, 3, 1 ) % n0, np, ft +NEWTON_KRYLOV_DPARAM= ( 0.0, 1e-20, -3, 1e-5 ) % r0, tp, rf, e +% +% ------------------------ CONVERGENCE CRITERIA ------------------------- % +% +ITER= 2000 +CONV_RESIDUAL_MINVAL= -11.5 +% +% --------------------------- INPUT / OUTPUT ---------------------------- % +% +MESH_FILENAME= mesh.su2 +MESH_FORMAT= SU2 +RESTART_SOL= YES +SOLUTION_FILENAME= solution_sst.dat +OUTPUT_WRT_FREQ= 500 +SCREEN_WRT_FREQ_INNER= 1 +HISTORY_OUTPUT= ( ITER, RMS_RES, AERO_COEFF ) +SCREEN_OUTPUT= ( INNER_ITER, RMS_RES, FORCE_X, LINSOL_RESIDUAL, LINSOL_RESIDUAL_TURB, AVG_CFL )