Table of Contents Author Guidelines Submit a Manuscript
The Scientific World Journal

Volume 2014 (2014), Article ID 671619, 7 pages
Research Article

Fast Transient Thermal Analysis of Non-Fourier Heat Conduction Using Tikhonov Well-Conditioned Asymptotic Waveform Evaluation

Department of Electrical Engineering, Faculty of Engineering, University of Malaya, 50603 Kuala Lumpur, Malaysia

Received 5 February 2014; Accepted 20 May 2014; Published 11 June 2014

Academic Editor: Navin K. Rastogi

Copyright © 2014 Sohel Rana et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Non-Fourier heat conduction model with dual phase lag wave-diffusion model was analyzed by using well-conditioned asymptotic wave evaluation (WCAWE) and finite element method (FEM). The non-Fourier heat conduction has been investigated where the maximum likelihood (ML) and Tikhonov regularization technique were used successfully to predict the accurate and stable temperature responses without the loss of initial nonlinear/high frequency response. To reduce the increased computational time by Tikhonov WCAWE using ML (TWCAWE-ML), another well-conditioned scheme, called mass effect (ME) T-WCAWE, is introduced. TWCAWE with ME (TWCAWE-ME) showed more stable and accurate temperature spectrum in comparison to asymptotic wave evaluation (AWE) and also partial Pade AWE without sacrificing the computational time. However, the TWCAWE-ML remains as the most stable and hence accurate model to analyze the fast transient thermal analysis of non-Fourier heat conduction model.

1. Introduction

The classical Fourier heat conduction law follows a property, which is not physical, that is, an infinite thermal wave propagation speed [13]. Fourier’s law is quite accurate for most common engineering problems. In some cases, the Fourier law cannot explain, such as near the absolute zero temperature, thermal gradient, which is extreme [3]. Non-Fourier conduction accounts for two-phase lag due to thermal wave propagation speed and another one due to relaxation time of electron. D. Y. Tzou proposed a non-Fourier model based on two-phase lag as follows [1, 3, 4]: where stands for finite relaxation time for electron phonon and represents finite thermal wave speed [5].

Conventional solver based on iterative, as fourth-order Runge-Kutta technique, is computationally expensive, though its result is accurate. In contrast, asymptotic waveform evaluation (AWE) technique is much faster solver [6]. Taylor’s series was used to expand the AWE transfer function [710]. Then the required moments are found from the coefficient of Taylor’s series [911]. AWE technique is able to solve up to three-dimensional models [12]. However, the limitation of AWE model is that this model cannot predict the temperature responses accurately as it is ill-conditioned moment matching [6, 13]. Then Loh et al. [6] developed technique, which is based on the AWE to remove the instability of AWE, called partial Pade AWE. In partial Pade AWE, selected poles and residues from any heat imposed boundary are used to calculate the temperature responses of the whole system. However, this technique is unable to predict the actual temperature. This prompted the present study to use another technique, which is based on well-conditioned asymptotic waveform evaluation (WCAWE).

WCAWE algorithm was introduced by Slone et al. [13, 14] to solve frequency domain electromagnetic problem based on the FEM. In this model, a correction term was introduced during the poles and residues calculation. Maximum likelihood (ML) was used successfully to calculate poles and residue to bring down the instability of AWE and also partial Pade AWE. Later, Liu et al. solved Fourier heat conduction by using the WCAWE technique [15]. The model proposed in the present work is based on the WCAWE algorithm to analyze the non-Fourier thermal problem. To avoid the singularity problem, our proposed model is embedded with Tikhonov regularization technique [16] to further enhance stable responses. Hence, the present model is recognized as TWCAWE. Later, we developed another well-conditioned scheme, called TWCAWE using the mass effect (TWCAWE-ME) that reduces the computational time compared to TWCAWE using the ML scheme (TWCAWE-ML).

2. Non-Fourier Model

