Skip to content

Commit 1c9a1fb

Browse files
bigfootedbigfooted
authored
Merge pull request #72 from su2code/TFC_python
TFC psi case
2 parents ab97d67 + 08f2255 commit 1c9a1fb

21 files changed

Lines changed: 131594 additions & 0 deletions

File tree

‎multiphysics/TFC_python/adiabatic/inlet.dat‎

Lines changed: 35 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,35 @@
1+
NMARK= 1
2+
MARKER_TAG= inlet
3+
NROW=30
4+
NCOL=9
5+
# COORD-X COORD-Y TEMPERATURE VELOCITY NORMAL-X NORMAL-Y TKE DISSIPATION SPECIES_0
6+
-0.1 0.0 673 52.113575 1 0 0.022036247 2249.2866 0
7+
-0.1 0.00072056672 673 52.113575 1 0 0.022036247 2249.2866 0
8+
-0.1 0.0014123108 673 52.130169 1 0 0.021853773 2227.1838 0
9+
-0.1 0.0020763851 673 52.155247 1 0 0.021966202 2223.8691 0
10+
-0.1 0.0027138963 673 52.185356 1 0 0.025456807 2238.5793 0
11+
-0.1 0.003325907 673 52.181305 1 0 0.11307145 2565.6211 0
12+
-0.1 0.0039134375 673 51.904926 1 0 0.91227305 3924.4668 0
13+
-0.1 0.0044774669 673 51.180588 1 0 3.0400403 5234.1846 0
14+
-0.1 0.0050189346 673 50.250782 1 0 5.6927605 5843.0571 0
15+
-0.1 0.0055387439 673 49.246353 1 0 8.354352 6278.2769 0
16+
-0.1 0.006037761 673 48.157909 1 0 11.022179 6693.3311 0
17+
-0.1 0.0065168175 673 46.984482 1 0 13.741647 7145.3198 0
18+
-0.1 0.0069767111 673 45.774849 1 0 16.35457 7609.7793 0
19+
-0.1 0.0074182097 673 44.565525 1 0 18.708464 8072.2192 0
20+
-0.1 0.0078420481 673 43.370502 1 0 20.779472 8543.3164 0
21+
-0.1 0.0082489327 673 42.1898 1 0 22.594801 9048.4883 0
22+
-0.1 0.0086395424 673 41.016678 1 0 24.199833 9620.7598 0
23+
-0.1 0.0090145273 673 39.840576 1 0 25.640738 10297.504 0
24+
-0.1 0.0093745124 673 38.646294 1 0 26.959324 11122.591 0
25+
-0.1 0.0097200992 673 37.411186 1 0 28.190733 12151.757 0
26+
-0.1 0.010051861 673 36.11153 1 0 29.362619 13455.807 0
27+
-0.1 0.010370354 673 34.725197 1 0 30.4884 15132.082 0
28+
-0.1 0.010676106 673 33.222515 1 0 31.5849 17335.744 0
29+
-0.1 0.010969629 673 31.561712 1 0 32.676102 20327.223 0
30+
-0.1 0.01125141 673 29.680511 1 0 33.802517 24576.211 0
31+
-0.1 0.011521921 673 27.479441 1 0 35.046871 31017.684 0
32+
-0.1 0.011781611 673 24.78553 1 0 36.619263 41792.531 0
33+
-0.1 0.012030913 673 21.265036 1 0 39.229969 62961.266 0
34+
-0.1 0.012270243 673 15.721132 1 0 44.032578 119438.84 0
35+
-0.1 0.0125 673 0 1 0 1e-10 211634.48 0

‎multiphysics/TFC_python/adiabatic/psi.cfg‎

