短脉冲激光加热分数阶导热及其热应力研究

短脉冲激光加热引起材料内部复杂的传热过程及热变形,现有的以Fourier定律或Cattaneo-Vernotte松弛方程结合弹性理论为框架建立起来热应力理论在刻画其热物理过程存在严重缺陷. 本文基于分数阶微积分理论, 以半空间为研究对象, 建立了分数阶Cattaneo热传导方程和相应的热应力方程, 给出了问题的初始条件和边界条件, 采用拉普拉斯变换方法, 给出了非高斯时间分布激光热源辐射下温度场和热应力场的解析解, 研究了短脉冲激光加热的温度场及热应力场的热物理行为. 数值计算中, 首先对理论解进行数值验证, 然后取分数阶变量$p=0.5$研究温度场和热应力场的变化特点及激光参数对温度和热应力的影响,最后数值计算分数阶参数对温度和热应力场的影响. 计算结果表明, 分数阶Cattaneo传热方程和热应力方程描述的温度和热应力任然具有波动特性,与经典的Fourier传热模型和标准的Cattaneo传热模型相比, 分数阶阶次越大, 热波波速越小, 热波波动性越明显; 反之, 则热波波速越大, 热扩散性越强.激光加热和冷却的速度越快, 温度上升和下降的速度越快, 压应力和拉应力交替变化越快, 温度变化幅值越小, 热应力幅值影响不明显.      

Performance enhancement of a DC-operated micropump with electroosmosis in a hybrid nanofluid: fractional Cattaneo heat flux problem

The purpose of this investigation is to theoretically shed some light on the effect of the unsteady electroosmotic flow (EOF) of an incompressible fractional second-grade fluid with low-dense mixtures of two spherical nanoparticles, copper, and titanium. The flow of the hybrid nanofluid takes place through a vertical micro-channel. A fractional Cattaneo model with heat conduction is considered. For the DC-operated micropump, the Lorentz force is responsible for the pressure difference through the microchannel. The Debye-Hükel approximation is utilized to linearize the charge density. The semi-analytical solutions for the velocity and heat equations are obtained with the Laplace and finite Fourier sine transforms and their numerical inverses. In addition to the analytical procedures, a numerical algorithm based on the finite difference method is introduced for the given domain. A comparison between the two solutions is presented. The variations of the velocity heat transfer against the enhancements in the pertinent parameters are thoroughly investigated graphically. It is noticed that the fractional-order parameter provides a crucial memory effect on the fluid and temperature fields. The present work has theoretical implications for biofluid-based microfluidic transport systems.

     

Wave propagation model of heat conduction and group speed

In view of the finite relaxation model of non-Fourier’s law, the Cattaneo and Vernotte (CV) model and Fourier’s law are presented in this work for comparing wave propagation modes. Independent variable translation is applied to solve the partial differential equation. Results show that the general form of the time spatial distribution of temperature for the three media comprises two solutions: those corresponding to the positive and negative logarithmic heating rates. The former shows that a group of heat waves whose spatial distribution follows the exponential function law propagates at a group speed; the speed of propagation is related to the logarithmic heating rate. The total speed of all the possible heat waves can be combined to form the group speed of the wave propagation. The latter indicates that the spatial distribution of temperature, which follows the exponential function law, decays with time. These features show that propagation accelerates when heated and decelerates when cooled. For the model media that follow Fourier’s law and correspond to the positive heat rate of heat conduction, the propagation mode is also considered the propagation of a group of heat waves because the group speed has no upper bound. For the finite relaxation model with non-Fourier media, the interval of group speed is bounded and the maximum speed can be obtained when the logarithmic heating rate is exactly the reciprocal of relaxation time. And for the CV model with a non-Fourier medium, the interval of group speed is also bounded and the maximum value can be obtained when the logarithmic heating rate is infinite.     

Non-fourier heat conduction in a plane slab with prescribed boundary heat flux

The non-stationary heat conduction in an infinitely wide plane slab with a prescribed boundary heat flux is studied. An arbitrary time dependent boundary heat flux is considered and a non-vanishing thermal relaxation time is assumed. The temperature and the heat flux density distributions are determined analytically by employing Cattaneo-Vernotte's constitutive equation for the heat flux density. It is proved that the temperature and the heat flux density distributions can be incompatible with the hypothesis of local thermodynamic equilibrium. A comparison with the solution which would be obtained by means of Fourier's law is performed by considering the limit of a vanishing thermal relaxation time.     

Heat wave phenomena in solids subjected to time dependent surface heat flux