The non-Fourier model with two-phase lag is shown in (2) in the case of fast thermal conduction. This model is converted into a normalized hyperbolic equation given by where where is the thermal diffusivity.

FEM meshing is performed to generate the rectangular element. The nodes of the element are denoted by , , , and . Applying FEM based on the Galerkin weighted residual technique of (2), we obtain

The above non-Fourier heat conduction equation shown in (4) is converted into a linear equation based on the Galerkin weighted residual technique, as shown in where mass matrix , damping matrix , stiffness matrix , and . The elemental matrices obtained from (4) are as follows: The Laplace transform of (5) is shown below: The system solution can be approximated by using polynomial equation in s-domain, as given by the following equation: The moments, , are the coefficients of Taylor series expansion about (Maclaurin series).

3. TWCAWE Algorithm

The TWCAWE algorithm is initiated by applying the Laplace transform on (5) where the moments obtained are used to estimate the zero state response (ZSR) and zero input response (ZIR). The moment computation for ZSR and ZSR is shown in (10) and (12), respectively.

3.1. Zero State Response

In ZSR, the initial condition is assumed to be zero, and : By equating the same powers of , the moments are generated from the equation below: where is the order of Pade approximation.

3.2. Zero Input Response

In ZIR, the input force is assumed to be zero, : By equating the same powers of , the moments are generated from the equation below:

The moment matching process in AWE is inherently ill conditioned. To overcome this instability, WCAWE was proposed [13]. In our proposed model, to reduce the instability, we removed the singularity problem by changing the stiffness matrix. Consider a well-conditioned approximation problem ; the residual becomes smaller for the proper choice of [16]. The singularity condition of is removed by adding term where, after modification, it becomes where is the regulation parameter, which depends on the order of the equation, whereas is the identical matrix. The approximate inverse family is defined by . As given in (10) and (12), moment calculation involves inverse of stiffness matrix . The resulting moments, especially higher order moments, are product of lower order moments and, hence, the singularity problem can accumulate the compounded error. Therefore, in TWCAWE, computing moments by the inverse of matrix is now replaced by inverse of . Subsequently, the ill-conditioned moments are converted to well-conditioned moments. The new moments, called TWCAWE moments, are obtained from ICAWE moments. The ICAWE moment subspace is and the TWCAWE moments subspace is denoted as follows: , where the relation between ICAWE and TWCAWE moments is where . Modified Gram-Schmidt technique is required to orthonormalize that gives a more accurate solution [11]. In TWCAWE algorithm, can be computed as shown in (16):

The nodal moment can be extracted from the global moment matrix for any arbitrary node as follows: The transient response for any arbitrary node can be approximated by using Pade approximation and then further simplified to partial fractions [17], as shown in

Poles and residues can be found by solving

4. T-WCAWE Using Maximum Likelihood Scheme

In T-WCAWE model, a correction term is introduced to obtain more stable responses. This correction term is calculated from matrix. The matrix is computed from ill-conditioned moments and well-conditioned moments, as shown in (16). This matrix is nonsingular [13, 14]. The correction term is obtained by using maximum likelihood (ML) from where is the number of nodes, is the order of the equation, and is the order of the product. Now, the new T-WCAWE poles and residues are calculated from (21) and (22), respectively:

5. TWCAWE with Mass Effect Scheme

In partitioning of FEM, it was reported that the resulting interpolation function can lead to singular or ill-conditioned. This singularity occurs because some interpolations are linearly dependent and ill-condition occurs because the functions are very close to each other. To minimize this singularity and ill-condition, the mass matrix is changed by introducing a factor “ ” [18]. For small value of , TWCAWE-ME scheme reduces the effect of mass, as given in where .

Here, where is the grid aspect ratio and is the length of each element in axis. and are consistent and row sum-lumped mass matrices, respectively. The lumped mass matrix is calculated using the total mass and allocating the lumped masses in the ratio of the diagonal element of the consistent mass matrix [19]. In this scheme, from (5) is replaced by new calculated .

5.1. Transient Response

