Skip to content

Commit bdf7b21

Browse files
oleburghardtoleburghardt
committed
Move updated config files for CHT tutorial to multiphysics/steady_cht.
1 parent 9ab93e4 commit bdf7b21

File tree

6 files changed

+99447
-0
lines changed

6 files changed

+99447
-0
lines changed

‎multiphysics/steady_cht/cht_2d_3cylinders.cfg‎

Lines changed: 118 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,118 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: 2D cylinder array with CHT couplings %
5+
% Author: O. Burghardt, T. Economon %
6+
% Institution: Chair for Scientific Computing, TU Kaiserslautern %
7+
% Date: August 8, 2019 %
8+
% File Version 6.0.1 "Falcon" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
12+
%
13+
% Physical governing equations (EULER, NAVIER_STOKES,
14+
% WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
15+
% POISSON_EQUATION)
16+
SOLVER= MULTIPHYSICS
17+
%
18+
% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT, DISCRETE_ADJOINT)
19+
MATH_PROBLEM= DIRECT
20+
%
21+
% Restart solution (NO, YES)
22+
RESTART_SOL= NO
23+
%
24+
% Configuration file list, one for each physical zone
25+
CONFIG_LIST= (flow_cylinder.cfg, solid_cylinder1.cfg, solid_cylinder2.cfg, solid_cylinder3.cfg)
26+
%
27+
% Definition of the interface
28+
MARKER_ZONE_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
29+
MARKER_CHT_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
30+
%
31+
% Objective function in gradient evaluation (DRAG, LIFT, SIDEFORCE, MOMENT_X,
32+
% MOMENT_Y, MOMENT_Z, EFFICIENCY,
33+
% EQUIVALENT_AREA, NEARFIELD_PRESSURE,
34+
% FORCE_X, FORCE_Y, FORCE_Z, THRUST,
35+
% TORQUE, TOTAL_HEATFLUX,
36+
% MAXIMUM_HEATFLUX, INVERSE_DESIGN_PRESSURE,
37+
% INVERSE_DESIGN_HEATFLUX, SURFACE_TOTAL_PRESSURE,
38+
% SURFACE_MASSFLOW, SURFACE_STATIC_PRESSURE, SURFACE_MACH)
39+
% For a weighted sum of objectives: separate by commas, add OBJECTIVE_WEIGHT and MARKER_MONITORING in matching order.
40+
OBJECTIVE_FUNCTION= TOTAL_HEATFLUX
41+
%
42+
% List of weighting values when using more than one OBJECTIVE_FUNCTION. Separate by commas and match with MARKER_MONITORING.
43+
OBJECTIVE_WEIGHT= 1.0
44+
%
45+
% Number of total iterations
46+
OUTER_ITER = 15000
47+
OUTPUT_WRT_FREQ = 15000
48+
%
49+
% Mesh input file
50+
MESH_FILENAME= mesh_cht_3cyl_ffd.su2
51+
MESH_OUT_FILENAME= mesh_cht_3cyl_out.su2
52+
%
53+
% Mesh input file format (SU2, CGNS, NETCDF_ASCII)
54+
MESH_FORMAT= SU2
55+
%
56+
% Output file format
57+
OUTPUT_FILES= (RESTART, TECPLOT, PARAVIEW, SURFACE_TECPLOT, SURFACE_PARAVIEW)
58+
%
59+
% Multizone convergence criteria
60+
CONV_RESIDUAL_MINVAL= -20
61+
62+
% -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------%
63+
%
64+
% Tolerance of the Free-Form Deformation point inversion
65+
FFD_TOLERANCE= 1E-12
66+
%
67+
% Maximum number of iterations in the Free-Form Deformation point inversion
68+
FFD_ITERATIONS= 500
69+
%
70+
% FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4,
71+
% X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8)
72+
% 2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0,
73+
% 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
74+
FFD_DEFINITION= (MAIN_BOX, -0.1, -0.6, 0.0, 1.1, -0.6, 0.0, 1.1, 0.6, 0.0, -0.1, 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
75+
%
76+
% FFD box degree: 3D case (x_degree, y_degree, z_degree)
77+
% 2D case (x_degree, y_degree, 0)
78+
FFD_DEGREE= ( 24, 1, 0)
79+
%
80+
% Surface continuity at the intersection with the FFD (1ST_DERIVATIVE, 2ND_DERIVATIVE)
81+
FFD_CONTINUITY= 2ND_DERIVATIVE
82+
83+
% ----------------------- DESIGN VARIABLE PARAMETERS --------------------------%
84+
%
85+
% Kind of deformation (NO_DEFORMATION, TRANSLATION, ROTATION, SCALE,
86+
% FFD_SETTING, FFD_NACELLE
87+
% FFD_CONTROL_POINT, FFD_CAMBER, FFD_THICKNESS, FFD_TWIST
88+
% FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D, FFD_TWIST_2D,
89+
% HICKS_HENNE, SURFACE_BUMP)
90+
DV_KIND= FFD_CONTROL_POINT_2D
91+
%
92+
% Marker of the surface in which we are going apply the shape deformation
93+
DV_MARKER= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
94+
%
95+
% Parameters of the shape deformation
96+
% - NO_DEFORMATION ( 1.0 )
97+
% - TRANSLATION ( x_Disp, y_Disp, z_Disp ), as a unit vector
98+
% - ROTATION ( x_Orig, y_Orig, z_Orig, x_End, y_End, z_End )
99+
% - SCALE ( 1.0 )
100+
% - ANGLE_OF_ATTACK ( 1.0 )
101+
% - FFD_SETTING ( 1.0 )
102+
% - FFD_CONTROL_POINT ( FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Disp, y_Disp, z_Disp )
103+
% - FFD_NACELLE ( FFD_BoxTag, rho_Ind, theta_Ind, phi_Ind, rho_Disp, phi_Disp )
104+
% - FFD_GULL ( FFD_BoxTag, j_Ind )
105+
% - FFD_ANGLE_OF_ATTACK ( FFD_BoxTag, 1.0 )
106+
% - FFD_CAMBER ( FFD_BoxTag, i_Ind, j_Ind )
107+
% - FFD_THICKNESS ( FFD_BoxTag, i_Ind, j_Ind )
108+
% - FFD_TWIST ( FFD_BoxTag, j_Ind, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End )
109+
% - FFD_CONTROL_POINT_2D ( FFD_BoxTag, i_Ind, j_Ind, x_Disp, y_Disp )
110+
% - FFD_CAMBER_2D ( FFD_BoxTag, i_Ind )
111+
% - FFD_THICKNESS_2D ( FFD_BoxTag, i_Ind )
112+
% - FFD_TWIST_2D ( FFD_BoxTag, x_Orig, y_Orig )
113+
% - HICKS_HENNE ( Lower Surface (0)/Upper Surface (1)/Only one Surface (2), x_Loc )
114+
% - SURFACE_BUMP ( x_Start, x_End, x_Loc )
115+
DV_PARAM= ( MAIN_BOX, 12, 1, 0.0, 1.0 )
116+
%
117+
% Value of the shape deformation
118+
DV_VALUE= 0.1

‎multiphysics/steady_cht/flow_cylinder.cfg‎

Lines changed: 234 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,234 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Steady incompressible laminar flow around heated cylinders %
5+
% %
6+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7+
8+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
9+
%
10+
% Physical governing equations (EULER, NAVIER_STOKES,
11+
% WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
12+
% POISSON_EQUATION)
13+
SOLVER= INC_NAVIER_STOKES
14+
%
15+
% If Navier-Stokes, kind of turbulent model (NONE, SA)
16+
KIND_TURB_MODEL= NONE
17+
%
18+
% Data written to history file
19+
WRT_ZONE_HIST= YES
20+
HISTORY_OUTPUT= (ITER, RMS_RES, HEAT)
21+
%
22+
% Number of inner iteration
23+
INNER_ITER=1
24+
25+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
26+
%
27+
% Min value of the residual (log10 of the residual)
28+
CONV_RESIDUAL_MINVAL= -20
29+
%
30+
% Start convergence criteria at iteration number
31+
CONV_STARTITER= 0
32+
33+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
34+
%
35+
% Farfield boundary marker(s) (NONE = no marker)
36+
MARKER_FAR= ( farfield )
37+
%
38+
% Marker(s) of the surface to be plotted or designed
39+
MARKER_PLOTTING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
40+
%
41+
% Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated
42+
MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
43+
44+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
45+
%
46+
% Density model within the incompressible flow solver.
47+
% Options are CONSTANT (default), BOUSSINESQ, or VARIABLE. If VARIABLE,
48+
% an appropriate fluid model must be selected.
49+
INC_DENSITY_MODEL= VARIABLE
50+
%
51+
% Solve the energy equation in the incompressible flow solver
52+
INC_ENERGY_EQUATION = YES
53+
%
54+
% Initial density for incompressible flows (1.2886 kg/m^3 by default)
55+
INC_DENSITY_INIT= 0.00042
56+
%
57+
% Initial velocity for incompressible flows (1.0,0,0 m/s by default)
58+
INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
59+
%
60+
% Initial temperature for incompressible flows that include the
61+
% energy equation (288.15 K by default). Value is ignored if
62+
% INC_ENERGY_EQUATION is false.
63+
INC_TEMPERATURE_INIT= 288.15
64+
%
65+
% Non-dimensionalization scheme for incompressible flows. Options are
66+
% INITIAL_VALUES (default), REFERENCE_VALUES, or DIMENSIONAL.
67+
% INC_*_REF values are ignored unless REFERENCE_VALUES is chosen.
68+
INC_NONDIM= DIMENSIONAL
69+
70+
% ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
71+
%
72+
% Fluid model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS,
73+
% CONSTANT_DENSITY, INC_IDEAL_GAS)
74+
FLUID_MODEL= INC_IDEAL_GAS
75+
%
76+
% Specific heat at constant pressure, Cp (1004.703 J/kg*K (air)).
77+
% Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
78+
SPECIFIC_HEAT_CP= 1004.703
79+
%
80+
% Molecular weight for an incompressible ideal gas (28.96 g/mol (air) default)
81+
% Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
82+
MOLECULAR_WEIGHT= 28.96
83+
84+
% --------------------------- VISCOSITY MODEL ---------------------------------%
85+
%
86+
% Viscosity model (SUTHERLAND, CONSTANT_VISCOSITY).
87+
VISCOSITY_MODEL= CONSTANT_VISCOSITY
88+
%
89+
% Molecular Viscosity that would be constant (1.716E-5 by default)
90+
MU_CONSTANT= 1.7893e-05
91+
%
92+
% Sutherland Viscosity Ref (1.716E-5 default value for AIR SI)
93+
MU_REF= 1.716E-5
94+
%
95+
% Sutherland Temperature Ref (273.15 K default value for AIR SI)
96+
MU_T_REF= 273.15
97+
%
98+
% Sutherland constant (110.4 default value for AIR SI)
99+
SUTHERLAND_CONSTANT= 110.4
100+
101+
% --------------------------- THERMAL CONDUCTIVITY MODEL ----------------------%
102+
%
103+
% Conductivity model (CONSTANT_CONDUCTIVITY, CONSTANT_PRANDTL).
104+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
105+
%
106+
% Molecular Thermal Conductivity that would be constant (0.0257 by default)
107+
KT_CONSTANT= 0.0257
108+
%
109+
% Laminar Prandtl number (0.72 (air), only for CONSTANT_PRANDTL)
110+
PRANDTL_LAM= 0.72
111+
%
112+
% Turbulent Prandtl number (0.9 (air), only for CONSTANT_PRANDTL)
113+
PRANDTL_TURB= 0.90
114+
115+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
116+
%
117+
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
118+
NUM_METHOD_GRAD= GREEN_GAUSS
119+
%
120+
% Courant-Friedrichs-Lewy condition of the finest grid
121+
CFL_NUMBER= 50.0
122+
%
123+
% Adaptive CFL number (NO, YES)
124+
CFL_ADAPT= NO
125+
%
126+
% Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
127+
% CFL max value )
128+
CFL_ADAPT_PARAM= ( 1.5, 0.5, 10.0, 10000.0 )
129+
%
130+
% Runge-Kutta alpha coefficients
131+
RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 )
132+
133+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
134+
%
135+
% Linear solver or smoother for implicit formulations (BCGSTAB, FGMRES, SMOOTHER_JACOBI,
136+
% SMOOTHER_ILU, SMOOTHER_LUSGS,
137+
% SMOOTHER_LINELET)
138+
LINEAR_SOLVER= FGMRES
139+
%
140+
% Preconditioner of the Krylov linear solver (ILU, LU_SGS, LINELET, JACOBI)
141+
LINEAR_SOLVER_PREC= ILU
142+
%
143+
% Linael solver ILU preconditioner fill-in level (0 by default)
144+
LINEAR_SOLVER_ILU_FILL_IN= 0
145+
%
146+
% Minimum error of the linear solver for implicit formulations
147+
LINEAR_SOLVER_ERROR= 1E-15
148+
%
149+
% Max number of iterations of the linear solver for the implicit formulation
150+
LINEAR_SOLVER_ITER= 10
151+
152+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
153+
%
154+
% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
155+
% TURKEL_PREC, MSW)
156+
CONV_NUM_METHOD_FLOW= FDS
157+
%
158+
% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
159+
% Required for 2nd order upwind schemes (NO, YES)
160+
MUSCL_FLOW= YES
161+
%
162+
% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG,
163+
% BARTH_JESPERSEN, VAN_ALBADA_EDGE)
164+
SLOPE_LIMITER_FLOW= NONE
165+
%
166+
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
167+
TIME_DISCRE_FLOW= EULER_IMPLICIT
168+
169+
% ----------- SLOPE LIMITER AND DISSIPATION SENSOR DEFINITION -----------------%
170+
%
171+
% Coefficient for the Venkat's limiter (upwind scheme). A larger values decrease
172+
% the extent of limiting, values approaching zero cause
173+
% lower-order approximation to the solution (0.05 by default)
174+
VENKAT_LIMITER_COEFF= 0.05
175+
%
176+
% Coefficient for the adjoint sharp edges limiter (3.0 by default).
177+
ADJ_SHARP_LIMITER_COEFF= 3.0
178+
%
179+
% Freeze the value of the limiter after a number of iterations
180+
LIMITER_ITER= 999999
181+
%
182+
% 1st order artificial dissipation coefficients for
183+
% the Lax–Friedrichs method ( 0.15 by default )
184+
LAX_SENSOR_COEFF= 0.15
185+
%
186+
% 2nd and 4th order artificial dissipation coefficients for
187+
% the JST method ( 0.5, 0.02 by default )
188+
JST_SENSOR_COEFF= ( 0.5, 0.05 )
189+
%
190+
% 1st order artificial dissipation coefficients for
191+
% the adjoint Lax–Friedrichs method ( 0.15 by default )
192+
ADJ_LAX_SENSOR_COEFF= 0.15
193+
%
194+
% 2nd, and 4th order artificial dissipation coefficients for
195+
% the adjoint JST method ( 0.5, 0.02 by default )
196+
ADJ_JST_SENSOR_COEFF= ( 0.5, 0.02 )
197+
198+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
199+
%
200+
% Restart flow input file
201+
SOLUTION_FILENAME= solution_flow.dat
202+
%
203+
% Restart adjoint input file
204+
SOLUTION_ADJ_FILENAME= solution_adj.dat
205+
%
206+
% Output file convergence history (w/o extension)
207+
CONV_FILENAME= history
208+
%
209+
% Output file with the forces breakdown
210+
BREAKDOWN_FILENAME= forces_breakdown.dat
211+
%
212+
% Output file restart flow
213+
RESTART_FILENAME= restart_flow.dat
214+
%
215+
% Output file restart adjoint
216+
RESTART_ADJ_FILENAME= restart_adj.dat
217+
%
218+
% Output file flow (w/o extension) variables
219+
VOLUME_FILENAME= flow
220+
%
221+
% Output file adjoint (w/o extension) variables
222+
VOLUME_ADJ_FILENAME= adjoint
223+
%
224+
% Output Objective function
225+
VALUE_OBJFUNC_FILENAME= of_eval.dat
226+
%
227+
% Output objective function gradient (using continuous adjoint)
228+
GRAD_OBJFUNC_FILENAME= of_grad.dat
229+
%
230+
% Output file surface flow coefficient (w/o extension)
231+
SURFACE_FILENAME= surface_flow
232+
%
233+
% Output file surface adjoint coefficient (w/o extension)
234+
SURFACE_ADJ_FILENAME= surface_adjoint

0 commit comments

Comments
 (0)