Skip to content

Commit e872149

Browse files
TobiKattmannTobiKattmann
committed
Add preliminary species transport files
1 parent 58ea309 commit e872149

17 files changed

Lines changed: 19618 additions & 0 deletions

‎incompressible_flow/Inc_Species_Transport/1__FFD-box-writing/mesh_out.su2‎

Lines changed: 7244 additions & 0 deletions
Large diffs are not rendered by default.

‎incompressible_flow/Inc_Species_Transport/1__FFD-box-writing/primitiveVenturi.su2‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
../primitiveVenturi.su2

‎incompressible_flow/Inc_Species_Transport/1__FFD-box-writing/species3_primitiveVenturi.cfg‎

Lines changed: 231 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,231 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Species mixing with 3 species, i.e. 2 transport equations %
5+
% Author: T. Kattmann %
6+
% Institution: Bosch Thermotechniek B.V. %
7+
% Date: 2021/10/14 %
8+
% File Version 7.2.1 "Blackbird" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
12+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
13+
%
14+
SOLVER= INC_RANS
15+
KIND_TURB_MODEL= SST
16+
%
17+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
18+
%
19+
INC_DENSITY_MODEL= CONSTANT
20+
INC_DENSITY_INIT= 1.1766
21+
%
22+
INC_VELOCITY_INIT= ( 1.00, 0.0, 0.0 )
23+
%
24+
INC_ENERGY_EQUATION= YES
25+
INC_TEMPERATURE_INIT= 300.0
26+
%
27+
INC_NONDIM= INITIAL_VALUES
28+
%
29+
% -------------------- FLUID PROPERTIES ------------------------------------- %
30+
%
31+
FLUID_MODEL= CONSTANT_DENSITY
32+
%
33+
CONDUCTIVITY_MODEL= CONSTANT_CONDUCTIVITY
34+
THERMAL_CONDUCTIVITY_CONSTANT= 0.0357
35+
%
36+
PRANDTL_LAM= 0.72
37+
TURBULENT_CONDUCTIVITY_MODEL= NONE
38+
PRANDTL_TURB= 0.90
39+
%
40+
VISCOSITY_MODEL= CONSTANT_VISCOSITY
41+
MU_CONSTANT= 1.716E-5
42+
%
43+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
44+
%
45+
MARKER_HEATFLUX= ( wall, 0.0 )
46+
MARKER_SYM= ( axis )
47+
%
48+
SPECIFIED_INLET_PROFILE= NO
49+
INLET_FILENAME= inlet_venturi.dat
50+
INC_INLET_TYPE= VELOCITY_INLET VELOCITY_INLET
51+
MARKER_INLET= ( gas_inlet, 300, 1.0, 1.0, 0.0, 0.0,\
52+
air_axial_inlet, 300, 1.0, 0.0, -1.0, 0.0 )
53+
SPECIES_USE_STRONG_BC= NO
54+
MARKER_INLET_SPECIES= (gas_inlet, 0.5, 0.5,\
55+
air_axial_inlet, 0.6, 0.0 )
56+
%
57+
INC_OUTLET_TYPE= PRESSURE_OUTLET
58+
MARKER_OUTLET= ( outlet, 0.0 )
59+
%
60+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
61+
%
62+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
63+
%
64+
CFL_NUMBER= 2000
65+
CFL_REDUCTION_SPECIES= 1.0
66+
CFL_REDUCTION_TURB= 1.0
67+
%
68+
% Run commented Iter for good results
69+
ITER= 1000
70+
%
71+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
72+
%
73+
LINEAR_SOLVER= FGMRES
74+
LINEAR_SOLVER_PREC= ILU
75+
LINEAR_SOLVER_ERROR= 1E-8
76+
LINEAR_SOLVER_ITER= 5
77+
78+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
79+
%
80+
CONV_NUM_METHOD_FLOW= FDS
81+
MUSCL_FLOW= YES
82+
SLOPE_LIMITER_FLOW = NONE
83+
TIME_DISCRE_FLOW= EULER_IMPLICIT
84+
%
85+
% -------------------- SCALAR TRANSPORT ---------------------------------------%
86+
%
87+
KIND_SCALAR_MODEL= PASSIVE_SCALAR
88+
DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
89+
DIFFUSIVITY_CONSTANT= 0.001
90+
%
91+
CONV_NUM_METHOD_SPECIES= SCALAR_UPWIND
92+
MUSCL_SPECIES= NO
93+
SLOPE_LIMITER_SPECIES = NONE
94+
%
95+
TIME_DISCRE_SPECIES= EULER_IMPLICIT
96+
%
97+
SPECIES_INIT= 1.0, 0.0
98+
SPECIES_CLIPPING= YES
99+
SPECIES_CLIPPING_MIN= 0.0, 0.0
100+
SPECIES_CLIPPING_MAX= 1.0, 1.0
101+
%
102+
% -------------------- TURBULENT TRANSPORT ---------------------------------------%
103+
%
104+
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
105+
MUSCL_TURB= NO
106+
%
107+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
108+
%
109+
CONV_FIELD= RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TKE, RMS_SPECIES
110+
CONV_RESIDUAL_MINVAL= -18
111+
CONV_STARTITER= 10
112+
%
113+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
114+
%
115+
MESH_FILENAME= primitiveVenturi.su2
116+
%
117+
SCREEN_OUTPUT= INNER_ITER WALL_TIME \
118+
RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION RMS_SPECIES_0 RMS_SPECIES_1 \
119+
LINSOL_ITER LINSOL_RESIDUAL \
120+
LINSOL_ITER_TURB LINSOL_RESIDUAL_TURB \
121+
LINSOL_ITER_SPECIES LINSOL_RESIDUAL_SPECIES \
122+
SURFACE_SPECIES_VARIANCE
123+
SCREEN_WRT_FREQ_INNER= 10
124+
%
125+
HISTORY_OUTPUT= ITER RMS_RES LINSOL SPECIES_COEFF SPECIES_COEFF_SURF
126+
CONV_FILENAME= history
127+
MARKER_ANALYZE= gas_inlet, air_axial_inlet, outlet
128+
MARKER_ANALYZE_AVERAGE= AREA
129+
%
130+
OUTPUT_FILES= RESTART_ASCII, PARAVIEW_MULTIBLOCK
131+
VOLUME_OUTPUT= RESIDUAL, PRIMITIVE
132+
OUTPUT_WRT_FREQ= 1000
133+
%
134+
RESTART_SOL= NO
135+
READ_BINARY_RESTART= NO
136+
RESTART_FILENAME= restart
137+
SOLUTION_FILENAME= solution
138+
%
139+
WRT_PERFORMANCE= YES
140+
%
141+
% -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------%
142+
%
143+
% Tolerance of the Free-Form Deformation point inversion
144+
FFD_TOLERANCE= 1E-10
145+
%
146+
% Maximum number of iterations in the Free-Form Deformation point inversion
147+
FFD_ITERATIONS= 500
148+
%
149+
% FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4,
150+
% X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8)
151+
% 2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0,
152+
% 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
153+
% Start at the lowest leftest corner and turn counter-clockwise
154+
FFD_DEFINITION= (BOX, \
155+
0.065, 0.01, 0.0, \
156+
0.15, 0.01, 0.0 \
157+
0.15, 0.02, 0.0, \
158+
0.065, 0.02, 0.0, \
159+
0.0, 0.0 ,0.0, \
160+
0.0 ,0.0, 0.0, \
161+
0.0, 0.0, 0.0, \
162+
0.0, 0.0, 0.0 )
163+
%
164+
% FFD box degree: 3D case (x_degree, y_degree, z_degree)
165+
% 2D case (x_degree, y_degree, 0)
166+
FFD_DEGREE= (6, 1, 0)
167+
%
168+
% Surface grid continuity at the intersection with the faces of the FFD boxes.
169+
% To keep a particular level of surface continuity, SU2 automatically freezes the right
170+
% number of control point planes (NO_DERIVATIVE, 1ST_DERIVATIVE, 2ND_DERIVATIVE, USER_INPUT)
171+
FFD_CONTINUITY= NO_DERIVATIVE
172+
%
173+
% Definition of the FFD planes to be frozen in the FFD (x,y,z).
174+
% Value from 0 FFD degree in that direction. Pick a value larger than degree if you don't want to fix any plane.
175+
%FFD_FIX_I= (0,2,3)
176+
%FFD_FIX_J= (0,2,3)
177+
%FFD_FIX_K= (0,2,3)
178+
179+
% ----------------------- DESIGN VARIABLE PARAMETERS --------------------------%
180+
%
181+
% Kind of deformation (NO_DEFORMATION, SCALE_GRID, TRANSLATE_GRID, ROTATE_GRID,
182+
% FFD_SETTING, FFD_NACELLE,
183+
% FFD_CONTROL_POINT, FFD_CAMBER, FFD_THICKNESS, FFD_TWIST
184+
% FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D,
185+
% FFD_TWIST_2D, HICKS_HENNE, SURFACE_BUMP, SURFACE_FILE)
186+
DV_KIND= FFD_SETTING
187+
%
188+
% Marker of the surface in which we are going apply the shape deformation
189+
% NOTE: for deformation the outlet should be a MARKER_SYM to hinder the mesh being ripped apart.
190+
DV_MARKER= ( wall )
191+
%
192+
% Parameters of the shape deformation
193+
% - NO_DEFORMATION ( 1.0 )
194+
% - FFD_SETTING ( 1.0 )
195+
% - FFD_CONTROL_POINT ( FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Disp, y_Disp, z_Disp )
196+
DV_PARAM= ( 1.0 )
197+
%
198+
% Value of the shape deformation
199+
DV_VALUE= 0.01
200+
201+
% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
202+
%
203+
% Linear solver or smoother for implicit formulations (FGMRES, RESTARTED_FGMRES, BCGSTAB)
204+
DEFORM_LINEAR_SOLVER= FGMRES
205+
%
206+
% Preconditioner of the Krylov linear solver (ILU, LU_SGS, JACOBI)
207+
DEFORM_LINEAR_SOLVER_PREC= ILU
208+
%
209+
% Number of smoothing iterations for mesh deformation
210+
DEFORM_LINEAR_SOLVER_ITER= 1000
211+
%
212+
% Number of nonlinear deformation iterations (surface deformation increments)
213+
DEFORM_NONLINEAR_ITER= 1
214+
%
215+
% Minimum residual criteria for the linear solver convergence of grid deformation
216+
DEFORM_LINEAR_SOLVER_ERROR= 1E-14
217+
%
218+
% Print the residuals during mesh deformation to the console (YES, NO)
219+
DEFORM_CONSOLE_OUTPUT= YES
220+
%
221+
% Deformation coefficient (linear elasticity limits from -1.0 to 0.5, a larger
222+
% value is also possible)
223+
DEFORM_COEFF = 1E6
224+
%
225+
% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME,
226+
% WALL_DISTANCE, CONSTANT_STIFFNESS)
227+
DEFORM_STIFFNESS_TYPE= WALL_DISTANCE
228+
%
229+
% Deform the grid only close to the surface. It is possible to specify how much
230+
% of the volumetric grid is going to be deformed in meters or inches (1E6 by default)
231+
DEFORM_LIMIT = 1E6

‎incompressible_flow/Inc_Species_Transport/2__mesh-deform-test/primitiveVenturi.su2‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
../1__FFD-box-writing/mesh_out.su2

0 commit comments

Comments
 (0)