Next Article in Journal
Block Generalized Locally Toeplitz Sequences: From the Theory to the Applications
Next Article in Special Issue
A New Type of Generalization on W—Asymptotically J λ—Statistical Equivalence with the Number of α
Previous Article in Journal
Neutrosophic Incidence Graphs With Application
Previous Article in Special Issue
Some Summation Theorems for Generalized Hypergeometric Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Exact Solutions to Non-Fourier Heat Equations with Substantial Derivative

by
Konstantin Zhukovsky
*,
Dmitrii Oskolkov
and
Nadezhda Gubina
Faculty of Physics, M.V. Lomonosov Moscow State University, Leninskie Gory, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Axioms 2018, 7(3), 48; https://doi.org/10.3390/axioms7030048
Submission received: 10 April 2018 / Revised: 3 July 2018 / Accepted: 13 July 2018 / Published: 18 July 2018
(This article belongs to the Special Issue Mathematical Analysis and Applications)

Abstract

:
One-dimensional equations of telegrapher’s-type (TE) and Guyer–Krumhansl-type (GK-type) with substantial derivative considered and operational solutions to them are given. The role of the exponential differential operators is discussed. The examples of their action on some initial functions are explored. Proper solutions are constructed in the integral form and some examples are studied with solutions in elementary functions. A system of hyperbolic-type inhomogeneous differential equations (DE), describing non-Fourier heat transfer with substantial derivative thin films, is considered. Exact harmonic solutions to these equations are obtained for the Cauchy and the Dirichlet conditions. The application to the ballistic heat transport in thin films is studied; the ballistic properties are accounted for by the Knudsen number. Two-speed heat propagation process is demonstrated—fast evolution of the ballistic quasi-temperature component in low-dimensional systems is elucidated and compared with slow diffusive heat-exchange process. The comparative analysis of the obtained solutions is performed.

1. Introduction

Recent progress in technology and science has driven interest to studies of heat conduction beyond common Fourier law [1]: t T = k x 2 T , where k is the thermal diffusivity, T is the temperature; Fourier law describes heat conduction in homogeneous matter at normal conditions well. New heat sources, such as lasers, microwaves etc, are employed in medicine, science and material processing for melting, welding, cutting, drilling, etc. Often heat source and treated media are in motion. Some of modern materials and media have one or two dimensions: ultra-thin films, layers and nano-wires. Highly inhomogeneous porous and multilayered media are also common in industry. First major deviations from Fourier law were found in liquid Helium and in some solid crystal dielectrics at low temperatures <25 ° K [2,3,4,5]. Proper phenomenon was called Second Sound [6]; its satisfactory qualitative mathematical description was proposed by Cattaneo and Vernotte [7], who supposed phonon heat transport in addition to Fourier heat diffusion. In terms of temperature it can be put as follows:
( τ t 2 + t ) T = k F 2 T
where kF is the Fourier thermal diffusivity, τ = k F / C 2 is the relaxation time of the heat waves propagation, which relates the moments of the temperature change and of the respective heat flux change. Cattaneo-Vernotte constitution implies a phase lag between the heat flux vector and the temperature gradient; in addition to heat diffusion the temperature perturbation propagates in matter like damped sound-wave at finite speed C = k F / τ . The relaxation time τ is associated with the phonon–phonon interaction time; at normal conditions it is very small: τ ≈ 10−13 s. The Cattaneo Equation (1) is the particular case of the telegrapher’s equation (TE).
( t 2 + ε t ) F ( x , t ) = ( α x 2 + κ ) F ( x , t ) ,   ε ,   α ,   κ = const
which takes its name because it describes the electric signal propagation in long electric lines without radiation [8]. For heat conduction we have τ = 1 / ε , k F = α / ε , μ = κ / ε and μ = 0 in Equation (1); the source term, κ 0 , describes heat exchange with the environment of low excess temperatures. However, precise quantitative description of the Second Sound with Equation (1) was not successful. There are reports on non-Fourier heat transport in highly inhomogeneous matter even at normal conditions [9,10,11,12,13,14,15,16,17], in fuel droplets [18], in biomaterials [19], in energy saving and insulating materials [20], in graphene, nanofibers, carbon nanotubes, silicon wires, etc. [21,22,23,24,25,26,27]. Heat transport in these cases was found close to that described by Guyer and Krumhansl (GK) in [28,29]. In the one-dimensional case in terms of temperature alone it has the following form:
( t 2 + ε t δ t , x , x 3 ) T ( x , t ) = ( α x 2 + κ ) T ( x , t ) ,   α ,   ε ,   δ ,   κ = const .
In its pure form the Guyer–Krumhansl law has κ = 0. The following formulation:
( τ t 2 + t ) T ( x , t ) = ( D b t , x , x 3 + k F x 2 + μ ) T ( x , t ) ,
involves the parameters τ = 1 / ε , μ = κ / ε , k F = α / ε , and D b = δ / ε , where the latter, kb, is the ballistic type heat conductivity. The parameters in the above Equations (3) and (4) have the following dimensions: [ τ ] = s , [ k F ] = m 2 s , [ D b ] = m 2 , [ μ ] = 1 s , [ κ ] = 1 s 2 , [ α ] = m 2 s 2 , [ δ ] = m 2 s , and [ ε ] = 1 s .
Despite GK equation is usually associated with ballistic properties, which manifest when the characteristic system scale L is comparable or less than the mean free path l of the phonons L < l, contradictory opinions on it exist [30,31,32]. In particular, the term “ballistic” is questioned in the context of GK type heat conduction in macroscopic inhomogeneous materials at room temperature [10,11]. Moreover, more complicated system of equations was proposed for the ballistic heat transport in [30], while heat propagation in highly inhomogeneous materials was reported to be close to GK law [9,11].
Solutions to Hyperbolic Heat Conduction Equation (HHCE) can be obtained both analytically and numerically [33,34,35,36,37,38,39,40], although numerical methods seem to be more commonly used [41,42,43,44]. Analytical study gives deeper insight in the problem; operational analytical approach and solutions to HHE were developed in [45,46,47,48,49,50]. This method easily handles also other linear DE of high order and fractional DE [51,52,53,54,55,56,57,58]. Use of the exponential differential operators, such as the heat operator S = e x 2 [59] allows operational solution of GK-type Equation (4) as demonstrated in [47,48,49,50].
In what follows we will obtain some exact solutions to non-Fourier heat transfer in GK-type equation with the substantial derivative, where the speed of the media is constant, and will obtain some particular operational solutions to modified GK-type equation in this case. Moreover, we will study the important case of periodic initial conditions, which occur, for example, in polymer electrolyte fuel cells (PEFCs) [60,61,62,63,64,65,66,67,68], as well as in thin membranes and highly inhomogeneous porous periodic structures, in printed wired heating boards [69] etc. We will analyze the solutions for the heat transport, described by the system of inhomogeneous partial derivative equation, and use the Knudsen number to account for the ballistic conditions in 1- and 2-dimensional structures [70].

2. Ballistic Heat Transport Equations with Substantial Derivative

