Skip to content

Commit 4e23045

Browse files
TobiKattmannTobiKattmann
committed
Add unsteady CHT tutorial files.
1 parent 3d73157 commit 4e23045

File tree

6 files changed

+98642
-0
lines changed

6 files changed

+98642
-0
lines changed

‎multiphysics/unsteady_cht/cht_2d_3cylinders.cfg‎

Lines changed: 61 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,61 @@
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+
%
22+
CONFIG_LIST = (flow_cylinder.cfg, solid_cylinder1.cfg, solid_cylinder2.cfg, solid_cylinder3.cfg)
23+
%
24+
%
25+
MARKER_ZONE_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
26+
%
27+
%
28+
MARKER_CHT_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
29+
%
30+
%
31+
TIME_DOMAIN = YES
32+
%
33+
%
34+
TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
35+
%
36+
%
37+
TIME_STEP= 0.05
38+
%
39+
%
40+
MAX_TIME= 100.0
41+
%
42+
%
43+
TIME_ITER= 2000
44+
%
45+
% Number of total iterations
46+
OUTER_ITER = 100
47+
%
48+
% Courant-Friedrichs-Lewy condition of the finest grid
49+
CFL_NUMBER= 100.0
50+
%
51+
% Mesh input file
52+
MESH_FILENAME= mesh_cht_3cyl.su2
53+
%
54+
% Mesh input file format (SU2, CGNS, NETCDF_ASCII)
55+
MESH_FORMAT= SU2
56+
%
57+
%
58+
HISTORY_OUTPUT= (ITER, BGS_RES[0], BGS_RES[1], BGS_RES[2], BGS_RES[3], HEAT[0], AERO_COEFF[0])
59+
%
60+
%
61+
OUTPUT_FILES=(RESTART, PARAVIEW_MULTIBLOCK)

‎multiphysics/unsteady_cht/flow_cylinder.cfg‎

