Skip to content

Commit 498393d

Browse files
bigfootedbigfooted
committed
TFC psi case
1 parent ee9d357 commit 498393d

6 files changed

Lines changed: 63721 additions & 0 deletions

File tree

‎multiphysics/TFC_python/inlet.dat‎

Lines changed: 45 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,45 @@
1+
NMARK= 1
2+
MARKER_TAG= inlet
3+
NROW=40
4+
NCOL=9
5+
# COORD-X COORD-Y TEMPERATURE VELOCITY NORMAL-X NORMAL-Y TKE DISSIPATION SPECIES_0
6+
-1.000000000000000e-01 1.154166394211445e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
7+
-1.000000000000000e-01 0.000000000000000e+00 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
8+
-1.000000000000000e-01 2.204457814215873e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
9+
-1.000000000000000e-01 3.160223011172386e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
10+
-1.000000000000000e-01 4.029969345405003e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
11+
-1.000000000000000e-01 4.821438503131977e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
12+
-1.000000000000000e-01 5.541675434666160e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
13+
-1.000000000000000e-01 6.197091047630005e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
14+
-1.000000000000000e-01 6.793519252514839e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
15+
-1.000000000000000e-01 7.336268913835854e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
16+
-1.000000000000000e-01 7.830171102836587e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
17+
-1.000000000000000e-01 8.279622096002559e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
18+
-1.000000000000000e-01 8.688622500268596e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
19+
-1.000000000000000e-01 9.060812867677235e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
20+
-1.000000000000000e-01 9.399506107798567e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
21+
-1.000000000000000e-01 9.707716957055516e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
22+
-1.000000000000000e-01 9.988188827727904e-03 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
23+
-1.000000000000000e-01 1.024341822719099e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
24+
-1.000000000000000e-01 1.047567697992562e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
25+
-1.000000000000000e-01 1.068703244700517e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
26+
-1.000000000000000e-01 1.087936591733933e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
27+
-1.000000000000000e-01 1.105438937094192e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
28+
-1.000000000000000e-01 1.121366071745546e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
29+
-1.000000000000000e-01 1.135859764414100e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
30+
-1.000000000000000e-01 1.149049024534486e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
31+
-1.000000000000000e-01 1.161051251507272e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
32+
-1.000000000000000e-01 1.171973278455936e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
33+
-1.000000000000000e-01 1.181912322770104e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
34+
-1.000000000000000e-01 1.190956853210104e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
35+
-1.000000000000000e-01 1.199187375572167e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
36+
-1.000000000000000e-01 1.206677150914347e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 1.0e+04 0.000000000000000e+00
37+
-1.000000000000000e-01 1.213492847168445e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 2.000000000000000e+00 1.0e+04 0.000000000000000e+00
38+
-1.000000000000000e-01 1.219695130380266e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e+00 1.1e+04 0.000000000000000e+00
39+
-1.000000000000000e-01 1.225339207572129e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.500000000000000e+01 1.5e+04 0.000000000000000e+00
40+
-1.000000000000000e-01 1.230475317865302e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 2.800000000000000e+01 2.5e+04 0.000000000000000e+00
41+
-1.000000000000000e-01 1.235149178369867e-02 6.730000000000000e+02 4.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 3.000000000000000e+01 5.0e+04 0.000000000000000e+00
42+
-1.000000000000000e-01 1.239402391748979e-02 6.730000000000000e+02 3.900000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 2.800000000000000e+01 1.0e+05 0.000000000000000e+00
43+
-1.000000000000000e-01 1.243272816372596e-02 6.730000000000000e+02 3.500000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+01 1.9e+05 0.000000000000000e+00
44+
-1.000000000000000e-01 1.246794902452563e-02 6.730000000000000e+02 0.500000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e+00 2.0e+05 0.000000000000000e+00
45+
-1.000000000000000e-01 1.250000000000000e-02 6.730000000000000e+02 0.000000000000000e+01 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e-10 2.0e+05 0.000000000000000e+00

‎multiphysics/TFC_python/psi.cfg‎

