FEATool Multiphysics
v1.16.1
Finite Element Analysis Toolbox

Transient heat diffusion in a 1 m rod where one end is kept at T = 25 degrees and the other subject to constant heat flux q_{n} = 1. This onedimensional time dependent problem is transformed to a static two dimensional one. The computed results can be compared with the analytic solution and Transient Heat Diffusion in a Rod example
\[ T_{ref} = (24+x) + \sum_{n=1}^\infty 8/(12n)^2/\pi^2cos((n1/2)\pi x)e^{((n1/2)^2\pi^2)t} \]
This model is available as an automated tutorial by selecting Model Examples and Tutorials... > Heat Transfer > SpaceTime Transformation of Heat Conduction from the File menu. Or alternatively, follow the stepbystep instructions below.
In the New Model dialog box, select 2D for the Space Dimensions, change the name of the space dimension coordinate names to x to represent the length dimension and yt to represent the time axis (t is a reserved variable name for time/test function and should not be used), and select Heat Transfer from the Select Physics dropdown menu. Leave the dependent variable name to the default and press OK to finish the physics mode selection.
The original problem examines a 1 m rod during 0.2 s, so the equivalent 2D geometry is a rectangle spanning 0 x 1 and 0 x 0.2 where the ydirection will represent the time axis.
0
into the x_{min} edit field.1
into the x_{max} edit field.0
into the y_{min} edit field.Enter 0.2
into the y_{max} edit field.
0.01
into the Grid Size edit field.Press the Generate button to call the grid generation algorithm.
The density, heat capacity, and thermal conductivity are all assumed equal to one.
1
into the Density edit field.1
into the Heat capacity edit field.1
into the Thermal conductivity edit field.Although the problem doesn't involve any convection, the resulting time discretization term is here represented by a convective term in the positive ydirection (time flowing from y = 0 upwards).
0
into the Convection velocity in xdirection edit field.1
into the Convection velocity in ytdirection edit field.0
into the Heat source edit field.Enter 25
into the Initial condition for T edit field.
The modified equation must be changed to remove the thermal conduction term from the time axis (ydirection).
Enter k_ht*Tx_x + rho_ht*cp_ht*u_ht*Tx_t + rho_ht*cp_ht*v_ht*Tyt_t = q_ht
into the Equation for T edit field.
In the original onedimensional the initial temperature is a constant 25 degrees, in the transformed problem this is equivalent to prescribing a constant temperature at the bottom boundary.
The left and right boundaries are prescribed fixed temperature and heat flux conditions just as in the original problem.
25
into the Temperature edit field.25
into the Temperature edit field.1
into the Inward heat flux edit field.Lastly, the top boundary where time unknown is represented by a free/or convective outflow boundary condition.
Select Convective flux/outflow from the Heat Transfer dropdown menu.
Press OK to plot and visualize the selected postprocessing options.
The resulting solution shows how the uniform temperature field starts at 25 degrees at the bottom, and with time smoothly loses more and more heat from the left boundary as expected. The solution can be compared with the one dimensional problem and analytic solution.
The spacetime transformation of heat conduction 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_heattransfer8 script file), or GUI script (.fes) file.