Qualitative properties of solutions in the time differential dual-phase-lag model of heat conduction

In this paper we study the time differential dual-phase-lag model of heat conduction incorporating the microstructural interaction effect in the fast-transient process of heat transport. We analyze the influence of the delay times upon some qualitative properties of the solutions of the initial boundary value problems associated to such a model. Thus, the uniqueness results are established under the assumption that the conductivity tensor is positive definite and the delay times τ_q and τ_T vary in the set {0 ≤ τ_q ≤ 2τ_T} ∪ {0 < 2τ_T < τ_q}. For the continuous dependence problem we establish two different estimates. The first one is obtained for the delay times with 0 ≤ τ_q ≤ 2τ_T, which agrees with the thermodynamic restrictions on the model in concern, and the solutions are stable. The second estimate is established for the delay times with 0 < 2τ_T < τ_q and it allows the solutions to have an exponential growth in time. The spatial behavior of the transient solutions and the steady-state vibrations is also addressed. For the transient solutions we establish a theorem of influence domain, under the assumption that the delay times are in {0 < τ_q ≤ 2τ_T} ∪ {0 < 2τ_T < τ_q}. While for the amplitude of the harmonic vibrations we obtain an exponential decay estimate of Saint–Venant type, provided the frequency of vibration is lower than a critical value and without any restrictions upon the delay times.     

Finite integral transform-based analytical solutions of dual phase lag bio-heat transfer equation

A finite integral transform (FIT)-based analytical solution to the dual phase lag (DPL) bio-heat transfer equation has been developed. One of the potential applications of this analytical approach is in the field of photo-thermal therapy, wherein the interest lies in determining the thermal response of laser-irradiated biological samples. In order to demonstrate the applicability of the generalized analytical solutions, three problems have been formulated: (1) time independent boundary conditions (constant surface temperature heating), (2) time dependent boundary conditions (medium subjected to sinusoidal surface heating), and (3) biological tissue phantoms subjected to short-pulse laser irradiation. In the context of the case study involving biological tissue phantoms, the FIT-based analytical solutions of Fourier, as well as non-Fourier, heat conduction equations have been coupled with a numerical solution of the transient form of the radiative transfer equation (RTE) to determine the resultant temperature distribution. Performance of the FIT-based approach has been assessed by comparing the results of the present study with those reported in the literature. A comparison of DPL-based analytical solutions with those obtained using the conventional Fourier and hyperbolic heat conduction models has been presented. The relative influence of relaxation times associated with the temperature gradients (τ_T) and heat flux (τ_q) on the resultant thermal profiles has also been discussed. To the best of the knowledge of the authors, the present study is the first successful attempt at developing complete FIT-based analytical solution(s) of non-Fourier heat conduction equation(s), which have subsequently been coupled with numerical solutions of the transient form of the RTE. The work finds its importance in a range of areas such as material processing, photo-thermal therapy, etc.     

Efficient and accurate spectral method for the time-fractional dual-phase-lag heat transfer model and its parameter estimation

In the current paper, a heat transfer model is suggested based on a time-fractional dual-phase-lag (DPL) model. We discuss the model in two parts, the direct problem and the inverse problem. Firstly, for solving it, the finite difference/Legendre spectral method is constructed. In the temporal direction, we employ the weighted and shifted Grünwald approximation, which can achieve second order convergence. In the spatial direction, the Legendre spectral method is used, it can obtain spectral accuracy. The stability and convergence are theoretically analyzed. For the inverse problem, the Bayesian method is used to construct an algorithm to estimate the four parameters for the model, namely, the time-fractional order α, the time-fractional order β, the delay time τ_T, and the relaxation time τ_q. Next, numerical experiments are provided to demonstrate the effectiveness of our scheme, with the values of τ_q and τ_T for processed meat employed. We also make a comparison with another method. The data obtained for the direct problem are used in the parameter estimation. The paper provides an accurate and efficient numerical method for the time-fractional DPL model.     

Stagnation-point flow and heat transfer over an exponentially shrinking sheet

An analysis is carried out to investigate the stagnation-point flow and heat transfer over an exponentially shrinking sheet. Using the boundary layer approximation and a similarity transformation in exponential form, the governing mathematical equations are transformed into coupled, nonlinear ordinary differential equations which are then solved numerically by a shooting method with fourth order Runge-Kutta integration scheme. The analysis reveals that a solution exists only when the velocity ratio parameter satisfies the inequality −1.487068 ? c/a. Also, the numerical calculations exhibit the existence of dual solutions for the velocity and the temperature fields; and it is observed that their boundary layers are thinner for the first solution (in comparison with the second). Moreover, the heat transfer from the sheet increases with an increase in c/a for the first solution, while the heat transfer decreases with increasing c/a for the second solution, and ultimately heat absorption occurs.     

