Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

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.

Efficient Inverse Analysis for Solving a Coupled Conduction,Convection and Radiation Problem Involving Non⁃gray Participating Media

The presence of non-gray radiative properties in a reheating furnace’s medium that absorbs,emits,and involves non-gray creates more complex radiative heat transfer problems.Furthermore,it adds difficulty to solving the coupled conduction,convection,and radiation problem,leading to suboptimal efficiency that fails to meet real-time control demands.To overcome this difficulty,comparable gray radiative properties of non-gray media are proposed and estimated by solving an inverse problem.However,the required iteration numbers by using a least-squares method are too many and resulted in a very low inverse efficiency.It is necessary to present an efficient method for the equivalence.The Levenberg-Marquardt algorithm is utilized to solve the inverse problem of coupled heat transfer,and the gray-equivalent radiative characteristics are successfully recovered.It is our intention that the issue of low inverse efficiency,which has been observed when the least-squares method is employed,will be resolved.To enhance the performance of the Levenberg-Marquardt algorithm,a modification is implemented for determining the damping factor.Detailed investigations are also conducted to evaluate its accuracy,stability of convergence,efficiency,and robustness of the algorithm.Subsequently,a comparison is made between the results achieved using each method.

NUMERICAL METHOD OF MIXED FINITE VOLUME-MODIFIED UPWIND FRACTIONAL STEP DIFFERENCE FOR THREE-DIMENSIONAL SEMICONDUCTOR DEVICE TRANSIENT BEHAVIOR PROBLEMS

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.

A new complex variable meshless method for transient heat conduction problems

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.

Meshless analysis of three-dimensional steady-state heat conduction problems

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.

A mixed Newton-Tikhonov method for nonlinear ill-posed problems

Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical problems are complicated. We propose a mixed Newton-Tikhonov method, i.e., one step Newton-Tikhonov method with several other steps of simplified Newton-Tikhonov method. Convergence and stability of this method are proved under some conditions. Numerical experiments show that the proposed method has obvious advantages over the classical Newton method in terms of computational costs.

Exact analytical solution to three-dimensional phase change heat transfer problems in biological tissues subject to freezing

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.

Newton-type methods and their modifications for inverse heat conduction problems

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.

A complex variable meshless local Petrov-Galerkin method for transient heat conduction problems

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.

Three-dimensional numerical manifold method for heat conduction problems with a simplex integral on the boundary

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.

A Numerical Method on Inverse Determination of Heat Transfer Coefficient Based on Thermographic Temperature Measurement

The heat transfer coefficient in a multidimensional heat conduction problem is obtained from the solution of the inverse heat conduction problem based on the thermographic temperature measurement. The modified one-dimensional correction method （MODCM）, along with the finite volume method, is employed for both two- and three-dimensional inverse problems. A series of numerical experiments are conducted in order to verify the effectiveness of the method. In addition, the effect of the temperature measurement error, the ending criterion of the iteration, etc. on the result of the inverse problem is investigated. It is proved that the method is a simple, stable and accurate one that can solve successfully the inverse heat conduction problem.

Spatial interfacial heat transfer and surface characteristics during gravity casting of A356 alloy

As one of the key boundary conditions during casting solidification process, the interfacial heat transfer coefficient (IHTC) affects the temperature variation and distribution. Based on the improved nonlinear estimation method (NEM), thermal measurements near both bottom and lateral metal-mold interfaces throughout A356 gravity casting process were carried out and applied to solving the inverse heat conduction problem (IHCP). Finite element method (FEM) is employed for modeling transient thermal fields implementing a developed NEM interface program to quantify transient IHTCs. It is found that IHTCs at the lateral interface become stable after the volumetric shrinkage of casting while those of the bottom interface reach the steady period once a surface layer has solidified. The stable value of bottom IHTCs is 750 W/(m^2·℃), which is approximately 3 times that at the lateral interface. Further analysis of the interplay between spatial IHTCs and observed surface morphology reveals that spatial heat transfer across casting-mold interfaces is the direct result of different interface evolution during solidification process.

An efficient inverse approach for reconstructing time-and space-dependent heat flux of participating medium

