Jong Keun Moon¹, Seung Jo Kim²

CMES-Computer Modeling in Engineering & Sciences, Vol.67, No.3, pp. 239-264, 2010, DOI:10.3970/cmes.2010.067.239

Abstract Large structural problems with high precision and complexity require a high-performance computation using the efficient parallel algorithm. The purpose of this paper is to present the parallel performance of thermal-structural coupled analysis tested on a parallel cluster system. In the coupled analysis, the heat transfer analysis is carried out, and then the structural analysis is performed based on temperature distribution. For the automatic and efficient connection of two parallel analysis modules, the several communication patterns were studied. The parallel performance was demonstrated for the sample and the real application problems, such as a laminated composite More >

C.H. Tsai¹, D.L. Young², J. Kolibal³

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.1, pp. 53-72, 2010, DOI:10.3970/cmes.2010.066.053

Abstract The time evolution method of fundamental solutions (MFS) is proposed to solve backward heat conduction problems (BHCPs). The time evolution MFS belongs to one of the mesh-free numerical methods and is essentially composed of a sequence of diffusion fundamental solutions which exactly satisfy the heat conduction equations. Through correct treatment of temporal evolution, the resulting system of the time evolution MFS is smaller, and effectively decreases the possibility of ill-conditioning induced by such strongly ill-posed problems. Both one-dimensional and two-dimensional BHCPs are examined in this study, and the numerical results demonstrate the accuracy and stability More >

B.I. Yun¹, W.T. Ang^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.3, pp. 217-244, 2010, DOI:10.3970/cmes.2010.065.217

Abstract A dual-reciprocity boundary element method is proposed for simulating numerically axisymmetric dual-phase-lag heat conduction in nonhomogeneous thermally isotropic media. The properties of the media, such as thermal conductivity and specific heat, are assumed to vary continuously in space. To check its validity and assess its accuracy, the proposed method is first applied to solve some specific test problems with known solutions. It is then used to simulate the axisymmetric dual-phase-lag heat conduction in a particular nonhomogeneous medium subject to a concentrated surface heating. The effects of the dual phase lags and the spatial variations of More >

Han-Taw Chen¹, Li-Shie Liu¹, Shin-Ku Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.2, pp. 155-178, 2010, DOI:10.3970/cmes.2010.065.155

Abstract The finite difference method in conjunction with the least-squares scheme and experimental measured temperatures is proposed to solve a two-dimensional steady-state inverse heat conduction problem in order to predict the natural-convection heat transfer coefficient under the isothermal situation h^−iso from a three fin array mounted on a horizontal plate and fin efficiency η_f for various values of the fin spacing and fin height. The measured fin temperatures and ambient temperature are obtained from the present experimental apparatus conducted in a small wind tunnel. The heat transfer coefficient on a fin is non-uniform for the present problem,… More >

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

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 101-116, 2010, DOI:10.3970/cmes.2010.063.101

Abstract A transient 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 transient heat source acted on the FGM plates. The temperature distribution of the FGM plates is obtained in different time and different position. Some useful results for transient heat conduction are shown in figures. Applications of HNM to transient More >

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

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 1-28, 2010, DOI:10.3970/cmes.2010.063.001

Abstract A general analytical solution for a reinforced elliptical hole embedded in an infinite matrix subjected to a point heat source is provided in this paper. 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 stresses in the reinforcement layer and the matrix are derived explicitly in a series form. Some numerical results are provided to investigate the effects of the material combinations and geometric configurations on the interfacial stresses. The solution obtained can be treated as Green's functions More >

Xianqiao Wang¹, James D. Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.2, pp. 150-170, 2010, DOI:10.3970/cmes.2010.062.150

Abstract The concurrent methods for coupling molecular dynamics with continuum thermodynamics offer a myriad of challenging problems, mostly related with energy transmission, wave reflection, and damage propagation at the interfaces between the continuum description and the discrete description. In this work, by virtue of the atomistic field theory (AFT), we present an analysis to reconcile the compatibility between atomic region and continuum region and to calculate the matching temperature field of a heat conduction problem in a concurrent atomistic/continuum system. First, formulation of AFT with finite temperature and its corresponding finite element implementation are briefly introduced. More >

A. Arefmanesh¹, M. Najafi², M. Nikfar³

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.2, pp. 113-149, 2010, DOI:10.3970/cmes.2010.062.113

Abstract As some new applications of the meshless local Petrov-Galerkin method (MLPG) with unity as the test function, a number of buoyancy-driven fluid flow natural convection heat transfer problems in cavities with differentially-heated wavy side walls were analyzed. Cavities with a single wavy wall on one side as well as two wavy walls erected on both sides were considered. For the cases of the double wavy walls, two different configurations in terms of the two walls facing each other on the two sides of the cavities symmetrically or non-symmetrically were investigated. All the simulations performed in… More >

Hao-Chueh Mai¹, Kuen-Hau Lin¹, Cheng-Hsiu Yang¹, Chao-An Lin^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.1, pp. 45-62, 2010, DOI:10.3970/cmes.2010.061.045

Abstract A thermal BGK lattice Boltzmann model for flows with viscous heat dissipation is proposed. In this model, the temperature is solved by a separate thermal distribution function, where the equilibrium distribution function is similar to its hydrodynamic counterpart, except that the leading quantity is temperature. The viscous dissipation rate is obtained by computing the second-order moments of non-equilibrium distribution function, which avoids the discretization of the complex gradient term, and can be easily implemented. The proposed thermal lattice Boltzmann model is scrutinized by computing two-dimensional thermal Poiseuille flow, thermal Couette flow, natural convection in a More >

David B. Carrington¹, Xiuling Wang², Darrell W. Pepper³

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.1, pp. 23-44, 2010, DOI:10.3970/cmes.2010.061.023

Abstract A two-equation turbulence closure model (k-ω) using an h-adaptive grid technique and finite element method (FEM) has been developed to simulate low Mach flow and heat transfer. These flows are applicable to many flows in engineering and environmental sciences. Of particular interest in the engineering modeling areas are: combustion, solidification, and heat exchanger design. Flows for indoor air quality modeling and atmospheric pollution transport are typical types of environmental flows modeled with this method. The numerical method is based on a hybrid finite element model using an equal-order projection process. The model includes thermal and More >