For ultra-thin films and wires neither pure GK nor Cattaneo laws describe exactly heat transport phenomena, but their combination with Fourier law; it forms the following inhomogeneous system of differential equations for the ballistic θ b ( x , t ) and diffusive component θ d ( x , t ) of the complete dimensionless quasi-temperature θ = θ d + θ b (see [30]):
{ t 2 + 2 t 10 K n b 2 3 x 2 3 K n b 2 t , x , x 3 + 1 } θ b ( x , t ) = 0 ,
( t 2 + K n b 2 K n d 2 t K n b 4 3 K n d 2 x 2 ) θ d ( x , t ) = ( t + K n b 2 K n d 2 ) θ b ( x , t ) ,
where indices b and d stand for “ballistic” and “diffusive”, respectively. The dimensionless quasi-temperature θ can be described as the non-dimensional energy, associated with the internal energy; the quantities θ d and θ b must therefore be understood to quasi-temperatures, defined as a measure of corresponding internal energy components u d and u b , respectively to which they are related by the simple relations θ d = u d / c θ b = u b / c (see [30] for details), where c is the heat capacity.
It should be understood though, that the ballistic heat transport is effectively described by both contributions in θ = θ d + θ b , θ d ( x , t ) and θ b ( x , t ) are distinguished for convenience. It is easy to see that the above Equations (5) and (6) are of GK-type (5) and of telegrapher’s type (6); the inhomogeneous right hand side (r.h.s.) of Equation (6) is given by the solution of (5). In the above equations and in all following equations the space- and time-coordinates are obviously dimensionless as well as the proper coefficients in the equations are. The dimensionless form of equations for the temperature is useful at least because of the space-time temperature distribution is expressed in our work explicitly in terms of x and t. Proper renormalization in SI units or else how is elementary (see, for example, [30] for details) and good, especially in view of heavy notations in the analytical solutions, which we will obtain in what follows.
Low-dimensional system is characterized by the dimensionless Knudsen number K n = l / L , which for the phonon heat transport describes the ratio of the free mean path of the phonons l to the characteristic scale of the system L. Knudsen number usually arises in problems, where the scale of the system and the characteristic scale of the processes in it compare with each other, such as for gas flow in ultra-narrow channels [71], etc. In distinct ballistic case in (5), Kn = 1, we have α = 3.333, δ = 3, ε = 2, κ = −1 in (3) and kF = 5/3, Db = 3/2, τ = 1/2, μ = −1/2 in (4). We assume dimensionless equations here and in what follows. In the weak ballistic case, Kn = 0.1, we get α = 0.03333, δ=0.03, ε = 2, κ = −1 in (3) and kF = 0.01666, kb = 0.015, τ = 1/2, μ = −1/2 in (4). In the distinct ballistic case, Kn = 1, all heat transport terms in GK-type equation contribute more or less equally, in the weak ballistic case, Kn = 0.1, the Cattaneo wave-term prevails. The inhomogeneous system of PDE (5) and (6) for ballistic heat transport thin film was studied numerically in [30]; based on the periodic analytical solutions to GK-type equation [72] some solutions to (5) and (6) were obtained in [73,74].
Above Equations (5) and (6), as well as Equations (2)–(4), describe non-Fourier heat transport in stationary media. In the case when the observer moves relatively to the media with constant speed v , the substantial derivative D / D t = / t + v should be used instead of the time-derivative / t to describe the variation of quantity F along a path. Generally speaking, both the material equations and the energy, mass and momentum conservation equations in case of moving observer should be considered with account for the substantial derivative, resulting in a proper heat conduction equation for the temperature alone. This is relatively simple in the case of Fourier heat diffusion and even for relativistic heat equation [75,76], some more complicated for GK-type equation. While not intending to rigorously derive a set of equations for the heat transport in case of moving observer, accounting for the relativistic heat flux, the energy balance, etc., we will solve the extended form of the known GK-type equation with the simple substantial derivative. Derivation of a physically comprehensive one-dimensional analogue of GK-equation in case of moving observer remains beyond the scope of the present study. The following modified form of GK-type equation arises upon the use of the substantial time derivative:
( t 2 + ε t + 2 v t , x 2 δ t , x , x 3 ) T ( x , t ) = ( ( α v 2 ) x 2 ε v x + δ v x 3 + κ ) T ( x , t ) .
This simply presumes that the temperature as a function of time is recorded by a floating instrument in a flow, such as a weather balloon in meteorology and oceanography, which implies the substantial derivative along the pathline traveled. With this said, we notice that the above Equation (7) differs from GK-type Equation (4) by some additional terms: the second mixed time-space derivative in the left hand side (l.h.s.) and the third- and first-order space-derivatives in the r.h.s. We need to emphasize that this equation is related to a rigid domain meanwhile the observer is moving with speed v . Moreover, the coefficient for the second-order space-derivative in the r.h.s. now contains the speed v , which can compensate α, eliminating the second-order space-derivative, if v = α , or invert the sign of this term, if v > α . Telegrapher’s equation accordingly modifies as follows:
( t 2 + ε t + 2 v t , x 2 ) T ( x , t ) = ( ( α v 2 ) x 2 ε v x + κ ) T ( x , t ) .
Differently from GK-type equation with substantial derivative (7), TE (8) has only second-order differential operators; both GK-type and TE can be solved by the operational method, employing exponential differential operators e D ^ and involving the heat operator S ^ = e x 2 [59]. The operational solution to the third-order PDE (7) apparently involves exponential differential operator of the third order e a x 3 . However, in what follows we will demonstrate as this solution effectively reduces to the shift of the solution in a stationary media. Considering thin films, we have obvious relations between the coefficients in (7) and (5): α b 10 K n b 2 / 3 , ε b 2 , δ 3 K n b 2 , κ b 1 , while for Equation (6) we see that α d K n b 4 / 3 K n d 2 , ε d K n b 2 / K n d 2 . For heat transport in moving thin film, the system of DE (5) and (6) for the complete quasi-temperature θ = θ d + θ b takes the following form:
( t 2 + 2 t + 2 v t , x 2 ( 10 3 K n b 2 v 2 ) x 2 3 K n b 2 t , x , x 3 + 2 v x 3 v K n b 2 x 3 + 1 ) θ b ( x , t ) = 0 ,
( t 2 + K n b 2 K n d 2 t + 2 v t , x 2 ( K n b 4 3 K n d 2 v 2 ) x 2 + v K n b 2 K n d 2 x ) θ d ( x , t ) = ( t + v x + K n b 2 K n d 2 ) θ b ( x , t ) .
In what follows we will approach the GK-type equation with substantial derivative (7) operationally, and will also provide exact analytical harmonic solution to the system of PDE (9) and (10), describing ballistic heat transfer with substantial derivative.

3. Operational Approach to Transport Equations