Local sensitivity analysis for the heat flux-temperature integral relationship in the half-space

This paper presents insight into the heat flux-temperature (q″ ? T) integral relationship based on constant thermophysical properties. This relationship is often used in one-dimensional, transient heat transfer studies involving null-point calorimetry and heat flux investigations. This study focuses on a short transient studies where energy has not fully penetrated the body as the result of an imposed surface heating condition. A full nonlinear heat transfer model is developed involving a half-space planar region. Temperature results are then introduced into the constant property integral relationship and a newly derived Kirchoff integral relationship for retrieving the local heat flux. Good agreement is observed between the fully nonlinear results and locally linearized system. Additionally, a sensitivity study is presented which involves perturbing the average thermophysical properties of thermal conductivity and heat capacity.     

On the uniqueness and reconstruction for an inverse problem of the fractional diffusion process

Consider an inverse problem for the time-fractional diffusion equation in one dimensional spatial space. The aim is to determine the initial status and heat flux on the boundary simultaneously from heat measurement data given on the other boundary. Using the Laplace transform and the unique extension technique, the uniqueness for this inverse problem is proven. Then we construct a regularizing scheme for the reconstruction of boundary flux for known initial status. The convergence rate of the regularizing solution is established under some a priori information about the exact solution. Moreover, the initial distribution can also be recovered approximately from our regularizing scheme. Finally we present some numerical examples, which show the validity of the proposed reconstruction scheme.     

Transient nonlinear optically-thick radiative–convective double-diffusive boundary layers in a Darcian porous medium adjacent to an impulsively started surface: Network simulation solutions

A boundary-layer model is described for the two-dimensional nonlinear transient thermal convection heat and mass transfer in an optically-thick fluid in a Darcian porous medium adjacent to an impulsively started vertical surface, in the presence of significant thermal radiation and buoyancy forces in an (X¹,Y¹,t¹) coordinate system. An algebraic approximation is employed to simplify the integro-differential equation of radiative transfer for unidirectional flux normal to the plate into the boundary-layer regime, by incorporating this flux term in the energy conservation equation. The conservation equations are non-dimensionalized into an (X,Y,T) coordinate system and solved using the Network Simulation Method (NSM), a robust numerical technique which demonstrates high efficiency and accuracy. The transient variation of non-dimensional streamwise velocity component (u) and temperature (T) and concentration (C) functions is computed for various selected values of Stark number (radiation–conduction interaction parameter) and Darcy number. Transient velocity (u) and steady-state local skin friction (τ_X) are also studied for various thermal Grashof number (Gr), species Grashof number (Gm), Schmidt number (Sc) and Stark number (N) values. These computations for the infinite permeability case (Da → ∞) are compared with previous finite difference solutions [Prasad et al. Int J Therm Sci 2007;46(12):1251–8] and shown to be in excellent agreement. An increase in Darcy number is seen to accelerate the flow and boost velocity. A decrease in Stark number (corresponding to an increase in thermal radiation heat transfer contribution) is shown to increase the velocity values. Temperature function is observed to fall in value with a rise in Da and increase with decrease in N (corresponding to an increase in thermal radiation heat transfer contribution). Applications of the study include rocket combustion chambers, astrophysical flows, spacecraft thermal fluid dynamics in debris-laden environments (cosmic dust), heat transfer in forest fire spread, geochemical contamination and ceramic materials processing.     

Characterization of the spectrum of regular boundary value problems for the Sturm-Liouville operator

On the interval (0, τ), we consider the spectral problem generated by the Sturm-Liouville equation with an arbitrary complex-valued potential q(x) ∈ L ₂(0, τ) and with regular (but not strengthened-regular) boundary conditions. Under certain additional assumptions, we establish necessary and sufficient conditions for a set of complex numbers to be the spectrum of such an operator.     

Atomic decomposition of Besov-type and Triebel-Lizorkin-type spaces

Abstract With the help of the maximal function caracterizations of the Besov-type space Bs,τp,q and the TriebelLizorkin-type space Fs,τp,q,we present the atomic decomposition of these function spaces.Our results cover the results on classical Besov and Triebel-Lizorkin spaces by taking τ=0.     

Effects of thermal radiation and magnetic field on unsteady mixed convection flow and heat transfer over an exponentially stretching surface with suction in the presence of internal heat generation/absorption

