In this paper,using the fractional Fourier law,we obtain the fractional heat conduction equation with a time-fractional derivative in the spherical coordinate system.The method of variable separation is used to solve ...In this paper,using the fractional Fourier law,we obtain the fractional heat conduction equation with a time-fractional derivative in the spherical coordinate system.The method of variable separation is used to solve the timefractional heat conduction equation.The Caputo fractional derivative of the order 0 〈 α≤ 1 is used.The solution is presented in terms of the Mittag-Leffler functions.Numerical results are illustrated graphically for various values of fractional derivative.展开更多
A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh rat...A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.展开更多
In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its ...In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its high ability to conduct heat and electrical conductivity, in addition to being a flexible and malleable metal that is easy to form without being broken, making it one of the basic minerals that humans have benefited from for thousands of years, it is one of the first minerals. That has been discovered and extracted, and still plays a major role in the development of societies. The obtained solutions are compared with the available exact solutions and the obtained solutions using the finite difference method. The results indicate that the finite difference method is a highly effective method for obtaining approximate solutions for the thermal conductivity equation for copper. It is also clear from the numerical results from copper in the high conductivity of heat and electricity.展开更多
It is known that Fourier’s heat equation, which is parabolic, implies an infinite velocity propagation, or, in other words, that the mechanism of heat conduction is established instantaneously under all conditions. T...It is known that Fourier’s heat equation, which is parabolic, implies an infinite velocity propagation, or, in other words, that the mechanism of heat conduction is established instantaneously under all conditions. This is unacceptable on physical grounds in spite of the fact that Fourier’s law agrees well with experiment. However, discrepancies are likely to occur when extremely short distances or extremely short time intervals are considered, as they must in some modern problems of aero-thermodynamics. Cattaneo and independently Vernotte proved that such process can be described by Heaviside’s telegraph equation. This paper shows that this fact can be derived using calculus of variations, by application of the Euler-Lagrange equation. So, we proved that the equation of heat conduction with finite velocity propagation of the thermal disturbance can be obtained as a solution to one variational problem.展开更多
The paper proposes an approximate solution to the classical (parabolic) multidimensional 2D and 3D heat conduction equation for a 5 × 5 cm aluminium plate and a 5 × 5 × 5 cm aluminum cube. An approximat...The paper proposes an approximate solution to the classical (parabolic) multidimensional 2D and 3D heat conduction equation for a 5 × 5 cm aluminium plate and a 5 × 5 × 5 cm aluminum cube. An approximate solution of the generalized (hyperbolic) 2D and 3D equation for the considered plate and cube is also proposed. Approximate solutions were obtained by applying calculus of variations and Euler-Lagrange equations. In order to verify the correctness of the proposed approximate solutions, they were compared with the exact solutions of parabolic and hyperbolic equations. The paper also presents the research on the influence of time parameters τ as well as the relaxation times τ ∗ to the variation of the profile of the temperature field for the considered aluminum plate and cube.展开更多
This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discr...This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.展开更多
To eliminate oscillation and overbounding of finite element solutions of classical heat conduction equation, the author and Xiao have put forward two new concepts of monotonies and have derived and proved several crit...To eliminate oscillation and overbounding of finite element solutions of classical heat conduction equation, the author and Xiao have put forward two new concepts of monotonies and have derived and proved several criteria. This idea is borrowed here to deal with generalized conduction equation and finite element criteria for eliminating oscillation and overbounding are also presented. Some new and useful conclusions are drawn.展开更多
In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat condu...In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.展开更多
This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem ...This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem is considered steady-state but solved at each time iteration of the transient conduction problem. The discrete ordinate method along with the decentered streamline-upwind Petrov-Galerkin method is developed. Since specular reflection is considered on borders, a very accurate algorithm has been developed for calculation of partition ratio coefficients of incident solid angles to the several reflected solid angles. The developed algorithms are tested on a paraboloid-shaped geometry used for example on concentrated solar power technologies.展开更多
Several researchers have dealt with the one-dimensional fractional heat conduction equation in the last decades,but as far as we know,no one has investigated such a problem from the perspective of developing suitable ...Several researchers have dealt with the one-dimensional fractional heat conduction equation in the last decades,but as far as we know,no one has investigated such a problem from the perspective of developing suitable fractional-order methods.This has actually motivated us to address this problem by the way of establishing a proper fractional approach that involves employing a combination of a novel fractional difference formula to approximate the Caputo differentiator of orderαcoupled with the modified three-point fractional formula to approximate the Caputo differentiator of order 2α,where 0<α≤1.As a result,the fractional heat conduction equation is then reexpressed numerically using the aforementioned formulas,and by dividing the considered mesh into multiple nodes,a system is generated and algebraically solved with the aid of MATLAB.This would allow us to obtain the desired approximate solution for the problem at hand.展开更多
In this paper the nonlinear heat-conduction equations rhoc partial derivativew/partial derivativet = partial derivative/partial derivativex (k partial derivativew/partial derivativex) with Dirichlet boundary condition...In this paper the nonlinear heat-conduction equations rhoc partial derivativew/partial derivativet = partial derivative/partial derivativex (k partial derivativew/partial derivativex) with Dirichlet boundary condition and the nonlinear boundary condition are studied. The asymptotic behavior of the global of solution are analyzed by using Lyapuunov function. As its application, the approximate solutions are constructed.展开更多
To better describe the phenomenon of non-Fourier heat conduction, the fractional Cattaneo heat equation is introduced from the generalized Cattaneo model with two fractional derivatives of different orders. The anomal...To better describe the phenomenon of non-Fourier heat conduction, the fractional Cattaneo heat equation is introduced from the generalized Cattaneo model with two fractional derivatives of different orders. The anomalous heat conduction under the Neumann boundary condition in a semi-infinity medium is investigated. Exact solutions are obtained in series form of the H-function by using the Laplace transform method. Finally, numerical examples are presented graphically when different kinds of surface temperature gradient are given. The effects of fractional parameters are also discussed.展开更多
The classical heat conduction equation is derived from the assumption that the temperature increases immediately after heat transfer, but the increase of temperature is a slow process, so the memory-dependent heat con...The classical heat conduction equation is derived from the assumption that the temperature increases immediately after heat transfer, but the increase of temperature is a slow process, so the memory-dependent heat conduction model has been reconstructed. Numerical results show that the solution of the initial boundary value problem of the new model is similar to that of the classical heat conduction equation, but its propagation speed is slower than that of the latter. In addition, the propagation speed of the former is also affected by time delay and kernel function.展开更多
In this paper we consider a non-standard inverse heat conduction problem for determining surface heat flux from an interior observation which appears in some applied subjects. This problem is ill-posed in the sense th...In this paper we consider a non-standard inverse heat conduction problem for determining surface heat flux from an interior observation which appears in some applied subjects. This problem is ill-posed in the sense that the solution (if it exists)does not depend continuously on the data. A Fourier method is applied to formulate a regularized approximation solution, and some sharp error estimates are also given.展开更多
In this paper,by me as of beundary element method,we try to deal with the initial -boundary value problem for a class of linear parunolic equations,which is a linear heat conduction equation. We tresent a boundary int...In this paper,by me as of beundary element method,we try to deal with the initial -boundary value problem for a class of linear parunolic equations,which is a linear heat conduction equation. We tresent a boundary integral equation for the solution to the problem and its variational formalation The well-posedness of the variational formulation is proved. And the error estimates for the approsutate solutions are provided. The results of this paper are more general than those of[1]展开更多
We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the lo...We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the local strong solutions in terms of the gradient of the velocity only, which coincides with the famous Beale-Kato-Majda criterion for ideal incompressible flows.展开更多
The effect of internal heating source on the film momentum and thermal transport characteristic of thin finite power-law liquids over an accelerating unsteady horizontal stretched interface is studied. Unlike most cla...The effect of internal heating source on the film momentum and thermal transport characteristic of thin finite power-law liquids over an accelerating unsteady horizontal stretched interface is studied. Unlike most classical works in this field, a general surface temperature distribution of the liquid film and the generalized Fourier's law for varying thermal conductivity are taken into consideration. Appropriate similarity transformations are used to convert the strongly nonlinear governing partial differential equations (PDEs) into a boundary value problem with a group of two-point ordinary differential equations (ODEs). The correspondence between the liquid film thickness and the unsteadiness parameter is derived with the BVP4C program in MATLAB. Numerical solutions to the self-similarity ODEs are obtained using the shooting technique combined with a Runge-Kutta iteration program and Newton's scheme. The effects of the involved physical parameters on the fluid's horizontal velocity and temperature distribution are presented and discussed.展开更多
The inverse heat conduction problem (IHCP) is a severely ill-posed problem in the sense that the solution ( if it exists) does not depend continuously on the data. But now the results on inverse heat conduction pr...The inverse heat conduction problem (IHCP) is a severely ill-posed problem in the sense that the solution ( if it exists) does not depend continuously on the data. But now the results on inverse heat conduction problem are mainly devoted to the standard inverse heat conduction problem. Some optimal error bounds in a Sobolev space of regularized approximation solutions for a sideways parabolic equation, i. e. , a non-standard inverse heat conduction problem with convection term which appears in some applied subject are given.展开更多
The present paper deals with the determination of quasi-static thermal stresses due to an instantaneous point heat source of strength g_(pi) situated at certain circle along the radial direction of the circular plate ...The present paper deals with the determination of quasi-static thermal stresses due to an instantaneous point heat source of strength g_(pi) situated at certain circle along the radial direction of the circular plate and releasing its heat spontaneously at time t=τ.A circular plate is considered having arbitrary initial temperature and subjected to time dependent heat flux at the fixed circular boundary of r=b.The governing heat conduction equation is solved by using the integral transform method,and results are obtained in series form in terms of Bessel functions.The mathematical model has been constructed for copper material and the thermal stresses are discussed graphically.展开更多
Heat induced in a finite Silver Selenide slab by an external pulsed laser source is studied in dimensionless scale according to hyperbolic heat conduction model (HHCE) using Laplace integral transform technique. The t...Heat induced in a finite Silver Selenide slab by an external pulsed laser source is studied in dimensionless scale according to hyperbolic heat conduction model (HHCE) using Laplace integral transform technique. The temperature profile, the critical time required to initiate phase transition and that to initiate damage at the front surface are obtained for different pulses and are illustrated graphically.展开更多
基金supported by the National Natural Science Foundation of China(11072134 and 11102102)
文摘In this paper,using the fractional Fourier law,we obtain the fractional heat conduction equation with a time-fractional derivative in the spherical coordinate system.The method of variable separation is used to solve the timefractional heat conduction equation.The Caputo fractional derivative of the order 0 〈 α≤ 1 is used.The solution is presented in terms of the Mittag-Leffler functions.Numerical results are illustrated graphically for various values of fractional derivative.
文摘A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.
文摘In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its high ability to conduct heat and electrical conductivity, in addition to being a flexible and malleable metal that is easy to form without being broken, making it one of the basic minerals that humans have benefited from for thousands of years, it is one of the first minerals. That has been discovered and extracted, and still plays a major role in the development of societies. The obtained solutions are compared with the available exact solutions and the obtained solutions using the finite difference method. The results indicate that the finite difference method is a highly effective method for obtaining approximate solutions for the thermal conductivity equation for copper. It is also clear from the numerical results from copper in the high conductivity of heat and electricity.
文摘It is known that Fourier’s heat equation, which is parabolic, implies an infinite velocity propagation, or, in other words, that the mechanism of heat conduction is established instantaneously under all conditions. This is unacceptable on physical grounds in spite of the fact that Fourier’s law agrees well with experiment. However, discrepancies are likely to occur when extremely short distances or extremely short time intervals are considered, as they must in some modern problems of aero-thermodynamics. Cattaneo and independently Vernotte proved that such process can be described by Heaviside’s telegraph equation. This paper shows that this fact can be derived using calculus of variations, by application of the Euler-Lagrange equation. So, we proved that the equation of heat conduction with finite velocity propagation of the thermal disturbance can be obtained as a solution to one variational problem.
文摘The paper proposes an approximate solution to the classical (parabolic) multidimensional 2D and 3D heat conduction equation for a 5 × 5 cm aluminium plate and a 5 × 5 × 5 cm aluminum cube. An approximate solution of the generalized (hyperbolic) 2D and 3D equation for the considered plate and cube is also proposed. Approximate solutions were obtained by applying calculus of variations and Euler-Lagrange equations. In order to verify the correctness of the proposed approximate solutions, they were compared with the exact solutions of parabolic and hyperbolic equations. The paper also presents the research on the influence of time parameters τ as well as the relaxation times τ ∗ to the variation of the profile of the temperature field for the considered aluminum plate and cube.
文摘This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.
文摘To eliminate oscillation and overbounding of finite element solutions of classical heat conduction equation, the author and Xiao have put forward two new concepts of monotonies and have derived and proved several criteria. This idea is borrowed here to deal with generalized conduction equation and finite element criteria for eliminating oscillation and overbounding are also presented. Some new and useful conclusions are drawn.
文摘In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.
文摘This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem is considered steady-state but solved at each time iteration of the transient conduction problem. The discrete ordinate method along with the decentered streamline-upwind Petrov-Galerkin method is developed. Since specular reflection is considered on borders, a very accurate algorithm has been developed for calculation of partition ratio coefficients of incident solid angles to the several reflected solid angles. The developed algorithms are tested on a paraboloid-shaped geometry used for example on concentrated solar power technologies.
文摘Several researchers have dealt with the one-dimensional fractional heat conduction equation in the last decades,but as far as we know,no one has investigated such a problem from the perspective of developing suitable fractional-order methods.This has actually motivated us to address this problem by the way of establishing a proper fractional approach that involves employing a combination of a novel fractional difference formula to approximate the Caputo differentiator of orderαcoupled with the modified three-point fractional formula to approximate the Caputo differentiator of order 2α,where 0<α≤1.As a result,the fractional heat conduction equation is then reexpressed numerically using the aforementioned formulas,and by dividing the considered mesh into multiple nodes,a system is generated and algebraically solved with the aid of MATLAB.This would allow us to obtain the desired approximate solution for the problem at hand.
文摘In this paper the nonlinear heat-conduction equations rhoc partial derivativew/partial derivativet = partial derivative/partial derivativex (k partial derivativew/partial derivativex) with Dirichlet boundary condition and the nonlinear boundary condition are studied. The asymptotic behavior of the global of solution are analyzed by using Lyapuunov function. As its application, the approximate solutions are constructed.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11102102, 11072134, and 91130017)the Natural Science Foundation of Shandong Province, China (Grant No. ZR2009AQ014)the Independent Innovation Foundation of Shandong University, China (Grant No. 2010ZRJQ002)
文摘To better describe the phenomenon of non-Fourier heat conduction, the fractional Cattaneo heat equation is introduced from the generalized Cattaneo model with two fractional derivatives of different orders. The anomalous heat conduction under the Neumann boundary condition in a semi-infinity medium is investigated. Exact solutions are obtained in series form of the H-function by using the Laplace transform method. Finally, numerical examples are presented graphically when different kinds of surface temperature gradient are given. The effects of fractional parameters are also discussed.
文摘The classical heat conduction equation is derived from the assumption that the temperature increases immediately after heat transfer, but the increase of temperature is a slow process, so the memory-dependent heat conduction model has been reconstructed. Numerical results show that the solution of the initial boundary value problem of the new model is similar to that of the classical heat conduction equation, but its propagation speed is slower than that of the latter. In addition, the propagation speed of the former is also affected by time delay and kernel function.
基金NNSF of China (No. 10271050)NSF of Gansu Province, China (ZS021-A25-001-Z)
文摘In this paper we consider a non-standard inverse heat conduction problem for determining surface heat flux from an interior observation which appears in some applied subjects. This problem is ill-posed in the sense that the solution (if it exists)does not depend continuously on the data. A Fourier method is applied to formulate a regularized approximation solution, and some sharp error estimates are also given.
文摘In this paper,by me as of beundary element method,we try to deal with the initial -boundary value problem for a class of linear parunolic equations,which is a linear heat conduction equation. We tresent a boundary integral equation for the solution to the problem and its variational formalation The well-posedness of the variational formulation is proved. And the error estimates for the approsutate solutions are provided. The results of this paper are more general than those of[1]
基金supported by the China Postdoctoral Science Foundation (20090450333)supported by the National Basic Research Program (2005CB321700)NSFC (40890154)
文摘We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the local strong solutions in terms of the gradient of the velocity only, which coincides with the famous Beale-Kato-Majda criterion for ideal incompressible flows.
基金Project supported by the Scientific Research Funds of Huaqiao University(No.14BS310)the National Natural Science Foundation of China(Nos.51276014 and 51476191)
文摘The effect of internal heating source on the film momentum and thermal transport characteristic of thin finite power-law liquids over an accelerating unsteady horizontal stretched interface is studied. Unlike most classical works in this field, a general surface temperature distribution of the liquid film and the generalized Fourier's law for varying thermal conductivity are taken into consideration. Appropriate similarity transformations are used to convert the strongly nonlinear governing partial differential equations (PDEs) into a boundary value problem with a group of two-point ordinary differential equations (ODEs). The correspondence between the liquid film thickness and the unsteadiness parameter is derived with the BVP4C program in MATLAB. Numerical solutions to the self-similarity ODEs are obtained using the shooting technique combined with a Runge-Kutta iteration program and Newton's scheme. The effects of the involved physical parameters on the fluid's horizontal velocity and temperature distribution are presented and discussed.
文摘The inverse heat conduction problem (IHCP) is a severely ill-posed problem in the sense that the solution ( if it exists) does not depend continuously on the data. But now the results on inverse heat conduction problem are mainly devoted to the standard inverse heat conduction problem. Some optimal error bounds in a Sobolev space of regularized approximation solutions for a sideways parabolic equation, i. e. , a non-standard inverse heat conduction problem with convection term which appears in some applied subject are given.
文摘The present paper deals with the determination of quasi-static thermal stresses due to an instantaneous point heat source of strength g_(pi) situated at certain circle along the radial direction of the circular plate and releasing its heat spontaneously at time t=τ.A circular plate is considered having arbitrary initial temperature and subjected to time dependent heat flux at the fixed circular boundary of r=b.The governing heat conduction equation is solved by using the integral transform method,and results are obtained in series form in terms of Bessel functions.The mathematical model has been constructed for copper material and the thermal stresses are discussed graphically.
文摘Heat induced in a finite Silver Selenide slab by an external pulsed laser source is studied in dimensionless scale according to hyperbolic heat conduction model (HHCE) using Laplace integral transform technique. The temperature profile, the critical time required to initiate phase transition and that to initiate damage at the front surface are obtained for different pulses and are illustrated graphically.