The final transient response for any arbitrary node is where

6. Results Analysis

In the present study, the non-Fourier heat conduction is analyzed by considering two-dimensional rectangular configurations using FEM and Tikhonov based WCAWE algorithm. The simulation work was carried out on a rectangular two-dimensional slab meshed to generate the rectangular elements, as shown in Figure 1 where the instantaneous boundary condition was applied at the left boundary of the slab. Boundary condition employed for this simulation work is given as where of temperature at the left boundary edge from the initial value . The applied normalized boundary conditions are stated as follows:

Figure 1: Rectangular two-dimensional slab.

Figures 2, 3, and 4 show the temperature distribution for the above rectangular slab considering different values of as 0.5, 0.05, and 0.0001, respectively. The temperature spectrum shown in Figures 24 is plotted by considering the node situated in the middle of the slab; all the responses are calculated by TWCAWE-ML scheme. It is clear from the figure that all temperature responses are continuous and accurate for different values of . In this scheme, the individual poles and residues are used to calculate the individual temperature responses. Hence, this scheme accurately predicts the temperature distribution across the system. The initial high frequency and delay due to relaxation time of electron are clear in Figure 4 for .

Figure 2: Normalized temperature response spectrum along the center of the slab in the case of .
Figure 3: Normalized temperature response spectrum along the center of the slab in the case of .
Figure 4: Normalized temperature response spectrum along the center of the slab for the case of .

Though the fourth-order Runge-Kutta technique is slow, the results obtained from this method are more accurate due to its stability. Therefore, our proposed model is compared with the results obtained by Runge-Kutta technique in Figures 5 and 6.

Figure 5: Comparison with TWCAWE-ML, TWCAWE-ME, Runge-Kutta technique, and AWE using partial Pade in the case of .
Figure 6: Comparison of TWCAWE-ML, TWCAWE-ME, Runge-Kutta technique, and ICAWE for the case of .

To show the accuracy of our proposed model, we compared both schemes, TWCAWE-ML and TWCAWE-ME, against AWE using partial Pade and Runge-Kutta technique, as shown in Figure 5. The comparison has been presented with . The temperature behavior obtained from both schemes of TWCAWE is similar to Runge-Kutta technique. Partial Pade AWE, however, is unable to predict the temperature behavior as Runge-Kutta technique. In partial Pade AWE, arbitrarily chosen poles and residues are used to approximate the temperature behavior for the whole system. Hence, this method was unable to calculate the temperature behavior for other nodes. On the other hand, both schemes of TWCAWE model use its own poles and residue to calculate the temperature behavior.

The comparison between our proposed two schemes is shown in Figure 6 where . TWCAWE-ME is able to calculate the delay and high frequency, but the temperature behavior deviates from Runge-Kutta results. On the other hand, the temperature behavior obtained from TWCAWE-ML is in close proximity to Runge-Kutta technique.

To bring out the effectiveness of employing Tikhonov technique, we have performed another comparison in Figure 7 in the case of . This figure shows that the temperature behavior obtained from AWE technique largely deviates from Runge-Kutta behavior. In AWE model, unstable poles make moment matching process ill conditioned.

Figure 7: Comparison of TWCAWE-ML, WCAWE-ML, and ICAWE for .

The result of the WCAWE-ML model is not similar to Runge-Kutta technique due to the singularity problem of the stiffness matrix as shown in Figure 7. The error is compounded during moment calculation and it becomes more evident for higher order moment calculation. This causes WCAWE-ML to be unstable, thus producing inaccurate initial high frequency responses, making it difficult to observe the delay. However, by embedding Tikhonov technique into existing WCAWE-ML, singularity problem was minimized. This is carried out by adding a new term based on Tikhonov regulation technique as given in (13). The new term reduces the symmetrical nature of stiffness matrix, thus minimizing the singularity problem with negligible effect on the temperature behavior. By observing temperature behavior on node 5 which is near the imposed boundary condition, Tikhonov based WCAWE-ML (TWCAWE-ML) is in close proximity to Runge-Kutta technique, as shown in Figure 7, compared to WCAWE-ML which is devoid of Tikhonov technique. The model proposed in the present work is able to reduce the effect of singularity as well as instability of the previous model.

