Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using t...Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.展开更多
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.展开更多
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.展开更多
To make heat conduction equation embody the essence of physical phenomenon under study, dimensionless factors were introduced and the transient heat conduction equation and its boundary conditions were transformed to ...To make heat conduction equation embody the essence of physical phenomenon under study, dimensionless factors were introduced and the transient heat conduction equation and its boundary conditions were transformed to dimensionless forms. Then, a theoretical solution model of transient heat conduction problem in one-dimensional double-layer composite medium was built utilizing the natural eigenfunction expansion method. In order to verify the validity of the model, the results of the above theoretical solution were compared with those of finite element method. The results by the two methods are in a good agreement. The maximum errors by the two methods appear when τ（τ is nondimensional time） equals 0.1 near the boundaries of ζ =1 （ζ is nondimensional space coordinate） and ζ =4. As τ increases, the error decreases gradually, and when τ =5 the results of both solutions have almost no change with the variation of coordinate 4.展开更多
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 a hybrid graded element model for the transient heat conduction problem in functionally graded materials （FGMs）. First, a Laplace transform approach is used to handle the time variable. Then, a f...This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials （FGMs）. First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest＇s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.展开更多
A meshless numerical model is developed for analyzing transient heat conductions in three-dimensional （3D） axisymmetric continuously nonhomogeneous functionally graded materials （FGMs）. Axial symmetry of geometry ...A meshless numerical model is developed for analyzing transient heat conductions in three-dimensional （3D） axisymmetric continuously nonhomogeneous functionally graded materials （FGMs）. Axial symmetry of geometry and boundary conditions reduces the original 3D initial-boundary value problem into a two-dimensional （2D） problem. Local weak forms are derived for small polygonal sub-domains which surround nodal points distributed over the cross section. In order to simplify the treatment of the essential boundary conditions, spatial variations of the temperature and heat flux at discrete time instants are interpolated by the natural neighbor interpolation. Moreover, the using of three-node triangular finite element method （FEM） shape functions as test functions reduces the orders of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with the traditional two-point difference technique in the time domain. Two numerical examples are investigated and excellent results are obtained, demonstrating the potential application of the proposed approach.展开更多
Uniform heating of complex surfaces,especially non-developable surfaces,is a crucial problem that traditional rigid heaters cannot solve.Inspired by flexible electronic devices,a novel design for the stretchable heati...Uniform heating of complex surfaces,especially non-developable surfaces,is a crucial problem that traditional rigid heaters cannot solve.Inspired by flexible electronic devices,a novel design for the stretchable heating film is proposed with the flexible serpentine wire embedded in the soft polymer film,which can be attached to non-developable surfaces conformally.It provides a new way for the stretchable heaters to realize uniform heating of complex surfaces.However,the thermal field of flexible serpentine heaters(FSHs)depends on the configurations of the embedded serpentine heating wire,which requires accurate theoretical prediction of real-time temperature distribution.Therefore,the analytical model for the transient heat conduction in FSHs is solved by the separation of variables method and validated by the finite element analysis(FEA)in this paper.Based on this model,the effects of the geometric parameters,such as the radius and the length of the serpentine heaters,on the thermal uniformity are systematically investigated.This study can help to design and fabricate flexible heaters with uniform heating in the future.展开更多
By introducing the radial basis functions(RBFs)into the reproducing kernel particle method(RKPM),the calculating accuracy and stability of the RKPM can be improved,and a novel meshfree method of the radial basis RKPM(...By introducing the radial basis functions(RBFs)into the reproducing kernel particle method(RKPM),the calculating accuracy and stability of the RKPM can be improved,and a novel meshfree method of the radial basis RKPM(meshfree RRKPM)is proposed.Meanwhile,the meshfree RRKPM is applied to transient heat conduction problems(THCP),and the corresponding equations of the meshfree RRKPM for the THCP are derived.The twopoint time difference scheme is selected to discretize the time of the THCP.Finally,the numerical results illustrate the effectiveness of the meshfree RRKPM for the THCP.展开更多
The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the ...The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.展开更多
On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is ...On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.展开更多
With the improved moving least-squares (IMLS) approximation, an orthogonal function system with a weight function is used as the basis function. The combination of the element-free Galerkin (EFG) method and the IMLS a...With the improved moving least-squares (IMLS) approximation, an orthogonal function system with a weight function is used as the basis function. The combination of the element-free Galerkin (EFG) method and the IMLS approximation leads to the development of the improved element-free Galerkin (IEFG) method. In this paper, the IEFG method is applied to study the partial differential equations that control the heat flow in three-dimensional space. With the IEFG technique, the Galerkin weak form is employed to develop the discretized system equations, and the penalty method is applied to impose the essential boundary conditions. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the transient heat conduction equations and the boundary and initial conditions are time dependent, the scaling parameter, number of nodes and time step length are considered in a convergence study.展开更多
This paper describes a new two-dimensional(2-D)control volume finite element method(CV-FEM)for transient heat conduction in multilayer functionally graded materials(FGMs).To deal with the mixed-grid problem,9-node qua...This paper describes a new two-dimensional(2-D)control volume finite element method(CV-FEM)for transient heat conduction in multilayer functionally graded materials(FGMs).To deal with the mixed-grid problem,9-node quadrilateral grids and 6-node triangular grids are used.The unknown temperature and material properties are stored at the node.By using quadratic triangular grids and quadratic quadrilateral grids,the present method offers greater geometric flexibility and the potential for higher accuracy than the linear CV-FEM.The properties of the FGMs are described by exponential,quadratic and trigonometric grading functions.Some numerical tests are studied to demonstrate the performance of the developed method.First,the present CV-FEM with mixed high-order girds provides a higher accuracy than the linear CV-FEM based on the same grid size.Second,the material properties defined location is proved to have a significant effect on the accuracy of the numerical results.Third,the present method provides better numerical solutions than the conventional FEM for the FGMs in conjunction with course high-order grids.Finally,the present method is also capable of analysis of transient heat conduction in multilayer FGM.展开更多
The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this pap...The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.展开更多
The aim of this paper is the formulation of the finite element method in polar coordinates to solve transient heat conduction problems. It is hard to find in the literature a formulation of the finite element method(F...The aim of this paper is the formulation of the finite element method in polar coordinates to solve transient heat conduction problems. It is hard to find in the literature a formulation of the finite element method(FEM) in polar or cylindrical coordinates for the solution of heat transfer problems. This document shows how to apply the most often used boundary conditions. The global equation system is solved by the Crank-Nicolson method. The proposed algorithm is verified in three numerical tests. In the first example, the obtained transient temperature distribution is compared with the temperature obtained from the presented analytical solution. In the second numerical example, the variable boundary condition is assumed. In the last numerical example the component with the shape different than cylindrical is used. All examples show that the introduction of the polar coordinate system gives better results than in the Cartesian coordinate system. The finite element method formulation in polar coordinates is valuable since it provides a higher accuracy of the calculations without compacting the mesh in cylindrical or similar to tubular components. The proposed method can be applied for circular elements such as boiler drums, outlet headers, flux tubes. This algorithm can be useful during the solution of inverse problems, which do not allow for high density grid. This method can calculate the temperature distribution in the bodies of different properties in the circumferential and the radial direction. The presented algorithm can be developed for other coordinate systems. The examples demonstrate a good accuracy and stability of the proposed method.展开更多
Accurately solving transient nonlinear inverse heat conduction problems in complex structures is of great importance to provide key parameters for modeling coupled heat transfer process and the structure’s optimizati...Accurately solving transient nonlinear inverse heat conduction problems in complex structures is of great importance to provide key parameters for modeling coupled heat transfer process and the structure’s optimization design.The finite element method in ABAQUS is employed to solve the direct transient nonlinear heat conduction problem.Improved particle swarm optimization(PSO)method is developed and used to solve the transient nonlinear inverse problem.To investigate the inverse performances,some numerical tests are provided.Boundary conditions at inaccessible surfaces of a scramjet combustor with the regenerative cooling system are inversely identified.The results show that the new methodology can accurately and efficiently determine the boundary conditions in the scramjet combustor with the regenerative cooling system.By solving the transient nonlinear inverse problem,the improved particle swarm optimization for solving the transient nonlinear inverse heat conduction problem in a complex structure is verified.展开更多
Accurate prediction of junction temperature is crucial for the efficient thermal design of silicon nano-devices. In nano-scale semiconductor devices, significant ballistic effects occur due to the mean free path of ph...Accurate prediction of junction temperature is crucial for the efficient thermal design of silicon nano-devices. In nano-scale semiconductor devices, significant ballistic effects occur due to the mean free path of phonons comparable to the heat source size and device scale. We employ a three-dimensional non-gray Monte Carlo simulation to investigate the transient heat conduction of silicon nanofilms with both single and multiple heat sources. The accuracy of the present method is first verified in the ballistic and diffusion limits. When a single local heat source is present, the width of the heat source has a significant impact on heat conduction in the domain. Notably, there is a substantial temperature jump at the boundary when the heat source width is 10 nm.With increasing heat source width, the boundary temperature jump weakens. Furthermore, we observe that the temperature excitation rate is independent of the heat source width, while the temperature influence range expands simultaneously with the increase in heat source width. Around 500 ps, the temperature and heat flux distribution in the domain stabilize. In the case of dual heat sources, the hot zone is broader than that of a single heat source, and the temperature of the hot spot decreases as the heat source spacing increases. However, the mean heat flux remains unaffected. Upon reaching a spacing of 200 nm between the heat sources, the peak temperature in the domain remains unchanged once a steady state is reached. These findings hold significant implications for the thermal design of silicon nano-devices with local heat sources.展开更多
This paper deals with the determination of temperature distribution and thermal deflection function of a thin circular plate with the stated conditions. The transient heat conduction equation is solved by using Marchi...This paper deals with the determination of temperature distribution and thermal deflection function of a thin circular plate with the stated conditions. The transient heat conduction equation is solved by using Marchi-Zgrablich transform and Laplace transform simultaneously and the results of temperature distribution and thermal deflection function are obtained in terms of infinite series of Bessel function and it is solved for special case by using Mathcad 2007 software and represented graphically by using Microsoft excel 2007.展开更多
Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditi...Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.展开更多
The cuckoo search algorithm(CS)is improved by using the conjugate gradient method(CGM),and the CS-CGM is proposed.The unknown inner boundary shapes are generated randomly and evolved by Lévy flights and eliminati...The cuckoo search algorithm(CS)is improved by using the conjugate gradient method(CGM),and the CS-CGM is proposed.The unknown inner boundary shapes are generated randomly and evolved by Lévy flights and elimination mechanism in the CS and CS-CGM.The CS,CGM and CS-CGM are examined for the prediction of a pipe’s inner surface.The direct problem is two-dimensional transient heat conduction in functionally graded materials(FGMs).Firstly,the radial integration boundary element method(RIBEM)is applied to solve the direct problem.Then the three methods are compared to identify the pipe’s inner surface with the information of measured temperatures.Finally,the influences of timepoints,measurement point number and random noise on the inverse results are investigated.It is found that the three algorithms are promising and can be used to identify the pipe’s inner surface.The CS-CGM has higher accuracy and faster convergence speed than the CS and CGM.The CS and CS-CGM are insensitive to the initial values.The CGM and CS-CGM are more insensitive to the measurement noises compared with the CS.With the increase of timepoints and measurement points,and with the decrease of measurement noises,the inverse results are more accurate.展开更多
基金This work was financially supported by the Key Science and Technology Project of Longmen Laboratory(No.LMYLKT-001)Innovation and Entrepreneurship Training Program for College Students of Henan Province(No.202310464050)。
文摘Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.
基金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.
基金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.
基金Projects(50576007,50876016) supported by the National Natural Science Foundation of ChinaProjects(20062180) supported by the National Natural Science Foundation of Liaoning Province,China
文摘To make heat conduction equation embody the essence of physical phenomenon under study, dimensionless factors were introduced and the transient heat conduction equation and its boundary conditions were transformed to dimensionless forms. Then, a theoretical solution model of transient heat conduction problem in one-dimensional double-layer composite medium was built utilizing the natural eigenfunction expansion method. In order to verify the validity of the model, the results of the above theoretical solution were compared with those of finite element method. The results by the two methods are in a good agreement. The maximum errors by the two methods appear when τ（τ is nondimensional time） equals 0.1 near the boundaries of ζ =1 （ζ is nondimensional space coordinate） and ζ =4. As τ increases, the error decreases gradually, and when τ =5 the results of both solutions have almost no change with the variation of coordinate 4.
基金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.
文摘This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials （FGMs）. First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest＇s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.
基金Project supported by the National Natural Science Foundation of China(Grant No.11002054)the Foundation of Hunan Educational Committee(Grant No.12C0059).
文摘A meshless numerical model is developed for analyzing transient heat conductions in three-dimensional （3D） axisymmetric continuously nonhomogeneous functionally graded materials （FGMs）. Axial symmetry of geometry and boundary conditions reduces the original 3D initial-boundary value problem into a two-dimensional （2D） problem. Local weak forms are derived for small polygonal sub-domains which surround nodal points distributed over the cross section. In order to simplify the treatment of the essential boundary conditions, spatial variations of the temperature and heat flux at discrete time instants are interpolated by the natural neighbor interpolation. Moreover, the using of three-node triangular finite element method （FEM） shape functions as test functions reduces the orders of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with the traditional two-point difference technique in the time domain. Two numerical examples are investigated and excellent results are obtained, demonstrating the potential application of the proposed approach.
基金the National Natural Science Foundation of China(No.11772030)the Aeronautical Science Foundation of China(No.2018ZC51030)+1 种基金the Natural Science Foundation of Zhejiang Province of China(No.Y21A020002)the Opening Fund of State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology(No.GZ19117)。
文摘Uniform heating of complex surfaces,especially non-developable surfaces,is a crucial problem that traditional rigid heaters cannot solve.Inspired by flexible electronic devices,a novel design for the stretchable heating film is proposed with the flexible serpentine wire embedded in the soft polymer film,which can be attached to non-developable surfaces conformally.It provides a new way for the stretchable heaters to realize uniform heating of complex surfaces.However,the thermal field of flexible serpentine heaters(FSHs)depends on the configurations of the embedded serpentine heating wire,which requires accurate theoretical prediction of real-time temperature distribution.Therefore,the analytical model for the transient heat conduction in FSHs is solved by the separation of variables method and validated by the finite element analysis(FEA)in this paper.Based on this model,the effects of the geometric parameters,such as the radius and the length of the serpentine heaters,on the thermal uniformity are systematically investigated.This study can help to design and fabricate flexible heaters with uniform heating in the future.
基金supported by the Natural Science Foundation of Shandong Province(Grant No.ZR2017MA028)supported by the Natural Science Foundation of Shandong Province(Grant No.ZR2020MA059).
文摘By introducing the radial basis functions(RBFs)into the reproducing kernel particle method(RKPM),the calculating accuracy and stability of the RKPM can be improved,and a novel meshfree method of the radial basis RKPM(meshfree RRKPM)is proposed.Meanwhile,the meshfree RRKPM is applied to transient heat conduction problems(THCP),and the corresponding equations of the meshfree RRKPM for the THCP are derived.The twopoint time difference scheme is selected to discretize the time of the THCP.Finally,the numerical results illustrate the effectiveness of the meshfree RRKPM for the THCP.
基金supported by the Key Program of the National Natural Science Foundation of China (Grand No. 51138001)the China-German Cooperation Project (Grand No. GZ566)+1 种基金the Innovative Research Groups Funded by the National Natural Science Foundation of China (Grand No. 51121005)the Special Funds for the Basic Scientific Research Expenses for the Central University (Grant No. DUT13LK16)
文摘The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.
基金supported by the National Natural Science Foundation of China(Grant No.51078250)the Research Project by Shanxi Scholarship Council of Shanxi Province,China(Grant No.2013-096)the Scientific&Technological Innovation Program for Postgraduates of Taiyuan University of Science and Technology,China(Grant No.20125026)
文摘On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.
基金the National Natural Science Foundation of China (Grant No. 11171208)Shanghai Leading Academic Discipline Project (Grant No. S30106)
文摘With the improved moving least-squares (IMLS) approximation, an orthogonal function system with a weight function is used as the basis function. The combination of the element-free Galerkin (EFG) method and the IMLS approximation leads to the development of the improved element-free Galerkin (IEFG) method. In this paper, the IEFG method is applied to study the partial differential equations that control the heat flow in three-dimensional space. With the IEFG technique, the Galerkin weak form is employed to develop the discretized system equations, and the penalty method is applied to impose the essential boundary conditions. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the transient heat conduction equations and the boundary and initial conditions are time dependent, the scaling parameter, number of nodes and time step length are considered in a convergence study.
基金The financial support from the Fundamental Research Funds for the Central Universities HEUCFP201711
文摘This paper describes a new two-dimensional(2-D)control volume finite element method(CV-FEM)for transient heat conduction in multilayer functionally graded materials(FGMs).To deal with the mixed-grid problem,9-node quadrilateral grids and 6-node triangular grids are used.The unknown temperature and material properties are stored at the node.By using quadratic triangular grids and quadratic quadrilateral grids,the present method offers greater geometric flexibility and the potential for higher accuracy than the linear CV-FEM.The properties of the FGMs are described by exponential,quadratic and trigonometric grading functions.Some numerical tests are studied to demonstrate the performance of the developed method.First,the present CV-FEM with mixed high-order girds provides a higher accuracy than the linear CV-FEM based on the same grid size.Second,the material properties defined location is proved to have a significant effect on the accuracy of the numerical results.Third,the present method provides better numerical solutions than the conventional FEM for the FGMs in conjunction with course high-order grids.Finally,the present method is also capable of analysis of transient heat conduction in multilayer FGM.
基金support by the National Natural Science Foundation of China(grant No.51779033,51409038)the National Key Research and Development Plan(grant No.2016YFB0201001)the National Natural Science Foundation of China(grant No.51421064)
文摘The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.
文摘The aim of this paper is the formulation of the finite element method in polar coordinates to solve transient heat conduction problems. It is hard to find in the literature a formulation of the finite element method(FEM) in polar or cylindrical coordinates for the solution of heat transfer problems. This document shows how to apply the most often used boundary conditions. The global equation system is solved by the Crank-Nicolson method. The proposed algorithm is verified in three numerical tests. In the first example, the obtained transient temperature distribution is compared with the temperature obtained from the presented analytical solution. In the second numerical example, the variable boundary condition is assumed. In the last numerical example the component with the shape different than cylindrical is used. All examples show that the introduction of the polar coordinate system gives better results than in the Cartesian coordinate system. The finite element method formulation in polar coordinates is valuable since it provides a higher accuracy of the calculations without compacting the mesh in cylindrical or similar to tubular components. The proposed method can be applied for circular elements such as boiler drums, outlet headers, flux tubes. This algorithm can be useful during the solution of inverse problems, which do not allow for high density grid. This method can calculate the temperature distribution in the bodies of different properties in the circumferential and the radial direction. The presented algorithm can be developed for other coordinate systems. The examples demonstrate a good accuracy and stability of the proposed method.
基金supported by the National Natural Science Foundation of China(Nos.12172078,51576026)Fundamental Research Funds for the Central Universities in China(No.DUT21LK04)。
文摘Accurately solving transient nonlinear inverse heat conduction problems in complex structures is of great importance to provide key parameters for modeling coupled heat transfer process and the structure’s optimization design.The finite element method in ABAQUS is employed to solve the direct transient nonlinear heat conduction problem.Improved particle swarm optimization(PSO)method is developed and used to solve the transient nonlinear inverse problem.To investigate the inverse performances,some numerical tests are provided.Boundary conditions at inaccessible surfaces of a scramjet combustor with the regenerative cooling system are inversely identified.The results show that the new methodology can accurately and efficiently determine the boundary conditions in the scramjet combustor with the regenerative cooling system.By solving the transient nonlinear inverse problem,the improved particle swarm optimization for solving the transient nonlinear inverse heat conduction problem in a complex structure is verified.
基金supported by the National Natural Science Foundation of China (Grant No. 52076088)the Core Technology Research Project of Shunde District, Foshan, China (Grant No. 2130218002932)。
文摘Accurate prediction of junction temperature is crucial for the efficient thermal design of silicon nano-devices. In nano-scale semiconductor devices, significant ballistic effects occur due to the mean free path of phonons comparable to the heat source size and device scale. We employ a three-dimensional non-gray Monte Carlo simulation to investigate the transient heat conduction of silicon nanofilms with both single and multiple heat sources. The accuracy of the present method is first verified in the ballistic and diffusion limits. When a single local heat source is present, the width of the heat source has a significant impact on heat conduction in the domain. Notably, there is a substantial temperature jump at the boundary when the heat source width is 10 nm.With increasing heat source width, the boundary temperature jump weakens. Furthermore, we observe that the temperature excitation rate is independent of the heat source width, while the temperature influence range expands simultaneously with the increase in heat source width. Around 500 ps, the temperature and heat flux distribution in the domain stabilize. In the case of dual heat sources, the hot zone is broader than that of a single heat source, and the temperature of the hot spot decreases as the heat source spacing increases. However, the mean heat flux remains unaffected. Upon reaching a spacing of 200 nm between the heat sources, the peak temperature in the domain remains unchanged once a steady state is reached. These findings hold significant implications for the thermal design of silicon nano-devices with local heat sources.
文摘This paper deals with the determination of temperature distribution and thermal deflection function of a thin circular plate with the stated conditions. The transient heat conduction equation is solved by using Marchi-Zgrablich transform and Laplace transform simultaneously and the results of temperature distribution and thermal deflection function are obtained in terms of infinite series of Bessel function and it is solved for special case by using Mathcad 2007 software and represented graphically by using Microsoft excel 2007.
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.
基金the National Natural Science Foundationof China(Nos.11672098,11502063)the Natural Science Foundation of Anhui Province(No.1608085QA07).
文摘The cuckoo search algorithm(CS)is improved by using the conjugate gradient method(CGM),and the CS-CGM is proposed.The unknown inner boundary shapes are generated randomly and evolved by Lévy flights and elimination mechanism in the CS and CS-CGM.The CS,CGM and CS-CGM are examined for the prediction of a pipe’s inner surface.The direct problem is two-dimensional transient heat conduction in functionally graded materials(FGMs).Firstly,the radial integration boundary element method(RIBEM)is applied to solve the direct problem.Then the three methods are compared to identify the pipe’s inner surface with the information of measured temperatures.Finally,the influences of timepoints,measurement point number and random noise on the inverse results are investigated.It is found that the three algorithms are promising and can be used to identify the pipe’s inner surface.The CS-CGM has higher accuracy and faster convergence speed than the CS and CGM.The CS and CS-CGM are insensitive to the initial values.The CGM and CS-CGM are more insensitive to the measurement noises compared with the CS.With the increase of timepoints and measurement points,and with the decrease of measurement noises,the inverse results are more accurate.