Hyperbolic heat conduction in a plane slab, infinitely long solid cylinder and solid sphere with a time dependent boundary heat flux is analytically studied. The solution is based on the separation of variables method and Duhamel’s principle. The temperature distribution, the propagation and reflection of the temperature wave and the effect of geometry on the shape of the wave front are studied for the case of a rectangular pulsed boundary heat flux. Comparisons with the solution obtained for Fourier heat conduction are performed by considering the limit of a vanishing thermal relaxation time.     

Transportation of heat through Cattaneo-Christov heat flux model in non-Newtonian fluid subject to internal resistance of particles

Thermal conduction which happens in all phases(liquid,solid,and gas) is the transportation of internal energy through minuscule collisions of particles and movement of electrons within a working body.The colliding particles comprise electrons,molecules,and atoms,and transfer disorganized microscopic potential and kinetic energy,mutually known as the internal energy.In engineering sciences,heat transfer comprises the processes of convection,thermal radiation,and sometimes mass transportation.Typically,more than one of these procedures may happen in a given circumstance.We use the Cattaneo-Christov(CC) heat flux model instead of the Fourier law of heat conduction to discuss the behavior of heat transportation.A mathematical model is presented for the Cattaneo-Christov double diffusion(CCDD) in the flow of a non-Newtonian nanofluid(the Jeffrey fluid) towards a stretched surface.The magnetohydrodynamic(MHD) fluid is considered.The behaviors of heat and mass transportation rates are discussed with the CCDD.These models are based on Fourier's and Fick's laws.The convective transportation in nanofluids is discussed,subject to thermophoresis and Brownian diffusions.The nonlinear governing flow expression is first altered into ordinary differential equations via appropriate transformations,and then numerical solutions are obtained through the built-in-shooting method.The impact of sundry flow parameters is discussed on the velocity,the skin friction coefficient,the temperature,and the concentration graphically.It is reported that the velocity of material particles decreases with higher values of the Deborah number and the ratio of the relaxation to retardation time parameter.The temperature distribution enhances when the Brownian motion and thermophoresis parameters increase.The concentration shows contrasting impact versus the Lewis number and the Brownian motion parameter.It is also noticed that the skin friction coefficient decreases when the ratio of the relaxation to retardation time parameter increases.     

The Cattaneo type space-time fractional heat conduction equation

The classical heat conduction equation is generalized using a generalized heat conduction law. In particular, we use the space-time Cattaneo heat conduction law that contains the Caputo symmetrized fractional derivative instead of gradient ${{\partial_x}}$ and fractional time derivative instead of the first order partial time derivative ${{\partial_t}}$ . The existence of the unique solution to the initial-boundary value problem corresponding to the generalized model is established in the space of distributions. We also obtain explicit form of the solution and compare it numerically with some limiting cases.     

A finite deformation fractional viscoplastic model for the glass transition behavior of amorphous polymers

The stress response of amorphous polymers exhibits tremendous change during the glass transition region, from soft viscoelastic response to stiff viscoplastic response. In order to describe the temperature-dependent and rate-dependent stress response of amorphous polymers, we extend the one-dimensional small strain fractional Zener model to the three-dimensional finite deformation model. The Eyring model is adopted to represent the stress-activated viscous flow. A phenomenological evolution equation of yield strength is used to describe the strain softening behaviors. We demonstrate that the stress response predicted by the three-dimensional model is consistent with that of one-dimensional model under uniaxial deformation, which confirms the validity of the extension. The model is then applied to describe the stress response of an amorphous thermoset at various temperatures and strain rates, which shows good agreement between experiments and simulation. We further perform a parameter study to investigate the influence of the model parameters on the stress response. The results show that a smaller fractional order results in a larger yield strain while has little effect on the yield stress when the temperature is below the glass transition temperature. For the stress relaxation tests, a smaller fractional order leads to a slower relaxation rate.     

Rewetting analysis of hot surfaces with internal heat source by the heat balance integral method

A two region conduction-controlled rewetting model of hot vertical surfaces with internal heat generation and boundary heat flux subjected to constant but different heat transfer coefficient in both wet and dry region is solved by the Heat Balance Integral Method (HBIM). The HBIM yields the temperature field and quench front temperature as a function of various model parameters such as Peclet number, Biot number and internal heat source parameter of the hot surface. Further, the critical (dry out) internal heat source parameter is obtained by setting Peclet number equal to zero, which yields the minimum internal heat source parameter to prevent the hot surface from being rewetted. Using this method, it has been possible to derive a unified relationship for a two-dimensional slab and tube with both internal heat generation and boundary heat flux. The solutions are found to be in good agreement with other analytical results reported in literature.     

