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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its ...In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its high ability to conduct heat and electrical conductivity, in addition to being a flexible and malleable metal that is easy to form without being broken, making it one of the basic minerals that humans have benefited from for thousands of years, it is one of the first minerals. That has been discovered and extracted, and still plays a major role in the development of societies. The obtained solutions are compared with the available exact solutions and the obtained solutions using the finite difference method. The results indicate that the finite difference method is a highly effective method for obtaining approximate solutions for the thermal conductivity equation for copper. It is also clear from the numerical results from copper in the high conductivity of heat and electricity.展开更多
In the current work, transient heat conduction in a semi-infinite medium is considered for its many applications in various heat fields. Here, the homotopy analysis method （HAM） is applied to solve this problem and ...In the current work, transient heat conduction in a semi-infinite medium is considered for its many applications in various heat fields. Here, the homotopy analysis method （HAM） is applied to solve this problem and analytical results are compared with those of the exact and integral methods results. The results show that the HAM can give much better approximations than the other approximate methods： Changes in heat fluxes and profiles of temperature are obtained at different times and positions for copper, iron and aluminum.展开更多
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.展开更多
This paper is devoted to the homogenization and statistical multiscale analysis of a transient heat conduction problem in random porous materials with a nonlinear radiation boundary condition.A novel statistical multi...This paper is devoted to the homogenization and statistical multiscale analysis of a transient heat conduction problem in random porous materials with a nonlinear radiation boundary condition.A novel statistical multiscale analysis method based on the two-scale asymptotic expansion is proposed.In the statistical multiscale formulations,a unified linear homogenization procedure is established and the second-order correctors are introduced for modeling the nonlinear radiative heat transfer in random perforations,which are our main contributions.Besides,a numerical algorithm based on the statistical multiscale method is given in details.Numerical results prove the accuracy and efficiency of our method for multiscale simulation of transient nonlinear conduction and radiation heat transfer problem in random porous materials.展开更多
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.展开更多
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.展开更多
The model of transient behavior of semiconductor with heat-conduction is an initial and boundary problem. Alternating-direction multistep preconditioned iterative methods and theory analyses are given in this paper. E...The model of transient behavior of semiconductor with heat-conduction is an initial and boundary problem. Alternating-direction multistep preconditioned iterative methods and theory analyses are given in this paper. Electric potential equation is approximated by mixed finite element method, concentration and heat-conduction equations are approximated by Galerkin alternating-direction multistep methods. Error estimates of optimal order in L2 are demonstrated.展开更多
The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue...The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue-discontinue problems and extend to THM coupling.In this study,we extended the 3D-NMM to simulate both steady-state and transient heat conduction problems.The modelling was carried out using the raster methods(RSM).For the system equation,a variational method was employed to drive the discrete equations,and the crucial boundary conditions were solved using the penalty method.To solve the boundary integral problem,the face integral of scalar fields and two-dimensional simplex integration were used to accurately describe the integral on polygonal boundaries.Several numerical examples were used to verify the results of 3D steady-state and transient heat-conduction problems.The numerical results indicated that the 3D-NMM is effective for handling 3D both steadystate and transient heat conduction problems with high solution accuracy.展开更多
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 transient plane source(TPS)method is developed recently to measure the thermal conductivity of materials.In the measurement,the heating power is influenced by the heat which is transferred via the probe electrical...The transient plane source(TPS)method is developed recently to measure the thermal conductivity of materials.In the measurement,the heating power is influenced by the heat which is transferred via the probe electrical leads.This fact further influences the measurement accuracy of thermal conductivity.To solve this problem,the influence of heat loss through the electrical leads on the heating power is studied theoretically.The mathematical formula of heat loss is deduced,and the corresponding correction model is presented.A series of measurement experiments on different materials have been conducted by using the hot disk thermal constant analyzer.The results show that the influence of the heat loss on the measurement is sensitive to different test materials and probes with different sizes.When the thermal conductivity of the material is greater than 0.2 W/(m·K),the influence of the heat loss is less than 0.16%,which can be ignored.As to the lower thermal conductivity materials,it is necessary to compensate the heat loss through the electrical leads,and the accuracy of thermal conductivity measurement can be effectively improved.展开更多
基金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.
基金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.
基金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.
文摘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.
文摘In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its high ability to conduct heat and electrical conductivity, in addition to being a flexible and malleable metal that is easy to form without being broken, making it one of the basic minerals that humans have benefited from for thousands of years, it is one of the first minerals. That has been discovered and extracted, and still plays a major role in the development of societies. The obtained solutions are compared with the available exact solutions and the obtained solutions using the finite difference method. The results indicate that the finite difference method is a highly effective method for obtaining approximate solutions for the thermal conductivity equation for copper. It is also clear from the numerical results from copper in the high conductivity of heat and electricity.
文摘In the current work, transient heat conduction in a semi-infinite medium is considered for its many applications in various heat fields. Here, the homotopy analysis method （HAM） is applied to solve this problem and analytical results are compared with those of the exact and integral methods results. The results show that the HAM can give much better approximations than the other approximate methods： Changes in heat fluxes and profiles of temperature are obtained at different times and positions for copper, iron and aluminum.
基金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.
基金This work was financially supported by the National Natural Science Foundation of China(11501449)the Fundamental Research Funds for the Central Universities(3102017zy043)+1 种基金the fund of the State Key Laboratory of Solidification Processing in NWPU(SKLSP201628)the National Key Research and Development Program of China(2016YFB1100602).
文摘This paper is devoted to the homogenization and statistical multiscale analysis of a transient heat conduction problem in random porous materials with a nonlinear radiation boundary condition.A novel statistical multiscale analysis method based on the two-scale asymptotic expansion is proposed.In the statistical multiscale formulations,a unified linear homogenization procedure is established and the second-order correctors are introduced for modeling the nonlinear radiative heat transfer in random perforations,which are our main contributions.Besides,a numerical algorithm based on the statistical multiscale method is given in details.Numerical results prove the accuracy and efficiency of our method for multiscale simulation of transient nonlinear conduction and radiation heat transfer problem in random porous materials.
基金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 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.
基金This research was surpported by the National Natural Science Foundation , Mathematical TY Foun-dation (TY10126029) of China and the Youth Foundation of Shandong University.
文摘The model of transient behavior of semiconductor with heat-conduction is an initial and boundary problem. Alternating-direction multistep preconditioned iterative methods and theory analyses are given in this paper. Electric potential equation is approximated by mixed finite element method, concentration and heat-conduction equations are approximated by Galerkin alternating-direction multistep methods. Error estimates of optimal order in L2 are demonstrated.
基金supported by the National Natural Science Foundation of China(Grant Nos.42277165,41920104007,and 41731284)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant Nos.CUGCJ1821 and CUGDCJJ202234)the National Overseas Study Fund(Grant No.202106410040)。
文摘The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue-discontinue problems and extend to THM coupling.In this study,we extended the 3D-NMM to simulate both steady-state and transient heat conduction problems.The modelling was carried out using the raster methods(RSM).For the system equation,a variational method was employed to drive the discrete equations,and the crucial boundary conditions were solved using the penalty method.To solve the boundary integral problem,the face integral of scalar fields and two-dimensional simplex integration were used to accurately describe the integral on polygonal boundaries.Several numerical examples were used to verify the results of 3D steady-state and transient heat-conduction problems.The numerical results indicated that the 3D-NMM is effective for handling 3D both steadystate and transient heat conduction problems with high solution accuracy.
基金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 transient plane source(TPS)method is developed recently to measure the thermal conductivity of materials.In the measurement,the heating power is influenced by the heat which is transferred via the probe electrical leads.This fact further influences the measurement accuracy of thermal conductivity.To solve this problem,the influence of heat loss through the electrical leads on the heating power is studied theoretically.The mathematical formula of heat loss is deduced,and the corresponding correction model is presented.A series of measurement experiments on different materials have been conducted by using the hot disk thermal constant analyzer.The results show that the influence of the heat loss on the measurement is sensitive to different test materials and probes with different sizes.When the thermal conductivity of the material is greater than 0.2 W/(m·K),the influence of the heat loss is less than 0.16%,which can be ignored.As to the lower thermal conductivity materials,it is necessary to compensate the heat loss through the electrical leads,and the accuracy of thermal conductivity measurement can be effectively improved.
