Skip to content

Commit 55b2723

Browse files
pcarruscagpcarruscag
authored
Merge pull request #1651 from su2code/swbli
Add V&V case, Mach 5 Shock-Wave Boundary-Layer Interaction
2 parents ead83fe + 4c7fe68 commit 55b2723

3 files changed

Lines changed: 214 additions & 0 deletions

File tree

‎TestCases/vandv.py‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -64,6 +64,22 @@ def main():
6464
bump_sst1994m.test_vals = [-13.069139, -10.325333, -5.558903, 0.024576, 0.004967]
6565
test_list.append(bump_sst1994m)
6666

67+
# SWBLI SA
68+
swbli_sa = TestCase('swbli_sa')
69+
swbli_sa.cfg_dir = "vandv/rans/swbli"
70+
swbli_sa.cfg_file = "config_sa.cfg"
71+
swbli_sa.test_iter = 20
72+
swbli_sa.test_vals = [-11.029255, -10.511982, -11.400926, -10.128471, -14.536798, 0.002233, -2.608466, 2.786]
73+
test_list.append(swbli_sa)
74+
75+
# SWBLI SST
76+
swbli_sst = TestCase('swbli_sst')
77+
swbli_sst.cfg_dir = "vandv/rans/swbli"
78+
swbli_sst.cfg_file = "config_sst.cfg"
79+
swbli_sst.test_iter = 20
80+
swbli_sst.test_vals = [-11.039738, -10.619752, -11.577767, -10.276345, -11.117432, -3.84103, 0.002315, -2.316445, -5.392453, 2.786]
81+
test_list.append(swbli_sst)
82+
6783
#################
6884
### RUN TESTS ###
6985
#################

‎TestCases/vandv/rans/swbli/config_sa.cfg‎