The above Equations (7) and (8) are the particular cases of the following DE with the coordinate-dependent operators ε ^ ( x ) and D ^ ( x ) :
( t 2 + ε ^ ( x ) t ) F ( x , t ) = D ^ ( x ) F ( x , t ) ,
which can be solved using the operational method. Equation (11) models a broad spectrum of physical phenomena. It becomes telegrapher’s equation for D ^ ( x ) = α x 2 + κ and constant ε term and further reduces to Cattaneo heat equation for D ^ ( x ) = x 2 ; GK-type equation appears when ε ^ = ε δ x 2 . Some of these equations were studied in [45,46,47,48,49,50]. Other second-order PDE and fractional DE were explored with the help of the operational approach in [51,52,53,54,55,56,57,58]. The formal converging particular operational solution to (11) reads as follows:
F ( x , t ) = e t ε ^ ( x ) 2 e t 2 ε ^ 2 ( x ) + 4 D ^ ( x ) C ( x ) ,
where C ( x ) can be obtained from the initial condition F ( x , 0 ) = f ( x ) . The particular form of the initial function will be chosen below. The other branch of the solution contains the positive argument in the exponential, e t 2 ε ^ 2 ( x ) + 4 D ^ ( x ) . No Laplace transforms exist for it. However, symmetry with respect to inversion, t t ,   ε ε ,   δ δ ,   v v , allows writing the other solution to Equation (7), based on the one we obtain with the help of the Laplace transforms
e t V = t 2 π 0 d ξ ξ ξ e t 2 4 ξ ξ V ,   t > 0 .
The solution has the following integral form, provided it converges:
F ( x , t ) = e t 2 ε ^ ( x ) t 4 π 0 d ξ ξ ξ e t 2 16 ξ e ξ ε ^ 2 ( x ) e 4 ξ D ^ ( x ) f ( x ) .
The ability to perform analytical integration in (14) depends on the explicit form of the operators ε ^ , D ^ and of the initial function f ( x ) ; the numerical calculation can be done though. Accounting for the explicit form of the operators D ^ = ( α v 2 ) x 2 ε v x + δ v x 3 + κ and ε ^ = ε + 2 v x δ x 2 in Equation (7), we obtain from (14) the following integral:
F ( x , t ) = e t 2 ( ε + 2 v x ) t 4 π 0 d ξ ξ ξ e t 2 16 ξ ξ ( ε 2 + 4 κ ) e ( t 2 δ + 2 ξ ( ε δ 2 α ) ) x 2 e ξ δ 2 x 4 f ( x ) ,
which benefits from the use of the operational identity for p ^ = a D ^ , D ^ = 2 (see [59,77]):
e p ^ 2 = 1 π exp ( ξ 2 + 2 ξ p ^ ) d ξ ,
and yields
e a D ^ 2 f ( x ) = 1 π exp ( ξ 2 + 2 ξ a D ^ ) f ( x ) d ξ .
Applying formula (17) to the fourth-order exponential differential operator e ξ δ 2 x 4 f ( x ) = e ζ 2 + 2 i ζ δ ξ x 2 f ( x ) d ζ / π , we get the heat operator [59] S ^ = e ν x 2 ; collecting the second-order derivative terms in the exponential, we get e ζ 2 + ( ( t δ / 2 ) 4 ξ α + 2 ξ ε δ + 2 i ζ δ ξ ) x 2 f ( x ) d ζ . This yields the following particular bounded solution to GK-type heat equation with substantial derivative (7):
F ( x , t ) = e t 2 ε t 4 π Θ ^ 0 d ξ ξ ξ e t 2 16 ξ ξ ( ε 2 + 4 κ ) e ζ 2 S ^ f ( x ) d ζ ,
which involves the shift operator Θ ^ = e y x , where y = v t , and the heat operator S ^ = e η x 2 , where η = 2 i ζ ξ δ ξ 2 ( 2 α + ε δ ) + t δ / 2 . The shift operator produces the translation along x e y ( x + α ) f ( x ) = e y α f ( x + y ) and the heat operator produces Gauss transforms (17) and operational relations [59], so that the solution with account for the motion contains the translation x v t of the stationary solution.
Let us consider the initial polynomial function f ( x ) = k x k . The particular solution (18) to GK-type equation with substantial derivative (7) arises upon the application of the operational rule
S ^ e γ x x k = e η x 2 x k e γ x = e γ x + γ 2 ν H k ( x + 2 γ ν , ν ) ,   k Integers ,   γ Reals ,
which is easy to prove [57,58]. For γ = 0 we immediately obtain the following integral:
F ( x , t ) | f ( x ) = k x k = k t e t 2 ε 4 π 0 d ξ ξ ξ e t 2 16 ξ ξ ( ε 2 + 4 κ ) e ζ 2 H k ( x v t , η ) d ζ ,
where
η = 2 ξ ε δ 4 ξ α + t δ / 2 + i ζ 2 ξ δ .
The above integral can be taken in elementary functions if we account for the explicit form of the Hermite polynomials:
H k ( x v t , η ) = k ! r = 0 [ k / 2 ] ( x v t ) k 2 r η r ( k 2 r ) ! r ! = ( i ) k η k 2 H k ( i ( x v t ) 2 η ) .
For example, for f ( x ) = x 2 the solution reads as follows:
F ( x , t ) | f ( x ) = x 2 = e t 2 ( V + ε ) ( ( x v t ) 2 + t δ + t V ( δ ε 2 α ) ) ,   V = ε 2 + 4 κ 2 .
More general case of the exponential-polynomial function f ( x ) = x k e γ x is cumbersome and we omit proper expressions for conciseness. However, we have performed the integration explicitly in Wolfram Mathematica program and below we give the example of the solution for F ( x , 0 ) = x 2 e x :
F ( x , t ) | f ( x ) = x 2 e x = e x t 2 ( r 2 v δ + ε ) ( x 2 + t r 3 / 2 ( a + b x ) + t 2 r 3 / 2 c ) , r = 4 α + δ 2 2 δ ε + ε 2 + 4 κ , a = r δ ( 4 α + ( ε δ ) 2 ) + δ 2 ( 6 α δ 2 3 ε ( ε δ ) ) + ε 2 ( 2 α + δ ε ) 4 κ ( 2 α + r δ 3 δ 2 + δ ε ) , b = 2 r ( 2 α + δ 2 r v δ ( r + ε ) ) , c = r 3 / 2 v 2 2 r v d + 2 p + δ 2 Δ d = 2 α + δ ( r + δ ε ) , p = 2 α 2 ( r 4 δ ) 2 α δ ( 3 δ 2 + r ε + ε 2 2 δ ( r + 2 ε ) + 4 κ ) , Δ = r ( ε δ ) 2 + ( ε δ ) 3 + 2 κ ( r + 2 ( ε δ ) ) .
Evidently, all the values in the above expressions are dimensionless as well as the coordinate x and the time t.
In the harmonic ansatz the evolution of the initial function f ( x ) = exp ( i n x ) can be easily obtained from (18) as follows:
F ( x , t ) | f ( x ) = e i n x = exp ( i n ( x v t ) t 2 ( ε ¯ + ε ¯ 2 + 4 ( κ 2 α n 2 ) ) ) ,   ε ¯ = ε + n 2 δ .
Other solutions and their detailed study will be performed in forthcoming publications. All the above solutions to equations with substantial derivative contain the shift, x v t , respectively to the solution in a stationary media, which depends on x. Further examples can be easily considered with the help of the operational approach. In the following chapter we will explore in details the particular case of the harmonic solution. Evidently, it can be easily generalized for any periodic solution, expandable in Fourier series. Similarly, the Fourier integral transforms technique apply.

4. Exact Periodic Solutions to GK-Type Equation with Substantial Derivative

