Jiongyang Wu¹, Inanc Senocak¹, Guoyu Wang², Yulin Wu³, Wei Shyy¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 679-690, 2003, DOI:10.3970/cmes.2003.004.679

Abstract Cavitation appears in a wide variety of fluid machinery, and can often cause negative impacts on performance and structural integrity. A main computational difficulty for cavitation is the large density ratio between liquid and vapor phases, around 1000 for water under normal temperature and pressure conditions. Moreover, cavitating flows are usually turbulent and the interfacial dynamics is complex. The fast time scales associated with turbulent cavitation also poses substantial challenges computationally and experimentally. In the present study, pressure-based algorithms are adopted to simulate three-dimensional turbulent cavitating flows in a hollow-jet valve. The Favre-averaged Navier-Stokes equations… More >

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

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 665-678, 2003, DOI:10.3970/cmes.2003.004.665

Abstract The numerical implementation of the truly Meshless Local Petrov-Galerkin (MLPG) type weak-forms of the displacement and traction boundary integral equations is presented, for solids undergoing small deformations. In the accompanying part I of this paper, the general MLPG/BIE weak-forms were presented [Atluri, Han and Shen (2003)]. The MLPG weak forms provide the most general basis for the numerical solution of the non-hyper-singular displacement and traction BIEs [given in Han, and Atluri (2003)], which are simply derived by using the gradients of the displacements of the fundamental solutions [Okada, Rajiyah, and Atluri (1989a,b)]. By employing the… More >

J. A. Nairn¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 649-664, 2003, DOI:10.3970/cmes.2003.004.649

Abstract A new algorithm is described which extends the material point method (MPM) to allow explicit cracks within the model material. Conventional MPM enforces velocity and displacement continuity through its background grid. This approach is incompatible with cracks which are displacement and velocity discontinuities. By allowing multiple velocity fields at special nodes near cracks, the new method (called CRAMP) can model cracks. The results provide an ``exact'' MPM analysis for cracks. Comparison to finite element analysis and to experiments show it gets good results for crack problems. The intersection of crack surfaces is prevented by implementing More >

Jan Sladek¹, Vladimir Sladek¹, Chuanzeng Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 637-648, 2003, DOI:10.3970/cmes.2003.004.637

Abstract A new computational method for solving transient elastodynamic initial-boundary value problems in continuously non-homogeneous solids, based on the meshless local Petrov-Galerkin (MLPG) method, is proposed in the present paper. The moving least squares (MLS) is used for interpolation and the modified fundamental solution as the test function. The local Petrov-Galerkin method for unsymmetric weak form in such a way is transformed to the local boundary integral equations (LBIE). The analyzed domain is divided into small subdomains, in which a weak solution is assumed to exist. Nodal points are randomly spread in the analyzed domain and More >

E. J. Sellountos¹, D. Polyzos²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 619-636, 2003, DOI:10.3970/cmes.2003.004.619

Abstract A new meshless local Petrov-Galerkin (MLPG) method for solving two dimensional frequency domain elastodynamic problems is proposed. Since the method utilizes, in its weak formulation, either the elastostatic or the frequency domain elastodynamic fundamental solution as test function, it is equivalent to the local boundary integral equation (LBIE) method. Nodal points spread over the analyzed domain are considered and the moving least squares (MLS) interpolation scheme for the approximation of the interior and boundary variables is employed. Two integral equations suitable for the integral representation of the displacement fields in the local sub- domains are… More >

Ch. Zhang², A. Savaidis³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 603-618, 2003, DOI:10.3970/cmes.2003.004.603

Abstract A novel non-hypersingular time-domain traction BEM is presented for three-dimensional (3-D) transient elastodynamic crack analysis. The initial-boundary value problem is formulated as a set of non-hypersingular time-domain traction boundary integral equations (BIEs). To solve the time-domain traction BIEs, a time-stepping scheme based on the convolution quadrature formula of Lubich (1988a,b; 1994) for temporal discretization and a collocation method for spatial discretization is adopted. Numerical examples are given for an unbounded solid with a penny-shaped crack under a tensile and shear impact loading. A comparison of the present time-domain BEM with the conventional one shows that More >

R. Mathur¹, S. G. Advani², B. K. Fink³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 587-602, 2003, DOI:10.3970/cmes.2003.004.587

Abstract Real number-coded hybrid genetic algorithms for optimal design of resin injection locations for the resin transfer molding process are evaluated in this paper. Resin transfer molding (RTM) is widely used to manufacture composite parts with material and geometric complexities, especially in automotive and aerospace sectors. The sub-optimal location of the resin injection locations (gates) can leads to the formation of resin starved regions and require long mold fill times, thus affecting the part quality and increasing manufacturing costs. There is a need for automated design algorithms and software that can determine the best gate and… More >

Q. Li¹, S. Shen¹, Z. D. Han¹, S. N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 571-586, 2003, DOI:10.3970/cmes.2003.004.571

Abstract In this paper, a truly meshless method, the Meshless Local Petrov-Galerkin (MLPG) Method, is developed for three-dimensional elasto-statics. The two simplest members of MLPG family of methods, the MLPG type 5 and MLPG type 2, are combined, in order to reduce the computational requirements and to obtain high efficiency. The MLPG5 method is applied at the nodes inside the 3-D domain, so that any domain integration is eliminated altogether, if no body forces are involved. The MLPG 2 method is applied at the nodes that are on the boundaries, and on the interfaces of material More >

M. Kröger¹, M. Müller², J. Nievergelt²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 559-570, 2003, DOI:10.3970/cmes.2003.004.559

Abstract We present a novel method for generating starting polymer structures for molecular simulations in the dense phase. The work describes the ingredients of an algorithm for the creation of large, dense or diluted amorphous polymeric systems close to equilibrium and provides measures for its quality. The model systems are made of semiflexible (wormlike) repulsive multibead chains. The key feature of the method is its efficiency, in particular for large systems, while approaching given local and global chain characteristics. Its output has been proven to serve as an excellent basis for subsequent off-lattice molecular dynamics computer More >

T. Iwasaki¹ and H. Miura¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 551-558, 2003, DOI:10.3970/cmes.2003.004.551

Abstract We have developed a molecular-dynamics technique for investigating migration-induced failures in interconnect films for ULSIs. This technique was used to simulate grain-boundary grooving in Al and Cu films. The simulations showed that the grain-boundary grooves are formed by atomic diffusion at the grain boundary. To clarify what kind of underlay material is effective in suppressing this diffusion, we calculated the dependence of groove depth on the kind of underlay material. The calculation showed that the groove depth of the Al film decreases in the order: Al/Ta, Al/W, and Al/TiN while that of the Cu film More >