Lines changed: 98 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,98 @@
1+
% See https://su2code.github.io/vandv/swbli/
2+
%
3+
% ------------------------------- SOLVER -------------------------------- %
4+
%
5+
SOLVER= RANS
6+
KIND_TURB_MODEL= SA
7+
SA_OPTIONS= WITHFT2, NEGATIVE
8+
REF_DIMENSIONALIZATION= FREESTREAM_VEL_EQ_MACH
9+
%
10+
% ----------------------------- FREESTREAM ------------------------------ %
11+
%
12+
MACH_NUMBER= 5
13+
AOA= 0
14+
INIT_OPTION= TD_CONDITIONS
15+
FREESTREAM_OPTION= TEMPERATURE_FS
16+
FREESTREAM_TEMPERATURE= 68.3333
17+
FREESTREAM_PRESSURE= 4006.88
18+
REYNOLDS_LENGTH= 0.523
19+
%
20+
% -------------------------- REFERENCE VALUES --------------------------- %
21+
%
22+
REF_AREA= 0.523
23+
REF_LENGTH= 0.523
24+
REF_ORIGIN_MOMENT_X= 0.0
25+
REF_ORIGIN_MOMENT_Y= 0.0
26+
REF_ORIGIN_MOMENT_Z= 0.0
27+
%
28+
% ------------------------- BOUNDARY CONDITIONS ------------------------- %
29+
%
30+
MARKER_HEATFLUX= ( bottom, 0.0, top, 0.0 )
31+
MARKER_SUPERSONIC_INLET= ( inlet, 68.3333, 4006.88, 828.58, 0, 0 )
32+
MARKER_OUTLET= ( outlet, 4006.88 )
33+
MARKER_SYM= ( sym )
34+
MARKER_PLOTTING= ( bottom )
35+
MARKER_MONITORING= ( bottom )
36+
%
37+
% ---------------------------- FLUID MODELS ----------------------------- %
38+
%
39+
FLUID_MODEL= STANDARD_AIR
40+
GAMMA_VALUE= 1.4
41+
GAS_CONSTANT= 287.058
42+
%
43+
VISCOSITY_MODEL= SUTHERLAND
44+
MU_REF= 1.716E-5
45+
MU_T_REF= 273.15
46+
SUTHERLAND_CONSTANT= 110.4
47+
%
48+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
49+
PRANDTL_LAM= 0.72
50+
PRANDTL_TURB= 0.90
51+
%
52+
% ----------------------- SPATIAL DISCRETIZATION ------------------------ %
53+
%
54+
NUM_METHOD_GRAD= GREEN_GAUSS
55+
CONV_NUM_METHOD_FLOW= ROE
56+
ENTROPY_FIX_COEFF= 1e-5
57+
MUSCL_FLOW= YES
58+
SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
59+
VENKAT_LIMITER_COEFF= 0.1
60+
%
61+
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
62+
MUSCL_TURB= NO
63+
%
64+
% ---------- PSEUDOTIME INTEGRATION / CONVERGENCE ACCELERATION ---------- %
65+
%
66+
USE_VECTORIZATION= YES
67+
TIME_DISCRE_FLOW= EULER_IMPLICIT
68+
TIME_DISCRE_TURB= EULER_IMPLICIT
69+
%
70+
CFL_NUMBER= 1
71+
CFL_REDUCTION_TURB= 1
72+
CFL_ADAPT= YES
73+
CFL_ADAPT_PARAM= ( 1.0, 1.05, 1, 100, 0.8 )
74+
%
75+
LINEAR_SOLVER= FGMRES
76+
LINEAR_SOLVER_PREC= ILU
77+
LINEAR_SOLVER_ERROR= 0.0001
78+
LINEAR_SOLVER_ITER= 3
79+
%
80+
NEWTON_KRYLOV= YES
81+
NEWTON_KRYLOV_IPARAM= ( 0, 3, 1 ) % n0, np, ft
82+
NEWTON_KRYLOV_DPARAM= ( 0.0, 1e-20, -3, 1e-5 ) % r0, tp, rf, e
83+
%
84+
% ------------------------ CONVERGENCE CRITERIA ------------------------- %
85+
%
86+
ITER= 2000
87+
CONV_RESIDUAL_MINVAL= -11.5
88+
%
89+
% --------------------------- INPUT / OUTPUT ---------------------------- %
90+
%
91+
MESH_FILENAME= mesh.su2
92+
MESH_FORMAT= SU2
93+
RESTART_SOL= YES
94+
SOLUTION_FILENAME= solution_sa.dat
95+
OUTPUT_WRT_FREQ= 500
96+
SCREEN_WRT_FREQ_INNER= 1
97+
HISTORY_OUTPUT= ( ITER, RMS_RES, AERO_COEFF )
98+
SCREEN_OUTPUT= ( INNER_ITER, RMS_RES, FORCE_X, LINSOL_RESIDUAL, AVG_CFL )

‎TestCases/vandv/rans/swbli/config_sst.cfg‎