An exact harmonic solution e i n x for the inhomogeneous system of PDE (9) and (10) can be obtained straight from the operational solution (18) (see (25)), or by separating the variables: T ( x , t ) = X ( x ) y ( t ) . GK-type equation with substantial derivative (7) in the harmonic ansatz, T ( x , t ) = e i n x y ( t ) F ( x , t ) reduces to the following ordinary differential equation for y ( t ) :
y ( t ) + ε e f f y ( t ) + λ e f f y ( t ) = 0 ,
where ε e f f = ε ¯ + 2 v i n , ε ¯ = ε + n 2 δ , and λ e f f = ( α v 2 ) n 2 + i v n ε ¯ κ . The function F ( x , t ) T ( x , t ) is introduced here for clarity of notations. The exact solution to GK equation with substantial derivative (7) then easily follows from the solution to (26):
F ( x , t ) = e i n x y ( t ) ,   y ( t ) = e t 2 ε e f f ( C 1 e t 2 U + C 2 e t 2 U ) ,   U = ε e f f 2 4 λ e f f = ε ¯ 2 + 4 ( κ α n 2 ) ,
where the constants C 1 , C 2 are determined either from the initial or boundary conditions. According to the theory of separation of variables, the whole solution has to satisfy initial or boundary conditions. Evidently, telegrapher’s equation with substantial derivative in the harmonic ansatz similarly reads
y ( t ) + ε ˜ y ( t ) + λ ˜ y ( t ) = 0 ,
and differs from (26) by the substitutions ε e f f ε ˜ = ε + 2 v i n , ε ¯ ε , λ e f f λ ˜ = ( α v 2 ) n 2 + i v n ε κ . It’s solution y ( t ) is evidently (27), where δ = 0.
Consider, for example, the Cauchy initial conditions
F ( x , 0 ) = A e i n x ,   t F ( x , 0 ) = B e i n x .
In the stationary case, v = 0, the Cauchy conditions (37), where initially θ ( x , t ) t = 0 , have the meaning of zero heat flux q if B = 0 for both ballistic and diffusive components of quasi-temperature [30]. If v ≠ 0, then the dimensionless equation θ t + v θ x + q = 0 describes the one-dimensional energy balance. The initially zero heat flux q then corresponds to B = i n v A .
Then, GK-type Equation (7) with substantial derivative in the harmonic ansatz has the solution (27), satisfying ODE (26), where the coefficients C 1 , C 2 are determined from the Cauchy conditions (29): C 1 + C 2 = A and C 1 ( ε e f f + U ) C 2 ( ε e f f U ) = 2 B , and read as follows:
C 1 = A 2 B + A ε e f f / 2 U ,   C 2 = A 2 + B + A ε e f f / 2 U .
TE with substantial derivative (8) has the harmonic solution, whose time-dependent part y ( t ) satisfies ODE (28):
F ( x , t ) = e i n x y ( t ) ,   y ( t ) = e t 2 ε ˜ ( B 1 e t 2 u + B 2 e t 2 u ) ,   u = ε ˜ 2 4 λ ˜ = ε 2 + 4 ( κ α n 2 ) ,
where the coefficients B 1 , B 2 are obtained from the Cauchy initial conditions (29):
B 1 = A 2 B + A ε ˜ / 2 u ,   B 2 = A 2 + B + A ε ˜ / 2 u .
Now, let us consider the Dirichlet boundary conditions in the moments of time t = 0 and t = T :
F ( x , 0 ) = A e i n x ,   F ( x , T ) = G e i n x .
The harmonic solutions (27) and (31) satisfy respectively GK Equation (7) and telegrapher’s Equation (8) with substantial derivative with the coefficients, determined by Equation (33):
C 1 = e T 2 ε e f f A e T 2 ( ε e f f U ) G e T 2 U e T 2 U ,   C 2 = A C 1 ,
and
B 1 = e T 2 ε ˜ A e T 2 ( ε ˜ u ) G e T 2 u e T 2 u ,   B 2 = A B 1 .
Similarly, we can consider evolution of any function, expandable in Fourier series.
Note, that for some values of the coefficients α and δ the quantities U and u can assume negative values. Albeit it is not very obvious in the form of the harmonic solution (27) to GK-type equation with substantial derivative (26) and the solution (31) to TE with substantial derivative (28), these solutions remains real at any moment of time for a real initial function, because of the complex exponential is compensated by proper complex parts in the coefficients C 1 , 2 and B 1 , 2 .

5. Exact Periodic Solutions to Ballistic Heat Transport in Thin Films

