Ching Kong Chao^1,2, Chin Kun Chen¹, Fu Mo Chen³

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.3, pp. 283-298, 2009, DOI:10.3970/cmes.2009.047.283

Abstract Analytical exact solutions of a fundamental heat conduction problem for a three-phase elliptical composite under a remote uniform heat flow are provided in this paper. The steady-state temperature and heat flux fields in each phase of an elliptical composite are analyzed in detail. Investigations on the present heat conduction problem are tedious due to the presence of material inhomogeneities and geometric discontinuities. Based on the technique of conformal mapping and the method of analytical continuation in conjunction with the alternating technique, the general expressions of the temperature and heat flux are derived explicitly in a More >

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 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… More >

W.F.Yuan¹, K.H.Tan ¹

CMC-Computers, Materials & Continua, Vol.13, No.1, pp. 49-62, 2009, DOI:10.3970/cmc.2009.013.049

Abstract A simple cellular automaton (CA) scheme is proposed to simulate heat conduction in anisotropic domains. The CA is built on random nodes rather than an irregular grid. The local rule used in the CA is defined by physical concepts instead of differential equations. The accuracy of the proposed approach is verified by classical examples. As an application of the proposed method, the CA approach is incorporated into fibre model which is widely used in finite element analysis to calculate the temperature distribution on the cross-section of composite beams. Numerical examples demonstrate that the proposed scheme More >

Chih-Ping Wu^1,2, Shao-En Huang²

CMC-Computers, Materials & Continua, Vol.12, No.3, pp. 251-282, 2009, DOI:10.3970/cmc.2009.012.251

Abstract A modified Pagano method is developed for the three-dimensional (3D) coupled analysis of simply-supported, doubly curved functionally graded (FG) piezo-thermo-elastic shells under thermal loads. Four different loading conditions, applied on the lateral surfaces of the shells, are considered. The material properties of FG shells are regarded as heterogeneous through the thickness coordinate, and then specified to obey an exponent-law dependent on this. The Pagano method, conventionally used for the analysis of multilayered composite elastic plates/shells, is modified to be feasible for the present analysis of FG piezo-thermo-elastic plates/shells. The modifications include that a displacement-based formulation More >

LiviuMarin ¹

CMC-Computers, Materials & Continua, Vol.12, No.1, pp. 71-100, 2009, DOI:10.3970/cmc.2009.012.071

Abstract In this paper, the alternating iterative algorithm originally proposed by Kozlov, Maz'ya and Fomin (1991) is numerically implemented for the Cauchy problem in anisotropic heat conduction using a meshless method. Every iteration of the numerical procedure consists of two mixed, well-posed and direct problems which are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point More >

J.H. Tian^1,2, X. Han², S.Y. Long³, G.Q. Xie⁴

CMC-Computers, Materials & Continua, Vol.10, No.3, pp. 229-242, 2009, DOI:10.3970/cmc.2009.010.229

Abstract A heat conduction analysis of the functionally graded material (FGM) plates has been investigated based on the hybrid numerical method (HNM). HNM combines the layer element method with the method of Fourier transforms and proves to be efficient and reliable. The FGM plates are infinite large and the material properties vary continuously through thickness. The heat source continually acted one the FGM plates. The temperature distribution of the FGM plates is obtained in different time and different position. Some useful results for heat conduction problems are shown in figures. This article applies HNM to heat 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 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 More >

Chih-Wen Chang¹, Chein-Shan Liu², Jiang-Ren Chang¹

The International Conference on Computational & Experimental Engineering and Sciences, Vol.3, No.2, pp. 69-80, 2007, DOI:10.3970/icces.2007.003.069

Abstract By using a quasi-boundary regularization we can formulate a two-point boundary value problem of the backward heat conduction equation. The ill-posed problem is analyzed by using the semi-discretization numerical schemes. Then, the resulting ordinary differential equations in the discretized space are numerically integrated towards the time direction by the Lie-group shooting method to find the unknown initial conditions. The key point is based on the erection of a one-step Lie group element G(T) and the formation of a generalized mid-point Lie group element G(r). Then, by imposing G(T) = G(r) we can seek the missing More >