Table 1 shows the computational time required by analyzing different methods in the present work. From the table, TWCAWE-ML is found to be almost two times faster than Runge-Kutta technique. The orthonormalization and maximum likelihood approximation process in T-WCAWE poses computational load that makes it slower compared to AWE or partial Pade AWE. As TWCAWE-ME does not require orthonormalization and maximum likelihood approximation, it performs as fast as AWE or AWE using partial Pade despite some sacrifice in temperature response quality in comparison to TWCAWE-ML. However, TWCAWE-ME showed more stable and hence accurate results compared to either AWE or AWE using partial Pade.

Table 1: Required time for different methods.

7. Conclusion

In the present work, the heat conduction based on non-Fourier model has been solved using FEM. The phase lag responsible for finite relaxation time is varied to analyze the capability of our proposed model for predicting the heat conduction in two-dimensional model. From the results, both schemes discussed in this work are able to reduce the instability of AWE. From the comparison, we can conclude that T-WCAWE using ML results is in close proximity to Runge-Kutta technique compared to TWCAWE-ME. However, TWCAWE-ME scheme is 3.3 times faster than the conventional solver, namely, Runge-Kutta technique (RK4), and needs 1.7 times less computational time than T-WCAWE with ML scheme. Despite this, the TWCAWE-ML is more accurate for predicting the initial frequency as well as delay based on the close proximity to Runge-Kutta technique, hence emerging as the most stable and accurate technique to analyze the fast transient thermal analysis of non-Fourier heat conduction model when compared to AWE, partial Pade AWE, and TWCAWE-ME.