Ballistic heat transport in thin films is described by the system of PDE (9) and (10). In order to obtain the exact solution to this problem we reduced GK-type Equation (9) in the harmonic ansatz to (26) and solved it (see (27)); the result naturally involves the particular solution (25). This solution now constitutes the r.h.s. of Equation (10) for θ d . Consider the Cauchy initial problem for θ b :
θ b ( x , 0 ) = A e i n x ,   t θ b ( x , 0 ) = B e i n x ,
and the Cauchy conditions for θ d :
θ d ( x , 0 ) = V e i n x ,   t θ d ( x , 0 ) = W e i n x ,   V , W = const .
The quasi-temperature component θ b ( x , t ) is in fact given by Equation (27):
θ b ( x , t ) = e i n x ( C 1 e t 2 E 1 + C 2 e t 2 E 2 ) ,   E 1 = ( ε e f f + U ) ,   E 2 = ( ε e f f U ) ,
where ε e f f = ε ¯ + 2 i v n , U = ε ¯ 2 + 4 ( κ α b n 2 ) , ε ¯ = ε b + n 2 δ , α b = 10 K n b 2 / 3 , ε b = 2 , δ = 3 K n b 2 , κ b = 1 , v is the speed of the media; for the Cauchy initial conditions (36) (see (29) with F θ b ) the coefficients C 1 , 2 are given by Equation (30). The solution (38) for θ b ( x , t ) contributes to the r.h.s. of the telegrapher’s Equation (6) for the component θ d ( x , t ) . The exponential differential operators e D ^ ( x ) do not bring new harmonics to the initial content, θ d ( x , t ) = Θ d ( t ) e i n x , and Equation (6) reduces to the following inhomogeneous ODE for Θ d ( t ) :
( d 2 d t 2 + ε ˜ d d d t + λ ˜ d ) Θ d ( t ) = + P e t 2 E 1 + Q e t 2 E 2 ,
where
ε ˜ d = ε d + i 2 v n ,   λ ˜ d = ( α d v 2 ) n 2 + i v n ε d ,   ε d = K n b 2 K n d 2 ,   α d = K n b 4 3 K n d 2 ,   P = C 1 ( ε d ε ¯ + U 2 + i v n ) ,   Q = C 2 ( ε d ε ¯ U 2 + i v n ) ,
where C 1 , 2 are given by Equation (30). The above Equation (39) for Θ d ( t ) possesses analytical solution, which, in turn, yields the evolving in time harmonic solution for the diffusive quasi-temperature component θ d ( x , t ) in the following form:
θ d ( x , t ) = e i n x ( 4 ( P S e t 2 E 1 + Q L e t 2 E 2 ) + D 1 e t 2 ( ε ˜ d + R ) + D 2 e t 2 ( ε ˜ d R ) ) ,
R = ε ˜ d 2 4 λ ˜ d = ε d 2 4 n 2 α d , S = E 1 2 + 4 λ ˜ d 2 E 1 ε ˜ d = ( ε ¯ ε d + U ) 2 + 4 α d n 2 ε d 2 , L = E 2 2 + 4 λ ˜ d 2 E 2 ε ˜ d = ( ε ¯ ε d U ) 2 + 4 α d n 2 ε d 2 ,
where the constants D 1 , 2 are determined from the Cauchy initial conditions (37) for θ d ( x , t ) :
D 1 = V 2 ( 1 ε ˜ d R ) 1 R ( W + 2 P S ( ε ¯ ε d + U + R ) + 2 Q L ( ε ¯ ε d U + R ) ) ,
D 2 = V 2 ( 1 + ε ˜ d R ) + 1 R ( W + 2 P S ( ε ¯ ε d + U R ) + 2 Q L ( ε ¯ ε d U R ) ) .
The exact analytical harmonic solutions for the θ b ( x , t ) is given by (38) and (30), and for θ d ( x , t ) , it is given by (41)–(44); their sum represents the complete solution for the ballistic heat transfer in thin films [30] with the Cauchy initial conditions (36), (37), valid for arbitrary values of the Knudsen numbers K n d , K n b . Note, that a real initial distribution of quasi-temperature θ ( x , 0 ) evolves in real domain even for imaginary exponentials in the solutions due to complex coefficients C 1 , C 2 and D 1 , D 2 .
Now let us consider the Dirichlet conditions at the boundaries
θ b ( x , 0 ) = A e i n x ,   θ b ( x , T ) = G e i n x ,   θ d ( x , 0 ) = V e i n x ,   θ d ( x , T ) = W e i n x ,
where A , G , V , W = const (see also Equation (33)). The evolution of θ b ( x , t ) is given by Equation (38), where C 1 , 2 are set by Equation (34). The evolution of the counterpart θ d ( x , t ) is governed by Equations (39) and (40). The r.h.s. of Equation (39) is set by Equation (38) for the ballistic component θ b ( x , t ) , where C 1 , 2 are given by (34). The explicit solution θ d ( x , t ) to inhomogeneous DE (39) is possible, although it is much more cumbersome than in the Cauchy case:
θ d ( x , t ) = e i n x ( b 2 w 1 e E 2 t + b 1 w 2 e E 1 t + q 1 w 1 w 2 e 1 2 t ( ε d + R ) q 2 w 1 w 2 e 1 2 t ( ε d R ) ) ,
where
w 1 = E 2 2 E 2 ε ˜ d + λ ˜ d ,   w 2 = E 1 2 E 1 ε ˜ d + λ ˜ d , q 1 = 1 e R T 1 ( e E 2 T w 1 b 1 d 2 + e E 2 T b 2 d 2 w 2 + d 1 w 1 w 2 ) , q 2 = 1 e R T 1 ( e E 2 T ( λ ˜ d w 2 ) ( b 2 d 6 w 1 d 4 ) + e E 1 T ( w 1 ( λ ˜ d d 3 b 1 d 5 ) + λ ˜ d b 2 d 7 ) ) , b 1 = C 1 ( ε ˜ d ε e f f + U 2 + i v n ) ,   b 2 = C 2 ( ε ˜ d ε e f f U 2 + i v n ) ,   d 1 = V e R T W e T 2 ( ε ˜ d + R ) , d 2 = e T 2 ( ε ˜ d + R ) e T ( E 1 + R ) ,   d 3 = V e E 1 T W e ( 2 E 1 + ε ˜ d + R ) T 2 ,   d 4 = V e E 2 T W e ( 2 E 2 + ε ˜ d + R ) T 2 , d 5 = e E 1 T e ( ε ˜ d + R ) T 2 ,   d 6 = e E 2 T e ( ε ˜ d + R ) T 2 ,   d 7 = e E 1 T + e ( ε ˜ d + R + 2 ( E 1 E 2 ) ) T 2 .
The solution for θ d ( x , t ) in the form (46), (47) explicitly involves the Dirichlet conditions and together with the solution for θ b ( x , t ) , given in the Dirichlet case by (38) and (34), they describe the evolution of the complete harmonic quasi-temperature θ = θ d + θ b .
Let us now consider some examples of the harmonic solutions to GK-type heat equation with substantial derivative. Consider, for example, the Cauchy problem with initial conditions θ b ( x , 0 ) = θ d ( x , 0 ) = e i x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 . Let’s first consider K n b = K n d = 1 . Proper solutions for θ b ( x , t ) , θ d ( x , t ) and θ ( x , t ) in the stationary media are presented in Figure 1.
It can be seen in the top left plot in Figure 1 as the ballistic component θ b ( x , t ) rapidly relaxes to the stationary state. The diffusive component θ d ( x , t ) behaves similarly, but the relaxation process takes much more time. The behavior of ballistic and diffusive components in the domain is more distinct for K n < 1 . The complete quasi-temperature θ ( x , t ) (see bottom plot in Figure 1) follows the behavior of both the ballistic θ b ( x , t ) and diffusive θ d ( x , t ) constituents. However, the latter, being the solution of Equation (39), is strongly influenced by the solution (38) for the ballistic component θ b ( x , t ) , which sets the inhomogeneous part in Equation (39). Thus, the evolution of the quasi-temperature θ ( x , t ) in thin films has two speeds: the ballistic constituent θ b ( x , t ) develops rapidly and sets the r.h.s. in Equation (39) for the diffusive constituent θ d ( x , t ) . The evolution of θ d ( x , t ) is much slower: after the first instances it proceeds under the influence of already relaxed true ballistic constituent θ b ( x , t ) .
In the case when the rate of change of temperature for a given position in a field depends both on the instantaneous rate of change of temperature at that location (∂/∂t) as well as on the rate at which the temperature is convected to that location by the fluid motion, the behaviors of the solutions change significantly. The solutions for Knd = Knb = 1 and v = 10, are presented in Figure 2, for v = −10 in Figure 3.
The ballistic constituent θ b rapidly vanishes with time; the direction of the colorful waves in the plot depends on the sign of the speed v (see top left plots in Figure 2 and Figure 3). Diffusive counterpart of the complete quasi-temperature θ d has waves of increasing amplitude as seen in top right plots in Figure 2 and Figure 3. The complete quasi-temperature θ = θ d + θ b has obvious non-Fourier behavior, which largely follows that of θ d .
For smaller values of Knudsen number, K n = 0.1 , the behavior of the harmonic solutions to ballistic heat propagation in thin films, obeying the Eqs. system (9) and (10) with Cauchy conditions (36) (37), has less distinctive wave interference with slower fade of the solution for θ b and weaker amplitude increase speed for θ d ; proper solutions for v = +10 are presented in Figure 4.
It is important to note, that the whole system is symmetric under the speed inversion, which becomes obvious with account for formulae (38)–(44).
Note, that for Knb = Knd = 0.1 θ b ( x , t ) and θ d ( x , t ) constituents behave differently (see Figure 4). The diffusive constituent slowly grows as the ballistic constituent shows gradual amplitude decrease, as shown in Figure 5. Moreover, for v = 10 the complete quasi-temperature θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) monotonously increases, following the behavior of the diffusive constituent as shown in Figure 5. The maximum principle, established though for parabolic equations, is violated for the diffusive component and complete quasi-temperature with Cauchy conditions for non-Fourier ballistic heat transport in thin films. In all of the cases, the evolution of θ b ( x , t ) occurs faster than that of θ d ( x , t ) ; the two-speed heat propagation process can be seen in every set of plots in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5.
The exploration of the solutions for Dirichlet conditions will be done elsewhere.
It should be noted that the one-dimensional energy balance equation yields zero value for the initial heat flux, if B = i n v A . Above we have obtained the solutions, which allow nonzero initial conditions on the first time derivative; we have investigated some examples while writing this paper. For example, for the case of K n b = K n d = 1 , v = 10 the maximum principle is violated for both zero and nonzero initial conditions imposed on the first time-derivative of the ballistic component; however, there was no significant qualitative difference in the behaviors of the diffusive component in this case.
Explicit study for the solution in thin films with zero initial heat flux and non-zero velocity v will be done in forthcoming publications.

6. Conclusions