The decentralized fuzzy inference method(DFIM)is employed as an optimization technique to reconstruct time-and space-dependent heat flux of two-dimensional(2D)participating medium.The forward coupled radiative and conductive heat transfer problem is solved by a combination of finite volume method and discrete ordinate method.The reconstruction task is formulated as an inverse problem,and the DFIM is used to reconstruct the unknown heat flux.No prior information on the heat flux distribution is required for the inverse analysis.All retrieval results illustrate that the time-and spacedependent heat flux of participating medium can be exactly recovered by the DFIM.The present method is proved to be more efficient and accurate than other optimization techniques.The effects of heat flux form,initial guess,medium property,and measurement error on reconstruction results are investigated.Simulated results indicate that the DFIM is robust to reconstruct different kinds of heat fluxes even with noisy data.

Improved Particle Swarm Optimization for Solving Transient Nonlinear Inverse Heat Conduction Problem in Complex Structure

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.

A meshless model for transient heat conduction analyses of 3D axisymmetric functionally graded solids

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.

OPTIMAL ERROR BOUND IN A SOBOLEV SPACE OF REGULARIZED APPROXIMATION SOLUTIONS FOR A SIDEWAYS PARABOLIC EQUATION

The inverse heat conduction problem （IHCP） is a severely ill-posed problem in the sense that the solution （ if it exists） does not depend continuously on the data. But now the results on inverse heat conduction problem are mainly devoted to the standard inverse heat conduction problem. Some optimal error bounds in a Sobolev space of regularized approximation solutions for a sideways parabolic equation, i. e. , a non-standard inverse heat conduction problem with convection term which appears in some applied subject are given.

Resolving double-sided inverse heat conduction problem using calibration integral equation method

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.

Reconstruction of Temperature Field in 3-D, Absorbing, Emitting, and Anisotropically Scattering Medium

The soft measurement technology of flame temperature field is an efficient method to learn the combustion status in furnace. Generally, it reconstructs the temperature field in furnace through the image of flame, which is a process to solve radiative inverse problem. In this paper, the flame of pulverized coal is considered as 3-D, absorbing, emitting, and anisotropically scattering non-gray medium. Through the study on inverse problem of radiative heat transfer, the temperature field in this kind of medium has been reconstructed. The mechanism of 3-D radiative heat transfer in a rectangular media, which is 2 m×3 m× 5 m and full of CO2, N2 and carbon particles, is studied with Monte Carlo method. The 3-D temperature field in this rectangular space is reconstructed and the influence of particles density profile is discussed.

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper,we present a unified finite volume method preserving discrete maximum principle(DMP)for the conjugate heat transfer problems with general interface conditions.We prove the existence of the numerical solution and the DMP-preserving property.Numerical experiments show that the nonlinear iteration numbers of the scheme in[24]increase rapidly when the interfacial coefficients decrease to zero.In contrast,the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero,which reveals that the unified scheme is more robust than the scheme in[24].The accuracy and DMP-preserving property of the scheme are also veri ed in the numerical experiments.

Multi-Modes Multiscale Approach of Heat Transfer Problems in Heterogeneous Solids with Uncertain Thermal Conductivity

Stochastic temperature distribution should be carefully inspected in the thermal-failure design of heterogeneous solids with unexpected random energy excitations.Stochastic multiscale modeling for these problems involve multiscale and highdimensional uncertain thermal parameters,which remains limitation of prohibitive computation.In this paper,we propose a multi-modes based constrained energy minimization generalized multiscale finite element method(MCEM-GMsFEM),which can transform the original stochastic multiscale model into a series of recursive multiscale models sharing the same deterministic material parameters by multiscale analysis.Thus,MCEM-GMsFEM reveals an inherent low-dimensional representation in random space,and is designed to effectively reduce the complexity of repeated computation of discretized multiscale systems.In addition,the convergence analysis is established,and the optimal error estimates are derived.Finally,several typical random fluctuations on multiscale thermal conductivity are considered to validate the theoretical results in the numerical examples.The numerical results indicate that the multi-modes multiscale approach is a robust integrated method with the excellent performance.