T. Furumura¹, L. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.2, pp. 153-168, 2004, DOI:10.3970/cmes.2004.006.153

Abstract Recent developments of the Earth Simulator, a high-performance parallel computer, has made it possible to realize realistic 3D simulations of seismic wave propagations on a regional scale including higher frequencies. Paralleling this development, the deployment of dense networks of strong ground motion instruments in Japan (K-NET and KiK-net) has now made it possible to directly visualize regional seismic wave propagation during large earthquakes. Our group has developed an efficient parallel finite difference method (FDM) code for modeling the seismic wavefield and a 3D visualization technique, both suitable for implementation on the Earth Simulator. Large-scale 3D… More >

Y. Bei¹, B.J. Fregly¹, W.G. Sawyer¹, S.A. Banks^1,2, N.H. Kim¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.2, pp. 145-152, 2004, DOI:10.3970/cmes.2004.006.145

Abstract Mild wear of ultra-high molecular weight polyethylene tibial inserts continues to affect the longevity of total knee replacements (TKRs). Using static finite element and elasticity analyses, previous studies have hypothesized that polyethylene wear can be reduced by using a thicker tibial insert to decrease contact pressures. To date, no study has taken this hypothesis to the next step by performing dynamic analyses under in vivo functional conditions to quantify the relationship between contact pressures, insert thickness, and mild wear. This study utilizes multibody dynamic simulations incorporating elastic contact to perform such analyses. \textit {In vivo} fluoroscopic… More >

Z.Y. Qian¹, Z.D. Han¹, P. Ufimtsev¹, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.2, pp. 133-144, 2004, DOI:10.3970/cmes.2004.006.133

Abstract The weak-form of Helmholtz differential equation, in conjunction with vector test-functions (which are gradients of the fundamental solutions to the Helmholtz differential equation in free space) is utilized as the basis in order to directly derive non-hyper-singular boundary integral equations for the velocity potential, as well as its gradients. Thereby, the presently proposed boundary integral equations, for the gradients of the acoustic velocity potential, involve only O(r⁻²) singularities at the surface of a 3-D body. Several basic identities governing the fundamental solution to the Helmholtz differential equation for velocity potential, are also derived for the further More >

S.Aoki¹, K.Amaya², M.Urago³, A.Nakayama⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.2, pp. 123-132, 2004, DOI:10.3970/cmes.2004.006.123

Abstract The Fast Multipole Boundary Element Method(FMBEM) which is suitable for a large scale computation is applied to corrosion analysis. Many techniques of the FMBEM on the potential problems can be usefully employed. Additionally, some procedures are developed for corrosion analysis. To cope with the non-linearity due to the polarization curve, the Bi-CGSTAB iterative method which is commonly used in the FMBEM is modified. To solve infinite domain problems, the M₀⁰ which is obtained naturally in the multipole expansion is conveniently used. A pipe element for the FMBEM is developed. A couple of example problems are More >

Z.Y. Qian¹, Z.D. Han¹, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.1, pp. 115-122, 2004, DOI:10.3970/cmes.2004.006.115

Abstract This article has no abstract. More >

Haiyi Liang¹, C.H. Woo^1,2, Hanchen Huang³, A.H.W. Ngan⁴, T.X. Yu⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.1, pp. 105-114, 2004, DOI:10.3970/cmes.2004.006.105

Abstract Molecular dynamics (MD) simulations are performed to study crystalline plasticity during nano-indentation by comparing the elastic-plastic response of three copper substrates with surfaces (001), (110), and (111) crystallographic planes. The effects of elastic anisotropy and crystallographic symmetry on the reduced modulus, dislocation nucleation, and subsequent microstructure evolution, are investigated. The reduced modulus of (111) surface is found to be the largest, while that of (001) surface is the smallest. Elastic stress distribution calculated from finite element method (FEM) is qualitatively consistent with the MD simulation results. Significant differences exist in the deformation behavior in the More >

Shengping Shen¹, S. N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.1, pp. 91-104, 2004, DOI:10.3970/cmes.2004.006.091

Abstract An atomistic level stress tensor is defined with physical clarity, based on the SPH method. This stress tensor rigorously satisfies the conservation of linear momentum, and is appropriate for both homogeneous and inhomogeneous deformations. The formulation is easier to implement than other stress tensors that have been widely used in atomistic analysis, and is validated by numerical examples. The present formulation is very robust and accurate, and will play an important role in the multiscale simulation, and in molecular dynamics. An equivalent continuum is also defined for the molecular dynamics system, based on the developed More >

E. Pan¹, X. Jiang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.1, pp. 77-90, 2004, DOI:10.3970/cmes.2004.006.077

Abstract It is of great importance to understand the factors that contribute to the strain and electrical distributions, which are induced by the misfit strain between a buried quantum wire (QWR) and its surrounding matrix. One of the important factors is the shape or geometry of cross section of the QWR. Utilizing a recent exact closed-form solution [Pan (2004)], we study the model system of QWRs with different shapes and calculate both the surface and internal elastic and piezoelectric fields induced by QWRs embedded in semiconductor GaAs substrates by properly setting the size and location of More >

N. Mai-Duy¹, T. Tran-Cong²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.1, pp. 59-76, 2004, DOI:10.3970/cmes.2004.006.059

Abstract This paper presents a new domain decomposition technique based on the use of Boundary Integral Equations (BIEs) for the analysis of viscous flow problems. The domain of interest is divided into a number of non-overlapping subdomains and an iterative procedure is then employed to update the boundary conditions at interfaces. The new feature in the present work is that at each iteration, the relevant two subdomains, together containing a particular interface, are assumed to satisfy the governing BI equations which they do at the end of a convergent iterative process. Hence the boundary conditions on More >

António Tadeu¹, Julieta António and Nuno Simões

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.1, pp. 43-58, 2004, DOI:10.3970/cmes.2004.006.043

Abstract This Analytical Green's functions for the steady-state response of homogeneous three-dimensional unbounded, half-space, slab and layered solid media subjected to a spatially sinusoidal harmonic heat line source are presented. In the literature, this problem is frequently referred to as the two-and-a-half dimensional fundamental solution or 2.5D Green's functions.
The proposed equations are theoretically interesting in themselves and they are also useful as benchmark results for validating numerical applications. They are also of great practical use in the formulation of three dimensional heat transfer problems in layered solid formations using integral transform methods and/or boundary elements.… More >