Darcy-Forchheimer flow of variable conductivity Jeffrey liquid with Cattaneo-Christov heat flux theory

The role of the Cattaneo-Christov heat flux theory in the two-dimensional laminar flow of the Jeffrey liquid is discussed with a vertical sheet. The salient feature in the energy equation is accounted due to the implementation of the Cattaneo-Christov heat flux. A liquid with variable thermal conductivity is considered in the Darcy-Forchheimer porous space. The mathematical expressions of momentum and energy are coupled due to the presence of mixed convection. A highly nonlinear coupled system of equations is tackled with the homotopic algorithm. The convergence of the homotopy expressions is calculated graphically and numerically. The solutions of the velocity and temperature are expressed for various values of the Deborah number, the ratio of the relaxation time to the retardation time, the porosity parameter, the mixed convective parameter, the Darcy-Forchheimer parameter, and the conductivity parameter. The results show that the velocity and temperature are higher in Fourier's law of heat conduction cases in comparison with the Cattaneo-Christov heat flux model.     

Mathematical model of heat transfer in a fluid with account for its relaxation properties

Using the terms that take account for the temporal and spatial nonlocality (time variation of the heat flux and the temperature gradient) in the formula of Fourier’s law for the heat flux a differential equation for a fluid in motion is derived that contains the second time derivative and themixed derivative with respect to the spatial and temporal variables. Numerical solution of the problem of heat transfer in the laminar fluid flow in a plane channel demonstrates that, in view of the lag in the time variation of the heat flux from zero to a certain maximum value, the boundary condition of the first kind (thermal shock) cannot be instantaneously realized. The process of its stabilization on the wall is characterized by a certain time interval, whose duration is determined by the relaxation properties of the fluid. At large values of the dimensionless coefficients of the heat flux relaxation and the temperature gradient the boundary condition of the first kind can be realized only as the steady state is attainted, as Fo→∞. In this case, the flow does not contain temperature jumps and negative temperature values.     

Quantification of ignition time uncertainty based on the classical ignition theory and Fourier analysis

This study aims at modeling the effect of incoming heat flux fluctuations, on solid material ignition. In order to propose a general methodology based on the classical ignition theory that can be applied to any kind of solid target, kernels accounting for the target temperature response regarding an incoming heat flux are considered for thermally thick and thin solids with low or high thermal inertia. A Fourier decomposition of the incoming heat flux is then used to calculate the target response to harmonic heat fluxes. Finally, effects of harmonic fluctuations on ignition are discussed based on the previous analytical results, allowing us to discriminate situations where ignition time is expected to be rather predictable from situations where ignition time is expected to be less predictable thanks to an uncertainty quantification of the ignition time.     

Application of Fractional Derivatives in Thermal Analysis of Disk Brakes

This paper presents a Fractional Derivative Approach for thermal analysis of disk brakes. In this research, the problem is idealized as one-dimensional. The formulation developed contains fractional semi integral and derivative expressions, which provide an easy approach to compute friction surface temperature and heat flux as functions of time. Given the heat flux, the formulation provides a means to compute the surface temperature, and given the surface temperature, it provides a means to compute surface heat flux. A least square method is presented to smooth out the temperature curve and eliminate/reduce the effect of statistical variations in temperature due to measurement errors. It is shown that the integer power series approach to consider simple polynomials for least square purposes can lead to significant error. In contrast, the polynomials considered here contain fractional power terms. The formulation is extended to account for convective heat loss from the side surfaces. Using a simulated experiment, it is also shown that the present formulation predicts accurate values for the surface heat flux. Results of this study compare well with analytical and experimental results.     

Jeffrey fluid flow due to curved stretching surface with Cattaneo-Christov heat flux

The two-dimensional (2D) motion of the Jeffrey fluid by the curved stretching sheet coiled in a circle is investigated. The non-Fourier heat flux model is used for the heat transfer analysis. Feasible similarity variables are used to transform the highly nonlinear ordinary equations to partial differential equations (PDEs). The homotopy technique is used for the convergence of the velocity and temperature equations. The effects of the involved parameters on the physical properties of the fluid are described graphically. The results show that the curvature parameter is an increasing function of velocity and temperature, and the temperature is a decreasing function of the thermal relaxation time. Besides, the Deborah number has a reverse effect on the pressure and surface drag force.     