Lines changed: 172 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,172 @@
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+
MGLEVEL= 0
21+
%
22+
23+
% temperature and density
24+
FREESTREAM_TEMPERATURE = 673
25+
FREESTREAM_DENSITY = 2.55
26+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
27+
%
28+
%INC_DENSITY_MODEL= CONSTANT
29+
INC_DENSITY_MODEL= VARIABLE
30+
INC_DENSITY_INIT= 2.55
31+
%
32+
INC_VELOCITY_INIT= (40.00, 0.0, 0.0 )
33+
%
34+
INC_ENERGY_EQUATION= NO
35+
INC_TEMPERATURE_INIT= 673.0
36+
%
37+
INC_NONDIM= DIMENSIONAL
38+
%
39+
% -------------------- FLUID PROPERTIES ------------------------------------- %
40+
%
41+
%FLUID_MODEL= CONSTANT_DENSITY
42+
FLUID_MODEL= INC_IDEAL_GAS
43+
%
44+
CONDUCTIVITY_MODEL= CONSTANT_CONDUCTIVITY
45+
THERMAL_CONDUCTIVITY_CONSTANT= 0.0357
46+
%
47+
PRANDTL_LAM= 0.72
48+
TURBULENT_CONDUCTIVITY_MODEL= NONE
49+
PRANDTL_TURB= 0.90
50+
%
51+
VISCOSITY_MODEL= SUTHERLAND
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 = 1000.0
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 )
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.01
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.01
78+
MARKER_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+
89+
%
90+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
91+
%
92+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
93+
%
94+
CFL_NUMBER= 10.0
95+
CFL_REDUCTION_SPECIES= 1.0
96+
CFL_REDUCTION_TURB= 1.0
97+
CFL_ADAPT= YES
98+
CFL_ADAPT_PARAM= ( 0.95, 1.01, 5.0, 250)
99+
%
100+
ITER= 1
101+
%
102+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
103+
%
104+
LINEAR_SOLVER= FGMRES
105+
LINEAR_SOLVER_PREC= ILU
106+
LINEAR_SOLVER_ERROR= 1E-3
107+
LINEAR_SOLVER_ITER= 5
108+
109+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
110+
%
111+
CONV_NUM_METHOD_FLOW= FDS
112+
MUSCL_FLOW= NO
113+
SLOPE_LIMITER_FLOW = NONE
114+
TIME_DISCRE_FLOW= EULER_IMPLICIT
115+
%
116+
% -------------------- SCALAR TRANSPORT ---------------------------------------%
117+
%
118+
KIND_SCALAR_MODEL= SPECIES_TRANSPORT
119+
DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
120+
SCHMIDT_NUMBER_LAMINAR= 1.0
121+
DIFFUSIVITY_CONSTANT= 7.56e-5
122+
123+
% according to the paper
124+
SCHMIDT_NUMBER_TURBULENT= 0.7
125+
126+
%
127+
CONV_NUM_METHOD_SPECIES= BOUNDED_SCALAR
128+
MUSCL_SPECIES= NO
129+
SLOPE_LIMITER_SPECIES = NONE
130+
%
131+
TIME_DISCRE_SPECIES= EULER_IMPLICIT
132+
%
133+
SPECIES_INIT= 0.0
134+
SPECIES_CLIPPING= YES
135+
SPECIES_CLIPPING_MIN= 0.0
136+
SPECIES_CLIPPING_MAX= 1.0
137+
%
138+
% -------------------- TURBULENT TRANSPORT ---------------------------------------%
139+
%
140+
CONV_NUM_METHOD_TURB= BOUNDED_SCALAR
141+
MUSCL_TURB= NO
142+
%
143+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
144+
%
145+
CONV_FIELD= RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TKE, RMS_SPECIES
146+
CONV_RESIDUAL_MINVAL= -12
147+
CONV_STARTITER= 10
148+
%
149+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
150+
%
151+
MESH_FILENAME= psi.su2
152+
%
153+
SCREEN_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
155+
156+
SCREEN_WRT_FREQ_INNER= 1
157+
%
158+
HISTORY_OUTPUT= ITER RMS_RES LINSOL SPECIES_COEFF SPECIES_COEFF_SURF
159+
CONV_FILENAME= history
160+
MARKER_ANALYZE= gas_inlet, air_axial_inlet, outlet
161+
MARKER_ANALYZE_AVERAGE= AREA
162+
%
163+
OUTPUT_FILES= RESTART,RESTART_ASCII, PARAVIEW_MULTIBLOCK
164+
VOLUME_OUTPUT= RESIDUAL, SOLUTION, PRIMITIVE, SPECIES_UDS_0, RANK, TIMESTEP
165+
OUTPUT_WRT_FREQ= 100
166+
%
167+
READ_BINARY_RESTART= YES
168+
RESTART_FILENAME= restart
169+
SOLUTION_FILENAME= solution
170+
%
171+
WRT_PERFORMANCE= YES
172+
PYTHON_CUSTOM_SOURCE= YES