Lines changed: 166 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,166 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Turbulent premixed high pressure combustion chamber. %
5+
% Author: N. Beishuizen %
6+
% Institution: Bosch Thermotechniek B.V. %
7+
% Date: 2025/01/01 %
8+
% File Version 8.0 "Harrier" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
12+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
13+
%
14+
SOLVER= INC_RANS
15+
KIND_TURB_MODEL= SST
16+
SST_OPTIONS= V1994m
17+
18+
RESTART_SOL= YES
19+
20+
% temperature and density
21+
FREESTREAM_TEMPERATURE = 673
22+
FREESTREAM_DENSITY = 2.55
23+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
24+
%
25+
INC_DENSITY_MODEL= VARIABLE
26+
INC_DENSITY_INIT= 2.55
27+
%
28+
INC_VELOCITY_INIT= (40.00, 0.0, 0.0 )
29+
%
30+
INC_ENERGY_EQUATION= NO
31+
INC_TEMPERATURE_INIT= 673.0
32+
%
33+
INC_NONDIM= DIMENSIONAL
34+
%
35+
% -------------------- FLUID PROPERTIES ------------------------------------- %
36+
%
37+
FLUID_MODEL= INC_IDEAL_GAS
38+
%
39+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
40+
%THERMAL_CONDUCTIVITY_CONSTANT= 0.0357
41+
%
42+
PRANDTL_LAM= 0.72
43+
TURBULENT_CONDUCTIVITY_MODEL= CONSTANT_PRANDTL_TURB
44+
%TURBULENT_CONDUCTIVITY_MODEL= NONE
45+
PRANDTL_TURB= 0.90
46+
%
47+
%VISCOSITY_MODEL= SUTHERLAND
48+
VISCOSITY_MODEL= POLYNOMIAL_VISCOSITY
49+
% polynomial viscosity for methane-air (phi=0.50)
50+
MU_POLYCOEFFS= (1.097e-5, 3.437e-8, 3.332e-12, 0, 0)
51+
52+
MU_CONSTANT= 1.716E-5
53+
MU_REF = 1.716e-5
54+
MU_T_REF= 273.15
55+
SUTHERLAND_CONSTANT = 110.4
56+
57+
SPECIFIC_HEAT_CP = 1350
58+
%
59+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
60+
%
61+
MARKER_HEATFLUX= ( wall_top, 0.0, wall_side, 0.0, wall_pipe, 0.0, wall_out,0.0, wall_side_out,0.0 )
62+
63+
% note, case is axisymmetric
64+
MARKER_SYM= ( symmetry )
65+
AXISYMMETRIC= YES
66+
%
67+
SPECIFIED_INLET_PROFILE= YES
68+
INLET_MATCHING_TOLERANCE=1e-4
69+
INLET_FILENAME= inlet.dat
70+
INC_INLET_TYPE= VELOCITY_INLET
71+
INC_INLET_DAMPING= 0.1
72+
MARKER_INLET= ( inlet, 673, 40.0, 1.0, 0.0, 0.0)
73+
MARKER_INLET_TURBULENT = (inlet, 0.10, 15)
74+
MARKER_INLET_SPECIES= (inlet, 0.0)
75+
%
76+
INC_OUTLET_TYPE= PRESSURE_OUTLET
77+
INC_OUTLET_DAMPING= 0.1
78+
MARKER_OUTLET= ( outlet, 0.0 )
79+
%
80+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
81+
%
82+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
83+
%
84+
CFL_NUMBER= 10.0
85+
CFL_REDUCTION_SPECIES= 1.0
86+
CFL_REDUCTION_TURB= 1.0
87+
CFL_ADAPT= YES
88+
CFL_ADAPT_PARAM= ( 0.95, 1.01, 1.0, 250, 1.0e-4, 0)
89+
90+
%
91+
ITER= 1
92+
%
93+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
94+
%
95+
LINEAR_SOLVER= FGMRES
96+
LINEAR_SOLVER_PREC= ILU
97+
LINEAR_SOLVER_ERROR= 1E-6
98+
LINEAR_SOLVER_ITER= 10
99+
100+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
101+
%
102+
CONV_NUM_METHOD_FLOW= FDS
103+
MUSCL_FLOW= NO
104+
SLOPE_LIMITER_FLOW = NONE
105+
TIME_DISCRE_FLOW= EULER_IMPLICIT
106+
%
107+
% -------------------- SCALAR TRANSPORT ---------------------------------------%
108+
%
109+
KIND_SCALAR_MODEL= SPECIES_TRANSPORT
110+
%DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
111+
DIFFUSIVITY_MODEL= CONSTANT_SCHMIDT
112+
%DIFFUSIVITY_MODEL= CONSTANT_LEWIS
113+
SCHMIDT_NUMBER_LAMINAR= 1.0
114+
CONSTANT_LEWIS_NUMBER= 1.0 1.0
115+
DIFFUSIVITY_CONSTANT= 5.0e-5
116+
117+
% according to the paper
118+
SCHMIDT_NUMBER_TURBULENT= 0.7
119+
120+
%
121+
CONV_NUM_METHOD_SPECIES= BOUNDED_SCALAR
122+
MUSCL_SPECIES= NO
123+
SLOPE_LIMITER_SPECIES = NONE
124+
%
125+
TIME_DISCRE_SPECIES= EULER_IMPLICIT
126+
%
127+
SPECIES_INIT= 0.0
128+
SPECIES_CLIPPING= YES
129+
SPECIES_CLIPPING_MIN= 0.0
130+
SPECIES_CLIPPING_MAX= 1.0
131+
%
132+
% -------------------- TURBULENT TRANSPORT ---------------------------------------%
133+
%
134+
CONV_NUM_METHOD_TURB= BOUNDED_SCALAR
135+
MUSCL_TURB= NO
136+
%
137+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
138+
%
139+
CONV_FIELD= RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TKE, RMS_SPECIES
140+
CONV_RESIDUAL_MINVAL= -12
141+
CONV_STARTITER= 10
142+
%
143+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
144+
%
145+
MESH_FILENAME= psi.su2
146+
%
147+
SCREEN_OUTPUT= INNER_ITER WALL_TIME \
148+
RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION RMS_SPECIES_0
149+
150+
SCREEN_WRT_FREQ_INNER= 1
151+
%
152+
HISTORY_OUTPUT= ITER RMS_RES LINSOL SPECIES_COEFF SPECIES_COEFF_SURF
153+
CONV_FILENAME= history
154+
MARKER_ANALYZE= gas_inlet, air_axial_inlet, outlet
155+
MARKER_ANALYZE_AVERAGE= AREA
156+
%
157+
OUTPUT_FILES= RESTART, PARAVIEW_MULTIBLOCK
158+
VOLUME_OUTPUT= RESIDUAL, PRIMITIVE, SPECIES_UDS_0
159+
OUTPUT_WRT_FREQ= 100
160+
%
161+
READ_BINARY_RESTART= YES
162+
RESTART_FILENAME= restart
163+
SOLUTION_FILENAME= solution
164+
%
165+
WRT_PERFORMANCE= YES
166+
PYTHON_CUSTOM_SOURCE= YES

0 commit comments

Comments
 (0)