Realization of Thermal Inertia in Frequency Domain

To realize the lagging behavior in heat conduction observed in these two decades, this paper firstly theoretically excludes the possibility that the underlying thermal inertia is a result of the time delay in heat diffusion. Instead, we verify in experiments the electro-thermal analogy, wherein the thermal inertial is parameterized by thermal inductance that formulates hyperbolic heat-conduction. The thermal hyperbolicity exhibits a special frequency response in Bode plot, wherein the amplitude ratios is kept flat after crossing some certain frequency, as opposed to Fourier heat-conduction. We apply this specialty to design an instrument that reliably identifies thermal inductances of some materials in frequency domain. The instrument is embedded with a DSP-based frequency synthesizer capable of modulating frequencies in utmost high-resolution. Thermal inertia implies a new possibility for energy storage in analogy to inductive energy storage in electricity or mechanics.


Introduction
Thermal inertia in heat transfer was observed in the lab or in practice during the past decades.As early as 1946, Bosworth postulated his experimental observation on the principles of electro-OPEN ACCESS thermal analogy, wherein the electric inductance has a thermal analogy in heat transfer [1,2].In the last decade of the twentieth century, Kaminski and Vedavarz et al. tabulated the relaxation times of heat conduction in metals, superconductors, semiconductors, organic materials and porous materials [3][4][5].In 2003, Roetzel et al. made another table for sand, NaHCO 3 and meat [6].Between these two quantitative measurements, in 1995, Mitra et al. found the conduction of heat in processed meat to be propagated like a wave, providing evidence that heat-transfer dynamics are hyperbolic [7].In 2002, Jiang et al. applied laser pulses to excite apparent non-Fourier phenomena in porous materials of small-scale dimensions [8].Sousa et al. claimed that the theory of hyperbolic heat transfer in agar gel improves the accuracy of radiofrequency ablation models and procedures [9].Currently, thermal inertia seriously affects nano-scaled or bio-heat treatments, say, in [10][11][12].The thermal inertia under these investigations has been mathematically modeled by time-delay in heat diffusion or thermal inductance with an electro-thermal analogy.
In fact, thermal-inductance modeling can be traced back to one century to when Maxwell [13] and Nernst [14] suggested through theoretical observations that, in properly chosen conductors, at low temperatures, heat may have sufficient inertia for oscillatory discharge to happen.In 1944, Peshkov measured the propagation velocity of heat flux in liquid helium at 1.4 K to be 19 m/s [15], which contradicts Fourier's law of heat conduction.Taitel hypothesized that the transient temperature at the middle location of a slab with constant-temperature heat source at two ends can be higher than that of the sources, a phenomenon which is called Taitel's paradox [16].
Along with those experimental observations mentioned above there are a number of numerical studies [17][18][19][20][21][22], wherein first-order approximation of time-delayed modeling results in thermal-inductance modeling.Therefore, it is hard to distinguish thermal-inductance modeling from time-delay modeling merely based on computer simulation.On the other hand, these numerical responses matching experimental observations indirectly support the physical realization of thermal inertia by thermal inductance.Firstly in this paper, we will give a theoretical proof that the physical realization of thermal inertia by time-delayed diffusion contradicts the first law of thermodynamics.
Although these two realizations can have indistinguishable responses of temperature or heat-flux in computer modeling, only the thermal-inductance realization involves the following aspects of physical significances: (1) There is a new "storage-form" since energy is stored in an inductive element though which a current of entropy flows; (2) Thermal inductance dominates the performance of high-speed thermal processes, such as the heating by short-time laser pulses, and (3) Even with small temperature differences, large energies can still be pumped in the form of heat-flux into the material, provided that we find how to increase its thermal inductance.
The literature documents three aspects of non-Fourier phenomena in heat transfer, which stimulated our efforts to understand more about thermal inertia: (nF1) The temperature at some locations of heat-conduction materials subjected to a constant temperature heat source can be higher than that of the source during the transience of heat transfer, which implies that heat flux can flow from the lower temperature level to a higher level; (nF2) A temperature gradient pulse induces temperature vibrations at fixed locations, in which temperature ought to wave-propagate before internal energy is exhausted, and (nF3) At any location of a heat-conduction material, an abrupt change of temperature will not result in an abrupt change of heat flux-some time lags have been observed, from which we infer that the speed of heat-transfer is comparatively smaller than that of Fourier's prediction.
For these non-Fourier phenomena, Cattaneo [23,24] and Vernotte [25,26] proposed to replace the zero-order diffusion equation with the first-order equation: where T and q stand for temperature and heat-flux, respectively, and τ indicates the observed time-lag.
Coupling the CV-Equation to thermal capacitance yields a hyperbolic equation, which can capture the phenomena described in nF1, nF2 and nF3.However, the CV-Equation seems not in line with the Clausius formulation of the second law of thermodynamics, since the heat-flux is allowed to flow into higher-temperature regions from a lower-temperature region in the transient state.The mainstream thus acknowledged that the CV-Equation ought to be a first-order approximation of the time-delayed heat-diffusion equation: which apparently still complies with the Clausius formulation.Coupling the TD-Equation to thermal capacitance, on the other hand, yields a parabolic equation with time-delayed Laplacian (TDP-Equation).
Unfortunately, we will perform a bifurcation analysis to show that the TDP-Equation is ill-posed, that is, given any initial condition, none of (bounded) solutions exist, no matter how small the timedelay τ is.Specifically, accompanied by the ill-posed dynamics is the fact that even a tiny heat-flux pulse input to the heat-transfer domain insulated from a heat source can trigger the temperature to be unbounded.This contradicts the first law of thermodynamics that excludes the possibility that thermal potential, part of the internal energy, becomes infinite due to a finite energy input to the domain.Thus, the hypothesis of time-delayed heat diffusion is untrue.
The bifurcation analysis adopts the composite of Laplace transform and Galerkin projection: the Laplace-Galerkin transform in space-time [27].With this newly developed transform, we represent the TDP-Equation by a transfer function served for pure algebra, which is a ratio of polynomials dependent on two variables: one is for space and the other is for time.With the capability of pure algebra, we construct the dynamics as a feedback interconnection of two 2D transfer-functions: one is for thermal capacitance, and the other is for the time-delayed heat diffusion.This makes it possible to extend the Nyquist criterion and Root Locus on stability analyses.
On the other hand, the thermal circuit resulting from the electro-thermo analogy captures the non-Fourier phenomena and can be formulated as a positive-damped, hyperbolic equation, compatible with the passivity of heat conduction.As mentioned, the hyperbolic nature of heat conduction seems incompatible with the Clausius formulation of the second law of thermodynamics.This incompatibility can be remedied by the inclusion of the thermal kinetics stored in the thermal inductor into the internal energy concept.Such a remedy is supported by the theory-extended irreversible thermodynamics [28,29].
One step further from the theoretical work [30], in this work we perform experimental measurements to prove that thermal inertial can be parameterized by thermal inductance within the electro-thermal analogy.In advance of these designing experiments, it is numerically found that thermal hyperbolicity exhibits a special frequency response, for example, in a Bode plot the amplitude ratios are no longer decreasing after some frequency, as opposed to parabolic Fourier heat-conduction, those of which are monotonically decreasing along with the increase of input frequencies.This frequency-domain specialty is then applied to reliably identify thermal inductances of some materials.In fact, measurements of thermal inertia have been documented in several papers, among which are [7] and [31].They identified the time delays in heat conduction in the time domain, but we employ the thermal hyperbolicity to reliably discover in the frequency domain the fact that thermal inertia is due to thermal inductance rather than the time-delay.An instrument is made for this purpose that is embedded with a DSP-based frequency synthesizer to obtain high-resolution Bode plot, especially in the concerned low-frequency region, which is beyond the flavors of current frequency generators.Therefore we design a limit-cycle generator and implement it into DSP-engine embedded microcontrollers (DSP-controller) to fulfill this task.
As a limit-cycle generator, the filtered van der Pol oscillator is implemented into a DSP-controller that provides sinusoidal control signals and is capable of modulating frequency in a real-time fashion.Its dynamics is transformed into iterated computation of addition and multiplication performed by the DSP engine inside.This kind of DSP-based implementation [32][33][34] is in line with the function of the DSP engine that is a master of fast calculation of addition and multiplication of floating numbers.Microcontrollers are really needed in practice for frequency synthesizers, since analog circuits lack communication ports, crystal oscillators have no reliable low-frequency outputs, and general-purpose computers transport signals to peripheral devices very slowly.Furthermore, DSP-based linear oscillators are not candidates, since they produce in the long run no sustained oscillations.
For the time being, there are two popular frequency synthesizers in industry: one is by data scheduling [35] and the other is by online sinusoidal function [36,37], both of which are widely used in wireless communication.In the data scheduling type, the working chip stores a table of data in code space that digitally shape a sinusoidal function.The targeted frequency is determined by the number of times the whole table is fetched per second at a fixed instruction clock, which is adjustable by the number of entries to be fetched in the entire table.This requires a large ROM size for long tables and produces wide frequency-span for radio modulation.Such ROM size is usually unavailable in DSP-controllers, and building special ICs for this purpose prevails [38][39][40][41].One characteristic of data scheduling for frequency synthesizers is low frequency-resolution, since the number of entries being fetched is always an integer.This situation is worse at lower frequency-range.It prevents this method from being applied to perform system identification of heat conduction in frequency domain.On the other hand, with the limit-cycle approach as developed by this paper, the frequency resolution can be pushed to the maximum allowed by hardware, which is usually more than 16 megabytes.
In the type of on-line sinusoidal function, the working chip has to compute sinusoidal functions one time within each sampling period, wherein the sinusoidal function is approximated by Taylor-series expansion [37] or by other complex algebra.Such a high computational burden is beyond the capability of DSP-controllers, and often requires FPGA chips.However, we really need DSP-controllers for system identification, since not just the frequency synthesizer is in the chip, but also an on-line window Fourier transform, and color and white noises filters to remove sensing contamination.All of these modules are implemented through DSP-based implementation.Otherwise, a special IC or a single-task FPGA cannot be integrated into other modules to fulfill control purposes.Moreover, there is usually a variety of peripherals in DSP-controllers served for control purposes.For example, almost every Microchip dsPIC has analog-to-digital (A/D) and pulse-width modulation (PWM) peripherals, for inputting analog signals from sensors and outputting PWM signals to switching power converters.