In the present work we have studied and analytically solved one-dimensional heat transport equations of Cattaneo- and Guyer–Krumhansl-type with substantial derivative. Exponential differential operators and the operational method were used to obtain integral forms of exact particular solutions to heat transport equations. The media speed is accounted for by the exponential differential operator e v t x , which produces a shift x v t . Together with the heat operator S ^ they transform the initial temperature profile and determine its evolution. Their action on the initial polynomial n x n yields sums of the Hermite polynomials H n ( x , y ) ; the action on the monomial-exponential function yields more complicated sums of Hermite polynomials. Proper exact operational solutions were obtained in the sums of converging integrals of elementary functions. Several examples were considered, the solutions for F ( x , 0 ) = x 2 , F ( x , 0 ) = x 2 e x and F ( x , 0 ) = e i n x were given. The solutions exactly satisfy the heat transport equations, which has been proven by direct substitution.
The equations with substantial derivative of Guyer-Krumhansl-type and of telegrapher’s-type are demonstrated to have similar structure in the harmonic ansatz. Their exact harmonic solutions have been obtained by both operational method and the separation of variables. The latter reduces these PDEs to ODEs. We have shown that the harmonic solutions may have local and global extremums in the domain due to the interference of the Cattaneo heat waves; proper second-order time-derivative term also insures finite speed of the temperature perturbation propagation. We have shown that in a GK-type equation with dominant ballistic term δ, the Cattaneo heat waves are damped.
We have obtained the exact analytical solution to the ballistic heat transport in thin films by solving the system of inhomogeneous PDEs [30] in the harmonic ansatz; both Cauchy and Dirichlet conditions were explored. We found that true ballistic constituent of quasi-temperature, θ b , evolves much faster than the diffusive counterpart of quasi-temperature, θ d . We have studied the evolution of the complete quasi-temperature θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) under the ballistic conditions; the latter apply when the phonon mean free path l is comparable with the system scale L. The quasi-temperature θ first senses fast evolving true ballistic part θ b . The contribution of θ b sets the r.h.s. in the inhomogeneous DE for the component θ d . After relatively fast relaxation of θ b , further evolution of the complete quasi-temperature θ follows the diffusive constituent θ d . Thus the initial fast-developing true ballistic component θ b determines the evolution of the complete quasi-temperature of the system θ by imposing distribution for the diffusive quasi-temperature θ d at the beginning of the process.
The obtained harmonic solutions for the Cauchy problem (9), (10), (36), (37) were applied for the study of the temperature patterns in thin films for the Knudsen number values Kn = 0.1 and Kn = 1. In the stationary case v = 0, the true ballistic component θ b monotonously relaxes (see Figure 1 and Figure 5); its diffusive counterpart θ d behavior shows noticeable dependence on the value of the Knudsen numbers: for Kn = 1 it slowly relaxes (see Figure 5) and for Kn = 0.1 it droningly grows (see Figure 1). The effect of the speed of the observer on the perceived heat transport in thin films is distinguished for the pure ballistic and diffusive quasi-temperature constituents. In the presence of the motion, the extremum of the ballistic component θ b occurs inside the domain after the initial moment t = 0 (see Figure 2, Figure 3 and Figure 4). The diffusive quasi-temperature component monotonously grows in the domain both in the weak ballistic (see Figure 4) and strong ballistic cases (see Figure 2 and Figure 3). This behavior contrasts common Fourier heat diffusion, for which the maximum principle holds. It is also noticeable that the system shows symmetric behavior while the value of the speed is reversed, which does not contradict the Second law. The Knudsen number variation from 0.1 to 1 in the presence of the speed v = 10, are practically not sensed by the system. Thus the speed of the observer v influences the registered temperature rather than the ballistic heat transport contribution.

Author Contributions

K.Z. has made the major work on the paper. N.G. has made a contribution in the Section 3. D.O. has made the contribution in Section 4 and Section 5.

Acknowledgments

