Skip to content

Commit 56aedab

Browse files
TobiKattmannTobiKattmann
authored
Merge pull request #12 from su2code/unsteady_cht_tutorial
Add unsteady CHT tutorial files.
2 parents 3d73157 + 507fc80 commit 56aedab

File tree

6 files changed

+98270
-0
lines changed

6 files changed

+98270
-0
lines changed

‎multiphysics/unsteady_cht/cht_2d_3cylinders.cfg‎

Lines changed: 50 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,50 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: 2D transient cylinder array with CHT couplings %
5+
% Author: O. Burghardt, T. Economon %
6+
% Institution: Chair for Scientific Computing, TU Kaiserslautern %
7+
% Date: March 18, 2021 %
8+
% File Version 7.1.1 "Blackbird" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
%
12+
SOLVER= MULTIPHYSICS
13+
%
14+
CONFIG_LIST = (flow_cylinder.cfg, solid_cylinder1.cfg, solid_cylinder2.cfg, solid_cylinder3.cfg)
15+
%
16+
MARKER_ZONE_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
17+
MARKER_CHT_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
18+
%
19+
TIME_DOMAIN = YES
20+
TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
21+
TIME_STEP= 0.05
22+
MAX_TIME= 100.0
23+
TIME_ITER= 2000
24+
%
25+
% Number of total outer iterations
26+
OUTER_ITER = 100
27+
%
28+
% Courant-Friedrichs-Lewy condition
29+
CFL_NUMBER= 100.0
30+
%
31+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
32+
%
33+
MESH_FILENAME= mesh_cht_3cyl.su2
34+
%
35+
SCREEN_OUTPUT= ( TIME_ITER, OUTER_ITER, \
36+
BGS_PRESSURE[0], BGS_TEMPERATURE[0], BGS_TEMPERATURE[1], BGS_TEMPERATURE[2], BGS_TEMPERATURE[3], \
37+
DRAG[0], SURFACE_STATIC_TEMPERATURE[0], AVG_TEMPERATURE[1], TOTAL_HEATFLUX[1] )
38+
SCREEN_WRT_FREQ_OUTER= 25
39+
%
40+
HISTORY_OUTPUT= (ITER, BGS_RES[0], BGS_RES[1], BGS_RES[2], BGS_RES[3], HEAT[0], AERO_COEFF[0])
41+
%
42+
OUTPUT_FILES=(RESTART, PARAVIEW_MULTIBLOCK)
43+
%
44+
RESTART_SOL= NO
45+
READ_BINARY_RESTART= YES
46+
SOLUTION_FILENAME= solution
47+
RESTART_FILENAME= solution
48+
%
49+
VOLUME_FILENAME= flow
50+

‎multiphysics/unsteady_cht/flow_cylinder.cfg‎

