Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.1, pp. 1-22, 2009, DOI:10.3970/cmes.2009.047.001

Abstract In this paper an analytical method for approximating the solution of backward heat conduction problem is presented. The Fourier series expansion technique is used to formulate a first-kind Fredholm integral equation for the temperature field u(x,t) at any time t < T, when the data are specified at a final time T. Then we consider a direct regularization, instead of the Tikhonov regularization, by adding the term αu(x,t) to obtain a second-kind Fredholm integral equation. The termwise separable property of kernel function allows us by transforming it to a two-point boundary value problem, and thus a closed-form solution is derived.… More >

Donghuan Liu¹, Xiaoping Zheng^1,2, Yinghua Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.3, pp. 263-300, 2009, DOI:10.3970/cmes.2009.039.263

Abstract A discontinuous Galerkin (DG) finite element method for the heat conduction problems with local high gradient and thermal contact resistance is presented. The DG formulation is constructed by employing the stabilization term and the Bassi-Rebay numerical flux term. The stabilization term is defined by a penalization of the temperature jump at the interface. By eliminating the penalization term of the temperature jump in the region of local high gradient and imperfect contact interfaces, the present DG method is applied to solve problems involving local high gradient and thermal contact resistance where the numerical flux is obtained from the definition of… More >

J. Sladek¹, V. Sladek¹, C.L. Tan², S.N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 161-174, 2008, DOI:10.3970/cmes.2008.032.161

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for the solution of steady-state and transient heat conduction problems in a continuously non-homogeneous anisotropic medium. The Laplace transform is used to treat the time dependence of the variables for transient problems. The analyzed domain is covered by small subdomains with a simple geometry. A weak formulation for the set of governing equations is transformed into local integral equations on local subdomains by using a unit test function. Nodal points are randomly distributed in the 3D analyzed domain and each node is surrounded by a spherical subdomain to which… More >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 99-122, 2008, DOI:10.3970/cmes.2008.030.099

Abstract The application of the method of fundamental solutions (MFS) to inverse boundary value problems associated with the steady-state heat conduction in isotropic media in the presence of sources, i.e. the Poisson equation, is investigated in this paper. Based on the approach of Alves and Chen (2005), these problems are solved in two steps, namely by finding first an approximate particular solution of the Poisson equation and then the numerical solution of the resulting inverse boundary value problem for the Laplace equation. The resulting MFS discretised system of equations is ill-conditioned and hence it is solved by employing the singular value… More >

WU XueHong¹, SHEN ShengPing², TAO WenQuan^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.1, pp. 65-76, 2007, DOI:10.3970/cmes.2007.022.065

Abstract Meshless local Petrov-Galerkin collocation method is applied to compute two-dimensional heat conduction problems in irregular domain. By taking the Dirac's Delta function as the test function, the local domain integration is avoided. The essential boundary conditions can be implemented easily in this method. A case that has analytical solution shows the present method can obtain desired accuracy and efficient. Two cases in engineering are computed to validate the approach by comparing the present method with the finite volume method (FVM) solutions obtained from a commercial CFD package FLUENT 6.3. The results show that the present method is in good agreement… More >

N. Osawa¹, K. Hashimoto¹, J. Sawamura¹, J. Kikuchi², Y. Deguchi², T. Yamaura²

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.1, pp. 43-54, 2007, DOI:10.3970/cmes.2007.020.043

Abstract A new hypothesis regarding heat transmission during line heating is proposed. It states that the distribution of the temperature of the gas adjacent to the plate, T_G, and the overall local heat transfer coefficient, α, depend only on the distance from the torch. An identification technique for T_G and α is developed. The validity of the employed hypothesis and the proposed technique is demonstrated by comparing the measured and identified T_G during a spot heating test. The plate temperature calculated by direct heat conduction analysis closely approximates the one measured for the spot and line heating tests, when T_G and… More >

S. Bargmann, P. Steinmann¹

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 133-150, 2005, DOI:10.3970/cmes.2005.009.133

Abstract The present contribution is concerned with the modeling and computation of non-classical heat conduction. In the 90s Green and Naghdi presented a new theory which is fully consistent. We suggest a solution method based on finite elements for the spatial as well as for the temporal discretization. A numerical example is compared to existing experimental results in order to illustrate the performance of the method. More >

C.-W. Chang¹, C.-S. Liu², J.-R. Chang^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 13-38, 2005, DOI:10.3970/cmes.2005.010.013

Abstract In this paper, the inverse heat conduction problem governed by sideways heat equation is investigated numerically. The problem is ill-posed because the solution, if it exists, does not depend continuously on the data. To begin with, this ill-posed problem is analyzed by considering the stability of the semi-discretization numerical schemes. Then the resulting ordinary differential equations at the discretized times are numerically integrated towards the spatial direction by the group preserving scheme, and the stable range of the index r = 1/2ν Δt is investigated. When the numerical results are compared with exact solutions, it is found that they are… More >

J. Sladek¹, V. Sladek¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 309-318, 2004, DOI:10.3970/cmes.2004.006.309

Abstract Meshless methods based on the local Petrov-Galerkin approach are proposed for solution of steady and transient heat conduction problem in a continuously nonhomogeneous anisotropic medium. Fundamental solution of the governing partial differential equations and the Heaviside step function are used as the test functions in the local weak form. It is leading to derive local boundary integral equations which are given in the Laplace transform domain. The analyzed domain is covered by small subdomains with a simple geometry. To eliminate the number of unknowns on artificial boundaries of subdomains the modified fundamental solution and/or the parametrix with a convenient cut-off… More >

Y. Ochiai

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.4, pp. 435-446, 2001, DOI:10.3970/cmes.2001.002.435

Abstract The boundary element method (BEM) is useful in solving the steady heat conduction problem of orthotropic bodies without heat generation. However, for cases with arbitrary heat generation, a number of internal cells are necessary. In this paper, it is shown that the problem of steady heat conduction in orthotropic bodies with heat generation can be solved without internal cells by the triple-reciprocity BEM. In this method, the distribution of heat generation is interpolated using integral equations. In order to solve the problem, the values of heat generation at internal points and on the boundary are used. Furthermore, a new interpolation… More >