On the nature of the relaxation time,the Maxwell–Cattaneo and Fourier law in the thermodynamics of a continuous medium,and the scale effects in thermal conductivity

We present the theory of space–time elasticity and demonstrate that it is the extended reversible thermodynamics and gives the coupled model of thermoelasticity and heat conductivity and involves traditional thermoelasticity. We formulate the generally covariant variational model’s dynamic thermoelasticity and heat conductivity in which the basic kinematic and static variables are unified tensor objects (subject, matter). Variation statement defines the whole set of the initial-boundary problems for the 4D vector governing equation (Euler equation), the spatial projections of which define motion equations and the time projection gives the heat conductivity equation. We show that space–time elasticity directly implies the Fourier and the Maxwell–Cattaneo laws of heat conduction. However, space–time elasticity is richer than classical thermoelasticity, and it advocates its own equations of motion for coupled thermoelasticity. Moreover, we establish that the Maxwell–Cattaneo law and Fourier law can be defined for the reversible processes as compatibility equations without introducing dissipation. We argue that the present framework of space–time elasticity should prove adequate to describe the thermoelastic phenomena at low temperatures for interpreting the results of molecular simulations of heat conduction in solids and for the optimal heat and stress management in the microelectronic components and the thermoelectric devices.     

A uniqueness result for the Cattaneo-Christov heat conduction model applied to incompressible fluids

In the context of heat conduction governed by the celebrated Cattaneo equation, Christov has recently proposed a modification of the time derivative term in order to satisfy the objectivity principle. For such a model applied to an incompressible fluid, the uniqueness of the solution is here proved.     

分数阶理论下实心球体受热冲击作用的渐进分析

基于分数阶广义热弹性理论,针对实心球体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征,借助于Laplace 正、反变换技术及柱函数的渐进性质,推导了热冲击作用周期内位移场、温度场和应力场的渐进表达式. 通过计算,得到了不同传热能力下受热冲击作用时热波、热弹性的传播规律以及位移场、温度场及应力场的分布规律. 结果表明:分数阶参数取值的不同,热波、热弹性波的传播以及各物理场的分布均有所不同,分数阶参数可视为延迟时间的影响因子,通过改变延迟效应对热弹性行为的影响来改变热冲击的作用效果.      

Non-fourier temperature field in a solid homogeneous finite hollow cylinder

Analytical solution of the non-Fourier Axisymmetric temperature field within a finite hollow cylinder is investigated considering the Cattaneo-Vernotte constitutive heat flux relation. The solution is found for the most general linear time-independent boundary conditions. The material is assumed to be homogeneous and isotropic with temperature-independent thermal properties. The standard method of separation of variables is used. The present solution can be reduced to special problems of interest by choosing appropriate boundary condition parameters. The solution is applied for two special cases including constant heat flux and the Gaussian distribution heating of a cylinder, and their respective non-Fourier thermal behavior is studied.     

非稳态分数阶Oldroyd-B流体绕楔形体流动研究

研究了非稳态分数阶Oldroyd-B流体在多孔介质中通过楔形拉伸板的驻点流动问题。基于分数阶Oldroyd-B流体的本构模型建立了动量方程,并在其中引入了浮升力和驻点流动特征。此外,考虑了具有热松弛延迟时间的修正的分数阶Fourier定律,并将其应用于能量方程和对流换热边界条件。接着,采用与L1算法相结合的有限差分法求解控制偏微分方程。最后,分析了相关物理参数对流动的影响。结果表明,随着楔角参数的增加,流体受到的浮升力增大,导致速度加快;达西数越大,介质的孔隙度变大,流体的流动越快;此外,温度分布先略有上升后明显下降,这表明Oldroyd-B流体具有热延迟特性。     

Transient hyperbolic heat conduction in thick-walled FGM cylinders and spheres with exponentially-varying properties

This paper focuses on non-Fourier hyperbolic heat conduction analysis for heterogeneous hollow cylinders and spheres made of functionally graded material (FGM). All the material properties vary exponentially across the thickness, except for the thermal relaxation parameter which is taken to be constant. The cylinder and sphere are considered to be cylindrically and spherically symmetric, respectively, leading to one-dimensional heat conduction problems. The problems are solved analytically in the Laplace domain, and the results obtained are transformed to the real-time space using the modified Durbin’s numerical inversion method. The transient responses of temperature and heat flux are investigated for different inhomogeneity parameters and relative temperature change values. The comparisons of temperature distribution and heat flux between various time and material properties are presented in the form of graphs.