Exact and Analytic-Numerical Solutions of Lagging Models of Heat Transfer in a Semi-Infinite Medium

and Applied Analysis 3


Introduction
Non-Fourier models of heat conduction have increasingly been considered in recent years to model microscale and ultrafast, transient, nonequilibrium responses in heat and mass transfer, where thermal lags and nonclassical phenomena are present (see, e.g., [1] and references therein).The growing area of applications of these models include, among other examples, the processing of thin-film engineering structures with ultrafast lasers [2,3], the transfer of heat in nanofluids [4,5], or the exchange of heat in biological tissues [6][7][8].
In the Dual-Phase-Lag (DPL) model [9][10][11], the equation relating the heat flux vector q and the temperature , for time  and spatial point r, q (r,  +   ) = −∇ (r,  +   ) , where  > 0 is the thermal conductivity, incorporates two lags,   for the heat flux and   for the temperature gradient.When both lags are zero, the Fourier law is recovered, while for   > 0 and   = 0, it reduces to the Single-Phase-Lag (SPL) model [12].
Combining (1) with the principle of energy conservation, where   is the volumetric heat capacity and , the volumetric heat source, in the absence of heat sources a partial differential equation with delay is obtained [13,14] as where  = /  is the thermal diffusivity.When both lags are zero, the diffusion equation, a parabolic partial differential equation which represents the classical model for heat conduction and other transport phenomena, is obtained.Using first-order approximations in (1), a hyperbolic equation is derived, commonly referred to as the DPL model [9], here denoted as DPL(1, 1), which for   = 0 reduces to the Cattaneo-Vernotte (CV) model [15][16][17].
In semi-infinite domains, some particular problems have also been considered.Solutions for heat propagation according to DPL(1, 1) model in a semi-infinite solid, produced by suddenly raising the temperature at the boundary, were obtained in [11,20], using Laplace and Fourier transforms.Relations between the local values of heat flux and temperature, in the form of integral equations, in a semi-infinite solid were considered in [13,28].
In this work, using Fourier transforms, explicit solutions for lagging models of heat conduction in a semi-infinite domain, with different types of boundary conditions, are obtained, allowing the construction of numerical solutions with bounded errors.
It should be noted that Fourier transforms can also be used in time-dependent problems (e.g., [29,30]), and the approach of this work could also be useful for timedependent DPL models, which have already been proposed [31].

Solutions of DPL Models in a Semi-Infinite Domain
Consider a plate of infinite thickness,  ∈ [0, ∞], that can be heated either at its surface,  = 0, or up to a certain depth,  ∈ [0, ].We will consider, for  ≥ 0, either Dirichlet, (0, ) = 0, or Neumann,   (0, ) = 0, boundary conditions and also that lim Appropriate initial conditions must be provided for the different models.Thus, for DPL(1, 1) initial values for temperature and its time derivative have to be specified, while for DPL(2, 1) and DPL(2, 2), its second derivative also has to be given, and for the DH model, the initial condition for the temperature has to be specified for a time interval of  amplitude, For a wide class of initial functions [32,33], the method of Fourier transform can be used to eliminate derivatives in the spatial domain and to obtain expressions for the exact solutions in the form of an infinite integral, either using Fourier sine transforms for Dirichlet conditions,  or cosine transforms for Neumann conditions, where the functions T(, ), which are the corresponding Fourier sine or cosine transforms of (, ), are obtained as solutions of the transformed temporal problems, depending on the continuous set of eigenvalues  2 .
Hence, these problems can be solved, obtaining expressions for T(, ) in terms of the roots of the corresponding characteristic equation, and thus explicit expressions for the exact solutions for these models, in the form of ( 13) or (14), can be obtained.
For the DH model, the transformed temporal problems are initial-value problems for delay differential equations with general initial functions, T  (, ) +  2 T (,  − ) = 0,  > , T (, ) =  (, ) , 0 ≤  ≤ , (19) where (, ) is the appropriate Fourier transform, according to the boundary conditions, of the initial function (, ).To obtain constructive solutions for this problem, a combination of the steps method and a convolution integral can be applied [34,35], producing the following expression, for  ∈ [, ( + 1)], T (, ) =  (, ) +  (, 0) The solutions obtained with the Fourier transforms, as given in (13) or (14), can be shown to converge and provide exact solutions under adequate integrability and regularity conditions on the initial functions.Numerical integration is required in general to compute numerical approximations of these solutions, with errors that can be bounded in finite spatial and temporal domains by controlling errors in the numerical integrators or by appropriately truncating the infinite integrals.However, for some particular initial functions, the solutions given by ( 13) or ( 14) may reduce to finite integrals.

Numerical Examples
Numerical examples are presented in the following figures, where, in order to properly compare DPL and DH models, the initial interval for DH, where the initial function (, ) is given, is displaced to [−, 0], and the initial functions for DPL models are set so that the values of temperature and its first derivative at  = 0, and also its second derivative for DPL(2, 1) and DPL(2, 2), are matched to those of the DH model.The classical diffusion model, whose solution is available and readily obtained [36], is also included as reference.
First, we consider models with   = 0, so that DPL(2, 1) and DPL(2, 2) are equal, and an initial function with damped temperature oscillations, thought to be the result of a modulated heat source that is switched off at  = 0, showing the transient behavior for the different DPL models for different values of  (Figure 1), as well as their differences from classical diffusion (Figure 2).
In Figure 3, a more detailed view of the spatiotemporal behavior of the DH model (Figure 3(a)) and differences from DH of DPL(1, 1) and DPL(2, 1) (Figure 3(b)) are presented.
In Figure 4, different values of   and   , such that  =   −   is kept constant, are used, so that variations in the temperature evolution are observed in the DPL approximate models, but not in the DH model, which only depends on the value of .

Figure 2 :
Figure 2: Differences from classical diffusion for models DPL and DH, for the data shown in Figure 1.