The authors are grateful to Anatolii Borisov, Andrey Lobanov for the discussions of the heat equations; to Laszlo I. Kiss (UQAC, Department of Applied Sciences, Quebec University) and Yasar Demirel (Department of Chemical and Biological Engineering, University of Nebraska, Lincoln) for discussions of the heat evolution patterns.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fourier, J.P.J. The Analytical Theory of Heat; Cambridge University Press: London, UK, 1878. [Google Scholar]
  2. Ackerman, C.C.; Guyer, R.A. Temperature pulses in dielectric solids. Ann. Phys. 1968, 50, 128–185. [Google Scholar] [CrossRef]
  3. Ackerman, C.C.; Overton, W.C. Second sound in solid helium-3. Phys. Rev. Lett. 1969, 22, 764. [Google Scholar] [CrossRef]
  4. McNelly, T.F.; Rogers, S.J.; Channin, D.J.; Rollefson, R.; Goubau, W.M.; Schmidt, G.E.; Krumhansl, J.A.; Pohl, R.O. Heat pulses in NaF: Onset of second sound. Phys. Rev. Lett. 1970, 24, 100. [Google Scholar] [CrossRef]
  5. Narayanamurti, V.; Dynes, R.D. Observation of second sound in Bismuth. Phys. Rev. Lett. 1972, 26, 1461–1465. [Google Scholar] [CrossRef]
  6. Peshkov, V. Second sound in Helium II. J. Phys. (Moscow) 1944, 8, 381. [Google Scholar]
  7. Cattaneo, C. Sur une forme de l’equation de la chaleur eliminant le paradoxe d’une propagation instantanee. C. R. l’Acad. Sci. Paris 1958, 247, 431–433. [Google Scholar]
  8. Terman, F.E. Radio Engineers’ Handbook, 1st ed.; McGraw-Hill: New York, NY, USA, 1943. [Google Scholar]
  9. Both, S.; Czél, B.; Fülöp, T.; Ván, P.; Verhás, J. Deviation from the Fourier law in room-temperature heat pulse experiments. J. Non-Equilib. Thermodyn. 2016, 41, 41–48. [Google Scholar] [CrossRef]
  10. Kovács, R.; Ván, P. Models of ballistic propagation of heat at low temperatures. Int. J. Thermophys. 2016, 37, 95. [Google Scholar] [CrossRef]
  11. Ván, P.; Berezovski, A.; Fülöp, T.; Gróf, G.; Kovács, R.; Lovas, Á.; Verhás, J. Guyer-Krumhansl-type heat conduction at room temperature. EPL 2017, 118, 50005. [Google Scholar] [CrossRef]
  12. Cahill, D.G. Thermal conductivity measurement from 30 to 750 K: The 3ω method. Rev. Sci. Instrum. 1990, 61, 802–808. [Google Scholar] [CrossRef]
  13. Tang, D.W.; Araki, N. Non-Fourier heat conduction behaviour in finite mediums under pulse surface heating. Mater. Sci. Eng. A 2000, 292, 173–178. [Google Scholar] [CrossRef]
  14. Kaminski, W. Hyperbolic heat conduction equations for materials with a nonhomogeneous inner structure. J. Heat Transf. 1990, 112, 555–560. [Google Scholar] [CrossRef]
  15. Mitra, K.; Kumar, S.; Vedavarz, A.; Moallemi, M.K. Experimental evidence of hyperbolic heat conduction in processed meat. J. Heat Transf. 1995, 117, 568–573. [Google Scholar] [CrossRef]
  16. Herwig, H.; Beckert, K. Fourier versus non-Fourier heat conduction in materials with a nonhomogeneous inner structure. J. Heat Transf. 2000, 122, 363–365. [Google Scholar] [CrossRef]
  17. Roetzel, W.; Putra, N.; Das, S.K. Experiment and analysis for non-Fourier conduction in materials with non-homogeneous inner structure. Int. J. Therm. Sci. 2003, 42, 541–552. [Google Scholar] [CrossRef]
  18. Sazhin, S.S. Modelling of fuel droplet heating and evaporation: Recent results and unsolved problems. Fuel 2017, 196, 69–101. [Google Scholar] [CrossRef]
  19. Scott, E.P.; Tilahun, M.; Vick, B. The question of thermal waves in heterogeneous and biological materials. J. Biomech. Eng. 2009, 131, 074518. [Google Scholar] [CrossRef] [PubMed]
  20. Ricciu, R.; Besalduch, L.A.; Galatioto, A.; Ciulla, G. Thermal characterization of insulating materials. Renew. Sustain. Energy Rev. 2018, 82, 1765–1773. [Google Scholar] [CrossRef]
  21. Porrà, J.M.; Masoliver, J.; Weiss, G.H. When the telegrapher’s equation furnishes a better approximation to the transport equation than the diffusion approximation. Phys. Rev. E 1997, 55, 7771–7774. [Google Scholar] [CrossRef]
  22. Shiomi, J.; Maruyama, S. Non-Fourier heat conduction in a single-walled carbon nanotube: Classical molecular dynamics simulations. Phys. Rev. B 2006, 73, 205420. [Google Scholar] [CrossRef]
  23. Baringhaus, J.; Ruan, M.; Edler, F.; Tejeda, A.; Sicot, M.; Taleb-Ibrahimi, A.; Li, A.P.; Jiang, Z.; Conrad, E.H.; Berger, C.; et al. Exceptional ballistic transport in epitaxial graphene nanoribbons. Nature 2014, 506, 349–354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hochbaum, A.I.; Chen, R.; Delgado, R.D.; Liang, W.; Garnett, E.C.; Najarian, M.; Majumdar, A.; Yang, P. Enhanced thermoelectric performance of rough silicon nanowires. Nature (London) 2008, 451, 163–167. [Google Scholar] [CrossRef] [PubMed]
  25. Boukai, A.I.; Bunimovich, Y.; Tahir-Kheli, J.; Yu, J.-K.; Goddard, W.A.; Heath, J.R. Silicon nanowires as efficient thermoelectric materials. Nature (London) 2008, 451, 168. [Google Scholar] [CrossRef] [PubMed]
  26. Paddock, C.A.; Eesley, G.L. Transient thermoreflectance from thin metal films. J. Appl. Phys. 1986, 60, 285–290. [Google Scholar] [CrossRef]
  27. Maldovan, M. Transition between ballistic and diffusive heat transport regimes in silicon materials. Appl. Phys. Lett. 2012, 101, 113110. [Google Scholar] [CrossRef]
  28. Guyer, R.A.; Krumhansl, J.A. Solution of the linearized phonon Boltzmann equation. Phys. Rev. 1966, 148, 766–778. [Google Scholar] [CrossRef]
  29. Guyer, R.A.; Krumhansl, J.A. Thermal conductivity, second sound and phonon hydrodynamic phenomena in non-metallic crystals. Phys. Rev. 1966, 148, 778–788. [Google Scholar] [CrossRef]
  30. Lebon, G.; Machrafi, H.; Gremela, M.; Dubois, C. An extended thermodynamic model of transient heat conduction at sub-continuum scales. Proc. R. Soc. A 2011, 467, 3241–3256. [Google Scholar] [CrossRef] [Green Version]
  31. Müller, I.; Ruggeri, T. Rational Extended Thermodynamics; Springer: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  32. Rogolino, P.; Kovács, R.; Ván, P.; Cimmelli, V.A. Generalized heat-transport equations: Parabolic and hyperbolic models. Contin. Mech. Thermodyn. 2018. [Google Scholar] [CrossRef]
  33. Moosaie, A. Non-Fourier heat conduction in a finite medium with insulated boundaries and arbitrary initial conditions. Int. Commun. Heat Mass Transf. 2008, 35, 103–111. [Google Scholar] [CrossRef]
  34. Ahmadikia1, H.; Rismanian, M. Analytical solution of non-Fourier heat conduction problem on a fin under periodic boundary conditions. J. Mech. Sci. Technol. 2011, 25, 2919–2926. [Google Scholar] [CrossRef]
  35. Yen, C.C.; Wu, C.Y. Modelling hyperbolic heat conduction in a finite medium with periodic thermal disturbance and surface radiation. Appl. Math. Model. 2003, 27, 397–408. [Google Scholar] [CrossRef]
  36. Lewandowska, M. Hyperbolic heat conduction in the semi-infinite body with a time-dependent laser heat source. Heat Mass Transf. 2001, 37, 333–342. [Google Scholar] [CrossRef]
  37. Lewandowska, M.; Malinowski, L. An analytical solution of the hyperbolic heat conduction equation for the case of a finite medium symmetrically heated on both sides. Int. Commun. Heat Mass Transf. 2006, 33, 61–69. [Google Scholar] [CrossRef]
  38. Saedodin, S.; Torabi, M. Analytical solution of non-Fourier heat conduction in cylindrical coordinates. Int. Rev. Mech. Eng. 2009, 3, 726–732. [Google Scholar]
  39. Challamel, N.; Grazide, C.; Picandet, V.; Perrot, A.; Zhang, Y. A nonlocal Fourier’s law and its application to the heat conduction of one-dimensional and two-dimensional thermal lattices. C. R. Mec. 2016, 344, 388–401. [Google Scholar] [CrossRef]
  40. Saedodin, S.; Torabi, M. Algebraically explicit analytical solution of three-dimensional hyperbolic heat conduction equation. Adv. Theor. Appl. Mech. 2010, 3, 369–383. [Google Scholar]
  41. Chen, G. Ballistic-diffusive heat-conduction equations. Phys. Rev. Lett. 2001, 86, 2297–2300. [Google Scholar] [CrossRef] [PubMed]
  42. Hsiao, T.K.; Chang, H.K.; Liou, S.C.; Chu, M.W.; Lee, S.C.; Chang, C.W. Observation of room-temperature ballistic thermal conduction persisting over 8.3 mm in SiGe nanowires. Nat. Nanotechnol. 2013, 8, 534–538. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, Y.; Ye, W. Modified ballistic–diffusive equations for transient non-continuum heat conduction. Int. J. Heat Mass Transf. 2015, 83, 51–63. [Google Scholar] [CrossRef]
  44. Kovacs, R.; Van, P. Generalized heat conduction in heat pulse experiments. Int. J. Heat Mass Transf. 2015, 83, 613–620. [Google Scholar] [CrossRef] [Green Version]
  45. KZhukovsky, V. Operational method of solution of linear non-integer ordinary and partial differential equations. SpringerPlus 2016, 5, 119. [Google Scholar] [CrossRef] [PubMed]
  46. Zhukovsky, K. Operational approach and solutions of hyperbolic heat conduction equations. Axioms 2016, 5, 28. [Google Scholar] [CrossRef]
  47. Zhukovsky, K.V.; Srivastava, H.M. Analytical solutions for heat diffusion beyond Fourier law. Appl. Math. Comput. 2017, 293, 423–437. [Google Scholar] [CrossRef]
  48. Zhukovsky, K.V. Violation of the maximum principle and negative solutions with pulse propagation in Guyer–Krumhansl model. Int. J. Heat Mass Transf. 2016, 98, 523–529. [Google Scholar] [CrossRef]
  49. Zhukovsky, K.V. Exact solution of Guyer–Krumhansl type heat equation by operational method. Int. J. Heat Mass Transf. 2016, 96, 132–144. [Google Scholar] [CrossRef]
  50. Zhukovsky, K. Exact negative solutions for Guyer–Krumhansl type equation and the violation of the maximum principle. Entropy 2017, 19, 440. [Google Scholar] [CrossRef]
  51. Zhukovsky, K. Operational solution for some types of second order differential equations and for relevant physical problems. J. Math. Anal. Appl. 2017, 446, 628–647. [Google Scholar] [CrossRef]
  52. Zhukovsky, K.V. A method of inverse differential operators using ortogonal polynomials and special functions for solving some types of differential equations and physical problems. Mosc. Univ. Phys. Bull. 2015, 70, 93–100. [Google Scholar] [CrossRef]
  53. Zhukovsky, K. Solution of some types of differential equations: Operational calculus and inverse differential operators. Sci. World J. 2014, 2014. [Google Scholar] [CrossRef] [PubMed]
  54. Zhukovsky, K.V. Solving evolutionary-type differential equations and physical problems using the operator method. Theor. Math. Phys. 2017, 190, 52–68. [Google Scholar] [CrossRef]
  55. Dattoli, G.; Srivastava, H.M.; Zhukovsky, K.V. Operational methods and Differential Equations with applications to initial-value problems. Appl. Math. Comput. 2007, 184, 979–1001. [Google Scholar] [CrossRef]
  56. Zhukovsky, K.V. Operational solution of differential equations with derivatives of non-integer order, Black–Scholes type and heat conduction. Mosc. Univ. Phys. Bull. 2016, 71, 237–244. [Google Scholar] [CrossRef]
  57. Dattoli, G.; Srivastava, H.M.; Zhukovsky, K.V. Orthogonality properties of the Hermite and related polynomials. J. Comput. Appl. Math. 2005, 182, 165–172. [Google Scholar] [CrossRef]
  58. Dattoli, G.; Srivastava, H.M.; Zhukovsky, K.V. A new family of integral transforms and their applications. Integral Transforms Spec. Funct. 2006, 17, 31–37. [Google Scholar] [CrossRef]
  59. Srivastava, H.M.; Manocha, H.L. A Treatise on Generating Functions; Halsted Press (Ellis Horwood Limited): Chichester, UK; John Wiley and Sons: New York, NY, USA; Chichester, UK; Brisbane, Australia; Toronto, ON, Canada, 1984. [Google Scholar]
  60. Boucetta, A.; Ghodbane, H.; Ayad, M.Y.; Bahri, M. A review on the performance and modelling of proton exchange membrane fuel cells. AIP Conf. Proc. 2016, 1758, 030019. [Google Scholar]
  61. Arato, E.; Pinna, M.; Mazzoccoli, M.; Bosio, B. Gas-phase mass-transfer resistances at polymeric electrolyte membrane fuel cells electrodes: Theoretical analysis on the effectiveness of interdigitated and serpentine flow arrangements. Energies 2016, 9, 229. [Google Scholar] [CrossRef]
  62. Maidhily, M.; Rajalakshmi, N.; Dhathathreyan, K.S. Electrochemical impedance spectroscopy as a diagnostic tool for the evaluation of flow field geometry in polymer electrolyte membrane fuel cells. Renew. Energy 2013, 51, 79–84. [Google Scholar] [CrossRef]
  63. St-Pierre, J. Hydrogen mass transport in fuel cell gas diffusion electrodes. Fuel Cells 2011, 11, 263–273. [Google Scholar] [CrossRef]
  64. Zhukovsky, K.; Pozio, A. Maximum current limitations of the PEM fuel cell with serpentine gas supply channels. J. Power Sources 2004, 130, 95–105. [Google Scholar] [CrossRef]
  65. Zhukovsky, K.V. Three Dimensional model of gas transport in a porous diffuser of a polymer electrolyte fuel cell. AIChE J. 2003, 49, 3029–3036. [Google Scholar] [CrossRef]
  66. Zhukovsky, K. Modeling of the Current Limitations of PEFC. AIChE J. 2006, 52, 2356–2366. [Google Scholar] [CrossRef]
  67. Weber, A.Z.; Newman, J. Modeling transport in polymer-electrolyte fuel cells. Chem. Rev. 2004, 104, 4679–4726. [Google Scholar] [CrossRef] [PubMed]
  68. Kawase, M.; Sato, K.; Mitsui, R.; Asonuma, H.; Kageyama, M.; Yamaguchi, K.; Inoue, G. Electrochemical reaction engineering of polymer electrolyte fuel cell. AIChE J. 2017, 63, 249–256. [Google Scholar] [CrossRef]
  69. Rogié, B.; Monier-Vinard, E.; Nguyen, M.-N.; Bissuel, V.; Laraqi, N. Practical analytical modeling of 3D multi-layer Printed Wired Board with buried volumetric heating sources. Int. J. Therm. Sci. 2018, 129, 404–415. [Google Scholar] [CrossRef]
  70. Calvo-Schwarzwälder, M.; Hennessy, M.G.; Torres, P.; Myers, T.G.; Alvarez, F.X. A slip-based model for the size-dependent effective thermal conductivity of nanowires. Int. Commun. Heat Mass Transf. 2018, 91, 57–63. [Google Scholar] [CrossRef] [Green Version]
  71. Zhukovskij, K.V. Gas flow in long microchannels. Vestn. Mosk. Univ. Ser. 3 Fiz. Astron. 2001, 3, 49–54. [Google Scholar]
  72. Zhukovsky, K.V. A harmonic solution for the hyperbolic heat conduction equation and its relationship to the Guyer–Krumhansl Equation. Mosc. Univ. Phys. Bull. 2018, 73, 45–52. [Google Scholar] [CrossRef]
  73. Zhukovsky, K. Exact harmonic solution to ballistic type heat propagation in thin films and wires. Int. J. Heat Mass Transf. 2018, 120, 944–955. [Google Scholar] [CrossRef]
  74. Zhukovsky, K.; Oskolkov, D. Exact harmonic solutions to Guyer–Krumhansl-type equation and application to heat transport in thin films. Contin. Mech. Thermodyn. 2018. [Google Scholar] [CrossRef]
  75. Ali, Y.M.; Zhang, L.C. Relativistic heat conduction. Int. J. Heat Mass Transf. 2005, 48, 2397–2406. [Google Scholar] [CrossRef]
  76. Al-Khairy, R.T.; Al-Ofey, Z.M. Analytical solution of the hyperbolic heat conduction equation for moving semi-infinite medium under the effect of time-dependent laser heat source. J. Appl. Math. 2009, 2009. [Google Scholar] [CrossRef]
  77. Wolf, K.B. Integral Transforms in Science and Engineering; New York Plenum Press: New York, NY, USA, 1979. [Google Scholar]