In this paper, the problem of unsteady laminar two-dimensional boundary layer flow and heat transfer of an incompressible viscous fluid in the presence of thermal radiation, internal heat generation or absorption, and magnetic field over an exponentially stretching surface subjected to suction with an exponential temperature distribution is discussed numerically. The governing boundary layer equations are reduced to a system of ordinary differential equations. New numerical method using Mathematica has been used to solve such system after obtaining the missed initial conditions. Comparison of obtained numerical results is made with previously published results in some special cases, and found to be in a good agreement.     

Trattazione rapida dei problemi al contorno relativi alle equazioni differenziali lineari a derivate parziali

Let Σ andS be two real Hilbert spaces and Σ₀ a subspace of Σ. Moreover, supposeT:S→Σ be a bounded linear operator whose rangeT (S) is contained in Σ₀, andE:S→Σ be a linear operator such that the productE T:S→S is a bounded operator with a closed range. In this framework we present an artifice from which the alternative theorem for the equationE u=q(u?Σ₀,q ?S) follows. It is worthwhile to note thatEu=q may represent a boundary value problem for elliptic equations.     

Effects of thermal radiation on the boundary layer flow of a Jeffrey fluid over an exponentially stretching surface

In the present investigation we have analyzed the boundary layer flow of a Jeffrey fluid over an exponentially stretching surface. The effects of thermal radiation are carried out for two cases of heat transfer analysis known as (1) Prescribed exponential order surface temperature (PEST) and (2) Prescribed exponential order heat flux (PEHF). The highly nonlinear coupled partial differential equations of Jeffrey fluid flow along with the energy equation are simplified by using similarity transformation techniques based on boundary layer assumptions. The reduced similarity equations are then solved analytically by the homotopy analysis method (HAM). The convergence of the HAM series solution is obtained by plotting (^h/_2p)\hbar-curves for velocity and temperature. The effects of physical parameters on the velocity and temperature profiles are examined by plotting graphs.     

Exponentially decaying boundary layers as limiting cases of families of algebraically decaying ones

The boundary value problem for the similar stream function f = f(η;λ) of the Cheng–Minkowycz free convection flow over a vertical plate with a power law temperature distribution T_w(x) = T_∞ + Ax^λ in a porous medium is revisited. It is shown that in the λ-range − 1/2 < λ < 0 , the well known exponentially decaying “first branch” solutions for the velocity and temperature fields are not some isolated solutions as one has believed until now, but limiting cases of families of algebraically decaying multiple solutions. For these multiple solutions well converging analytical series expressions are given. This result yields a bridging to the historical quarreling concerning the feasibility of exponentially and algebraically decaying boundary layers. Owing to a mathematical analogy, our results also hold for the similar boundary layer flows induced by continuous surfaces stretched in viscous fluids with power-law velocities u_w(x)∼ x^λ. (Received: June 7, 2005)     

Optimal paths for minimizing entransy dissipation during heat transfer processes with generalized radiative heat transfer law

A common of finite-time heat transfer processes between high- and low-temperature sides with generalized radiative heat transfer law [q ∝ Δ(Tⁿ)] is studied in this paper. In general, the minimization of entropy generation in heat transfer processes is taken as the optimization objective. A new physical quantity, entransy, has been identified as a basis for optimizing heat transfer processes in terms of the analogy between heat and electrical conduction recently. Heat transfer analyses show that the entransy of an object describes its heat transfer ability, as the electrical energy in a capacitor describes its charge transfer ability. Entransy dissipation occurs during heat transfer processes, as a measure of the heat transfer irreversibility with the dissipation related thermal resistance. Under the condition of fixed heat load, the optimal configurations of hot and cold fluid temperatures for minimizing entransy dissipation are derived by using optimal control theory. The condition corresponding to the minimum entransy dissipation strategy with Newtonian heat transfer law (n = 1) is that corresponding to a constant heat flux rate, while the condition corresponding to the minimum entransy dissipation strategy with the linear phenomenological heat transfer law (n = −1) is that corresponding to a constant ratio of hot to cold fluid temperatures. Numerical examples for special cases with Newtonian, linear phenomenological and radiative heat transfer law (n = 4) are provided, and the obtained results are also compared with the conventional strategies of constant heat flux rate and constant hot fluid (reservoir) temperature operations and optimal strategies for minimizing entropy generation. Moreover, the effects of heat load changes on the optimal hot and fluid temperature configurations are also analyzed.     

The forced convection flow of a uniform stream over a flat surface with a convective surface boundary condition