‎multiphysics/TFC_python/psi.geo‎

Lines changed: 85 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,85 @@
1+
//+
2+
Point(1) = {0, 0, 0, 1.0};
3+
Point(2) = {0.18, 0, 0, 1.0};
4+
Point(3) = {0.18, 0.0375, 0, 1.0};
5+
Point(4) = {0, 0.0375, 0, 1.0};
6+
Point(5) = {0, 0.0125, 0, 1.0};
7+
Point(6) = {-0.1, 0.0125, 0, 1.0};
8+
Point(7) = {-0.1, 0.0, 0, 1.0};
9+
Point(8) = {0.18, 0.008, -0, 1.0};
10+
Point(9) = {-0.05, 0.0, 0, 1.0};
11+
Point(10) = {-0.05, 0.0125, 0, 1.0};
12+
Point(11) = {0.40, 0.0, 0, 1.0};
13+
Point(12) = {0.40, 0.0375, 0, 1.0};
14+
Point(13) = {0.40, 0.008, 0, 1.0};
15+
//+
16+
Line(1) = {7, 6};
17+
Line(2) = {1, 9};
18+
Line(3) = {1, 5};
19+
Line(4) = {5,10};
20+
Line(5) = {1, 2};
21+
Line(6) = {2, 8};
22+
Line(7) = {5, 8};
23+
Line(8) = {5, 4};
24+
Line(9) = {4, 3};
25+
Line(10) = {3, 8};
26+
Line(11) = {9,10};
27+
Line(12) = {7,9};
28+
Line(13) = {6,10};
29+
30+
Line(14) = {11,13};
31+
Line(15) = {13,12};
32+
Line(16) = {2,11};
33+
Line(17) = {8,13};
34+
Line(18) = {3,12};
35+
//+
36+
Curve Loop(1) = {13, -11, -12, 1};
37+
Plane Surface(1) = {1};
38+
//+
39+
Curve Loop(2) = {5, 6, -7, -3};
40+
Plane Surface(2) = {2};
41+
//+
42+
Curve Loop(3) = {-7, 8, 9, 10};
43+
Plane Surface(3) = {3};
44+
//+
45+
Curve Loop(4) = {4, -11, -2, 3};
46+
Plane Surface(4) = {4};
47+
//+
48+
Curve Loop(5) = {17, 15, -18, 10};
49+
Plane Surface(5) = {5};
50+
//+
51+
Curve Loop(6) = {16, 14, -17, -6};
52+
Plane Surface(6) = {6};
53+
//+
54+
Physical Curve("inlet", 11) = {1};
55+
Physical Curve("outlet", 12) = {14, 15};
56+
Physical Curve("wall_top", 13) = {9,18};
57+
Physical Curve("wall_side", 14) = {8};
58+
Physical Curve("wall_pipe", 15) = {4,13};
59+
Physical Curve("symmetry", 16) = {2, 5, 12,16};
60+
//+
61+
Physical Surface("interior", 17) = {1, 2, 3, 4,5,6};
62+
//+
63+
Transfinite Curve {2, 4} = 80 Using Progression 1.03;
64+
//+
65+
Transfinite Curve {5, 7, 9} = 200 Using Progression 1.01;
66+
//+
67+
Ny = 40;
68+
Transfinite Curve {1, 3, 11} = Ny Using Progression 0.91;
69+
Transfinite Curve {6} = Ny Using Progression 1.0;
70+
Transfinite Curve {14} = Ny Using Progression 1.0;
71+
//+
72+
Transfinite Curve { 8} = 50 Using Progression 1.08;
73+
Transfinite Curve {10} = 50 Using Progression 1;
74+
Transfinite Curve {15} = 50 Using Progression 1;
75+
Transfinite Curve {12, 13} = 40 Using Progression 1;
76+
Transfinite Curve {16, 17, 18} = 100 Using Progression 1.005;
77+
//+
78+
Transfinite Surface {1};
79+
Transfinite Surface {2};
80+
Transfinite Surface {3};
81+
Transfinite Surface {4};
82+
Transfinite Surface {5};
83+
Transfinite Surface {6};
84+
//+
85+
Recombine Surface {1, 3, 2, 4, 5, 6};

0 commit comments

Comments
 (0)