Lines changed: 156 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,156 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Transient incomp 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+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
19+
%
20+
% Farfield boundary marker(s) (NONE = no marker)
21+
MARKER_FAR= ( farfield )
22+
%
23+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
24+
%
25+
% Density model within the incompressible flow solver.
26+
% Options are CONSTANT (default), BOUSSINESQ, or VARIABLE. If VARIABLE,
27+
% an appropriate fluid model must be selected.
28+
INC_DENSITY_MODEL= VARIABLE
29+
%
30+
% Solve the energy equation in the incompressible flow solver
31+
INC_ENERGY_EQUATION = YES
32+
%
33+
% Initial density for incompressible flows (1.2886 kg/m^3 by default)
34+
INC_DENSITY_INIT= 0.0210322
35+
%
36+
% Initial velocity for incompressible flows (1.0,0,0 m/s by default)
37+
INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
38+
%
39+
% Initial temperature for incompressible flows that include the
40+
% energy equation (288.15 K by default). Value is ignored if
41+
% INC_ENERGY_EQUATION is false.
42+
INC_TEMPERATURE_INIT= 288.15
43+
%
44+
% Non-dimensionalization scheme for incompressible flows. Options are
45+
% INITIAL_VALUES (default), REFERENCE_VALUES, or DIMENSIONAL.
46+
% INC_*_REF values are ignored unless REFERENCE_VALUES is chosen.
47+
INC_NONDIM= DIMENSIONAL
48+
%
49+
% ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
50+
%
51+
% Fluid model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS,
52+
% CONSTANT_DENSITY, INC_IDEAL_GAS)
53+
FLUID_MODEL= INC_IDEAL_GAS
54+
%
55+
% Specific heat at constant pressure, Cp (1004.703 J/kg*K (air)).
56+
% Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
57+
SPECIFIC_HEAT_CP= 1004.703
58+
%
59+
% Molecular weight for an incompressible ideal gas (28.96 g/mol (air) default)
60+
% Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
61+
MOLECULAR_WEIGHT= 28.96
62+
%
63+
% --------------------------- VISCOSITY MODEL ---------------------------------%
64+
%
65+
% Viscosity model (SUTHERLAND, CONSTANT_VISCOSITY).
66+
VISCOSITY_MODEL= CONSTANT_VISCOSITY
67+
%
68+
% Molecular Viscosity that would be constant (1.716E-5 by default)
69+
MU_CONSTANT= 1.7893e-05
70+
%
71+
% Sutherland Viscosity Ref (1.716E-5 default value for AIR SI)
72+
MU_REF= 1.716E-5
73+
%
74+
% Sutherland Temperature Ref (273.15 K default value for AIR SI)
75+
MU_T_REF= 273.15
76+
%
77+
% Sutherland constant (110.4 default value for AIR SI)
78+
SUTHERLAND_CONSTANT= 110.4
79+
80+
% --------------------------- THERMAL CONDUCTIVITY MODEL ----------------------%
81+
%
82+
% Conductivity model (CONSTANT_CONDUCTIVITY, CONSTANT_PRANDTL).
83+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
84+
%
85+
% Molecular Thermal Conductivity that would be constant (0.0257 by default)
86+
KT_CONSTANT= 0.0257
87+
%
88+
% Laminar Prandtl number (0.72 (air), only for CONSTANT_PRANDTL)
89+
PRANDTL_LAM= 0.72
90+
%
91+
% Turbulent Prandtl number (0.9 (air), only for CONSTANT_PRANDTL)
92+
PRANDTL_TURB= 0.90
93+
%
94+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
95+
%
96+
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
97+
NUM_METHOD_GRAD= GREEN_GAUSS
98+
%
99+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
100+
%
101+
% Linear solver or smoother for implicit formulations (BCGSTAB, FGMRES, SMOOTHER_JACOBI,
102+
% SMOOTHER_ILU, SMOOTHER_LUSGS,
103+
% SMOOTHER_LINELET)
104+
LINEAR_SOLVER= FGMRES
105+
%
106+
% Preconditioner of the Krylov linear solver (ILU, LU_SGS, LINELET, JACOBI)
107+
LINEAR_SOLVER_PREC= ILU
108+
%
109+
% Linael solver ILU preconditioner fill-in level (0 by default)
110+
LINEAR_SOLVER_ILU_FILL_IN= 0
111+
%
112+
% Minimum error of the linear solver for implicit formulations
113+
LINEAR_SOLVER_ERROR= 1E-15
114+
%
115+
% Max number of iterations of the linear solver for the implicit formulation
116+
LINEAR_SOLVER_ITER= 20
117+
%
118+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
119+
%
120+
% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
121+
% TURKEL_PREC, MSW)
122+
CONV_NUM_METHOD_FLOW= FDS
123+
%
124+
% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
125+
% Required for 2nd order upwind schemes (NO, YES)
126+
MUSCL_FLOW= YES
127+
%
128+
% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG,
129+
% BARTH_JESPERSEN, VAN_ALBADA_EDGE)
130+
SLOPE_LIMITER_FLOW= NONE
131+
%
132+
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
133+
TIME_DISCRE_FLOW= EULER_IMPLICIT
134+
%
135+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
136+
%
137+
% Min value of the residual (log10 of the residual)
138+
CONV_RESIDUAL_MINVAL= -19
139+
%
140+
% Start convergence criteria at iteration number
141+
CONV_STARTITER= 10
142+
%
143+
% Number of elements to apply the criteria
144+
CONV_CAUCHY_ELEMS= 100
145+
%
146+
% Epsilon to control the series convergence
147+
CONV_CAUCHY_EPS= 1E-6
148+
%
149+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
150+
%
151+
MARKER_PLOTTING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
152+
%
153+
MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
154+
%
155+
MARKER_ANALYZE= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
156+

0 commit comments

Comments
 (0)