Figure 1. Evolution of quasi-temperature components in a thin stationary film with Knb = Knd = 1 for v = 0: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Figure 1. Evolution of quasi-temperature components in a thin stationary film with Knb = Knd = 1 for v = 0: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Axioms 07 00048 g001
Figure 2. Evolution of quasi-temperature components in a thin film with Knb = Knd = 1 for v = +10: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Figure 2. Evolution of quasi-temperature components in a thin film with Knb = Knd = 1 for v = +10: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Axioms 07 00048 g002
Figure 3. Evolution of quasi-temperature components in a thin film with Knb = Knd = 1 for v = −10: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Figure 3. Evolution of quasi-temperature components in a thin film with Knb = Knd = 1 for v = −10: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Axioms 07 00048 g003
Figure 4. Evolution of quasi-temperature components in a thin film with Knb = Knd = 0.1 for v = +10: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Figure 4. Evolution of quasi-temperature components in a thin film with Knb = Knd = 0.1 for v = +10: θ b ( x , t ) —top left plot, θ d ( x , t ) —top right plot and θ ( x , t ) = θ b ( x , t ) + θ d ( x , t ) —bottom plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Axioms 07 00048 g004
Figure 5. Evolution of quasi-temperature components in a stationary thin film with Knb = Knd = 0.1 for v = 0: θ b ( x , t ) —left plot, θ d ( x , t ) —right plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Figure 5. Evolution of quasi-temperature components in a stationary thin film with Knb = Knd = 0.1 for v = 0: θ b ( x , t ) —left plot, θ d ( x , t ) —right plot. The Cauchy conditions for PDEs (9) and (10) system are θ b ( x , 0 ) = θ d ( x , 0 ) = e i n x , t θ b ( x , 0 ) = t θ d ( x , 0 ) = 0 .
Axioms 07 00048 g005

Share and Cite

MDPI and ACS Style

Zhukovsky, K.; Oskolkov, D.; Gubina, N. Some Exact Solutions to Non-Fourier Heat Equations with Substantial Derivative. Axioms 2018, 7, 48. https://doi.org/10.3390/axioms7030048

AMA Style

Zhukovsky K, Oskolkov D, Gubina N. Some Exact Solutions to Non-Fourier Heat Equations with Substantial Derivative. Axioms. 2018; 7(3):48. https://doi.org/10.3390/axioms7030048

Chicago/Turabian Style

Zhukovsky, Konstantin, Dmitrii Oskolkov, and Nadezhda Gubina. 2018. "Some Exact Solutions to Non-Fourier Heat Equations with Substantial Derivative" Axioms 7, no. 3: 48. https://doi.org/10.3390/axioms7030048

APA Style

Zhukovsky, K., Oskolkov, D., & Gubina, N. (2018). Some Exact Solutions to Non-Fourier Heat Equations with Substantial Derivative. Axioms, 7(3), 48. https://doi.org/10.3390/axioms7030048

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop