Skip to content

Commit 0b01500

Browse files
frx-wintermutefrx-wintermute
committed
restore QuickStart subdirectory
...with the mesh and config file.
1 parent 46b42e3 commit 0b01500

4 files changed

Lines changed: 16054 additions & 0 deletions

File tree

‎QuickStart/.gitignore‎

Lines changed: 17 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,17 @@
1+
# In the QuickStart/ directory, we would like to
2+
# ignore any of the data needed to run the case
3+
# or their output data, i.e., meshes, solution
4+
# files, restart files, CSVs, etc.
5+
6+
# Things appearing in QuickStart/ to ignore:
7+
*.dat
8+
*.csv
9+
*.plt
10+
*.szplt
11+
*.vtu
12+
*.log
13+
*.out
14+
*.aux
15+
*.pdf
16+
*.png
17+
config_*.cfg

‎QuickStart/inv_NACA0012.cfg‎

Lines changed: 327 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,327 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Transonic inviscid flow around a NACA0012 airfoil %
5+
% Author: Thomas D. Economon %
6+
% Institution: Stanford University %
7+
% Date: 2014.06.11 %
8+
% File Version 7.5.1 "Blackbird" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
12+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
13+
%
14+
% Physical governing equations (EULER, NAVIER_STOKES,
15+
% WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
16+
% POISSON_EQUATION)
17+
SOLVER= EULER
18+
%
19+
% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT)
20+
MATH_PROBLEM= DIRECT
21+
%
22+
% Restart solution (NO, YES)
23+
RESTART_SOL= NO
24+
25+
% ----------- COMPRESSIBLE AND INCOMPRESSIBLE FREE-STREAM DEFINITION ----------%
26+
%
27+
% Mach number (non-dimensional, based on the free-stream values)
28+
MACH_NUMBER= 0.8
29+
%
30+
% Angle of attack (degrees)
31+
AOA= 1.25
32+
%
33+
% Free-stream pressure (101325.0 N/m^2 by default, only Euler flows)
34+
FREESTREAM_PRESSURE= 101325.0
35+
%
36+
% Free-stream temperature (273.15 K by default)
37+
FREESTREAM_TEMPERATURE= 273.15
38+
39+
% -------------- COMPRESSIBLE AND INCOMPRESSIBLE FLUID CONSTANTS --------------%
40+
%
41+
% Ratio of specific heats (1.4 (air), only for compressible flows)
42+
GAMMA_VALUE= 1.4
43+
%
44+
% Specific gas constant (287.87 J/kg*K (air), only for compressible flows)
45+
GAS_CONSTANT= 287.87
46+
47+
% ---------------------- REFERENCE VALUE DEFINITION ---------------------------%
48+
%
49+
% Reference origin for moment computation
50+
REF_ORIGIN_MOMENT_X = 0.25
51+
REF_ORIGIN_MOMENT_Y = 0.00
52+
REF_ORIGIN_MOMENT_Z = 0.00
53+
%
54+
% Reference length for pitching, rolling, and yawing non-dimensional moment
55+
REF_LENGTH= 1.0
56+
%
57+
% Reference area for force coefficients (0 implies automatic calculation)
58+
REF_AREA= 1.0
59+
%
60+
% Flow non-dimensionalization (DIMENSIONAL, FREESTREAM_PRESS_EQ_ONE,
61+
% FREESTREAM_VEL_EQ_MACH, FREESTREAM_VEL_EQ_ONE)
62+
REF_DIMENSIONALIZATION= DIMENSIONAL
63+
64+
% ----------------------- BOUNDARY CONDITION DEFINITION -----------------------%
65+
%
66+
% Marker of the Euler boundary (NONE = no marker)
67+
MARKER_EULER= ( airfoil )
68+
%
69+
% Marker of the far field (NONE = no marker)
70+
MARKER_FAR= ( farfield )
71+
72+
% ------------------------ SURFACES IDENTIFICATION ----------------------------%
73+
%
74+
% Marker(s) of the surface in the surface flow solution file
75+
MARKER_PLOTTING = ( airfoil )
76+
%
77+
% Marker(s) of the surface where the non-dimensional coefficients are evaluated.
78+
MARKER_MONITORING = ( airfoil )
79+
%
80+
% Marker(s) of the surface where obj. func. (design problem) will be evaluated
81+
MARKER_DESIGNING = ( airfoil )
82+
83+
% ------------- COMMON PARAMETERS TO DEFINE THE NUMERICAL METHOD --------------%
84+
%
85+
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
86+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
87+
%
88+
% Objective function in optimization problem (DRAG, LIFT, SIDEFORCE, MOMENT_X,
89+
% MOMENT_Y, MOMENT_Z, EFFICIENCY,
90+
% EQUIVALENT_AREA, NEARFIELD_PRESSURE,
91+
% FORCE_X, FORCE_Y, FORCE_Z, THRUST,
92+
% TORQUE, FREE_SURFACE, TOTAL_HEATFLUX,
93+
% MAXIMUM_HEATFLUX, INVERSE_DESIGN_PRESSURE,
94+
% INVERSE_DESIGN_HEATFLUX)
95+
OBJECTIVE_FUNCTION= DRAG
96+
%
97+
% Courant-Friedrichs-Lewy condition of the finest grid
98+
CFL_NUMBER= 1e3
99+
%
100+
% Adaptive CFL number (NO, YES)
101+
CFL_ADAPT= NO
102+
%
103+
% Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
104+
% CFL max value )
105+
CFL_ADAPT_PARAM= ( 0.1, 2.0, 10.0, 1e10 )
106+
%
107+
% Number of total iterations
108+
ITER= 250
109+
110+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
111+
%
112+
% Linear solver for implicit formulations (BCGSTAB, FGMRES)
113+
LINEAR_SOLVER= FGMRES
114+
%
115+
% Preconditioner of the Krylov linear solver (JACOBI, LINELET, LU_SGS)
116+
LINEAR_SOLVER_PREC= ILU
117+
%
118+
% Minimum error of the linear solver for implicit formulations
119+
LINEAR_SOLVER_ERROR= 1E-10
120+
%
121+
% Max number of iterations of the linear solver for the implicit formulation
122+
LINEAR_SOLVER_ITER= 10
123+
124+
% -------------------------- MULTIGRID PARAMETERS -----------------------------%
125+
%
126+
% Multi-Grid Levels (0 = no multi-grid)
127+
MGLEVEL= 3
128+
%
129+
% Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE)
130+
MGCYCLE= W_CYCLE
131+
%
132+
% Multi-Grid PreSmoothing Level
133+
MG_PRE_SMOOTH= ( 1, 2, 3, 3 )
134+
%
135+
% Multi-Grid PostSmoothing Level
136+
MG_POST_SMOOTH= ( 0, 0, 0, 0 )
137+
%
138+
% Jacobi implicit smoothing of the correction
139+
MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 )
140+
%
141+
% Damping factor for the residual restriction
142+
MG_DAMP_RESTRICTION= 1.0
143+
%
144+
% Damping factor for the correction prolongation
145+
MG_DAMP_PROLONGATION= 1.0
146+
147+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
148+
%
149+
% Convective numerical method
150+
%
151+
CONV_NUM_METHOD_FLOW= JST
152+
%
153+
% 2nd and 4th order artificial dissipation coefficients
154+
JST_SENSOR_COEFF= ( 0.5, 0.02 )
155+
%
156+
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
157+
TIME_DISCRE_FLOW= EULER_IMPLICIT
158+
159+
% ---------------- ADJOINT-FLOW NUMERICAL METHOD DEFINITION -------------------%
160+
%
161+
% Convective numerical method (JST, LAX-FRIEDRICH, ROE)
162+
CONV_NUM_METHOD_ADJFLOW= JST
163+
%
164+
% Reduction factor of the CFL coefficient in the adjoint problem
165+
CFL_REDUCTION_ADJFLOW= 0.01
166+
%
167+
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT)
168+
TIME_DISCRE_ADJFLOW= EULER_IMPLICIT
169+
170+
% ----------------------- DESIGN VARIABLE PARAMETERS --------------------------%
171+
%
172+
% Kind of deformation (NO_DEFORMATION, TRANSLATION, ROTATION, SCALE,
173+
% FFD_SETTING, FFD_NACELLE
174+
% FFD_CONTROL_POINT, FFD_CAMBER, FFD_THICKNESS, FFD_TWIST
175+
% FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D, FFD_TWIST_2D,
176+
% HICKS_HENNE, SURFACE_BUMP)
177+
DV_KIND= HICKS_HENNE
178+
%
179+
% Marker of the surface in which we are going apply the shape deformation
180+
DV_MARKER= ( airfoil )
181+
%
182+
% Parameters of the shape deformation
183+
% - NO_DEFORMATION ( 1.0 )
184+
% - TRANSLATION ( x_Disp, y_Disp, z_Disp ), as a unit vector
185+
% - ROTATION ( x_Orig, y_Orig, z_Orig, x_End, y_End, z_End )
186+
% - SCALE ( 1.0 )
187+
% - ANGLE_OF_ATTACK ( 1.0 )
188+
% - FFD_SETTING ( 1.0 )
189+
% - FFD_CONTROL_POINT ( FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Disp, y_Disp, z_Disp )
190+
% - FFD_NACELLE ( FFD_BoxTag, rho_Ind, theta_Ind, phi_Ind, rho_Disp, phi_Disp )
191+
% - FFD_GULL ( FFD_BoxTag, j_Ind )
192+
% - FFD_ANGLE_OF_ATTACK ( FFD_BoxTag, 1.0 )
193+
% - FFD_CAMBER ( FFD_BoxTag, i_Ind, j_Ind )
194+
% - FFD_THICKNESS ( FFD_BoxTag, i_Ind, j_Ind )
195+
% - FFD_TWIST ( FFD_BoxTag, j_Ind, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End )
196+
% - FFD_CONTROL_POINT_2D ( FFD_BoxTag, i_Ind, j_Ind, x_Disp, y_Disp )
197+
% - FFD_CAMBER_2D ( FFD_BoxTag, i_Ind )
198+
% - FFD_THICKNESS_2D ( FFD_BoxTag, i_Ind )
199+
% - FFD_TWIST_2D ( FFD_BoxTag, x_Orig, y_Orig )
200+
% - HICKS_HENNE ( Lower Surface (0)/Upper Surface (1)/Only one Surface (2), x_Loc )
201+
% - SURFACE_BUMP ( x_Start, x_End, x_Loc )
202+
DV_PARAM= ( 1, 0.5 )
203+
%
204+
% Value of the shape deformation
205+
DV_VALUE= 0.01
206+
207+
% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
208+
%
209+
% Number of smoothing iterations for FEA mesh deformation
210+
DEFORM_LINEAR_SOLVER_ITER= 500
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+
% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME,
222+
% WALL_DISTANCE, CONSTANT_STIFFNESS)
223+
DEFORM_STIFFNESS_TYPE= INVERSE_VOLUME
224+
225+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
226+
%
227+
CONV_FIELD= RMS_DENSITY
228+
%
229+
% Min value of the residual (log10 of the residual)
230+
CONV_RESIDUAL_MINVAL= -8
231+
%
232+
% Start Cauchy criteria at iteration number
233+
CONV_STARTITER= 10
234+
%
235+
% Number of elements to apply the criteria
236+
CONV_CAUCHY_ELEMS= 100
237+
%
238+
% Epsilon to control the series convergence
239+
CONV_CAUCHY_EPS= 1E-6
240+
%
241+
SCREEN_OUTPUT=(INNER_ITER, WALL_TIME, RMS_RES, LIFT, DRAG, CAUCHY_SENS_PRESS, CAUCHY_DRAG RMS_ADJ_DENSITY RMS_ADJ_ENERGY)
242+
243+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
244+
% Mesh input file
245+
MESH_FILENAME= mesh_NACA0012_inv.su2
246+
%
247+
% Mesh input file format (SU2, CGNS, NETCDF_ASCII)
248+
MESH_FORMAT= SU2
249+
%
250+
% Mesh output file
251+
MESH_OUT_FILENAME= mesh_out.su2
252+
%
253+
% Restart flow input file
254+
SOLUTION_FILENAME= solution_flow.dat
255+
%
256+
% Restart adjoint input file
257+
SOLUTION_ADJ_FILENAME= solution_adj.dat
258+
%
259+
% Output file format (TECPLOT, CSV)
260+
TABULAR_FORMAT= CSV
261+
%
262+
% Output file convergence history (w/o extension)
263+
CONV_FILENAME= history
264+
%
265+
% Output file restart flow
266+
RESTART_FILENAME= restart_flow.dat
267+
%
268+
% Output file restart adjoint
269+
RESTART_ADJ_FILENAME= restart_adj.dat
270+
%
271+
% Output file flow (w/o extension) variables
272+
VOLUME_FILENAME= flow
273+
%
274+
% Output file adjoint (w/o extension) variables
275+
VOLUME_ADJ_FILENAME= adjoint
276+
%
277+
% Output Objective function gradient (using continuous adjoint)
278+
GRAD_OBJFUNC_FILENAME= of_grad.dat
279+
%
280+
% Output file surface flow coefficient (w/o extension)
281+
SURFACE_FILENAME= surface_flow
282+
%
283+
% Output file surface adjoint coefficient (w/o extension)
284+
SURFACE_ADJ_FILENAME= surface_adjoint
285+
%
286+
% Writing solution file frequency
287+
OUTPUT_WRT_FREQ= 250
288+
%
289+
% Output file format
290+
OUTPUT_FILES= (RESTART, PARAVIEW, SURFACE_CSV)
291+
292+
% --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------%
293+
% Available flow based objective functions or constraint functions
294+
% DRAG, LIFT, SIDEFORCE, EFFICIENCY,
295+
% FORCE_X, FORCE_Y, FORCE_Z,
296+
% MOMENT_X, MOMENT_Y, MOMENT_Z,
297+
% THRUST, TORQUE, FIGURE_OF_MERIT,
298+
% EQUIVALENT_AREA, NEARFIELD_PRESSURE,
299+
% TOTAL_HEATFLUX, MAXIMUM_HEATFLUX,
300+
% INVERSE_DESIGN_PRESSURE, INVERSE_DESIGN_HEATFLUX,
301+
%
302+
% Available geometrical based objective functions or constraint functions
303+
% AIRFOIL_AREA, AIRFOIL_THICKNESS, AIRFOIL_CHORD, AIRFOIL_TOC, AIRFOIL_AOA,
304+
% WING_VOLUME, WING_MIN_THICKNESS, WING_MAX_THICKNESS, WING_MAX_CHORD, WING_MIN_TOC, WING_MAX_TWIST, WING_MAX_CURVATURE, WING_MAX_DIHEDRAL
305+
% STATION#_WIDTH, STATION#_AREA, STATION#_THICKNESS, STATION#_CHORD, STATION#_TOC,
306+
% STATION#_TWIST (where # is the index of the station defined in GEO_LOCATION_STATIONS)
307+
%
308+
% Available design variables
309+
% HICKS_HENNE ( 1, Scale | Mark. List | Lower(0)/Upper(1) side, x_Loc )
310+
% NACA_4DIGITS ( 4, Scale | Mark. List | 1st digit, 2nd digit, 3rd and 4th digit )
311+
% TRANSLATION ( 5, Scale | Mark. List | x_Disp, y_Disp, z_Disp )
312+
% ROTATION ( 6, Scale | Mark. List | x_Axis, y_Axis, z_Axis, x_Turn, y_Turn, z_Turn )
313+
% FFD_CONTROL_POINT_2D ( 15, Scale | Mark. List | FFD_Box_ID, i_Ind, j_Ind, x_Mov, y_Mov )
314+
% FFD_CAMBER_2D ( 16, Scale | Mark. List | FFD_Box_ID, i_Ind )
315+
% FFD_THICKNESS_2D ( 17, Scale | Mark. List | FFD_Box_ID, i_Ind )
316+
%
317+
% Optimization objective function with scaling factor
318+
% ex= Objective * Scale
319+
OPT_OBJECTIVE= DRAG * 0.001
320+
%
321+
% Optimization constraint functions with scaling factors, separated by semicolons
322+
% ex= (Objective = Value ) * Scale, use '>','<','='
323+
OPT_CONSTRAINT= ( LIFT > 0.328188 ) * 0.001; ( MOMENT_Z > 0.034068 ) * 0.001; ( AIRFOIL_THICKNESS > 0.11 ) * 0.001
324+
%
325+
326+
% Optimization design variables, separated by semicolons
327+
DEFINITION_DV= ( 30, 1.0 | airfoil | 0, 0.05 ); ( 30, 1.0 | airfoil | 0, 0.10 ); ( 30, 1.0 | airfoil | 0, 0.15 ); ( 30, 1.0 | airfoil | 0, 0.20 ); ( 30, 1.0 | airfoil | 0, 0.25 ); ( 30, 1.0 | airfoil | 0, 0.30 ); ( 30, 1.0 | airfoil | 0, 0.35 ); ( 30, 1.0 | airfoil | 0, 0.40 ); ( 30, 1.0 | airfoil | 0, 0.45 ); ( 30, 1.0 | airfoil | 0, 0.50 ); ( 30, 1.0 | airfoil | 0, 0.55 ); ( 30, 1.0 | airfoil | 0, 0.60 ); ( 30, 1.0 | airfoil | 0, 0.65 ); ( 30, 1.0 | airfoil | 0, 0.70 ); ( 30, 1.0 | airfoil | 0, 0.75 ); ( 30, 1.0 | airfoil | 0, 0.80 ); ( 30, 1.0 | airfoil | 0, 0.85 ); ( 30, 1.0 | airfoil | 0, 0.90 ); ( 30, 1.0 | airfoil | 0, 0.95 ); ( 30, 1.0 | airfoil | 1, 0.05 ); ( 30, 1.0 | airfoil | 1, 0.10 ); ( 30, 1.0 | airfoil | 1, 0.15 ); ( 30, 1.0 | airfoil | 1, 0.20 ); ( 30, 1.0 | airfoil | 1, 0.25 ); ( 30, 1.0 | airfoil | 1, 0.30 ); ( 30, 1.0 | airfoil | 1, 0.35 ); ( 30, 1.0 | airfoil | 1, 0.40 ); ( 30, 1.0 | airfoil | 1, 0.45 ); ( 30, 1.0 | airfoil | 1, 0.50 ); ( 30, 1.0 | airfoil | 1, 0.55 ); ( 30, 1.0 | airfoil | 1, 0.60 ); ( 30, 1.0 | airfoil | 1, 0.65 ); ( 30, 1.0 | airfoil | 1, 0.70 ); ( 30, 1.0 | airfoil | 1, 0.75 ); ( 30, 1.0 | airfoil | 1, 0.80 ); ( 30, 1.0 | airfoil | 1, 0.85 ); ( 30, 1.0 | airfoil | 1, 0.90 ); ( 30, 1.0 | airfoil | 1, 0.95 )

0 commit comments

Comments
 (0)