:Dimensionless temperature
:Dimensionless time
:Distance along “
:Distance along “
:Shape function matrix
:Element area
:Initial temperature
:Raised temperature
:Phase lag for temperature gradient
:Phase lag for heat flux.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This research work is supported by the University of Malaya High Impact Research (HIR) Grant (no. 51) sponsored by the Ministry of Higher Education (MOHE), Malaysia, as well as Science Fund (SF010-2013) sponsored by Ministry of Science, Technology and Innovation (MOSTI), Malaysia.


  1. G. Atefi, A. Bahrami, and M. R. Talaee, “Analytical solution of dual phase lagging heat conduction in a hollow sphere with time-dependent heat flux,” in New Aspects of Fluid Mechanics, Heat Transfer and Environment, vol. 1, pp. 114–126, 2010. View at Google Scholar
  2. B. Abdel-Hamid, “Modelling non-Fourier heat conduction with periodic thermal oscillation using the finite integral transform,” Applied Mathematical Modelling, vol. 23, no. 12, pp. 899–914, 1999. View at Publisher · View at Google Scholar · View at Scopus
  3. R. Quintanilla and R. Racke, “A note on stability in dual-phase-lag heat conduction,” International Journal of Heat and Mass Transfer, vol. 49, no. 7-8, pp. 1209–1213, 2006. View at Publisher · View at Google Scholar · View at Scopus
  4. J. Ordóñez-Miranda and J. J. Alvarado-Gil, “Exact solution of the dual-phase-lag heat conduction model for a one-dimensional system excited with a periodic heat source,” Mechanics Research Communications, vol. 37, no. 3, pp. 276–281, 2010. View at Publisher · View at Google Scholar · View at Scopus
  5. T. S. Cheah, K. N. Seetharamu, A. Q. Ghulam, Z. Z. Zainal, and T. Sundararajan, “Numerical modeling of microscale heat conduction effects in electronic package for different thermal boundary conditions,” in Proceedings of the 3rd Electronics Packaging Conference, vol. 3, pp. 53–59, 2000.
  6. J. S. Loh, I. A. Azid, K. N. Seetharamu, and G. A. Quadir, “Fast transient thermal analysis of Fourier and non-Fourier heat conduction,” International Journal of Heat and Mass Transfer, vol. 50, no. 21-22, pp. 4400–4408, 2007. View at Publisher · View at Google Scholar · View at Scopus
  7. S. Rana, K. Jeevan, R. Hari-krishnan, and A. W. Reza, “A well-condition asymptotic waveform evaluation method for heat conduction problems,” Advanced Materials Research, vol. 845, pp. 209–215, 2014. View at Google Scholar
  8. R. D. Slone, R. Lee, and J.-F. Lee, “Multipoint Galerkin asymptotic waveform evaluation for model order reduction of frequency domain FEM electromagnetic radiation problems,” IEEE Transactions on Antennas and Propagation, vol. 49, no. 10, pp. 1504–1513, 2001. View at Publisher · View at Google Scholar · View at Scopus
  9. C. R. Cockrell and F. B. Beck, “Asymptotic Waveform Evaluation (AWE) technique for frequency domain electromagnetic analysis,” NASA Technical Memorandum 110292, 1996. View at Google Scholar
  10. K. Gallivan, E. Grimme, and P. van Dooren, “Asymptotic waveform evaluation via a Lanczos method,” Applied Mathematics Letters, vol. 7, no. 5, pp. 75–80, 1994. View at Google Scholar · View at Scopus
  11. P. K. Chan, “Comments on “Asymptotic waveform evaluation for timing analysis”,” IEEE Transaction on Computer Aided Design, vol. 10, no. 8, pp. 1078–1079, 1991. View at Publisher · View at Google Scholar
  12. Y. Sun and I. S. Wichman, “On transient heat conduction in a one-dimensional composite slab,” International Journal of Heat and Mass Transfer, vol. 47, no. 6-7, pp. 1555–1559, 2004. View at Publisher · View at Google Scholar · View at Scopus
  13. R. D. Slone, R. Lee, and J.-F. Lee, “Well-conditioned asymptotic waveform evaluation for finite elements,” IEEE Transactions on Antennas and Propagation, vol. 51, no. 9, pp. 2442–2447, 2003. View at Publisher · View at Google Scholar · View at Scopus
  14. R. D. Slone, R. Lee, and J.-F. Lee, “Broadband model order reduction of polynomial matrix equations using single-point well-conditioned asymptotic waveform evaluation: derivations and theory,” International Journal for Numerical Methods in Engineering, vol. 58, no. 15, pp. 2325–2342, 2003. View at Publisher · View at Google Scholar · View at Scopus
  15. P. Liu, H. Li, L. Jin, W. Wu, S. X.-D. Tan, and J. Yang, “Fast thermal simulation for runtime temperature tracking and management,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 25, no. 12, pp. 2882–2892, 2006. View at Publisher · View at Google Scholar · View at Scopus
  16. A. Neumaier, “Solving ill-conditioned and singular linear systems: a tutorial on regularization,” SIAM Review, vol. 40, no. 3, pp. 636–666, 1998. View at Google Scholar · View at Scopus
  17. C. K. Ooi, K. N. Seetharamu, Z. A. Z. Alauddin, G. A. Quadir, K. S. Sim, and T. J. Goh, “Fast transient solutions for heat transfer,” in Proceedings of the Confernce on Covergent Technologies for the Asia-Pacific Region (IEEE TENCON '03), vol. 1, pp. 469–473, October 2003. View at Scopus
  18. S. Ham and K.-J. Bathe, “A finite element method enriched for wave propagation problems,” Computers and Structures, vol. 94-95, pp. 1–12, 2012. View at Publisher · View at Google Scholar · View at Scopus
  19. M. A. Christon, “The influence of the mass matrix on the dispersive nature of the semi-discrete, second-order wave equation,” Computer Methods in Applied Mechanics and Engineering, vol. 173, no. 1-2, pp. 147–166, 1999. View at Google Scholar · View at Scopus