The forced convection heat transfer resulting from the flow of a uniform stream over a flat surface on which there is a convective boundary condition is considered. In previous papers [5], [6], [7], [8] it was assumed that the convective heat transfer parameter h_f associated with the hot surface depended on x, where x measures distance along the surface, so that problem could be reduced to similarity form. Here it is assumed that this heat transfer parameter h_f is a constant, with the result that the temperature profiles and overall heat transfer characteristics evolve as the solution develops from the leading edge. The heat transfer near the leading edge (small x), which we find to be dominated by the surface heat flux, the solution at large distances along the surface (large x), which dominated by the surface temperature, are discussed. A numerical solution to the full problem is then obtained for a range of values of the Prandtl number to join these two solution regimes.     

Initially and eventually dominating nonoscillatory solutions in integro-differential difference equations

For the system T′(t) + p(t) T(t) + q(t) T(t ? τ) = ∝₀^tK(t ? μ) T(μ) dμ for t ? 0, T(t) = g(t) for t ∈ [?τ, 0], conditions have been obtained which ensure that a solution of this system is dominated by a nonoscillatory solution in the interval ¦τ, ∞).     

Long-time behavior of the difference solutions to the 2D GKS equation

We study the long-time behavior of the finite difference solution to the generalized Kuramoto-Sivashinsky equation in two space dimensions with periodic boundary conditions. The unique solvability of numerical solution is shown. It is proved that there exists a global attractor of the discrete dynamical system and the upper semicontinuity d(A_h,τ,A)→0. Finally, we obtain the long-time stability and convergence of the difference scheme. Our results show that the difference scheme can effectively simulate the infinite dimensional dynamical systems.     

Efficient long-time computations of time-domain boundary integrals for 2D and dissipative wave equation

Linear hyperbolic partial differential equations in a homogeneous medium, e.g., the wave equation describing the propagation and scattering of acoustic waves, can be reformulated as time-domain boundary integral equations. We propose an efficient implementation of a numerical discretization of such equations when the strong Huygens’ principle does not hold.For the numerical discretization, we make use of convolution quadrature in time and standard Galerkin boundary element method in space. The quadrature in time results in a discrete convolution of weights W_j with the boundary density evaluated at equally spaced time points. If the strong Huygens’ principle holds, W_j converge to 0 exponentially quickly for large enough j. If the strong Huygens’ principle does not hold, e.g., in even space dimensions or when some damping is present, the weights are never zero, thereby presenting a difficulty for efficient numerical computation.In this paper we prove that the kernels of the convolution weights approximate in a certain sense the time domain fundamental solution and that the same holds if both are differentiated in space. The tails of the fundamental solution being very smooth, this implies that the tails of the weights are smooth and can efficiently be interpolated. Further, we hint on the possibility to apply the fast and oblivious convolution quadrature algorithm of Schädle et al. to further reduce memory requirements for long-time computation. We discuss the efficient implementation of the whole numerical scheme and present numerical experiments.     

Exponentially decaying boundary layers as limiting cases of families of algebraically decaying ones

The boundary value problem for the similar stream function f  =  f(η;λ) of the Cheng–Minkowycz free convection flow over a vertical plate with a power law temperature distribution T_w(x)  =  T_∞ + Ax^λ in a porous medium is revisited. It is shown that in the λ-range  − 1/2  < λ  <  0 , the well known exponentially decaying “first branch” solutions for the velocity and temperature fields are not some isolated solutions as one has believed until now, but limiting cases of families of algebraically decaying multiple solutions. For these multiple solutions well converging analytical series expressions are given. This result yields a bridging to the historical quarreling concerning the feasibility of exponentially and algebraically decaying boundary layers. Owing to a mathematical analogy, our results also hold for the similar boundary layer flows induced by continuous surfaces stretched in viscous fluids with power-law velocities u_w(x)∼ x^λ.     

A fundamental region for Hecke's modular group

Hecke proved analytically that when λ ≥ 2 or when $\text{λ = 2}\text{cos}\text{(}\text{π}\text{q}\text{)}$, q ∈ Z, q ≥ 3, then $\text{B(λ) = {τ:}\text{Im}\text{τ > 0, |}\text{Re}\text{τ| <}\text{λ}\text{2}\text{, |τ| > 1}}$ is a fundamental region for the group G(λ) = 〈S_λ, T〉, where S_λ: τ → τ + λ and $\text{T: τ →}\text{?1}\text{τ}$. He also showed that B(λ) fails to be a fundamental region for all other λ > 0 by proving that G(λ) is not discontinuous. We give an elementary proof of these facts and prove a related result concerning the distribution of G(λ)-equivalent points.