2D Transfer Function
This section briefs the mathematical tool named by 2D Transfer Function, which will appear in the remaining sections.
Consider a class of Laplacian operators ) defined over the set of spatial functions: is a bounded domain.Based on the Sturm-Liouville theory [42], such an operator A is nonnegative-definite and self-adjoint, which implies that all of its eigenvalues are nonnegative real numbers, and its normalized eigenfunctions constitute an orthonormal, complete basis of . Denote its countable set of eigenvalues by  and the eigenfunction associated with , from spatial functions to modal functions is defined by: Completeness and orthonormality of the eigenfunctions basis imply that the Galerkin transform G has a unique inverse which is the composite of the Galerkin transform and the Laplace transform, and explicitly: Accordingly, the inverse Laplace-Galerkin transform is the composite of the inverse Laplace transform and the inverse Galerkin transform, that is: and explicitly: where C is a closed contour that encircles poles of where h is a ratio of two expressions of finite or some infinite length constructed from two independent variables, one standing for space and the other for time, allowing for the operations of addition, subtraction, multiplication, integer exponents in time, and fraction-order exponents in space.For example: Linear dynamics .
Take as an explanatory example the longitudinal wave dynamics G ˆ: We firstly check that the negative Laplacian Taking Laplace-Galerkin transform H on both sides of the differential equation yields: That is, the 2D transfer function of the dynamics G ˆ is: As the input u is that is the spatial-temporal unit-pulse since its Laplace-Galerkin transform is one, the corresponding output defines the impulse response g of the dynamics G ˆ. Taking the inverse Laplace-Galerkin transform of the 2D transfer function G yields the impulse response: To verify this solution, let us integrate the differential equation from . Thereby, the impulse response becomes the solution of the following homogeneous differential equation: which can be solved via the conventional Separation-of-Variables method.It is easy to find that both arrive at the same solution.
Representation of partial differential equations by 2D transfer functions renders its stability analysis to be in essence the same as 1D stability analysis.For example: (1) For a 2D transfer function: where M and N are polynomials of s and allows for fraction-orders of  , the set: contains all of its poles.Its coherent dynamics G ˆ is asymptotically stable if all poles are on the left-half plane.(2) For a dynamics that is represented by feedback-interconnection of G and H that have no poles in the right half-plane the assembly of loop transfer functions is defined by the set . Then, this dynamics is unstable if any member of L has its Nyquist plot encircling ) 0 , 1 ( in the complex plane.The instability means that the output of the loop will grow exponentially up to be unbound even the input is as tiny as infinitesimal.

Time-Delay Heat Conduction
Assume that there really exists a time-delay  in heat diffusion, i.e., Coupling it with the equation of energy conservation: leads to the time-delayed parabolic equation: over a bounded domain  , where Q stands for the heat generation in the interior of  .Consider the class of boundary conditions: , where 0 T and 0 q stands for inputs of temperature and heat-flux, respectively, on the boundary   .
Denote the equilibrium point of Equation ( 26) by T , i.e., and the perturbation around the equilibrium T by  .That is, the temperature T and the heat generation are decomposed, respectively, by: Note that Equation ( 26) has only the equilibrium point T that indicates the steady-state temperature observed in Fourier and non-Fourier heat-transfers.Substituting Equations ( 29) and (30) into Equation ( 26) yields the equation governing the perturbation dynamics: with homogeneous boundary conditions: 0   or 0   on   .Let Equation (31) be further rewritten as: where is a positive number.Equation ( 32) is one of parabolic equations with time-delayed

Laplacian of homogeneous boundary conditions (TDP-Equation).
In this section it will be shown that the TDP-Equation of Equation ( 12) is unstable, no matter how small the time delay  is.We will prove this with the 2D transfer-function in Section 2.
As mentioned in Section 2, all eigenvalues  's of the Laplacian operator 2    in Equation (32) are nonnegative real, and its eigenfunctions  's can constitute an orthonormal, complete basis of . Based on Equation ( 12), the 2D transfer-function of the TDP-Equation in Equation ( 32) is known as: This transfer function can further be represented by the feedback-interconnection of thermal capacitance s / 1 and time-delayed heat diffusion s e    , as shown in Figure 1.Nyquist criterion [43]   can then be applied to do stability analysis for this loop, as discussed in Section for all delay times 0   , so the feedback loop in Figure 1 is unstable.Therefore, we come to the result that the perturbation dynamics of Equation ( 32) around the unique equilibrium point T is unstable, no matter how small the delay time  is.Another view on the instability of Equation ( 32) is provided as follows.The 2D transfer-function in Equation ( 33) has its poles parameterized by each eigenvalue  : Since the transfer function s e s /   has infinite zeros in the right half-plane (RHP), which can be known from Pades' approximation of s e   by a ratio of polynomials [43], the Root Locus regarding to Equation (34) will go across the imaginary axis into the RHP as the positive gain  is large enough.
Therefore, there are always poles of the 2D transfer function in Equation (33) being in the RHP.That is, the perturbation dynamics in Equation ( 32) is unstable, as proved in the preceding paragraph.33) starts to have RHP poles when the eigenvalue  is larger than 100 , and are getting more RHP poles along with the increase of  .This yields the same message as that from the  As of such, even the time-delay in heat diffusion is infinitesimally small, no steady-states of temperature will be possibly observed.The inexistence of observably steady states arises from the ill-posedness of the mathematical modeling, instead of the size of parameter-the time delay in heat diffusion.This implies that the realization of non-Fourier phenomena by the time-delay in heat diffusion is contrary to the truth.Specifically, consider the heat-transfer domain  is insulated from heat source across the boundary, i.e., on the boundary of which gradients of temperature are zero.With the same analysis, this heat-transfer dynamics can be proved to be internally unstable.That is, even a tiny heat-flux pulse input to the domain triggers the temperature distributions to be unbounded.This implies that thermal potential (or the internal energy) is infinite due to a finite energy input to the domain, which obviously contradicts the first law of thermodynamics.

