Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attr...Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson＇s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method （RKPM）. A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.展开更多
This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional ...This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.展开更多
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.展开更多
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 article,a meshless method using the spacetime collocation for solving the two-dimensional backward heat conduction problem(BHCP)is proposed.The spacetime collocation meshless method(SCMM)is to derive the gener...In this article,a meshless method using the spacetime collocation for solving the two-dimensional backward heat conduction problem(BHCP)is proposed.The spacetime collocation meshless method(SCMM)is to derive the general solutions as the basis functions for the two-dimensional transient heat equation using the separation of variables.Numerical solutions of the heat conduction problem are expressed as a series using the addition theorem.Because the basis functions are the general solutions of the governing equation,the boundary points may be collocated on the spacetime boundary of the domain.The proposed method is verified by conducting several heat conduction problems.We also carry out numerical applications to compare the SCMM with other meshless methods.The results show that the SCMM is accurate and efficient.Furthermore,it is found that the recovered boundary data on inaccessible boundary can be obtained with high accuracy even though the over specified data are provided only at a 1/6 portion of the spacetime boundary.展开更多
By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variable...By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variables. Firstly, the spatial space and temporal domain are discretized by FEM and precise integral algorithm respectively. Then, the high accuracy semi-analytical solution of direct problem can be got. Finally, based on the solution, the computing model of inverse problem and expression of sensitivity analysis are established. Single variable and variables combined identifications including thermal parameters, boundary conditions and source-related terms etc. are given to validate the approach proposed in 1-D and 2-D cases. The effects of noise data and initial guess on the results are investigated. The numerical examples show the effectiveness of this approach.展开更多
This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several mo...This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.展开更多
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.展开更多
In this paper, we present an effective meshless method for solving the inverse heat conduction problems, with the Neumann boundary condition. A PDE-constrained optimization method is developed to get a global approxim...In this paper, we present an effective meshless method for solving the inverse heat conduction problems, with the Neumann boundary condition. A PDE-constrained optimization method is developed to get a global approximation scheme in both spatial and temporal domains, by using the fundamental solution of the governing equation as the basis function.Since the initial measured data contain some noises, and the resulting systems of equations are usually ill-conditioned, the Tikhonov regularization technique with the generalized crossvalidation criterion is applied to obtain more stable numerical solutions. It is shown that the proposed schemes are effective by some numerical tests.展开更多
Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green func...Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.展开更多
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.展开更多
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.展开更多
Non-standard backward heat conduction problem is ill-posed in the sense that the solution(if it exists) does not depend continuously on the data. In this paper, we propose a regularization strategy-quasi-reversibility...Non-standard backward heat conduction problem is ill-posed in the sense that the solution(if it exists) does not depend continuously on the data. In this paper, we propose a regularization strategy-quasi-reversibility method to analysis the stability of the problem. Meanwhile, we investigate the roles of regularization parameter in this method. Numerical result show that our algorithm is effective and stable.展开更多
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.展开更多
In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat condu...In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.展开更多
Under consideration is a nonclassical stationary problem on heat conduction in a body with the pre-set surface temperature and heat flow. The body contains inclusions at unknown locations and with unknown boundaries. ...Under consideration is a nonclassical stationary problem on heat conduction in a body with the pre-set surface temperature and heat flow. The body contains inclusions at unknown locations and with unknown boundaries. The body and inclusions have different constant thermal conductivities. The author explores the possibility of locating inclusions. The article presents an integral criterion based on which a few statements on identification of inclusions in a body are proved.展开更多
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.展开更多
A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh rat...A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.展开更多
The classical heat conduction equation is derived from the assumption that the temperature increases immediately after heat transfer, but the increase of temperature is a slow process, so the memory-dependent heat con...The classical heat conduction equation is derived from the assumption that the temperature increases immediately after heat transfer, but the increase of temperature is a slow process, so the memory-dependent heat conduction model has been reconstructed. Numerical results show that the solution of the initial boundary value problem of the new model is similar to that of the classical heat conduction equation, but its propagation speed is slower than that of the latter. In addition, the propagation speed of the former is also affected by time delay and kernel function.展开更多
The results of studies by solving the inverse thermal conductivity problem of the heat capacity of evaporator of the short linear heat pipes (HP’s) with a Laval nozzle-liked vapour channel and intended for cooling sp...The results of studies by solving the inverse thermal conductivity problem of the heat capacity of evaporator of the short linear heat pipes (HP’s) with a Laval nozzle-liked vapour channel and intended for cooling spacecraft and satellites with strict take-off mass regulation are presented. Mathematical formulation of the inverse problem for the HP’s thermal conductivity in one-dimensional coordinate system is accompanied by the measurement results using the monotonic heating method in a vacuum adiabatic calorimeter the HP’s surface temperatures along the longitudinal axis over the entire temperature load range, thermal resistance, and arrays of thermal power data on the evaporator Q<sub>ev</sub> and vortex flow calorimeter Q<sub>cond</sub> for the condensation surface allow us to estimate the average value of the evaporator heat capacity C<sub>ev</sub> by solving the inverse thermal conductivity problem in the HP’s evaporator region. Since at the beginning of working fluid boiling for a certain time interval, the temperature of the capillary-porous evaporator remains close to constant, and with the continuation of heating and by solving the inverse thermal conductivity problem, it becomes possible to calculate the heat capacity of the working evaporator and the evaporation specific heat of the boiling working fluid and compare it with the table values.展开更多
基金supported by the Natural Science Foundation of Ningbo,China (Grant Nos.2009A610014 and 2009A610154)the Natural Science Foundation of Zhejiang Province,China (Grant No.Y6090131)
文摘Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson＇s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method （RKPM）. A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.
基金This study was supported by“the Fundamental Research Funds for the Central Universities”(Grant No.2015B37814)the Postgraduate Research and Practice Innovation Program of Jiangsu Province(Grant No.KYLX15_0489)+1 种基金the National Natural Science Foundation of China(Grant No.51679081)“the Fundamental Research Funds for the Central Universities”(Grant No.2018B48514).
文摘This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.
基金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.
基金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.
文摘In this article,a meshless method using the spacetime collocation for solving the two-dimensional backward heat conduction problem(BHCP)is proposed.The spacetime collocation meshless method(SCMM)is to derive the general solutions as the basis functions for the two-dimensional transient heat equation using the separation of variables.Numerical solutions of the heat conduction problem are expressed as a series using the addition theorem.Because the basis functions are the general solutions of the governing equation,the boundary points may be collocated on the spacetime boundary of the domain.The proposed method is verified by conducting several heat conduction problems.We also carry out numerical applications to compare the SCMM with other meshless methods.The results show that the SCMM is accurate and efficient.Furthermore,it is found that the recovered boundary data on inaccessible boundary can be obtained with high accuracy even though the over specified data are provided only at a 1/6 portion of the spacetime boundary.
文摘By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variables. Firstly, the spatial space and temporal domain are discretized by FEM and precise integral algorithm respectively. Then, the high accuracy semi-analytical solution of direct problem can be got. Finally, based on the solution, the computing model of inverse problem and expression of sensitivity analysis are established. Single variable and variables combined identifications including thermal parameters, boundary conditions and source-related terms etc. are given to validate the approach proposed in 1-D and 2-D cases. The effects of noise data and initial guess on the results are investigated. The numerical examples show the effectiveness of this approach.
基金Project supported by the Key Disciplines of Shanghai Municipality (Grant No.S30104)the Shanghai Leading Academic Discipline Project (Grant No.J50101)
文摘This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.
基金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.
基金Supported by the National Natural Science Foundation of China(Grant Nos.1129014311471066+3 种基金11572081)the Fundamental Research of Civil Aircraft(Grant No.MJ-F-2012-04)the Fundamental Research Funds for the Central Universities(Grant No.DUT15LK44)the Scientific Research Funds of Inner Mongolia University for the Nationalities(Grant No.NMD1304)
文摘In this paper, we present an effective meshless method for solving the inverse heat conduction problems, with the Neumann boundary condition. A PDE-constrained optimization method is developed to get a global approximation scheme in both spatial and temporal domains, by using the fundamental solution of the governing equation as the basis function.Since the initial measured data contain some noises, and the resulting systems of equations are usually ill-conditioned, the Tikhonov regularization technique with the generalized crossvalidation criterion is applied to obtain more stable numerical solutions. It is shown that the proposed schemes are effective by some numerical tests.
基金Project supported by the National Natural Science Foundation of China (No. 50776097)
文摘Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.
基金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 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.
文摘Non-standard backward heat conduction problem is ill-posed in the sense that the solution(if it exists) does not depend continuously on the data. In this paper, we propose a regularization strategy-quasi-reversibility method to analysis the stability of the problem. Meanwhile, we investigate the roles of regularization parameter in this method. Numerical result show that our algorithm is effective and stable.
基金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.
文摘In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.
文摘Under consideration is a nonclassical stationary problem on heat conduction in a body with the pre-set surface temperature and heat flow. The body contains inclusions at unknown locations and with unknown boundaries. The body and inclusions have different constant thermal conductivities. The author explores the possibility of locating inclusions. The article presents an integral criterion based on which a few statements on identification of inclusions in a body are proved.
基金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.
文摘A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.
文摘The classical heat conduction equation is derived from the assumption that the temperature increases immediately after heat transfer, but the increase of temperature is a slow process, so the memory-dependent heat conduction model has been reconstructed. Numerical results show that the solution of the initial boundary value problem of the new model is similar to that of the classical heat conduction equation, but its propagation speed is slower than that of the latter. In addition, the propagation speed of the former is also affected by time delay and kernel function.
文摘The results of studies by solving the inverse thermal conductivity problem of the heat capacity of evaporator of the short linear heat pipes (HP’s) with a Laval nozzle-liked vapour channel and intended for cooling spacecraft and satellites with strict take-off mass regulation are presented. Mathematical formulation of the inverse problem for the HP’s thermal conductivity in one-dimensional coordinate system is accompanied by the measurement results using the monotonic heating method in a vacuum adiabatic calorimeter the HP’s surface temperatures along the longitudinal axis over the entire temperature load range, thermal resistance, and arrays of thermal power data on the evaporator Q<sub>ev</sub> and vortex flow calorimeter Q<sub>cond</sub> for the condensation surface allow us to estimate the average value of the evaporator heat capacity C<sub>ev</sub> by solving the inverse thermal conductivity problem in the HP’s evaporator region. Since at the beginning of working fluid boiling for a certain time interval, the temperature of the capillary-porous evaporator remains close to constant, and with the continuation of heating and by solving the inverse thermal conductivity problem, it becomes possible to calculate the heat capacity of the working evaporator and the evaporation specific heat of the boiling working fluid and compare it with the table values.
