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.展开更多
The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element ...The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.展开更多
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.展开更多
Guilin rice noodles, a unique cuisine from Guilin, Guangxi, is renowned both domestically and internationally as one of the top ten “Guilin Classics”. Utilizing a heat conduction model, this study explores the effec...Guilin rice noodles, a unique cuisine from Guilin, Guangxi, is renowned both domestically and internationally as one of the top ten “Guilin Classics”. Utilizing a heat conduction model, this study explores the effectiveness of the cooking process in sterilizing Guilin rice noodles before consumption. The model assumes that a large pot is filled with boiling water which is maintained at a constant high temperature heat resource through continuous gentle heating. And the room temperature is set as the initial temperature for the preheating process and the final temperature for the cooling process. The objective is to assess whether the cooking process achieves satisfactory sterilization results. The temperature distribution function of rice noodle with time is analytically obtained using the separation of variables method in the three-dimensional cylindrical coordinate system. Meanwhile, the thermal diffusion coefficient of Guilin rice noodles is obtained in terms of Riedel’ theory. By analyzing the elimination characteristics of Pseudomonas cocovenenans subsp. farinofermentans, this study obtains the optimal time required for effective sterilization at the core of Guilin rice noodles. The results show that the potential Pseudomonas cocovenenans subsp. farinofermentans will be completely eliminated through continuously preheating more than 31 seconds during the cooking process before consumption. This study provides a valuable reference of food safety standards in the cooking process of Guilin rice noodles, particularly in ensuring the complete inactivation of potentially harmful strains such as Pseudomonas cocovenenans subsp. farinofermentans.展开更多
With rapid developments in the field of very large-scale integrated circuits,heat dissipation has emerged as a significant factor that restricts the high-density integration of chips.Due to their high thermal conducti...With rapid developments in the field of very large-scale integrated circuits,heat dissipation has emerged as a significant factor that restricts the high-density integration of chips.Due to their high thermal conductivity and low thermal expansion coefficient,diamond/Cu composites have attracted considerable attention as a promising thermal management material.In this study,a surface tungsten carbide gradient layer coating of diamond particles has been realized using comprehensive magnetron sputtering technology and a heat treatment process.Diamond/Cu composites were prepared using high-temperature and high-pressure technology.The results show that,by adjusting the heat treatment process,tungsten carbide and di-tungsten carbide are generated by an in situ reaction at the tungsten–diamond interface,and W–WC–W_(2)C gradient layer-coated diamond particles were obtained.The diamond/Cu composites were sintered by high-temperature and high-pressure technology,and the density of surface-modified diamond/Cu composites was less than 4 g cm^(-3).The W–WC–W_(2)C@diamond/Cu composites have a thermal diffusivity as high as 331 mm^(2)s^(-1),and their thermal expansion coefficient is as low as 1.76×10^(-6)K^(-1).The interface coherent structure of the gradient layer-coated diamond/copper composite can effectively improve the interface heat transport efficiency.展开更多
A 1D radially self-consistent model in helicon plasmas has been established to investigate the influence of radial heat conduction on plasma transport and wave propagation.Two kinds of 1D radial fluid models,with and ...A 1D radially self-consistent model in helicon plasmas has been established to investigate the influence of radial heat conduction on plasma transport and wave propagation.Two kinds of 1D radial fluid models,with and without considering heat conduction,have been developed to couple the 1D plasma-wave interaction model,and self-consistent solutions have been obtained.It is concluded that in the low magnetic field range the radial heat conduction plays a moderate role in the transport of helicon plasmas and the importance depends on the application of the helicon source.It influences the local energy balance leading to enhancement of the electron temperature in the bulk region and a decrease in plasma density.The power deposition in the plasma is mainly balanced by collisional processes and axial diffusion,whereas it is compensated by heat conduction in the bulk region and consumed near the boundary.The role of radial heat conduction in the large magnetic field regime becomes negligible and the two fluid models show consistency.The local power balance,especially near the wall,is improved when conductive heat is taken into account.展开更多
This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by...This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.展开更多
Induction brazecoating technology is an important means to improve the surface properties of materials.In this paper,copper plate and corundum are selected as substrates for induction brazecoating respectively.The tem...Induction brazecoating technology is an important means to improve the surface properties of materials.In this paper,copper plate and corundum are selected as substrates for induction brazecoating respectively.The temperature variation of powder and paste coating is systematically studied,and the heat conduction mode and path in the brazecoating process are analyzed.The results show that rise of the coating temperature mainly depends on the heat absorption from the substrate.The liquid-solid interface conducts heat violently and advances step by step,which promotes the melting spread of metal filler metal.The powdery brazecoating material is in a free state,and there is a gas insulation film between the powder particles and the diamond,making it difficult for substrate to conduct heat to the coating.The binder not only assists forming,shell making and oxygen isolation,but also plays an important role during melting.展开更多
A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problem...A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problems in single and multi-layer materials. In the proposed approach, a new variational functional is firstly constructed for the proposed HFS-FE model and the related existence of extremum is presented. Then, the assumed internal potential field constructed by the linear combination of fundamental solutions at points outside the elemental domain under consideration is used as the internal interpolation function, which analytically satisfies the governing equation within each element. As a result, the domain integrals in the variational functional formulation can be converted into the boundary integrals which can significantly simplify the calculation of the element stiffness matrix. The independent frame field is also introduced to guarantee the inter-element continuity and the stationary condition of the new variational functional is used to obtain the final stiffness equations. The proposed method inherits the advantages of the hybrid Trefftz finite element method (HT-FEM) over the conventional finite element method (FEM) and boundary element method (BEM), and avoids the difficulty in selecting appropriate terms of T-complete functions used in HT-FEM, as the fundamental solutions contain usually one term only, rather than a series containing infinitely many terms. Further, the fundamental solutions of a problem are, in general, easier to derive than the T-complete functions of that problem. Finally, several examples are presented to assess the performance of the proposed method, and the obtained numerical results show good numerical accuracy and remarkable insensitivity to mesh distortion.展开更多
In order to simulate thermal strains, thermal stresses, residual stresses and microstructure of the steel during gas quenching by means of the numerical method, it is necessary to obtain an accurate boundary condition...In order to simulate thermal strains, thermal stresses, residual stresses and microstructure of the steel during gas quenching by means of the numerical method, it is necessary to obtain an accurate boundary condition of temperature field. The surface heat transfer coefficient is a key parameter. The explicit finite difference method, nonlinear estimation method and the experimental relation between temperature and time during gas quenching have been used to solve the inverse problem of heat conduction. The relationship between surface temperature and surface heat transfer coefficient of a cylinder has been given. The nonlinear surface heat transfer coefficients include the coupled effects between martensitic phase transformation and temperature.展开更多
A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is pres...In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.展开更多
A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue...A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue of the level-cut method, the difference discrete equations with fuzzy parameters are equivalently transformed into groups of interval equations. New stability analysis theory suited to fuzzy difference schemes is developed. Based on the parameter perturbation method, the interval ranges of the uncertain temperature field can be approximately predicted. Subsequently, fuzzy solutions to the original difference equations are obtained by the fuzzy resolution theorem. Two numerical examples are given to demonstrate the feasibility and efficiency of the presented method for solving both steady-state and transient heat conduction problems.展开更多
This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structu...This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.展开更多
Phase-field modeling approach has been used to study the oxidation behavior of pure Ni when considering heat conduction. In this calculation, the dependence of the coefficient of the Cahn–Hilliard equation Lc on the ...Phase-field modeling approach has been used to study the oxidation behavior of pure Ni when considering heat conduction. In this calculation, the dependence of the coefficient of the Cahn–Hilliard equation Lc on the temperature T was considered. To this end, high-temperature oxidation experiments and phase-field modeling for pure Ni were performed in air under atmospheric pressure at 600,700, and 800?C. The oxidation rate was measured by thermogravimetry and Lc at these temperatures was determined via interactive algorithm. With the Lc-T relationship constructed, oxidation behavior of Ni when considering heat conduction was investigated. The influence of temperature boundaries on the oxidation degree, oxide film thickness, and specific weight gain were discussed. The phase-field modeling approach proposed in this study will give some highlights of the oxidation resistance analysis and cooling measures design of thermal protection materials.展开更多
Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attr...Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson's equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.展开更多
We build a fractional dual-phase-lag model and the corresponding bioheat transfer equation, which we use to interpret the experiment results for processed meat that have been explained by applying the hyperbolic condu...We build a fractional dual-phase-lag model and the corresponding bioheat transfer equation, which we use to interpret the experiment results for processed meat that have been explained by applying the hyperbolic conduction. Analytical solutions expressed by H-functions are obtained by using the Laplace and Fourier transforms method. The inverse fractional dual-phase-lag heat conduction problem for the simultaneous estimation of two relaxation times and orders of fractionality is solved by applying the nonlinear least-square method. The estimated model parameters are given. Finally, the measured and the calculated temperatures versus time are compared and discussed. Some numerical examples are also given and discussed.展开更多
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.展开更多
In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE metho...In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.展开更多
This paper presents two exact explicit solutions for the three dimensional dual-phase lag (DLP) heat conduction equation, during the derivation of which the method of trial and error and the authors' previous exper...This paper presents two exact explicit solutions for the three dimensional dual-phase lag (DLP) heat conduction equation, during the derivation of which the method of trial and error and the authors' previous experiences are utilized. To the authors' knowledge, most solutions of 2D or 3D DPL models available in the literature are obtained by numerical methods, and there are few exact solutions up to now. The exact solutions in this paper can be used as benchmarks to validate numerical solutions and to develop numerical schemes, grid generation methods and so forth. In addition, they are of theoretical significance since they correspond to physically possible situations. The main goal of this paper is to obtain some possible exact explicit solutions of the dual-phase lag heat conduction equation as the benchmark solutions for computational heat transfer, rather than specific solutions for some given initial and boundary conditions. Therefore, the initial and boundary conditions are indeterminate before derivation and can be deduced from the solutions afterwards. Actually, all solutions given in this paper can be easily proven by substituting them into the governing equation.展开更多
文摘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.
基金Project supported by the National Natural Science Foundation of China (Nos. 12102043, 12072375U2241240)the Natural Science Foundation of Hunan Province of China (Nos. 2023JJ40698 and 2021JJ40710)。
文摘The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.
文摘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.
文摘Guilin rice noodles, a unique cuisine from Guilin, Guangxi, is renowned both domestically and internationally as one of the top ten “Guilin Classics”. Utilizing a heat conduction model, this study explores the effectiveness of the cooking process in sterilizing Guilin rice noodles before consumption. The model assumes that a large pot is filled with boiling water which is maintained at a constant high temperature heat resource through continuous gentle heating. And the room temperature is set as the initial temperature for the preheating process and the final temperature for the cooling process. The objective is to assess whether the cooking process achieves satisfactory sterilization results. The temperature distribution function of rice noodle with time is analytically obtained using the separation of variables method in the three-dimensional cylindrical coordinate system. Meanwhile, the thermal diffusion coefficient of Guilin rice noodles is obtained in terms of Riedel’ theory. By analyzing the elimination characteristics of Pseudomonas cocovenenans subsp. farinofermentans, this study obtains the optimal time required for effective sterilization at the core of Guilin rice noodles. The results show that the potential Pseudomonas cocovenenans subsp. farinofermentans will be completely eliminated through continuously preheating more than 31 seconds during the cooking process before consumption. This study provides a valuable reference of food safety standards in the cooking process of Guilin rice noodles, particularly in ensuring the complete inactivation of potentially harmful strains such as Pseudomonas cocovenenans subsp. farinofermentans.
基金National Natural Science Foundation of China(Grant No.52072327)the China National Key R&D Program(2021YFB3701802)+6 种基金Scientific and Technological Projects of Henan Province(No.232102231050)the Higher Education and Teaching Reformation Project(2014SJGLX064)the Project for Work-station of Zhongyuan scholars of Henan Province(Nos.214400510002,224400510023)the Science and Technology Major Project of Henan Province(No.221100230300)the Postgraduate Education Reform and QualityAcademic Degrees&Graduate Education Reform Project of Henan Province(No.2021SJGLX060Y)the Postgraduate Education Reform and Quality Improvement Project of Henan Province(No.YJS2022JD34)the Science and Technology on Plasma Physics Laboratory(Grant No.JCKYS2021212010).
文摘With rapid developments in the field of very large-scale integrated circuits,heat dissipation has emerged as a significant factor that restricts the high-density integration of chips.Due to their high thermal conductivity and low thermal expansion coefficient,diamond/Cu composites have attracted considerable attention as a promising thermal management material.In this study,a surface tungsten carbide gradient layer coating of diamond particles has been realized using comprehensive magnetron sputtering technology and a heat treatment process.Diamond/Cu composites were prepared using high-temperature and high-pressure technology.The results show that,by adjusting the heat treatment process,tungsten carbide and di-tungsten carbide are generated by an in situ reaction at the tungsten–diamond interface,and W–WC–W_(2)C gradient layer-coated diamond particles were obtained.The diamond/Cu composites were sintered by high-temperature and high-pressure technology,and the density of surface-modified diamond/Cu composites was less than 4 g cm^(-3).The W–WC–W_(2)C@diamond/Cu composites have a thermal diffusivity as high as 331 mm^(2)s^(-1),and their thermal expansion coefficient is as low as 1.76×10^(-6)K^(-1).The interface coherent structure of the gradient layer-coated diamond/copper composite can effectively improve the interface heat transport efficiency.
基金National Natural Science Foundation of China(No.51907039)Shenzhen Technology Project(Nos.JCYJ20190806142603534 and ZDSYS201707280904031)+1 种基金ESPEOS project(No.PID2019108034RB-I00/AEI/10.13039/501100011033)funded by the Agencia Estatal de Investigacion(Spanish National Research Agency)。
文摘A 1D radially self-consistent model in helicon plasmas has been established to investigate the influence of radial heat conduction on plasma transport and wave propagation.Two kinds of 1D radial fluid models,with and without considering heat conduction,have been developed to couple the 1D plasma-wave interaction model,and self-consistent solutions have been obtained.It is concluded that in the low magnetic field range the radial heat conduction plays a moderate role in the transport of helicon plasmas and the importance depends on the application of the helicon source.It influences the local energy balance leading to enhancement of the electron temperature in the bulk region and a decrease in plasma density.The power deposition in the plasma is mainly balanced by collisional processes and axial diffusion,whereas it is compensated by heat conduction in the bulk region and consumed near the boundary.The role of radial heat conduction in the large magnetic field regime becomes negligible and the two fluid models show consistency.The local power balance,especially near the wall,is improved when conductive heat is taken into account.
基金supported by the NationalNatural Science Foundation of China (No.11802151)the Natural Science Foundation of Shandong Province of China (No.ZR2019BA008)the China Postdoctoral Science Foundation (No.2019M652315).
文摘This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.
基金supported by the 2020 Ningbo “3315 Talent Introduction Plan” Innovative Team (C-Class)the major project of Ningbo “Scientific and Technological Innovation 2025”(Grant No.2020Z111)Science and Technology Major Project of Zhejiang Province (No.203ZP20220161)
文摘Induction brazecoating technology is an important means to improve the surface properties of materials.In this paper,copper plate and corundum are selected as substrates for induction brazecoating respectively.The temperature variation of powder and paste coating is systematically studied,and the heat conduction mode and path in the brazecoating process are analyzed.The results show that rise of the coating temperature mainly depends on the heat absorption from the substrate.The liquid-solid interface conducts heat violently and advances step by step,which promotes the melting spread of metal filler metal.The powdery brazecoating material is in a free state,and there is a gas insulation film between the powder particles and the diamond,making it difficult for substrate to conduct heat to the coating.The binder not only assists forming,shell making and oxygen isolation,but also plays an important role during melting.
文摘A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problems in single and multi-layer materials. In the proposed approach, a new variational functional is firstly constructed for the proposed HFS-FE model and the related existence of extremum is presented. Then, the assumed internal potential field constructed by the linear combination of fundamental solutions at points outside the elemental domain under consideration is used as the internal interpolation function, which analytically satisfies the governing equation within each element. As a result, the domain integrals in the variational functional formulation can be converted into the boundary integrals which can significantly simplify the calculation of the element stiffness matrix. The independent frame field is also introduced to guarantee the inter-element continuity and the stationary condition of the new variational functional is used to obtain the final stiffness equations. The proposed method inherits the advantages of the hybrid Trefftz finite element method (HT-FEM) over the conventional finite element method (FEM) and boundary element method (BEM), and avoids the difficulty in selecting appropriate terms of T-complete functions used in HT-FEM, as the fundamental solutions contain usually one term only, rather than a series containing infinitely many terms. Further, the fundamental solutions of a problem are, in general, easier to derive than the T-complete functions of that problem. Finally, several examples are presented to assess the performance of the proposed method, and the obtained numerical results show good numerical accuracy and remarkable insensitivity to mesh distortion.
基金This work has been supported by the National Natural Science Foundation of China (10162002) and Foundation for University Key Teacher by the Ministry of Education and The Yunnan Foundation of Natural Science (1999A0023M).
文摘In order to simulate thermal strains, thermal stresses, residual stresses and microstructure of the steel during gas quenching by means of the numerical method, it is necessary to obtain an accurate boundary condition of temperature field. The surface heat transfer coefficient is a key parameter. The explicit finite difference method, nonlinear estimation method and the experimental relation between temperature and time during gas quenching have been used to solve the inverse problem of heat conduction. The relationship between surface temperature and surface heat transfer coefficient of a cylinder has been given. The nonlinear surface heat transfer coefficients include the coupled effects between martensitic phase transformation and temperature.
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)the Innovation Fund for Graduate Student of Shanghai University of China (Grant No.SHUCX120125)
文摘In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.
基金supported by the National Special Fund for Major Research Instrument Development(2011YQ140145)111 Project(B07009)+1 种基金the National Natural Science Foundation of China(11002013)Defense Industrial Technology Development Program(A2120110001 and B2120110011)
文摘A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue of the level-cut method, the difference discrete equations with fuzzy parameters are equivalently transformed into groups of interval equations. New stability analysis theory suited to fuzzy difference schemes is developed. Based on the parameter perturbation method, the interval ranges of the uncertain temperature field can be approximately predicted. Subsequently, fuzzy solutions to the original difference equations are obtained by the fuzzy resolution theorem. Two numerical examples are given to demonstrate the feasibility and efficiency of the presented method for solving both steady-state and transient heat conduction problems.
基金funded by National Natural Science Foundation of China(NSFC)under Grant Nos.11702238,51904202,and 11902212Nanhu Scholars Program for Young Scholars of XYNU.
文摘This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.
基金supported by the Beijing Jiaotong University (Grant C15JB00080)
文摘Phase-field modeling approach has been used to study the oxidation behavior of pure Ni when considering heat conduction. In this calculation, the dependence of the coefficient of the Cahn–Hilliard equation Lc on the temperature T was considered. To this end, high-temperature oxidation experiments and phase-field modeling for pure Ni were performed in air under atmospheric pressure at 600,700, and 800?C. The oxidation rate was measured by thermogravimetry and Lc at these temperatures was determined via interactive algorithm. With the Lc-T relationship constructed, oxidation behavior of Ni when considering heat conduction was investigated. The influence of temperature boundaries on the oxidation degree, oxide film thickness, and specific weight gain were discussed. The phase-field modeling approach proposed in this study will give some highlights of the oxidation resistance analysis and cooling measures design of thermal protection materials.
基金supported by the Natural Science Foundation of Ningbo,China (Grant Nos.2009A610014 and 2009A610154)the Natural Science Foundation of Zhejiang Province,China (Grant No.Y6090131)
文摘Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson's equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.
基金supported by the National Natural Science Foundation of China(Grant Nos.11102102,11472161,and 91130017)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2014AQ015)the Independent Innovation Foundation of Shandong University,China(Grant No.2013ZRYQ002)
文摘We build a fractional dual-phase-lag model and the corresponding bioheat transfer equation, which we use to interpret the experiment results for processed meat that have been explained by applying the hyperbolic conduction. Analytical solutions expressed by H-functions are obtained by using the Laplace and Fourier transforms method. The inverse fractional dual-phase-lag heat conduction problem for the simultaneous estimation of two relaxation times and orders of fractionality is solved by applying the nonlinear least-square method. The estimated model parameters are given. Finally, the measured and the calculated temperatures versus time are compared and discussed. Some numerical examples are also given and discussed.
基金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.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Special Fund for Basic Scientific Research of Central Colleges of Chang’an University, China (Grant No. CHD2011JC080)
文摘In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.
基金supported by the National Natural Science Foundation of China (50576097) the National Defense Basic Research Program of China (DEDP 1003)
文摘This paper presents two exact explicit solutions for the three dimensional dual-phase lag (DLP) heat conduction equation, during the derivation of which the method of trial and error and the authors' previous experiences are utilized. To the authors' knowledge, most solutions of 2D or 3D DPL models available in the literature are obtained by numerical methods, and there are few exact solutions up to now. The exact solutions in this paper can be used as benchmarks to validate numerical solutions and to develop numerical schemes, grid generation methods and so forth. In addition, they are of theoretical significance since they correspond to physically possible situations. The main goal of this paper is to obtain some possible exact explicit solutions of the dual-phase lag heat conduction equation as the benchmark solutions for computational heat transfer, rather than specific solutions for some given initial and boundary conditions. Therefore, the initial and boundary conditions are indeterminate before derivation and can be deduced from the solutions afterwards. Actually, all solutions given in this paper can be easily proven by substituting them into the governing equation.