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.展开更多
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.展开更多
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.展开更多
EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, t...EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2) one order higher than its interpolation error O(h), the superclose results of order O(h2) in broken Hi-norm are obtained. At the same time, the global superconvergence in broken Hi-norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQrot element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.展开更多
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.展开更多
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.展开更多
This paper has solved the Chester modified heat conduction equation of the different relaxation time r value under different temperature conditions, different boundary conditions and the different initial conditions b...This paper has solved the Chester modified heat conduction equation of the different relaxation time r value under different temperature conditions, different boundary conditions and the different initial conditions by different means of methods. These solutions can help to obtain temperature field of laser thermal effects.展开更多
Soil temperature is a key variable in the control of underground hydro-thermal processes. To estimate soil temperature more accurately, this study proposed a solution method of the heat conduction equation of soil tem...Soil temperature is a key variable in the control of underground hydro-thermal processes. To estimate soil temperature more accurately, this study proposed a solution method of the heat conduction equation of soil temperature (improved heat conduction model) by applying boundary conditions that incorporate the annual and diurnal variations of soil surface temperature and the temporal variation of daily temperature amplitude, as well as the temperature difference between two soil layers in the Tanggula observation site of the Qinghai-Tibet Plateau of China. We employed both the improved heat conduction model and the classical heat conduction model to fit soil temperature by using the 5 cm soil layer as the upper boundary for soil depth. The results indicated that the daily soil temperature amplitude can be better described by the sinusoidal function in the improved model, which then yielded more accurate soil temperature simulating effect at the depth of 5 cm. The simulated soil temperature values generated by the improved model and classical heat conduction model were then compared to the observed soil temperature values at different soil depths. Statistical analyses of the root mean square error (RMSE), the normalized standard error (NSEE) and the bias demonstrated that the improved model showed higher accuracy, and the average values of RMSE, bias and NSEE at the soil depth of 10-105 cm were 1.41℃, 1.15℃ and 22.40%, respectively. These results indicated that the improved heat conduction model can better estimate soil temperature profiles compared to the traditional model.展开更多
This paper deals with a parabolic system in a multi dimentional bounded domain ΩR n with the smooth boundary Ω. We discuss an inverse parabolic problem of determining the indirectly measurable internal heat distri...This paper deals with a parabolic system in a multi dimentional bounded domain ΩR n with the smooth boundary Ω. We discuss an inverse parabolic problem of determining the indirectly measurable internal heat distribution at any intermediate moment from the heat distribution measurements in arbitrary accessible subdomain ωΩ at some time interval. Our main result is the Hlder stability estimate in the inverse problem and the proof is completed with a Carleman estimate and a eigenfunction expansion for parabolic equations.展开更多
The automated energy generating complex allows obtaining heat energy from waste coal-water slurry fuel (WCF) that is a mixture of fine coal particles from coal enrichment wastes with water. The mixture is blown into...The automated energy generating complex allows obtaining heat energy from waste coal-water slurry fuel (WCF) that is a mixture of fine coal particles from coal enrichment wastes with water. The mixture is blown into the swirl chamber under the pressure through the special sprayers. The received heat energy is used in different ways. One of the important issues is to estimate the heat losses through the walls of this chamber. In this paper we solved the boundary problem of mathematical physics to estimate the temperature fields in the walls of the swirl chamber. The obtained solution allows us to estimate the heat losses through the waUs of the swid chamber. The task of the mathematical physics has been solved by a numerical finite-difference method. The method for solving this prob- lem can be used in the calculation of temperature fields and evaluation of heat losses in other thermal power units.展开更多
Lunar equatorial regolith temperature profiles were simulated using the half-limited solid heat conduction model. Based on the infrared data measured using the Diviner radiometer on the Lunar Reconnaissance Orbiter la...Lunar equatorial regolith temperature profiles were simulated using the half-limited solid heat conduction model. Based on the infrared data measured using the Diviner radiometer on the Lunar Reconnaissance Orbiter launched by the United Sates in June 2009, three factors influencing temperature profiles were analyzed. The infrared brightness temperature data from Diviner channel 7 were used to retrieve surface temperature. In simulating regolith temperature profiles, the retrieved temperature, rather than temperatures calculated from solar radiance at the lunar surface, were used as the input for surface temperature in solving the heat-conductive equation. The results showed that the bottom-layer temperature at depths of 6 m approached almost 246 K after 10000 iterations. The temperature was different to the temperature of 250 K at the same depth encountered in simulations using solar radiance. Simulations from both methods of surface temperatures over a lunar day gave similar variations. At lunar night, the temperature difference between the two was about 2 K; the main differences occurred when the solar elevation angle was very low when surface temperatures are largely affected by terrain topography. With no certainty in lunar temperature profiles at present, the advantage of the retrieval method using infrared sensor data as input to the boundary conditions in solving the lunar heat conduction equation is that simulations of surface temperature variations are more accurate. This is especially true in areas with large variations in terrain topography, where surface temperatures vary greatly because of shading from the sunlight.展开更多
基金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.
文摘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.
文摘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.
基金Supported by the National Natural Science Foundation of China (Nos. 10971203 11101381)+3 种基金Tianyuan Mathe-matics Foundation of National Natural Science Foundation of China (No. 11026154)Natural Science Foundation of Henan Province (No. 112300410026)Natural Science Foundation of the Education Department of Henan Province (Nos. 2011A110020 12A110021)
文摘EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2) one order higher than its interpolation error O(h), the superclose results of order O(h2) in broken Hi-norm are obtained. At the same time, the global superconvergence in broken Hi-norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQrot element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.
文摘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.
文摘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.
基金This work was supported by the National Natural Science Foundation of China(No.60068001)and the Natural Science Foundation of Yunnan Province(No.2000A0021M)and ESF of Yunnan(No.0111054).
文摘This paper has solved the Chester modified heat conduction equation of the different relaxation time r value under different temperature conditions, different boundary conditions and the different initial conditions by different means of methods. These solutions can help to obtain temperature field of laser thermal effects.
基金financially supported by the National Basic Research Program of China(2013CBA01803)the key project of the Chinese Academy of Sciences(KJZD-EW-G03-02)+4 种基金the National Natural Science Foundation of China(4127108141271086)the One Hundred Talent Program of the Chinese Academy of Sciences(51Y551831)the Natural Science Foundation of Gansu Province(1308RJZA309)the support of the West Light Foundation of the Chinese Academy of Sciences
文摘Soil temperature is a key variable in the control of underground hydro-thermal processes. To estimate soil temperature more accurately, this study proposed a solution method of the heat conduction equation of soil temperature (improved heat conduction model) by applying boundary conditions that incorporate the annual and diurnal variations of soil surface temperature and the temporal variation of daily temperature amplitude, as well as the temperature difference between two soil layers in the Tanggula observation site of the Qinghai-Tibet Plateau of China. We employed both the improved heat conduction model and the classical heat conduction model to fit soil temperature by using the 5 cm soil layer as the upper boundary for soil depth. The results indicated that the daily soil temperature amplitude can be better described by the sinusoidal function in the improved model, which then yielded more accurate soil temperature simulating effect at the depth of 5 cm. The simulated soil temperature values generated by the improved model and classical heat conduction model were then compared to the observed soil temperature values at different soil depths. Statistical analyses of the root mean square error (RMSE), the normalized standard error (NSEE) and the bias demonstrated that the improved model showed higher accuracy, and the average values of RMSE, bias and NSEE at the soil depth of 10-105 cm were 1.41℃, 1.15℃ and 22.40%, respectively. These results indicated that the improved heat conduction model can better estimate soil temperature profiles compared to the traditional model.
文摘This paper deals with a parabolic system in a multi dimentional bounded domain ΩR n with the smooth boundary Ω. We discuss an inverse parabolic problem of determining the indirectly measurable internal heat distribution at any intermediate moment from the heat distribution measurements in arbitrary accessible subdomain ωΩ at some time interval. Our main result is the Hlder stability estimate in the inverse problem and the proof is completed with a Carleman estimate and a eigenfunction expansion for parabolic equations.
基金the project No.2010-218-02-174 by the Governmental Order of the Russian Federation of April 9,2010 No.218'On measures of federal support of cooperation between higher educational institutions and enterprises which realize complex projects of hi-tech production development
文摘The automated energy generating complex allows obtaining heat energy from waste coal-water slurry fuel (WCF) that is a mixture of fine coal particles from coal enrichment wastes with water. The mixture is blown into the swirl chamber under the pressure through the special sprayers. The received heat energy is used in different ways. One of the important issues is to estimate the heat losses through the walls of this chamber. In this paper we solved the boundary problem of mathematical physics to estimate the temperature fields in the walls of the swirl chamber. The obtained solution allows us to estimate the heat losses through the waUs of the swid chamber. The task of the mathematical physics has been solved by a numerical finite-difference method. The method for solving this prob- lem can be used in the calculation of temperature fields and evaluation of heat losses in other thermal power units.
文摘Lunar equatorial regolith temperature profiles were simulated using the half-limited solid heat conduction model. Based on the infrared data measured using the Diviner radiometer on the Lunar Reconnaissance Orbiter launched by the United Sates in June 2009, three factors influencing temperature profiles were analyzed. The infrared brightness temperature data from Diviner channel 7 were used to retrieve surface temperature. In simulating regolith temperature profiles, the retrieved temperature, rather than temperatures calculated from solar radiance at the lunar surface, were used as the input for surface temperature in solving the heat-conductive equation. The results showed that the bottom-layer temperature at depths of 6 m approached almost 246 K after 10000 iterations. The temperature was different to the temperature of 250 K at the same depth encountered in simulations using solar radiance. Simulations from both methods of surface temperatures over a lunar day gave similar variations. At lunar night, the temperature difference between the two was about 2 K; the main differences occurred when the solar elevation angle was very low when surface temperatures are largely affected by terrain topography. With no certainty in lunar temperature profiles at present, the advantage of the retrieval method using infrared sensor data as input to the boundary conditions in solving the lunar heat conduction equation is that simulations of surface temperature variations are more accurate. This is especially true in areas with large variations in terrain topography, where surface temperatures vary greatly because of shading from the sunlight.