Lines changed: 245 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,245 @@
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+
% Restart solution (NO, YES)
19+
RESTART_SOL= NO
20+
%
21+
% Objective function in gradient evaluation (DRAG, LIFT, SIDEFORCE, MOMENT_X,
22+
% MOMENT_Y, MOMENT_Z, EFFICIENCY,
23+
% EQUIVALENT_AREA, NEARFIELD_PRESSURE,
24+
% FORCE_X, FORCE_Y, FORCE_Z, THRUST,
25+
% TORQUE, TOTAL_HEATFLUX,
26+
% MAXIMUM_HEATFLUX, INVERSE_DESIGN_PRESSURE,
27+
% INVERSE_DESIGN_HEATFLUX, SURFACE_TOTAL_PRESSURE,
28+
% SURFACE_MASSFLOW, SURFACE_STATIC_PRESSURE, SURFACE_MACH)
29+
% For a weighted sum of objectives: separate by commas, add OBJECTIVE_WEIGHT and MARKER_MONITORING in matching order.
30+
OBJECTIVE_FUNCTION= TOTAL_HEATFLUX
31+
%
32+
% List of weighting values when using more than one OBJECTIVE_FUNCTION. Separate by commas and match with MARKER_MONITORING.
33+
OBJECTIVE_WEIGHT = 1.0
34+
%
35+
% Read binary restart files (YES, NO)
36+
READ_BINARY_RESTART = YES
37+
%
38+
% Data written to history file
39+
HISTORY_OUTPUT=(ITER, RMS_RES, HEAT)
40+
41+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
42+
%
43+
% Farfield boundary marker(s) (NONE = no marker)
44+
MARKER_FAR= ( farfield )
45+
%
46+
% Marker(s) of the surface to be plotted or designed
47+
MARKER_PLOTTING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
48+
%
49+
% Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated
50+
MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
51+
52+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
53+
%
54+
% Density model within the incompressible flow solver.
55+
% Options are CONSTANT (default), BOUSSINESQ, or VARIABLE. If VARIABLE,
56+
% an appropriate fluid model must be selected.
57+
INC_DENSITY_MODEL= VARIABLE
58+
%
59+
% Solve the energy equation in the incompressible flow solver
60+
INC_ENERGY_EQUATION = YES
61+
%
62+
% Initial density for incompressible flows (1.2886 kg/m^3 by default)
63+
INC_DENSITY_INIT= 0.0210322
64+
%
65+
% Initial velocity for incompressible flows (1.0,0,0 m/s by default)
66+
INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
67+
%
68+
% Initial temperature for incompressible flows that include the
69+
% energy equation (288.15 K by default). Value is ignored if
70+
% INC_ENERGY_EQUATION is false.
71+
INC_TEMPERATURE_INIT= 288.15
72+
%
73+
% Non-dimensionalization scheme for incompressible flows. Options are
74+
% INITIAL_VALUES (default), REFERENCE_VALUES, or DIMENSIONAL.
75+
% INC_*_REF values are ignored unless REFERENCE_VALUES is chosen.
76+
INC_NONDIM= DIMENSIONAL
77+
78+
% ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
79+
%
80+
% Fluid model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS,
81+
% CONSTANT_DENSITY, INC_IDEAL_GAS)
82+
FLUID_MODEL= INC_IDEAL_GAS
83+
%
84+
% Specific heat at constant pressure, Cp (1004.703 J/kg*K (air)).
85+
% Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
86+
SPECIFIC_HEAT_CP= 1004.703
87+
%
88+
% Molecular weight for an incompressible ideal gas (28.96 g/mol (air) default)
89+
% Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
90+
MOLECULAR_WEIGHT= 28.96
91+
92+
% --------------------------- VISCOSITY MODEL ---------------------------------%
93+
%
94+
% Viscosity model (SUTHERLAND, CONSTANT_VISCOSITY).
95+
VISCOSITY_MODEL= CONSTANT_VISCOSITY
96+
%
97+
% Molecular Viscosity that would be constant (1.716E-5 by default)
98+
MU_CONSTANT= 1.7893e-05
99+
%
100+
% Sutherland Viscosity Ref (1.716E-5 default value for AIR SI)
101+
MU_REF= 1.716E-5
102+
%
103+
% Sutherland Temperature Ref (273.15 K default value for AIR SI)
104+
MU_T_REF= 273.15
105+
%
106+
% Sutherland constant (110.4 default value for AIR SI)
107+
SUTHERLAND_CONSTANT= 110.4
108+
109+
% --------------------------- THERMAL CONDUCTIVITY MODEL ----------------------%
110+
%
111+
% Conductivity model (CONSTANT_CONDUCTIVITY, CONSTANT_PRANDTL).
112+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
113+
%
114+
% Molecular Thermal Conductivity that would be constant (0.0257 by default)
115+
KT_CONSTANT= 0.0257
116+
%
117+
% Laminar Prandtl number (0.72 (air), only for CONSTANT_PRANDTL)
118+
PRANDTL_LAM= 0.72
119+
%
120+
% Turbulent Prandtl number (0.9 (air), only for CONSTANT_PRANDTL)
121+
PRANDTL_TURB= 0.90
122+
123+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
124+
%
125+
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
126+
NUM_METHOD_GRAD= GREEN_GAUSS
127+
128+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
129+
%
130+
% Linear solver or smoother for implicit formulations (BCGSTAB, FGMRES, SMOOTHER_JACOBI,
131+
% SMOOTHER_ILU, SMOOTHER_LUSGS,
132+
% SMOOTHER_LINELET)
133+
LINEAR_SOLVER= FGMRES
134+
%
135+
% Preconditioner of the Krylov linear solver (ILU, LU_SGS, LINELET, JACOBI)
136+
LINEAR_SOLVER_PREC= ILU
137+
%
138+
% Linael solver ILU preconditioner fill-in level (0 by default)
139+
LINEAR_SOLVER_ILU_FILL_IN= 0
140+
%
141+
% Minimum error of the linear solver for implicit formulations
142+
LINEAR_SOLVER_ERROR= 1E-15
143+
%
144+
% Max number of iterations of the linear solver for the implicit formulation
145+
LINEAR_SOLVER_ITER= 5
146+
147+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
148+
%
149+
% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
150+
% TURKEL_PREC, MSW)
151+
CONV_NUM_METHOD_FLOW= FDS
152+
%
153+
% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
154+
% Required for 2nd order upwind schemes (NO, YES)
155+
MUSCL_FLOW= YES
156+
%
157+
% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG,
158+
% BARTH_JESPERSEN, VAN_ALBADA_EDGE)
159+
SLOPE_LIMITER_FLOW= NONE
160+
%
161+
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
162+
TIME_DISCRE_FLOW= EULER_IMPLICIT
163+
164+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
165+
%
166+
% Min value of the residual (log10 of the residual)
167+
CONV_RESIDUAL_MINVAL= -19
168+
%
169+
% Start convergence criteria at iteration number
170+
CONV_STARTITER= 10
171+
%
172+
% Number of elements to apply the criteria
173+
CONV_CAUCHY_ELEMS= 100
174+
%
175+
% Epsilon to control the series convergence
176+
CONV_CAUCHY_EPS= 1E-6
177+
178+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
179+
%
180+
% Restart flow input file
181+
SOLUTION_FILENAME= solution_flow.dat
182+
%
183+
% Restart adjoint input file
184+
SOLUTION_ADJ_FILENAME= solution_adj.dat
185+
%
186+
% Output file format (TECPLOT, TECPLOT_BINARY, PARAVIEW,
187+
% FIELDVIEW, FIELDVIEW_BINARY)
188+
TABULAR_FORMAT= TECPLOT
189+
%
190+
% Output file convergence history (w/o extension)
191+
CONV_FILENAME= history
192+
%
193+
% Output file with the forces breakdown
194+
BREAKDOWN_FILENAME= forces_breakdown.dat
195+
%
196+
% Output file restart flow
197+
RESTART_FILENAME= restart_flow.dat
198+
%
199+
% Output file restart adjoint
200+
RESTART_ADJ_FILENAME= restart_adj.dat
201+
%
202+
% Output file flow (w/o extension) variables
203+
VOLUME_FILENAME= flow
204+
%
205+
% Output file adjoint (w/o extension) variables
206+
VOLUME_ADJ_FILENAME= adjoint
207+
%
208+
% Output Objective function
209+
VALUE_OBJFUNC_FILENAME= of_eval.dat
210+
%
211+
% Output objective function gradient (using continuous adjoint)
212+
GRAD_OBJFUNC_FILENAME= of_grad.dat
213+
%
214+
% Output file surface flow coefficient (w/o extension)
215+
SURFACE_FILENAME= surface_flow
216+
%
217+
% Output file surface adjoint coefficient (w/o extension)
218+
SURFACE_ADJ_FILENAME= surface_adjoint
219+
%
220+
% Writing solution file frequency
221+
WRT_SOL_FREQ= 250
222+
%
223+
% Writing convergence history frequency
224+
WRT_CON_FREQ= 1
225+
226+
% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
227+
%
228+
% Linear solver or smoother for implicit formulations (FGMRES, RESTARTED_FGMRES, BCGSTAB)
229+
DEFORM_LINEAR_SOLVER= FGMRES
230+
%
231+
% Number of smoothing iterations for mesh deformation
232+
DEFORM_LINEAR_SOLVER_ITER= 200
233+
%
234+
% Number of nonlinear deformation iterations (surface deformation increments)
235+
DEFORM_NONLINEAR_ITER= 1
236+
%
237+
% Print the residuals during mesh deformation to the console (YES, NO)
238+
DEFORM_CONSOLE_OUTPUT= YES
239+
%
240+
% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME,
241+
% WALL_DISTANCE, CONSTANT_STIFFNESS)
242+
DEFORM_STIFFNESS_TYPE= INVERSE_VOLUME
243+
%
244+
% Visualize the deformation (NO, YES)
245+
VISUALIZE_VOLUME_DEF= YES

0 commit comments

Comments
 (0)