Skip to content

Commit 4b60859

Browse files
bigfootedbigfooted
committed
Merge branch 'tut_unst_cylinder' of https://github.com/su2code/Tutorials into tut_unst_cylinder
2 parents 475f10a + c66cf8b commit 4b60859

File tree

12 files changed

+386962
-0
lines changed

12 files changed

+386962
-0
lines changed

‎design/Inc_Turbulent_Bend_Wallfunctions/config.cfg‎

Lines changed: 222 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,222 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: turbulent 90 degree bend flow of Sudo (2001) %
5+
% Reynolds = 1.185*8.7*0.104/1.785e-5 ~ 60,000 %
6+
% Author: Nijso Beishuizen %
7+
% Date: 2024.01.13 %
8+
% File Version 8.0.0 "Harrier" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
12+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
13+
%
14+
SOLVER= INC_RANS
15+
INC_NONDIM= DIMENSIONAL
16+
KIND_TURB_MODEL= SST
17+
SST_OPTIONS= V2003m
18+
19+
% switch to RESTART_SOL= YES if you do not have an initial solution
20+
% RESTART_SOL= YES
21+
RESTART_SOL= YES
22+
%
23+
% ITER= 4000
24+
ITER= 100
25+
26+
%
27+
OBJECTIVE_FUNCTION= __OBJ_FUNC__
28+
OBJECTIVE_WEIGHT= 1.0
29+
%
30+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
31+
%
32+
INC_DENSITY_INIT= 1.185
33+
INC_DENSITY_REF = 1.185
34+
INC_VELOCITY_INIT= ( 0.0, 0.0, 1.0 )
35+
INC_VELOCITY_REF = 8.0
36+
INC_INLET_TYPE= VELOCITY_INLET
37+
INC_OUTLET_TYPE= PRESSURE_OUTLET
38+
%
39+
% --------------------------- VISCOSITY MODEL ---------------------------------%
40+
%
41+
VISCOSITY_MODEL= CONSTANT_VISCOSITY
42+
MU_CONSTANT= 1.785e-5
43+
%
44+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
45+
%
46+
47+
% ------------------------------- %
48+
SPECIFIED_INLET_PROFILE= YES
49+
INLET_FILENAME= inlet.dat
50+
% ------------------------------- %
51+
52+
%MARKER_HEATFLUX= ( wall_1, 0.0, wall_2,0.0, wall_bend,0.0)
53+
MARKER_HEATFLUX= ( wall,0.0)
54+
MARKER_INLET= ( inlet, 300.0, 8.7, 0.0, 0.0, 1.0 )
55+
MARKER_OUTLET= ( outlet, 0.0)
56+
MARKER_SYM= ( symmetry )
57+
%
58+
MARKER_MONITORING= ( outlet, inlet )
59+
MARKER_PLOTTING= ( wall, symmetry )
60+
% # for pressure drop we need outlet and inlet
61+
MARKER_ANALYZE= ( outlet, inlet )
62+
MARKER_ANALYZE_AVERAGE= MASSFLUX
63+
%
64+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
65+
%
66+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
67+
%__DIRECT__CFL_NUMBER= 300
68+
%__ADJOINT__CFL_NUMBER=200
69+
70+
CFL_ADAPT= NO
71+
72+
%
73+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
74+
%
75+
LINEAR_SOLVER= FGMRES
76+
LINEAR_SOLVER_PREC= ILU
77+
LINEAR_SOLVER_ERROR= 1E-06
78+
%__DIRECT__LINEAR_SOLVER_ITER= 5
79+
%__ADJOINT__LINEAR_SOLVER_ITER= 10
80+
81+
%
82+
% -------------------------- MULTIGRID PARAMETERS -----------------------------%
83+
%
84+
MGLEVEL= 0
85+
%
86+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
87+
%
88+
CONV_NUM_METHOD_FLOW= FDS
89+
MUSCL_FLOW= YES
90+
SLOPE_LIMITER_FLOW= NONE
91+
TIME_DISCRE_FLOW= EULER_IMPLICIT
92+
%
93+
% -------------------- TURBULENT NUMERICAL METHOD DEFINITION ------------------%
94+
%
95+
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
96+
MUSCL_TURB= NO
97+
%SLOPE_LIMITER_TURB= VENKATAKRISHNAN
98+
TIME_DISCRE_TURB= EULER_IMPLICIT
99+
FREESTREAM_TURBULENCEINTENSITY= 0.10
100+
FREESTREAM_TURB2LAMVISCRATIO= 100.0
101+
%
102+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
103+
%
104+
%__ADJOINT__CONV_FIELD= RMS_ADJ_PRESSURE
105+
%__DIRECT__CONV_FIELD= RMS_PRESSURE
106+
CONV_RESIDUAL_MINVAL= -12
107+
CONV_STARTITER= 10
108+
CONV_CAUCHY_ELEMS= 100
109+
CONV_CAUCHY_EPS= 1E-6
110+
%
111+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
112+
%
113+
MESH_FILENAME= sudo_coarse_FFD.su2
114+
SOLUTION_FILENAME= solution
115+
CONV_FILENAME= history
116+
RESTART_FILENAME= restart
117+
WRT_RESTART_OVERWRITE= YES
118+
WRT_VOLUME_OVERWRITE= YES
119+
WRT_SURFACE_OVERWRITE= YES
120+
WRT_PERFORMANCE= YES
121+
VOLUME_FILENAME= flow
122+
OUTPUT_WRT_FREQ= 100
123+
SCREEN_WRT_FREQ_INNER= 1
124+
SCREEN_WRT_FREQ_OUTER= 1
125+
% Note that this goes together with RESTART_ASCII
126+
READ_BINARY_RESTART= NO
127+
%
128+
OUTPUT_FILES= PARAVIEW_MULTIBLOCK RESTART_ASCII
129+
GRAD_OBJFUNC_FILENAME= of_grad.csv
130+
%__DIRECT__SCREEN_OUTPUT= INNER_ITER WALL_TIME RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION LINSOL_ITER LINSOL_RESIDUAL LINSOL_ITER_TURB LINSOL_RESIDUAL_TURB
131+
%__ADJOINT__SCREEN_OUTPUT= INNER_ITER WALL_TIME RMS_ADJ_PRESSURE RMS_ADJ_VELOCITY-X RMS_ADJ_VELOCITY-Y RMS_ADJ_TKE RMS_ADJ_DISSIPATION
132+
133+
% !!! Note that this determines the output of the objectives !!!
134+
HISTORY_OUTPUT= INNER_ITER RMS_RES FLOW_COEFF FLOW_COEFF_SURF
135+
136+
OUTPUT_PRECISION= 16
137+
138+
WRT_RESTART_COMPACT = YES
139+
140+
MESH_OUT_FILENAME= mesh_out
141+
%
142+
% -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------%
143+
%
144+
% Tolerance of the Free-Form Deformation point inversion
145+
FFD_TOLERANCE= 1E-10
146+
%
147+
% Maximum number of iterations in the Free-Form Deformation point inversion
148+
FFD_ITERATIONS= 500
149+
%
150+
% FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4,
151+
% X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8)
152+
% 2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0,
153+
% 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
154+
% Start at the lowest leftest corner and turn counter-clockwise
155+
%FFD_DEFINITION= (BOX, \
156+
%-0.06, 0.000, 0.26, \
157+
% 0.20, 0.000, 0.26 \
158+
% 0.20, 0.055, 0.26, \
159+
%-0.06, 0.055, 0.26, \
160+
%-0.06, 0.000 ,0.00, \
161+
% 0.20 ,0.000, 0.00, \
162+
% 0.20, 0.055, 0.00, \
163+
%-0.06, 0.055, 0.00 )
164+
%
165+
% FFD box degree: 3D case (i_degree, j_degree, k_degree)
166+
% 2D case (i_degree, j_degree, 0)
167+
% fixes the nodes on the surface plane that cross the domain
168+
FFD_CONTINUITY= 1ST_DERIVATIVE
169+
170+
%
171+
% # DV_KIND= FFD_SETTING
172+
% # 6 x 6 x 6 = 216 nodes
173+
% # 216 x 3 = 648 d.o.f. for all 3 directions
174+
%
175+
DV_KIND= __FFD_CTRL_PTS__
176+
%
177+
% Marker of the surface in which we are going apply the shape deformation
178+
%
179+
DV_MARKER= ( wall, symmetry)
180+
%
181+
% Parameters of the shape deformation
182+
% - FFD_SETTING ( 1.0 )
183+
% - FFD_CONTROL_POINT ( FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Disp, y_Disp, z_Disp )
184+
%
185+
DV_PARAM= __FFD_PARAM__
186+
%
187+
% 216 x 3 = 648 points
188+
%
189+
DV_VALUE= __FFD_PTS__
190+
%
191+
% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
192+
%
193+
% Linear solver or smoother for implicit formulations (FGMRES, RESTARTED_FGMRES, BCGSTAB)
194+
DEFORM_LINEAR_SOLVER= FGMRES
195+
%
196+
% Preconditioner of the Krylov linear solver (ILU, LU_SGS, JACOBI)
197+
DEFORM_LINEAR_SOLVER_PREC= ILU
198+
%
199+
% Number of smoothing iterations for mesh deformation
200+
DEFORM_LINEAR_SOLVER_ITER= 1000
201+
%
202+
% Number of nonlinear deformation iterations (surface deformation increments)
203+
DEFORM_NONLINEAR_ITER= 1
204+
%
205+
% Minimum residual criteria for the linear solver convergence of grid deformation
206+
DEFORM_LINEAR_SOLVER_ERROR= 1E-14
207+
%
208+
% Print the residuals during mesh deformation to the console (YES, NO)
209+
DEFORM_CONSOLE_OUTPUT= YES
210+
%
211+
% Deformation coefficient (linear elasticity limits from -1.0 to 0.5, a larger
212+
% value is also possible)
213+
DEFORM_COEFF = 1.0e6
214+
%
215+
% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME,
216+
% WALL_DISTANCE, CONSTANT_STIFFNESS)
217+
%DEFORM_STIFFNESS_TYPE= WALL_DISTANCE
218+
DEFORM_STIFFNESS_TYPE= CONSTANT_STIFFNESS
219+
%
220+
% Deform the grid only close to the surface. It is possible to specify how much
221+
% of the volumetric grid is going to be deformed in meters or inches (1E6 by default)
222+
DEFORM_LIMIT = 1E6

0 commit comments

Comments
 (0)