Thermal Inertia Parameterization
Once thermal inertia does not come from time-delayed diffusion of heat, let us consider the electro-thermal analogy in this section.
The Fourier heat-conduction is a coupling of the heat diffusion: and the energy conservation: where k stands for the thermal conductivity,  for the mass density, v C for specific heat (at constant volume), and Q for the heat generation.On grounds of electro-thermal analogy, where we choose the temperature T as a thermal voltage and the heat-flux q as a thermal current, it is inferred from Equations ( 35) and ( 36) that a differential element is of the thermal resistance being To fulfill the electro-thermal analogy, we include thermal inductance L into Equation (35) to capture the thermal inertial observed in Non-Fourier heat transfer: and write Equation (36) as: Equations ( 37) and ( 38) represent the heat-conduction in a differential element as an infinitesimal RLC circuit.Application of Kirchhoff's theorems or careful combination of Equations ( 37) and ( 38) yields a hyperbolic equation: The equilibrium point T is the same as that of the parabolic equation governing Fourier heat-conduction, and the perturbation  imposed on T is governed by: with a homogeneous boundary condition: Let us choose a characteristic length  indicating the size of the heat-conduction domain  and a characteristic time  indicating the time lag in heat diffusion R L /   , together with a referenced temperature r T , to make independent and dependent variables in Equations ( 41) and (42) dimensionless through: Substituting those dimensionless variables in Equation ( 42) into Equations ( 40) and ( 41) yields a non-dimensional equation: where 0  stands for the normalized domain of heat conduction, for the dimensionless entropy flow-rate into the heat-transfer domain, and Performing Laplace-Galerkin transform H on both sides of Equation ( 43) yields the 2D transfer-function G of the hyperbolic heat-conduction dynamics of Equations ( 43) and (44) to be: are in the left-half plane (LHP), or possibly, there is a single pole at the origin and the others are in LHP.This implies that the perturbation dynamics in Equation ( 40) around the unique equilibrium T is always internally stable, and therefore the hyperopic formulation of Non-Fourier heat transfer is of passivity and well-posedness.Given the heat generation u and the initial conditions: , one can perform the Laplace-Galerkin transform to yield: Then he/she can perform the inverse Laplace-Galerkin transform The transfer function of Equation ( 45) reveals that the  -mode response possesses a damping ratio   . As this dimensionless quantity  is getting larger, thermal inertia is getting easy to be observed.We give it a name: 2nd-sound number.

