Solved with COMSOL Multiphysics 5.
Axisymmetric Transient Heat Transfer
This example shows an axisymmetric transient thermal analysis with a step change to
1000 C at time 0. The example is taken from a NAFEMS benchmark collection (Ref.
1).
Model Definition
This model domain is 0.3 m-by- 0.4 m. For the boundary conditions, assume the
following:
The left boundary is the symmetry axis.
The other boundaries have a temperature of 1000 C. The entire domain is at 0 C
at the start, which represents a step change in temperature at the boundaries.
In the domain use the following material properties:
The density, , is 7850 kg/m3
The heat capacity is 460 J/(kgC)
The thermal conductivity is 52 W/(mC)
1 |
A X I S Y M M E T R I C TR A N S I E N T H E A T TR A N S F E R
Solved with COMSOL Multiphysics 5.0
Results
The following revolved surface plot shows the temperature distribution inside the
cylinder after 190 seconds:
Figure 1: Temperature distribution after 190 seconds.
The benchmark result for the target location (r = 0.1 m and z = 0.3 m) is a temperature
of 186.5 C. The COMSOL Multiphysics model, using a default mesh with about 720
elements, gives a temperature of roughly 186.4 C.
Reference
1. A.D. Cameron, J.A. Casey, and G.B. Simpson, NAFEMS Benchmark Tests for
Thermal Analysis (Summary), NAFEMS, Glasgow, 1986.
Model Library path: COMSOL_Multiphysics/Heat_Transfer/
heat_transient_axi
2 |
A X I S Y M M E T R I C TR A N S I E N T H E A T TR A N S F E R
Solved with COMSOL Multiphysics 5.0
Modeling Instructions
From the File menu, choose New.
NEW
1 In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D Axisymmetric.
2 In the Select physics tree, select Heat Transfer>Heat Transfer in Solids (ht).
3 Click Add.
4 Click Study.
5 In the Select study tree, select Preset Studies>Time Dependent.
6 Click Done.
GEOMETRY 1
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size section.
3 In the Width text field, type 0.3.
4 In the Height text field, type 0.4.
5 Click the Build All Objects button.
H E A T TR A N S F E R I N S O L I D S ( H T )
Temperature 1
1 On the Physics toolbar, click Boundaries and choose Temperature.
2 In the Settings window for Temperature, locate the Boundary Selection section.
3 From the Selection list, choose All boundaries.
4 Locate the Temperature section. In the T0 text field, type 1000[degC].
Heat Transfer in Solids 1
1 In the Model Builder window, under Component 1 (comp1)>Heat Transfer in Solids (ht)
click Heat Transfer in Solids 1.
2 In the Settings window for Heat Transfer in Solids, locate the Heat Conduction, Solid
section.
3 |
A X I S Y M M E T R I C TR A N S I E N T H E A T TR A N S F E R
Solved with COMSOL Multiphysics 5.0
3 From the k list, choose User defined. In the associated text field, type 52.
4 Locate the Thermodynamics, Solid section. From the Cp list, choose User defined. In
the associated text field, type 460.
5 From the list, choose User defined. In the associated text field, type 7850.
Initial Values 1
1 In the Model Builder window, under Component 1 (comp1)>Heat Transfer in Solids (ht)
click Initial Values 1.
2 In the Settings window for Initial Values, locate the Initial Values section.
3 In the T text field, type 0[degC].
MESH 1
In the Model Builder window, under Component 1 (comp1) right-click Mesh 1 and choose
Build All.
STUDY 1
Step 1: Time Dependent
1 In the Model Builder window, expand the Study 1 node, then click Step 1: Time
Dependent.
2 In the Settings window for Time Dependent, locate the Study Settings section.
3 In the Times text field, type range(0,10,190).
4 On the Model toolbar, click Compute.
RESULTS
Temperature, 3D (ht)
To get the plot shown in Figure 1, just change the unit as follows:
1 In the Model Builder window, expand the Temperature, 3D (ht) node, then click
Surface 1.
2 In the Settings window for Surface, locate the Expression section.
3 From the Unit list, choose degC.
4 On the 3D plot group toolbar, click Plot.
5 Click the Zoom Extents button on the Graphics toolbar.
4 |
A X I S Y M M E T R I C TR A N S I E N T H E A T TR A N S F E R
Solved with COMSOL Multiphysics 5.0
Isothermal Contours (ht)
The second default plot group visualizes the temperature field using a contour plot.
Data Sets
The benchmark value for the temperature at r = 0.1 and z = 0.3 is 186.5 Co. To
compare the value from the simulation, evaluate the temperature in that position.
1 On the Results toolbar, click Cut Point 2D.
2 In the Settings window for Cut Point 2D, locate the Point Data section.
3 In the r text field, type 0.1.
4 In the z text field, type 0.3.
Derived Values
1 On the Results toolbar, click Point Evaluation.
2 In the Settings window for Point Evaluation, locate the Data section.
3 From the Data set list, choose Cut Point 2D 1.
4 From the Time selection list, choose Last.
5 Locate the Expression section. From the Unit list, choose degC.
6 Click the Evaluate button.
5 |
A X I S Y M M E T R I C TR A N S I E N T H E A T TR A N S F E R
Solved with COMSOL Multiphysics 5.0
TABLE
Go to the Table window.
6 |
A X I S Y M M E T R I C TR A N S I E N T H E A T TR A N S F E R
