FEATool Multiphysics
v1.14
Finite Element Analysis Toolbox

FEATool supports modeling heat transfer through both conduction, that is heat transported by a diffusion process, and also convection, which is heat transported through a fluid through convection by a velocity field. The heat transfer physics mode supports both these processes, and defines the following equation
\[ \rho C_p\frac{\partial T}{\partial t} + \nabla\cdot(k\nabla c) = Q  \rho C_p\mathbf{u}\cdot\nabla T \]
where \(\rho\) is the density, \(C_p\) the heat capacity, \(k\) is the thermal conductivity, \(Q\) heat source term, and \(\mathbf{u}\) a vector valued convective velocity field.
This example models heat conduction in the form of transient cooling for shrink fitting of a two part assembly. A tungsten rod heated to 84 °C is inserted into a 10 °C chilled steel frame part. The time when the maximum temperature has cooled to 70 °C should be determined. The assembly is cooled due to convection through a surrounding medium kept at T_{inf} = 17 °C and a heat transfer coefficient of h = 50 W/(m^{2}K). The surrounding cooling medium is not modeled directly, thus the convective term is omitted, but the effects are incorporated into the model by the use of natural convection boundary conditions [1].
This section describes how to set up and solve the thermal shrink fitting example with the FEATool graphical user interface (GUI).
The model is available as an automated tutorial by selecting Model Examples and Tutorials... > Heat Transfer > Shrink Fitting of an Assembly from the File menu. Or alternatively, follow the stepbystep instructions below. Note that the CFDTool interface differ slightly from the FEATool Multiphysics instructions described in the following.
Select the Heat Transfer physics mode from the Select Physics dropdown menu. (Note that for CFDTool the physics selection is done in the Equation settings dialog box.)
0
into the x_{min} edit field.0.11
into the x_{max} edit field.0
into the y_{min} edit field.Enter 0.12
into the y_{max} edit field.
0.065 0
into the center edit field.0.015
into the x_{radius} edit field.Enter 0.015
into the y_{radius} edit field.
0.11 0.12
into the center edit field.0.035
into the x_{radius} edit field.Enter 0.035
into the y_{radius} edit field.
0 0.06
into the center edit field.0.025
into the x_{radius} edit field.Enter 0.025
into the y_{radius} edit field.
Press OK to finish and close the dialog box.
Enter R1  E1  E2  E3
into the Geometry Formula edit field.
Press OK to finish and close the dialog box.
0.065
into the x_{min} edit field.0.16
into the x_{max} edit field.0.05
into the y_{min} edit field.Enter 0.07
into the y_{max} edit field.
0.065 0.06
into the center edit field.0.01
into the x_{radius} edit field.Enter 0.01
into the y_{radius} edit field.
Select R2 and E4 in the geometry object Selection list box.
Select the Make copy of geometry object check box.
Select CS1 and CJ1 in the geometry object Selection list box.
Press the  / Subtract geometry objects Toolbar button.
The default grid may be too coarse ensure an accurate solution. Decreasing the grid size and generating a finer grid can resolve curved boundaries better.
Enter 0.0025
into the Grid Size edit field and press the Generate button to call the grid generation algorithm.
Equation and material coefficients are specified in Equation/Subdomain mode. In the Equation Settings dialog box enter the coefficients for density, heat capacity, thermal conductivity, and initial temperature for each material. (Note that the Equation Settings dialog box may look different for CFDTool.)
rho_tungsten
into the Density edit field.cp_tungsten
into the Heat capacity edit field.k_tungsten
into the Thermal conductivity edit field.Enter 84
into the Initial condition for T edit field.
rho_steel
into the Density edit field.cp_steel
into the Heat capacity edit field.k_steel
into the Thermal conductivity edit field.Enter 10
into the Initial condition for T edit field.
The Model Constants and Expressions functionality can be used to define and store convenient expressions which then are available in the point, equation, boundary coefficients, and as postprocessing expressions. Here it is used to define the material parameters.
Name  Expression 

rho_tungsten  19000 
cp_tungsten  134 
k_tungsten  163 
rho_steel  7500 
cp_steel  470 
k_steel  44 
h_coef  50 
In the Boundary Settings dialog box, select the Heat flux boundary condition for all the boundaries. Enter h_coef
in the edit field for the convective transfer coefficient h, and also enter 17
for the surrounding reference temperature T_{inf}.
Now that the problem is specified, press the Solve mode button to switch to solve mode. Since this is a time dependent study, open the solver settings and select the TimeDependent solver. Set the Time step to 0.25
, Simulation time to 16
, Time stopping criteria to 0
, then press Solve to start the solution process.
Enter 20
into the Number or specified vector of contour levels to plot edit field.
Press Apply to plot and visualize the selected postprocessing options.
Look back through the solutions and verify that the assembly has cooled to a temperature of 70 degrees around t = 13 s. Note that both the colorbar and Limits field will show the minimum and maximum surface plot value.
The shrink fitting of an assembly heat transfer model has now been completed and can be saved as a binary (.fea) model file, or exported as a programmable MATLAB mscript text file (available as the example ex_heattransfer5 script file), or GUI script (.fes) file.
[1] Krysl P. A Pragmatic Introduction to the Finite Element Method for Thermal and Stress Analysis. Pressure Cooker Press, USA, 2005.