Measurement Instrumentation Design
The parameterization of thermal inertia into thermal inductance makes possible to design a measurement instrument to perform parametric identification.As shown in Figure 4, the instrument to measure thermal inductance comprises a short-pulse laser and its driver for boundary actuation, low-frequency limit-cycle generator, USB4711 data acquisition, Matlab-to-dsPIC networking, distributed temperature sensor and cooler/heater, instrumentation amplifiers, constant-temperature heat sinker, and some supported mechanism.Meat is firstly chosen as the heat-conduction material.The short-pulse laser inputs pulse-width modulated (PWM) thermal current to the right boundary of the conduction material, where the temperature is converted by K-type thermocouples into voltage signals that are tremendously amplified by the instrumentation amplifier and then acquisitioned into a generous-purpose computer through USB4711 communication.The Matlab-to-PIC software is made for integration of dynamics synthesis in the generous-purpose computer and implementation of dynamics into DSP-based controllers.The DSP-based limit-cycle generator plays as the frequency synthesizer of the PWM duty of heat flow from the short-pulse laser.
We collect as the frequency response y of the heat conduction the steady-state sinusoidal responses across a span of input frequencies of heat flow u imposed on a preset DC-level.Temporal responses can further be sent into spectrum-decomposition to calculate out reliable amplitude ratios, figuring out a Bode magnitude plot.Based on Equation (45), the dimensionless frequency response in this experimental setup is where the modulus indicates the amplitude ratio and phase difference, respectively, of the output to the input as the input frequency is  .Figure 5 shows the Bode magnitude plot according to Equation (48), wherein Fourier heat-conduction is also plotted in dotted lines for comparison.The Bode plot in Figure 5 shows special frequency responses in a non-dimensional fashion, dubbed thermal hyperbolicity, listed as follows: (1) As opposed to the parabolic Fourier heat-conduction, the amplitude ratio of which is monotonically decreasing with the input frequency, the amplitude ratio of Non-Fourier heat-conduction is no longer decreasing in average after crossing some certain frequency named crossover frequency.
(2) After the crossover frequency, the amplitude ratio exhibits a series of resonant ripples at natural frequencies of all modes.The amplitudes of ripples are getting larger along with the increasing 2nd-sound number  .(3) Before the crossover frequency, the decay rate of amplitude ratio is reduced along with increase of the 2nd-sound number  .Any Non-Fourier heat-conduction has smaller decay rate than that of Fourier heat-conduction, which is nearly 10 dB per  log .
(4) Since larger amplitude ratios follows larger 2nd-sound number  , identification of thermal inductance in frequency domain is more accurate in a smaller scaled heat-conduction domain.
With these special properties, we are able to accurately and reliably identify thermal inductance in the frequency domain by adjusting the length of heat-conducting materials.Note that the thermal vibration implied in the non-dimensional modeling of Equation ( 43) involves a constant damping across all modes.This means that the damping ratio of each mode is inversely proportional to the square-root of its stiffness.In the Bode plot, both higher stiffness and damping ratio yield smaller amplitude ratios and both effects arrive at a balance after the crossover frequency.
To fulfill frequency-domain identification, a limit-cycle generator is programmed into a microcontroller embedded with a DSP engine (named by DSP-controller), frequencies of which are adjustable in a real-time fashion.In this work, we generalize the van der Po oscillator as the nominated limit-cycle generator.The van der Pol oscillator is a non-conservative oscillator with non-linear damping, which evolves in time according to the second order differential equation.Let us consider a more general class of van der Pol oscillators as: where y is the output, and three parameters-a, ω and µ are to be adjusted.The dynamics of Equation (49) generates a sustained oscillation with amplitude a and angular frequency ω in steady state, and the transience period is getting shorter along with increasing value of µ.
In our design, the frequency synthesizer is fulfilled by series of such a van der Pol oscillator with a first-order low-pass filter to shape outputs into a pure sinusoidal function even when µ is large at will for ultra short transience.Accordingly, the frequency synthesizer is the van der Pol dynamics: which has been represented by internal feedback-interconnection, in series with the first-order filter: By choosing three state variables: , the state-space realization of this limit-cycle generator becomes: It appears as a feedback-interconnected LPV form: , where the command (angular) frequency  , amplitude a and converging speed  constitute the slow-time parameter that will be on-line scheduled in a slow-time fashion.The frequency synthesizer in Equation ( 53) is then implemented into a DSP-engine embedded microcontroller with the Linear Parameter-varying Matrix Iteration Time-fixed (LPV-MIT) method as follows.Continuous-to-digital conversion of Equation (53) arrives at: where k and j are fast-time and slow-time indexes, respectively, and the system matrices ) , (   will be calculated on-line in a slow-time fashion since the resonant frequency  is slow-time updated in practice.Explicitly: where T is the fast-time sampling period that is short enough to make the above Taylor-series expansion accurate even with few terms, say, 2  n or 3 as T is 001 .0 second.At the present time t k the DSP-controller merely stores the current state x k , and the update of state from x k to x k+1 at next instant is fulfilled by DSP engine performing the addition and multiplication of floating numbers according to Equation (54).In that sense, any time can be treated as the initial time, which makes real-time processing most efficient.Moreover, the state update intervals are held identical to those in computer simulation, so that the real-timed operation matches the dynamics that has been verified by offline calculation, thus achieving robust implementation.At any instance, a pulse-width modular (PWM) signal in line with the output y is sent to the gate-driving circuit of switched power converters that drive the short-pulse laser.

Thermal Inductance Identification
With the instrument and instrumentation designed above, we start to measure temperature responses of meat excited by the short-pulse laser.
Figure 6 shows the typical responses as the laser delivers a series of pulses to the right boundary.Compared with the calculated response of Fourier heat-conduction, the measured non-Fourier response exhibits diffusion lagging initially, but quickly they are indistinguishable to each other until the thermal wave is reflected from the left end back to the right end as a small jump to be measured.This standing wave behavior is getting obvious as the length of meat is cut shorter.However, we find that the identification of thermal inductance through the measurement of standing-wave behavior is unreliable.Similar problems also happen at other measurements in time domain, wherein catching the jump with high-bandwidth sensing is always severely contaminated by white-noise, and filtering out white-noise often seriously degrades the jump-signal.Figure 7 shows the step response as the control signal to laser is a series of steps.Although the time constant therein reveals the information of thermal inductance, however, it is very close to that of Fourier heat-conduction.Moreover, as observed, the corresponding thermal current into the heat-conduction domain is not a pure step, so identification of thermal inductance with step responses is inaccurate or unreliable.Identification of thermal inductance in frequency domain is the most reliable because: (1) The noisy contamination of amplitude ratio can be perfectly removed with spectrum decomposition since merely the component at the input frequency is concerned, as shown in Figure 8. (2) The crossover frequency quantifies the thermal inductance and it is against any uncertainties that can happen, such as the real values of amplitude ratios.The finding of thermal inductance gives us fresh impetus to manufacture storages of thermal kinetics, and to control the thermal vibrations in small-scaled or fast-timed scaled heat-transfer dynamics.

Figure 1 .
Figure 1.Feedback realization of time-delay parabolic heat transfer.

Figure 2
Figure 2 shows the Nyquist plot s e s /  

Figure 2 .
Figure 2. Nyquist plot of s e s /  

Figure 3 .
It shows that the transfer function in Equation (

Figure 2 ,
Figure 2, in which the Nyquist plot of the transfer-function in Equation (33) begin to encircle ) 0 , 1 ( as the eigenvalue  is larger than 100 , and the number of encirclements is increased with larger  .

Figure 9
Figure 9 records the measured Bode plot upon 0.25 mm-length meat, wherein there are 130 points ranged from 01 .0 Hz to 10 Hz .This Bode plot evidently reveals the thermal hyperbolicity resulted from parameterization of thermal inertia into thermal inductance, and the crossover frequency being 4 Hz .This verifies the electro-thermal analogy and identify the thermal inductance with 83 .3  L W K m s   .In advance, we have measure the following quantities for the processed meat at 25 °C: the density 980  

Figure 9 .
Figure 9. Identification of thermal inductance in frequency domain.

Figure 10 .
Figure 10.The effect of heat-conduction domain of frequency-domain identification.
Table 1 lists the other three materials at room temperature.

Table 1 .
Thermal parameters measured at 25 degrees centigrade.