Determination of a time-dependent heat transfer coefficient from non-standard boundary measurements

In this paper the determination of the time-dependent heat transfer coefficient in one-dimensional transient heat conduction from a non-standard boundary measurement is investigated. For this inverse nonlinear ill-posed problem the uniqueness of the solution holds. Numerical results obtained using the boundary element method are presented and discussed.     

基于结合扩展精度技术的基本解方法的非线性功能梯度材料热传导问题求解

采用基本解方法结合扩展精度技术和Kirchhoff变换求解功能梯度材料的二维热传导问题.在求解瞬态热传导问题时运用Laplace变换处理时间变量,将时域问题转化为频域问题求解;采用基本解方法计算得到高精度的频域数值解,再分别采用Stehfest和Talbot这2种数值Laplace逆变换恢复原瞬态热传导问题的计算结果.通过3个非线性功能梯度材料的稳态和瞬态热传导基准算例,分析结合扩展精度技术的基本解方法的计算精度与扩展精度位数、边界布点数和虚拟边界参数三者之间的关系.比较Stehfest和Talbot这2种数值Laplace逆变换算法的优劣.采用结合扩展精度技术的基本解方法数值研究热传导系数随位置剧烈变化的功能梯度材料热传导行为.数值结果表明该方法具有求解精度高、适用性好等特点,能高效模拟非线性功能梯度材料的二维稳态与瞬态热传导行为.     

Reconstruction of heat transfer coefficients using the boundary element method

This paper describes the reconstruction of the heat transfer coefficient (space, Problem I, or time dependent, Problem II) in one-dimensional transient inverse heat conduction problems from surface temperature or average temperature measurements. Since the inverse problem posed does not involve internal temperature measurements, this means that non-destructive testing of materials can be performed. In the formulation, convective boundary conditions relate the boundary temperature to the heat flux. Numerical results obtained using the boundary element method are presented and discussed.     

2D and 3D transient heat conduction analysis by BEM via particular integrals

Particular integral formulations are presented for 2D and 3D transient potential flow (heat conduction) analysis. The results of the analysis are compared with an alternative formulation developed using the volume integral conversion approach. Although the mathematical foundation of the two methods are different both formulations are shown to produce almost identical results.For the particular integral formulation, the steady-state heat conduction equation is used as the complementary solution and two global shape functions (GSFs) are considered to approximate the transient term of the heat conduction equation.The numerical results for three example problems are given and compared with their analytical solutions.     

Numerical solution of integrodifferential heat conduction equation for a nonlocal medium

A heat conduction model is proposed based on the relations of rational thermodynamics of irreversible processes, taking into account the nonlocal nature of the medium and the finite speed of heat propagation. In the one-dimensional case, the numerical solution of the integral-differential equation of heat conduction is obtained by the finite-element method.     

Hybrid and other computer treatments for a given heat transfer problem

A series of mechanical engineering projects has been run in an endeavour to give final year students direct programming and computational experience. This was designed to permit hybrid, other computer and analytical treatments of a one-dimensional transient heat conduction problem. The scheme was ultimately successful, although the hybrid solution proved very difficult to achieve. A comparison of the errors of the various methods is given, and evaluations are made of the technical and educational aspects of the work. Results are also given of a similar study involving finite element techniques.     

Uncertainty in Thermal Basin Modeling: An Interval Finite Element Approach

Uncertainty assessment in basin modeling and reservoir characterization is traditionally treated by geostatistical methods which are normally based on stochastic probabilistic approaches. In this paper, we present an alternative approach which is based on interval arithmetic. Here, we discuss a fnite element formulation which uses interval numbers rather than real numbers to solve the transient heat conduction in sedimentary basins. For this purpose, a novel formulation was developed to deal with both the special interval arithmetic properties and the transient term in the differential Equation governing heat transfer. In this formulation, the “stiffness” matrix resulting from the discretization of the heat conduction equation is assembled with an element-by-element technique in which the elements are globally independent and the continuity is enforced by Lagrange multipliers. This formulation is an alternative to traditional Monte Carlo method, where it is necessary to run a simulation several times to estimate the uncertainty in the results.We have applied the newly developed techniques to a one-dimensional thermal basin simulation to assess their potential and limitations.We also compared the quality of our formulation with other solution methods for interval linear systems of equations.     

