@@ -25,7 +25,6 @@ FREESTREAM_TEMPERATURE = 673
2525FREESTREAM_DENSITY = 2.55
2626% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
2727%
28- %INC_DENSITY_MODEL = CONSTANT
2928INC_DENSITY_MODEL = VARIABLE
3029INC_DENSITY_INIT = 2.55
3130%
@@ -41,24 +40,30 @@ INC_NONDIM= DIMENSIONAL
4140%FLUID_MODEL = CONSTANT_DENSITY
4241FLUID_MODEL = INC_IDEAL_GAS
4342%
44- CONDUCTIVITY_MODEL = CONSTANT_CONDUCTIVITY
45- THERMAL_CONDUCTIVITY_CONSTANT = 0.0357
43+ CONDUCTIVITY_MODEL = CONSTANT_PRANDTL
44+ % THERMAL_CONDUCTIVITY_CONSTANT = 0.0357
4645%
4746PRANDTL_LAM = 0.72
48- TURBULENT_CONDUCTIVITY_MODEL = NONE
47+ TURBULENT_CONDUCTIVITY_MODEL = CONSTANT_PRANDTL_TURB
48+ %TURBULENT_CONDUCTIVITY_MODEL = NONE
4949PRANDTL_TURB = 0.90
5050%
51- VISCOSITY_MODEL = SUTHERLAND
51+ %VISCOSITY_MODEL = SUTHERLAND
52+ VISCOSITY_MODEL = POLYNOMIAL_VISCOSITY
53+ % polynomial viscosity for methane-air (phi =0.50)
54+ MU_POLYCOEFFS = (1.097e-5, 3.437e-8, 3.332e-12, 0, 0)
55+
56+
5257MU_CONSTANT = 1.716E-5
5358MU_REF = 1.716e-5
5459MU_T_REF = 273.15
5560SUTHERLAND_CONSTANT = 110.4
5661
57- SPECIFIC_HEAT_CP = 1000.0
62+ SPECIFIC_HEAT_CP = 1350
5863%
5964% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
6065%
61- MARKER_HEATFLUX = ( wall_top, 0.0,wall_side,0.0,wall_pipe,0.0, wall_out , 0.0 )
66+ MARKER_HEATFLUX = ( wall_top, 0.0, wall_side, 0.0, wall_pipe , 0.0 )
6267
6368% note, case is axisymmetric
6469MARKER_SYM = ( symmetry )
@@ -68,43 +73,34 @@ SPECIFIED_INLET_PROFILE= YES
6873INLET_MATCHING_TOLERANCE =1e-4
6974INLET_FILENAME = inlet.dat
7075INC_INLET_TYPE = VELOCITY_INLET
71- INC_INLET_DAMPING = 0.01
76+ INC_INLET_DAMPING = 0.1
7277MARKER_INLET = ( inlet, 673, 40.0, 1.0, 0.0, 0.0)
7378MARKER_INLET_TURBULENT = (inlet, 0.10, 15)
7479MARKER_INLET_SPECIES = (inlet, 0.0)
7580%
7681INC_OUTLET_TYPE = PRESSURE_OUTLET
77- INC_OUTLET_DAMPING = 0.01
82+ INC_OUTLET_DAMPING = 0.1
7883MARKER_OUTLET = ( outlet, 0.0 )
79-
80- % ------------------------ WALL FUNCTION DEFINITION --------------------------%
81- %
82- %MARKER_WALL_FUNCTIONS = ( wall_top, STANDARD_WALL_FUNCTION, wall_side,STANDARD_WALL_FUNCTION )
83- %WALLMODEL_KAPPA = 0.41
84- %WALLMODEL_B = 5.5
85- %WALLMODEL_MINYPLUS = 5.0
86- %WALLMODEL_MAXITER = 200
87- %WALLMODEL_RELFAC = 0.5
88-
8984%
9085% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
9186%
9287NUM_METHOD_GRAD = WEIGHTED_LEAST_SQUARES
9388%
94- CFL_NUMBER = 10.0
89+ CFL_NUMBER = 10.00
9590CFL_REDUCTION_SPECIES = 1.0
9691CFL_REDUCTION_TURB = 1.0
9792CFL_ADAPT = YES
98- CFL_ADAPT_PARAM = ( 0.95, 1.01, 5.0, 250)
93+ CFL_ADAPT_PARAM = ( 0.95, 1.01, 1.0, 250, 1.0e-4, 0)
94+
9995%
10096ITER = 1
10197%
10298% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
10399%
104100LINEAR_SOLVER = FGMRES
105101LINEAR_SOLVER_PREC = ILU
106- LINEAR_SOLVER_ERROR = 1E-3
107- LINEAR_SOLVER_ITER = 5
102+ LINEAR_SOLVER_ERROR = 1E-12
103+ LINEAR_SOLVER_ITER = 10
108104
109105% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
110106%
@@ -116,9 +112,12 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
116112% -------------------- SCALAR TRANSPORT ---------------------------------------%
117113%
118114KIND_SCALAR_MODEL = SPECIES_TRANSPORT
119- DIFFUSIVITY_MODEL = CONSTANT_DIFFUSIVITY
115+ %DIFFUSIVITY_MODEL = CONSTANT_DIFFUSIVITY
116+ DIFFUSIVITY_MODEL = CONSTANT_SCHMIDT
117+ %DIFFUSIVITY_MODEL = CONSTANT_LEWIS
120118SCHMIDT_NUMBER_LAMINAR = 1.0
121- DIFFUSIVITY_CONSTANT = 7.56e-5
119+ CONSTANT_LEWIS_NUMBER = 1.0 1.0
120+ DIFFUSIVITY_CONSTANT = 5.0e-5
122121
123122% according to the paper
124123SCHMIDT_NUMBER_TURBULENT = 0.7
@@ -143,15 +142,15 @@ MUSCL_TURB= NO
143142% --------------------------- CONVERGENCE PARAMETERS --------------------------%
144143%
145144CONV_FIELD = RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TKE, RMS_SPECIES
146- CONV_RESIDUAL_MINVAL = -12
145+ CONV_RESIDUAL_MINVAL = -18
147146CONV_STARTITER = 10
148147%
149148% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
150149%
151150MESH_FILENAME = psi.su2
152151%
153152SCREEN_OUTPUT = INNER_ITER WALL_TIME \
154- RMS_ENTHALPY RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION RMS_SPECIES_0, LINSOL_ITER, LINSOL_RESIDUAL, AVG_CFL
153+ RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION RMS_SPECIES_0
155154
156155SCREEN_WRT_FREQ_INNER = 1
157156%
@@ -160,8 +159,8 @@ CONV_FILENAME= history
160159MARKER_ANALYZE = gas_inlet, air_axial_inlet, outlet
161160MARKER_ANALYZE_AVERAGE = AREA
162161%
163- OUTPUT_FILES = RESTART,RESTART_ASCII, PARAVIEW_MULTIBLOCK
164- VOLUME_OUTPUT = RESIDUAL, SOLUTION, PRIMITIVE, SPECIES_UDS_0, RANK, TIMESTEP
162+ OUTPUT_FILES = RESTART, RESTART_ASCII, PARAVIEW_MULTIBLOCK
163+ VOLUME_OUTPUT = RESIDUAL, PRIMITIVE, SPECIES_UDS_0
165164OUTPUT_WRT_FREQ = 100
166165%
167166READ_BINARY_RESTART = YES
0 commit comments