C. Moulinec¹, R. Issa², J.-C. Marongiu³, D. Violeau⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.3, pp. 133-148, 2008, DOI:10.3970/cmes.2008.025.133

Abstract The gridless Smoothed Particle Hydrodynamics (SPH) numerical method is preferably used in Computational Fluid Dynamics (CFD) to simulate complex flows with one or several convoluted free surfaces. This type of flows requires distorted meshes with classical Eulerian mesh-based methods or very fine meshes with Volume of Fluid method. Few 3-D SPH simulations have been carried out to our knowlegde so far, mainly due to prohibitive computational investment since the number of particles required in 3-D is usually too large to be handled by a single processor. In this paper, a parallel 3-D SPH code is More >

Atsuya Oishi¹, Shinobu Yoshimura²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 115-132, 2008, DOI:10.3970/cmes.2008.025.115

Abstract This paper describes the finite element analyses of dynamic problems using graphics hardware. The graphics hardware, known as GPU that is an acronym of Graphics Processing Unit, was first developed only for processing 3D computer graphics. However it has obtained both flexible programmability using a high-level shader programming language such as OpenGL Shading Language (GLSL), and has also obtained fast numerical processing ability of over 160 GFLOPS that is much faster than CPU. In this paper, GPU is utilized for the finite element analyses of dynamic problems. Two different computational tasks in the dynamic finite More >

Y.Y. Shan¹, C. Shu^1,2, Z.L. Lu³

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 99-114, 2008, DOI:10.3970/cmes.2008.025.099

Abstract The local multiquadric-based differential quadrature (MQ-DQ) method proposed by [Shu, Ding, and Yeo (2003)] is a natural mesh-free approach for derivative approximation, which is easy to be implemented to solve problems with curved boundary. Previously, it has been well tested for the two-dimensional (2D) case. In this work, this mesh-free method was extended to simulate fluid flow problems with curved boundary in three-dimensional (3D) space. The main concern of this work is to numerically study the performance of the 3D local MQ-DQ method and demonstrate its capability and flexibility for simulation of 3D incompressible fluid More >

J.R. Jimenez-Octavio¹, O. Lopez-Garcia², E. Pilo¹, A. Carnicero²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 81-98, 2008, DOI:10.3970/cmes.2008.025.081

Abstract This paper presents a multidisciplinary design and optimization method of power transmission lines. This optimization method solves both mechanical and electrical problems by a new strongly coupled method that also optimizes the potential designs using a genetic algorithm. A multi-objective function is formulated to simplify a constrained typical optimization problem into an unconstrained one. The scope of this work is the sizing and configuration optimization problem with fixed topology. The method is applied to a railway overhead transmission line. The genetic algorithm is applied to mechanical, electrical and electromechanical optimization problems obtaining good results. Finally, More >

S. D. Harris¹, R. Mustata², L. Elliott², D. B. Ingham², D. Lesnic²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 69-80, 2008, DOI:10.3970/cmes.2008.025.069

Abstract Two homogeneous anisotropic materials are butted together to form a contact surface within a single composite material (the specimen). An inverse boundary element method (BEM) is developed to determine the components of the hydraulic conductivity tensor of each material and the position of the contact surface. A steady state flow is forced through the specimen by the application of a constant pressure differential on its opposite faces. Experimental measurements (simulated) of pressure and average hydraulic flux at exposed boundaries are then used in a modified least squares functional. This functional minimises the gap between the More >

P. Le¹, N. Mai-Duy², T. Tran-Cong³, G. Baker⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 43-68, 2008, DOI:10.3970/cmes.2008.025.043

Abstract This paper describes an integrated radial basis function network (IRBFN) method for the numerical modelling of the dynamics of strain localization due to strain softening in quasi-brittle materials. The IRBFN method is a truly meshless method that is based on an unstructured point collocation procedure. We introduce a new and effective regularization method to enhance the performance of the IRBFN method and alleviate the numerical oscillations associated with weak discontinuity at the elastic wave front. The dynamic response of a one dimensional bar is investigated using both local and non-local continuum models. Numerical results, which More >

Arezoo Emdadi¹, Edward J. Kansa², Nicolas Ali Libre^1,3, Mohammad Rahimian¹, Mohammad Shekarchi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 23-42, 2008, DOI:10.3970/cmes.2008.025.023

Abstract We present a new method based upon the paper of Volokh and Vilney (2000) that produces highly accurate and stable solutions to very ill-conditioned multiquadric (MQ) radial basis function (RBF) asymmetric collocation methods for partial differential equations (PDEs). We demonstrate that the modified Volokh-Vilney algorithm that we name the improved truncated singular value decomposition (IT-SVD) produces highly accurate and stable numerical solutions for large values of a constant MQ shape parameter, c, that exceeds the critical value of c based upon Gaussian elimination. More >

D. Soares Jr.¹, M. P. Vinagre²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 1-8, 2008, DOI:10.3970/cmes.2008.025.001

Abstract This work presents an alternative procedure to compute time-domain electromagnetic fields. The Boundary Element Method is here adopted to numerically analyze wave propagation problems, computing just a so-called primary field (either the electric or the magnetic field can be selected as primary field; the complementary field is here named secondary field). The secondary field is obtained following Maxwell's equations, i.e., considering space derivatives of the primary field (computed by the Complex Variable Method) and time integration procedures. This methodology is more efficient and flexible since fewer systems of equations must be solved at each time-step. More >

Arun K. Nair¹, Diana Farkas², Ronald D. Kriz¹

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 239-248, 2008, DOI:10.3970/cmes.2008.024.239

Abstract The indentation response of Ni thin films of thicknesses in the nano scale was studied using molecular dynamics simulations with embedded atom method (EAM) interatomic potentials. Simulations were performed in single crystal films in the [111] orientation with thicknesses of 7nm and 33nm. In the elastic regime, the loading curves observed start deviating from the Hertzian predictions for indentation depths greater than 2.5% of the film thickness. The observed loading curves are therefore dependent on the film thickness. The simulation results also show that the contact stress necessary to emit the first dislocation under the More >

Mustafa Uludoğan¹, D. Paula Guarin¹, Zully E. Gomez¹, Tahir Cagin¹, William A. Goddard III²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 215-238, 2008, DOI:10.3970/cmes.2008.024.215

Abstract Aiming at a presentation of the utility of the state of art of first-principles methods (PBE flavor of Density Functional Theory, DFT) in the area of materials science and engineering, we present our studies of the equation of state and ferroelectric-paraelectricphase transition in several ferroelectric systems, including BaTiO₃, PbTiO₃, SrTiO₃, AgNbO₃, Pb_xBa_1-xTiO₃ and Sr_xBa_1-xTiO₃. We also report the Born effective charges, optical dielectric constant, and phonon dispersion relation properties from Density Functional Perturbation Theory. Computed results are compared with other theoretical results (which were mostly on BaTiO₃, PbTiO₃, cubic SrTiO₃ using various approaches, as well as experiments. The studies on More >