A method of fundamental solutions for two-dimensional heat conduction

We investigate an application of the method of fundamental solutions (MFS) to heat conduction in two-dimensional bodies, where the thermal diffusivity is piecewise constant. We extend the MFS proposed in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction, Eng. Anal. Bound. Elem. 32 (2008), pp. 697–703] for one-dimensional heat conduction with the sources placed outside the space domain of interest, to the two-dimensional setting. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be obtained efficiently with small computational cost.     

Controller design for parabolic distributed parameter systems using finite integral transform techniques

A procedure is presented for designing a control system for distributed parameter systems of parabolic type based on the reduced-order decoupled state-space model obtained by a finite integral-transform technique. A Kalman filter is used as an observer to estimate the state variables, and state feedback control is performed. The method was applied to a one-dimensional heat conduction process and a moving bed adsorber. Both state estimation and control performances were satisfactory in spite of the model and parameter uncertainties. Following this controller design approach, the searching algorithms for the optimal sensors' and the optimal actuators' allocation problems were solved. These algorithms were applied to a one-dimensional heat conduction process in order to confirm the state estimator and controller performance. The fastest state estimation could be achieved by assigning the sensors at the optimal locations and the desired state distribution was realized with a few actuators located at the optimal positions.     

CAL in the heat-transfer laboratory

The paper describes the computer-linked laboratory experiment in one-dimensional transient heat-conduction used by the authors in a laboratory session accompanying a second-year B.Sc. (Mechanical Engineering) course in heat transfer. The physical experiment is described. and the accompanying computer program explained. An outline is presented of the basic mathematical theory underlying the computer prediction procedure, and the finite-difference analysis is described. The interactive program is used via a terminal in the laboratory first to simulate the physical experiment before it is performed, then to explore the properties of the solution procedure, next to record the data obtained by the students and finally to generate alphanumeric graphs comparing experiment and prediction. Specimens of the dialogue between the student and the computer in these various modes are given in Appendices. The results of evaluation of the innovation, both subjective and objective, are presented.     

Development of a sixth-order two-dimensional convection–diffusion scheme via Cole–Hopf transformation

In this paper, we develop a two-dimensional finite-difference scheme for solving the time-dependent convection–diffusion equation. The numerical method exploits Cole–Hopf equation to transform the nonlinear scalar transport equation into the linear heat conduction equation. Within the semi-discretization context, the time derivative term in the transformed parabolic equation is approximated by a second-order accurate time-stepping scheme, resulting in an inhomogeneous Helmholtz equation. We apply the alternating direction implicit scheme of Polezhaev to solve the Helmholtz equation. As the key to success in the present simulation, we develop a Helmholtz scheme with sixth-order spatial accuracy. As is standard practice, we validated the code against test problems which were amenable to exact solutions. Results show excellent agreement for the one-dimensional test problems and good agreement with the analytical solution for the two-dimensional problem.     

热色液晶瞬态测量全表面换热系数的技术研究

建立了基于瞬态实验技术的风洞和测试系统,用HSV色彩模型对颜色进行描述,得到颜色-温度关系;用半无限大导热理论对瞬态导热过程进行换热系数求解.测量了平板表面的换热系数分布,并与理论解进行了对比,两者相吻合.     

On the numerical implementation of continuous adjoint sensitivity for transient heat conduction problems using an isogeometric approach

A crucial problem of continuous adjoint shape sensitivity analysis is the numerical implementation of its lengthy formulations. In this paper, the numerical implementation of continuous adjoint shape sensitivity analysis is presented for transient heat conduction problems using isogeometric analysis, which can serve as a tutorial guide for beginners. Using the adjoint boundary and loading conditions derived from the design objective and the primary state variable fields, the numerical analysis procedure of the adjoint problem, which is solved backward in time, is demonstrated. Following that, the numerical integration algorithm of the shape sensitivity using a boundary approach is provided. Adjoint shape sensitivity is studied with detailed explanations for two transient heat conduction problems to illustrate the numerical implementation aspects of the continuous adjoint method. These two problems can be used as benchmark problems for future studies.     

Application of the method of fundamental solutions and radial basis functions for inverse transient heat source problem

