Skip to content

Commit f0d4fb9

Browse files
TobiKattmannTobiKattmann
committed
Add new solid heat conduction Testcase with periodic BC's.
1 parent a0d296c commit f0d4fb9

3 files changed

Lines changed: 108 additions & 0 deletions

File tree

‎TestCases/parallel_regression.py‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1258,6 +1258,22 @@ def main():
12581258
p1rad.tol = 0.00001
12591259
test_list.append(p1rad)
12601260

1261+
1262+
# #############################
1263+
# ### Solid Heat Conduction ###
1264+
# #############################
1265+
1266+
# 2D pins, periodically connected
1267+
solid_periodic_pins = TestCase('solid_periodic_pins')
1268+
solid_periodic_pins.cfg_dir = "solid_heat_conduction/periodic_pins"
1269+
solid_periodic_pins.cfg_file = "configSolid.cfg"
1270+
solid_periodic_pins.test_iter = 750
1271+
solid_periodic_pins.test_vals = [-15.739745, -14.448665, 300.900000, 425.320000, 0.000000, 5.000000, -1.448445] #last 7 lines
1272+
solid_periodic_pins.su2_exec = "mpirun -n 2 SU2_CFD"
1273+
solid_periodic_pins.timeout = 1600
1274+
solid_periodic_pins.tol = 0.00001
1275+
test_list.append(solid_periodic_pins)
1276+
12611277
# ###############################
12621278
# ### Conjugate heat transfer ###
12631279
# ###############################

‎TestCases/solid_heat_conduction/periodic_pins/README.md‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,10 @@
1+
# Periodic Pins
2+
3+
This a solid heat conduction testcase in which the `HEAT_EQUATION` solver runs standalone (i.e. not as CHT).
4+
The simulation domain is the solid domain of the `incomp_navierstokes/streamwise_periodic/chtPinArray_2d`-Testcase.
5+
Therefore the provided gmsh `.geo` file contains the full CHT mesh but only writes out the solid zone when called.
6+
7+
Note that using periodic boundary conditions for the solid zone made the solution take ~10x more iterations to converge , compared to the same setup using adiabatic walls.
8+
This was found for solid only as well as CHT cases.
9+
10+
Expected results are perfectly matched Temperatures at the periodic interface. Compare e.g. using Paraview's `Transform`-Filter with domain length 0.0111544m.

‎TestCases/solid_heat_conduction/periodic_pins/configSolid.cfg‎

Lines changed: 82 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,82 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Solid-only Heated Pins, periodically connected %
5+
% Author: T. Kattmann %
6+
% Institution: Bosch Thermotechniek B.V. %
7+
% Date: 2021.09.27 %
8+
% File Version 7.2.0 "Blackbird" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
%
12+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
13+
%
14+
SOLVER= HEAT_EQUATION
15+
%
16+
OBJECTIVE_FUNCTION= AVG_TEMPERATURE
17+
OBJECTIVE_WEIGHT= 1.0
18+
%
19+
OPT_OBJECTIVE= AVG_TOTALTEMP
20+
%
21+
% ---------------- (SOLIDS) CONDUCTION CONDITION DEFINITION -------------------%
22+
%
23+
INC_NONDIM= DIMENSIONAL
24+
FREESTREAM_TEMPERATURE= 345.0
25+
MATERIAL_DENSITY= 2719
26+
SPECIFIC_HEAT_CP = 871.0
27+
THERMAL_CONDUCTIVITY_CONSTANT= 200
28+
%
29+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
30+
%
31+
MARKER_HEATFLUX= (solid_pin1_inner, 5e5, \
32+
solid_pin2_inner, 5e5, \
33+
solid_pin3_inner, 0.0, \
34+
solid_pin1_walls, 0.0, \
35+
solid_pin2_walls, 0.0, \
36+
solid_pin3_walls, 0.0)
37+
%
38+
MARKER_PERIODIC= ( solid_pin1_periodic, solid_pin3_periodic, 0.0,0.0,0.0, 0.0,0.0,0.0, 0.0111544,0.0,0.0 )
39+
%
40+
MARKER_ISOTHERMAL= ( solid_pin1_interface, 300, solid_pin2_interface, 300, solid_pin3_interface, 300 )
41+
%
42+
% ------------------------ SURFACES IDENTIFICATION ----------------------------%
43+
%
44+
MARKER_MONITORING = ( solid_pin1_inner, solid_pin3_inner )
45+
%
46+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
47+
%
48+
%NUM_METHOD_GRAD= GREEN_GAUSS
49+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
50+
CFL_NUMBER= 1e8
51+
%
52+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
53+
%
54+
LINEAR_SOLVER= FGMRES
55+
LINEAR_SOLVER_PREC= ILU
56+
LINEAR_SOLVER_ERROR= 1E-15
57+
LINEAR_SOLVER_ITER= 5
58+
%
59+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
60+
%
61+
ITER= 1000
62+
%
63+
CONV_FIELD= RMS_TEMPERATURE
64+
CONV_RESIDUAL_MINVAL= -16
65+
CONV_STARTITER= 10
66+
%
67+
% -------------------- HEAT NUMERICAL METHOD DEFINITION -----------------------%
68+
%
69+
TIME_DISCRE_HEAT= EULER_IMPLICIT
70+
%
71+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
72+
%
73+
MESH_FILENAME= solid.su2
74+
%
75+
SCREEN_OUTPUT= INNER_ITER, RMS_TEMPERATURE, MAX_TEMPERATURE, AVG_TEMPERATURE, TOTAL_HEATFLUX, MAXIMUM_HEATFLUX, LINSOL_ITER, LINSOL_RESIDUAL
76+
SCREEN_WRT_FREQ_INNER= 50
77+
%
78+
HISTORY_OUTPUT= (ITER, RMS_RES, HEAT, LINSOL)
79+
%
80+
OUTPUT_FILES= RESTART, PARAVIEW_MULTIBLOCK
81+
VOLUME_OUTPUT= RESIDUAL
82+
OUTPUT_WRT_FREQ= 1000

0 commit comments

Comments
 (0)