TY - JOUR
T1 - A finite element method for non-Fourier heat conduction in strong thermal shock environments
AU - Wang, Bao-Lin
AU - Han, Jie-Cai
PY - 2010
Y1 - 2010
N2 - Non-Fourier effect is important in heat conduction in strong thermal environments. Currently, generally-purposed commercial finite element code for non-Fourier heat conduction is not available. In this paper, we develop a finite element code based on a hyperbolic heat conduction equation, which includes the non-Fourier effect in heat conduction. The finite element space discretization is used to obtain a system of differential equations for the time. The transient responses are obtained by solving the system of differential equations, based on the finite difference, mode superposition, or exact time integral. The code is validated by comparing the numerical results with exact solutions for some special cases. The stability analysis is conducted and it shows that the finite difference scheme is an ideal method for the transient solution of the temperature field. It is found that with mesh refining (decreasing mesh size) and/or high-order elements, the oscillation in the vicinity of sharp change vanishes, and can be essentially suppressed by the finite difference scheme. A relationship between the time step and the space length of the element was identified to ensure that numerical oscillation vanishes.
AB - Non-Fourier effect is important in heat conduction in strong thermal environments. Currently, generally-purposed commercial finite element code for non-Fourier heat conduction is not available. In this paper, we develop a finite element code based on a hyperbolic heat conduction equation, which includes the non-Fourier effect in heat conduction. The finite element space discretization is used to obtain a system of differential equations for the time. The transient responses are obtained by solving the system of differential equations, based on the finite difference, mode superposition, or exact time integral. The code is validated by comparing the numerical results with exact solutions for some special cases. The stability analysis is conducted and it shows that the finite difference scheme is an ideal method for the transient solution of the temperature field. It is found that with mesh refining (decreasing mesh size) and/or high-order elements, the oscillation in the vicinity of sharp change vanishes, and can be essentially suppressed by the finite difference scheme. A relationship between the time step and the space length of the element was identified to ensure that numerical oscillation vanishes.
UR - http://handle.uws.edu.au:8081/1959.7/548271
U2 - 10.1007/s11706-010-0090-4
DO - 10.1007/s11706-010-0090-4
M3 - Article
SN - 1673-7377
VL - 4
SP - 226
EP - 233
JO - Frontiers of Materials Science in China
JF - Frontiers of Materials Science in China
IS - 3
ER -