The paper considers the determination of heat sources in unsteady 2-D heat conduction problem. The determination of the strength of a heat source is achieved by using the boundary condition, initial condition and a known value of temperature in chosen points placed inside the domain. For the solution of the inverse problem of identification of the heat source the θ-method with the method of fundamental solution and radial basis functions is proposed. Due to ill conditioning of the inverse transient heat conduction problem the Tikhonov regularization method based on SVD decomposition was used. In order to determine the optimum value of the regularization parameter the L-curve criterion was used. For testing purposes of the proposed algorithm the 2-D inverse boundary-initial-value problems in square region Ω with the known analytical solutions are considered. The numerical results show that the proposed method is easy to implement and pretty accurate. Moreover the accuracy of the results does not depend on the value of the θ parameter and is greater in the case of the identification of the temperature field than in the case of the identification of the heat sources function.     

Source term identification in 1-D IHCP

We propose stable numerical solutions for the simultaneous identification of temperature, temperature gradient, and general source terms in the one-dimensional inverse heat conduction problem (IHCP).

The numerical solution consists of a regularization procedure, based on the mollification method,and a marching scheme for the solution of the stabilized problem. The stability, error analysis and implementation of the algorithm are presented together with a set of numerical results.     

基于无网格自然邻接点Petrov-Galerkin法求解带源参数瞬态热传导问题

基于无网格自然邻接点Petrov-Galerkin法,本文建立了一种求解带源参数瞬态热传导问题的新方法.为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻接点插值构造试函数.在局部多边形子域上采用局部Petrov-Galerkin方法建立瞬态热传导问题的积分弱形式.这些多边形子域可由Delaunay三角形创建.时间域则通过传统的两点差分法进行离散.最后通过算例验证了该数值算法的有效性和正确性.     

Unconditionally stable algorithms for nonlinear heat conduction

It is shown that many commonly used one-step algorithms which are unconditionally stable for linear transient heat conduction problems become conditionally stable in the nonlinear regime. Alternative algorithms are proposed which for linear problems are identical to those commonly used, whereas for nonlinear problems the unconditional stability behavior of the linear case is retained.     

Boundary element analysis of nonlinear transient heat conduction problems

This paper is concerned with a boundary element formulation and its numerical implementation for the nonlinear transient heat conduction problems with temperature-dependent material properties. By using the Kirchhoff transformation for the material properties a set of pseudo-linear integral equations is obtained in space and time for the fully three-dimensional nonlinear problems under consideration. The resulting boundary integral equations are solved by means of the usual boundary element method. Emphasis is placed on the numerical solution procedure employing constant elements with respect to time. It is shown that in this case there is no need to evaluate the domain integrals resulting from the nonlinearity of the problem. Finally, the powerful usefulness of the proposed method is demonstrated through the numerical computation of several sample problems.     

Rapid heating vibrations of FGM annular sector plates

In this paper, thermally induced vibration of annular sector plate made of functionally graded materials is analyzed. All of the thermomechanical properties of the FGM media are considered to be temperature dependent. Based on the uncoupled linear thermoelasticity theory, the one-dimensional transient Fourier type of heat conduction equation is established. The top and bottom surfaces of the plate are under various types of rapid heating boundary conditions. Due to the temperature dependency of the material properties, heat conduction equation becomes nonlinear. Therefore, a numerical method should be adopted. First, the generalized differential quadrature method (GDQM) is implemented to discretize the heat conduction equation across the plate thickness. Next, the governing system of time-dependent ordinary differential equations is solved using the successive Crank–Nicolson time marching technique. The obtained thermal force and thermal moment resultants at each time step from temperature profile are applied to the equations of motion. The equations of motion, based on the first-order shear deformation theory (FSDT), are derived with the aid of the Hamilton principle. Using the GDQM, two-dimensional domain of the sector plate and suitable boundary conditions are divided into a number of nodal points and differential equations are turned into a system of ordinary differential equations. To obtain the unknown displacement vector at any time, a direct integration method based on the Newmark time marching scheme is utilized. Comparison investigations are performed to validate the formulation and solution method of the present research. Various examples are demonstrated to discuss the influences of effective parameters such as power law index in the FGM formulation, thickness of the plate, temperature dependency, sector opening angle, values of the radius, in-plane boundary conditions, and type of rapid heating boundary conditions on thermally induced response of the FGM plate under thermal shock.