Peter Lucas¹, Alexander H. van Zuijlen¹, Hester Bijl¹

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.2, pp. 147-176, 2009, DOI:10.3970/cmes.2009.041.147

Abstract Mesh adaptation is a fairly established tool to obtain numerically accurate solutions for flow problems. Computational efficiency is, however, not always guaranteed for the adaptation strategies found in literature. Typically excessive mesh growth diminishes the potential efficiency gain. This paper, therefore, extends the strategy proposed by [Aftosmis and Berger (2002)] to compute the refinement threshold. The extended strategy computes the refinement threshold based on a user desired number of grid cells and adaptations, thereby ensuring high computational efficiency. Because our main interest is flow around wind turbines, the adaptation strategy has been optimized for flow More >

Hung-Chang Lee¹, Chein-Shan Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.1, pp. 1-26, 2009, DOI:10.3970/cmes.2009.041.001

Abstract The group-preserving schemes developed by Liu (2001) for integrating ordinary differential equations system were adopted the Cayley transform and Padé approximants to formulate the Lie group from its Lie algebra. However, the accuracy of those schemes is not better than second-order. In order to increase the accuracy by employing the group-preserving schemes on ordinary differential equations, according to an efficient technique developed by Runge and Kutta to raise the order of accuracy from the Euler method, we combine the Runge-Kutta method on the group-preserving schemes to obtain the higher-order numerical methods of group-preserving type. They More >

S.H. Pyo¹, H.K. Lee^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.3, pp. 271-306, 2009, DOI:10.3970/cmes.2009.040.271

Abstract This paper presents a micromechanical approach for predicting the elastic and multi-level damage response of aligned and randomly oriented whisker-/ short fiber-reinforced composites. Based on a combination of Eshelby's micromechanics and the evolutionary imperfect interface approach, the effective elastic moduli of the composites are derived explicitly. The modified Eshelby's tensor for spheroidal inclusions with slightly weakened interface [Qu (1993b)] is extended in the present study to model whiskers or short fibers having mild or severe imperfect interfaces. Aligned and random orientations of spheroidal reinforcements are considered. A multi-level damage model in accordance with the Weibull's More >

D. L. Young^1,3, K. H. Chen², T. Y. Liu³, L. H. Shen³, C. S. Wu³

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.3, pp. 225-270, 2009, DOI:10.3970/cmes.2009.040.225

Abstract In this article, a hypersingular meshless method (HMM) is extended to solve 3D potential problems for arbitrary domains after a 2D model was successfully developed (Young et al. 2005a). The solutions are represented by a distribution of the double layer potentials instead of the single layer potentials as generally used in the conventional method of fundamental solutions (MFS). By using the desingularization technique to regularize the singularity and hypersingularity of the double layer potentials, the source points can be located exactly on the real boundary to avoid the sensitivity of locating fictitious boundary for putting… More >

A. Takei¹, S. Yoshimura¹, H. Kanayama²

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.1, pp. 63-82, 2009, DOI:10.3970/cmes.2009.040.063

Abstract This paper describes a large-scale finite element analysis (FEA) for a high-frequency electromagnetic field of Maxwell equations including the displacement current. A stationary Helmholtz equation for the high-frequency electromagnetic field analysis is solved by considering an electric field and an electric scalar potential as unknown functions. To speed up the analysis, the hierarchical domain decomposition method (HDDM) is employed as a parallel solver. In this study, the Parent-Only type (Parallel processor mode: P-mode) of the HDDM is employed. In the P-mode, Parent processors perform the entire FEA. In this mode, all CPUs can be used More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.1, pp. 1-28, 2009, DOI:10.3970/cmes.2009.040.001

Abstract We consider an inverse problem of Laplace equation by recoverning boundary value on the inner circle of a two-dimensional annulus from the overdetermined data on the outer circle. The numerical results can be used to determine the Robin coefficient or crack's position inside a disk from the measurements of Cauchy data on the outer boundary. The Fourier series is used to formulate the first kind Fredholm integral equation for the unknown data f(θ) on the inner circle. Then we consider a Lavrentiev regularization, by adding an extra term αf(θ) to obtain the second kind Fredholm integral More >

W. Attaporn¹, H. Koguchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.3, pp. 237-262, 2009, DOI:10.3970/cmes.2009.039.237

Abstract In the present study, a stress singularity field along free edges meeting at a corner in a three-dimensional joint structure is investigated. The order of stress singularity is determined using an eigen analysis based on a finite element method. Intensities of stress singularity not only at the corner but also along the free edge of interface are determined directly without any post-processing by a new FEM formulation referred to as a three-dimensional enriched FEM. Result in the present analysis is also compared with that in another numerical method. It was slightly larger than the intensity More >

D. S. Liu¹, C.Y. Tsai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.1, pp. 29-48, 2009, DOI:10.3970/cmes.2009.039.029

Abstract Utilizing a thin copper substrate for illustration purposes, this study presents a novel numerical method for extracting the thermo-mechanical properties of a thin-film. In the proposed approach, molecular dynamics (MD) simulations are performed to establish the load-displacement response of a thin copper substrate nanoindented at temperatures ranging from 300~1400 K. The load data are then input to an artificial neural network (ANN), trained using a finite element model (FEM), in order to extract the material constants of the copper substrate. The material constants are then used to construct the corresponding stress-strain curve, from which the… More >

A. B. Mahmud Hasan^1,2, S.M. Guo¹, M.A. Wahab¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.4, pp. 399-410, 2009, DOI:10.3970/fdmp.2009.005.399

Abstract The effects of gas flow channel design were studied experimentally for increasing fuel cell performance and fuel optimization. Three types of gas flow channels (serpentine, straight and interdigitated) were designed on the basis of water flooding due to electrochemical reactions, electro-osmotic drag, etc. Experimental results indicate that the best cell performance can be obtained by arranging interdigitated gas flow channel at the anode side and serpentine gas flow channel at the cathode side. Detailed analysis on complex two phase water generation and electrochemical phenomena behind those results were analyzed in this work to find out More >

Shishir Sinha¹, Vinay Kumar²

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.3, pp. 297-312, 2009, DOI:10.3970/fdmp.2009.005.297

Abstract High-capacity ion exchange resins based on Poly Methyl Methacrylate (PMMA) was developed during this study having higher mechanical strength, stable and, hydrophilic in nature compared to Poly styrene (PS) based systems. PMMA was cross-linked by Di Vinyl Benzene (DVB) to get macro porous PMMA-DVB resins by suspension polymerization. PMMA-DVB resins were surface modified by nitration followed by amination under appropriate temperature conditions. Change in the color of resin from white to pale yellow and pale yellow to golden yellow as well as characteristic peaks of NO2at 700 cm-1and 1400 cm-1and 3100-3200 cm-1for NH2group in FTIR More >