Michael J. Rodgers¹, Leon M. Keer, Herbert S. Cheng

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.4, pp. 483-496, 2002, DOI:10.3970/cmes.2002.003.483

Abstract The steady-state temperature rise due to frictional heating on the surface of coated halfspaces and halfplanes is described by closed form expressions in the Fourier transformed frequency domain. These frequency response functions (FRFs) include the effects of the coating and the speed of the moving heat source and apply for all Peclet number regimes. Analytical inversion of these expressions for several special cases shows the Green's functions as infinite series of images, which may be costly and slowly convergent. Also, the influence coefficients integrated from these Green's functions are not available in closed form. Applying fast Fourier transform (FFT) methods… More >

Shuangbiao Liu, Michael J. Rodgers, Qian Wang, Leon M. Keer, Herbert S. Cheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.4, pp. 465-482, 2002, DOI:10.3970/cmes.2002.003.465

Abstract Computing the temperature rise and thermoelastic displacement of a material subjected to frictional heating is essential for the realistic modeling of the performance of mechanical components. This paper presents a novel set of frequency-domain expressions for the surface temperature rise and the surface normal thermoelastic displacement of a moving three-dimensional elastic halfspace subjected to arbitrary transient frictional heating, where the velocity of the body and its direction can be an arbitrary function of time. Frequency response functions are derived by using the Carslaw-Jaeger theory, the Seo-Mura result, and the Fourier transform. General formulas are expressed in the form of time… More >

Yuanzhong Hu¹, Yuchuan Liu, Hui Wang

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.4, pp. 447-454, 2002, DOI:10.3970/cmes.2002.003.447

Abstract Scuffing, a major cause of failure in automobile engines, is considered as a dynamic process in this study. Local adhesions may occur randomly in lubricated contacts due to the existence of asperity contact and breakdown of lubricating films. Scuffing would take place if the local events develop rapidly into a large-scale plastic deformation and catastrophic failure. A system dynamic model established in the present paper allows one to predict dynamic behavior of a tribological system through numerical solutions of a group of differential equations. Results show that a transition to adhesion begins when the surface temperature goes beyond a critical… More >

J.H. Wu¹, L.Y. Zhu², J.T. Tang³

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.3, pp. 225-240, 2012, DOI:10.3970/cmes.2012.086.225

Abstract Magnetic induction hyperthermia is one of hopeful methods for tumor therapy. In this method, several ferromagnetic seeds are needed to be implanted into the tumor. The seeds would produce energy, and cause the nearby tumor to die. Temperature prediction is significant before treatment. In addition, in clinical treatment, the tumor temperature has to be monitored in realtime. However, using as few thermometers as possible is the basic principle. Fortunately, the numerical simulation can contribute to realtime measurement. The seed temperature is modeled based on the Haider's method, which is treated as the thermal boundary in numerical simulation. We employ the… More >

Wei-Chung Tien¹, Yue-Tzu Yang¹, Cha’o-Kuang Chen^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.2, pp. 155-170, 2012, DOI:10.3970/cmes.2012.084.155

Abstract In this paper, the homotopy analysis method is applied to analyze the heat transfer of the oscillating base temperature processes occurring in a convective rectangular fin with variable thermal conductivity. This method is a powerful and easy-to-use tool for non-linear problems and it provides us with a simple way to adjust and control the convergence region of solution series. Without the need of iteration, the obtained solution is in the form of an infinite power series and the results indicated that the series has high accuracy by comparing it with those generated by the complex combination method. More >

Cao Guohua, Li Kai, Zhu Zhencai, Peng Weihong, Mao Xianbiao

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.4, pp. 403-424, 2012, DOI:10.3970/cmes.2012.083.403

Abstract In order to master the geometric and mechanical behavior of helical structure under complicated condition such as the hoisting rope in mine shaft and the transmitting cable in electric power, the thermal expansion characteristic of single helical structure is systematically investigated under temperature effect in different layer. Linearly explicit expressions of axial strain and increment of helical angle for the helical unit of the ith layer are developed. Based on theory of curve by Love and theory of wire rope by Costello, the linearly explicit expressions of tension, torsion and bending moment of the helical unit are presented. After that,… More >

T. Nakamura¹, R. Kobyashi¹, S. Ogata¹

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.4, pp. 357-386, 2011, DOI:10.3970/cmes.2011.073.357

Abstract The coarse-grained particle (CGP) method has been proposed to coarse-grain a crystalline system of atoms to meso-scale. In the method, virtual particles are distributed in the system, and the inter-particle interaction is calculated through the constrained statistical ensemble average of the atomic Hamiltonian at a given temperature. For simplicity, however, the harmonic approximation has been used for the inter-atomic interaction and hence anharmonicity at finite temperatures has been ignored. We improve the former CGP method to incorporate the anharmonicity of atomic system at finite temperatures into the inter-particle interaction. Also the divide-and-conquer strategy is applied to calculate the inter-particle interaction… 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. Then we develop a new… More >

S. Brischetto¹, E. Carrera²

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.3, pp. 249-302, 2010, DOI:10.3970/cmes.2010.056.249

Abstract This work considers the fully coupled thermo-mechanical analysis of one-layered and multilayered isotropic and composite shells. The temperature is assumed a primary variable as the displacement; it is therefore directly obtained from the model and this feature permits the temperature field to be evaluated through the thickness direction. Three problems are analyzed: - static analysis of shells with imposed temperature on the external surfaces; - static analysis of shells subjected to a mechanical load, with the possibility of considering the temperature field effects; - a free vibration problem, with the evaluation of the temperature field effects. In the first problem,… More >

Tai-Ming Chang¹, Chien-Chou Weng¹, Mei-Jiau Huang^1,2, Chun-KaiLiu², Chih-Kuang Yu²

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

Abstract We employ the non-equilibrium molecular dynamics (NEMD) simulation to calculate the in-plane thermal conductivity of silicon thin films of thickness 2.2nm and 11nm. To eliminate the finite-size effect, samples of various lengths are simulated and an extrapolation technique is applied. To perform the quantum correction which is necessary as the MD simulation temperature is lower than Debye temperature, the confined phonon spectra are obtained in advance via the EMD simulations. The investigation shows the thermal conductivities corrected based on the bulk and thin-film phonon densities of states are very close and they agree excellently with the theoretical predictions of a… More >