Skip to content

Add V&V case, Mach 5 Shock-Wave Boundary-Layer Interaction #1651

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 4 commits into from
May 31, 2022
Merged

Add V&V case, Mach 5 Shock-Wave Boundary-Layer Interaction #1651

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
881729f
add configs
pcarruscag May 29, 2022
bb22c9b
update vandv.py
pcarruscag May 29, 2022
8e11d38
Update TestCases/vandv.py
pcarruscag May 30, 2022
4c7fe68
Apply suggestions from code review
pcarruscag May 31, 2022
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
16 changes: 16 additions & 0 deletions TestCases/vandv.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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 ###
#################
Expand Down
98 changes: 98 additions & 0 deletions TestCases/vandv/rans/swbli/config_sa.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 )
100 changes: 100 additions & 0 deletions TestCases/vandv/rans/swbli/config_sst.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 )