Lines changed: 100 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,100 @@
1+
% See https://su2code.github.io/vandv/swbli/
2+
%
3+
% ------------------------------- SOLVER -------------------------------- %
4+
%
5+
SOLVER= RANS
6+
KIND_TURB_MODEL= SST
7+
SST_OPTIONS= V2003m
8+
REF_DIMENSIONALIZATION= FREESTREAM_VEL_EQ_MACH
9+
%
10+
% ----------------------------- FREESTREAM ------------------------------ %
11+
%
12+
MACH_NUMBER= 5
13+
AOA= 0
14+
INIT_OPTION= TD_CONDITIONS
15+
FREESTREAM_OPTION= TEMPERATURE_FS
16+
FREESTREAM_TEMPERATURE= 68.3333
17+
FREESTREAM_PRESSURE= 4006.88
18+
REYNOLDS_LENGTH= 0.523
19+
FREESTREAM_TURBULENCEINTENSITY= 5e-4
20+
FREESTREAM_TURB2LAMVISCRATIO= 0.01
21+
%
22+
% -------------------------- REFERENCE VALUES --------------------------- %
23+
%
24+
REF_AREA= 0.523
25+
REF_LENGTH= 0.523
26+
REF_ORIGIN_MOMENT_X= 0.0
27+
REF_ORIGIN_MOMENT_Y= 0.0
28+
REF_ORIGIN_MOMENT_Z= 0.0
29+
%
30+
% ------------------------- BOUNDARY CONDITIONS ------------------------- %
31+
%
32+
MARKER_HEATFLUX= ( bottom, 0.0, top, 0.0 )
33+
MARKER_SUPERSONIC_INLET= ( inlet, 68.3333, 4006.88, 828.58, 0, 0 )
34+
MARKER_OUTLET= ( outlet, 4006.88 )
35+
MARKER_SYM= ( sym )
36+
MARKER_PLOTTING= ( bottom )
37+
MARKER_MONITORING= ( bottom )
38+
%
39+
% ---------------------------- FLUID MODELS ----------------------------- %
40+
%
41+
FLUID_MODEL= STANDARD_AIR
42+
GAMMA_VALUE= 1.4
43+
GAS_CONSTANT= 287.058
44+
%
45+
VISCOSITY_MODEL= SUTHERLAND
46+
MU_REF= 1.716E-5
47+
MU_T_REF= 273.15
48+
SUTHERLAND_CONSTANT= 110.4
49+
%
50+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
51+
PRANDTL_LAM= 0.72
52+
PRANDTL_TURB= 0.90
53+
%
54+
% ----------------------- SPATIAL DISCRETIZATION ------------------------ %
55+
%
56+
NUM_METHOD_GRAD= GREEN_GAUSS
57+
CONV_NUM_METHOD_FLOW= ROE
58+
ENTROPY_FIX_COEFF= 1e-5
59+
MUSCL_FLOW= YES
60+
SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
61+
VENKAT_LIMITER_COEFF= 0.1
62+
%
63+
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
64+
MUSCL_TURB= NO
65+
%
66+
% ---------- PSEUDOTIME INTEGRATION / CONVERGENCE ACCELERATION ---------- %
67+
%
68+
USE_VECTORIZATION= YES
69+
TIME_DISCRE_FLOW= EULER_IMPLICIT
70+
TIME_DISCRE_TURB= EULER_IMPLICIT
71+
%
72+
CFL_NUMBER= 1
73+
CFL_REDUCTION_TURB= 1
74+
CFL_ADAPT= YES
75+
CFL_ADAPT_PARAM= ( 1.0, 1.05, 1, 100, 0.8 )
76+
%
77+
LINEAR_SOLVER= FGMRES
78+
LINEAR_SOLVER_PREC= ILU
79+
LINEAR_SOLVER_ERROR= 0.0001
80+
LINEAR_SOLVER_ITER= 3
81+
%
82+
NEWTON_KRYLOV= YES
83+
NEWTON_KRYLOV_IPARAM= ( 0, 3, 1 ) % n0, np, ft
84+
NEWTON_KRYLOV_DPARAM= ( 0.0, 1e-20, -3, 1e-5 ) % r0, tp, rf, e
85+
%
86+
% ------------------------ CONVERGENCE CRITERIA ------------------------- %
87+
%
88+
ITER= 2000
89+
CONV_RESIDUAL_MINVAL= -11.5
90+
%
91+
% --------------------------- INPUT / OUTPUT ---------------------------- %
92+
%
93+
MESH_FILENAME= mesh.su2
94+
MESH_FORMAT= SU2
95+
RESTART_SOL= YES
96+
SOLUTION_FILENAME= solution_sst.dat
97+
OUTPUT_WRT_FREQ= 500
98+
SCREEN_WRT_FREQ_INNER= 1
99+
HISTORY_OUTPUT= ( ITER, RMS_RES, AERO_COEFF )
100+
SCREEN_OUTPUT= ( INNER_ITER, RMS_RES, FORCE_X, LINSOL_RESIDUAL, LINSOL_RESIDUAL_TURB, AVG_CFL )

0 commit comments

Comments
 (0)