L. Mai-Cao¹, T. Tran-Cong²

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 157-188, 2008, DOI:10.3970/cmes.2008.031.157

Abstract This paper presents a new meshless numerical approach to solving a special class of moving interface problems known as the passive transport where an ambient flow characterized by its velocity field causes the interfaces to move and deform without any influences back on the flow. In the present approach, the moving interface is captured by the level set method at all time as the zero contour of a smooth function known as the level set function whereas one of the two new meshless schemes, namely the SL-IRBFN based on the semi-Lagrangian method and the Taylor-IRBFN More >

Y. H. Liu^1,2, S. S. Chen¹, J. Li¹, Z. Z. Cen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 145-156, 2008, DOI:10.3970/cmes.2008.031.145

Abstract A novel meshless method for structural dynamic analysis is presented and discussed in this paper. It is called meshless local natural neighbour interpolation (MLNNI) method, which uses a meshless spatial approximation based only on nodes. The MLNNI is derived from the generalized meshless local Petrov-Galerkin (MLPG) method as a special case. Local weak forms are developed using weighted residual method locally from the dynamic partial differential equation. In the construction of trial functions, the natural neighbour interpolation (NNI) is employed to simplify the treatment of the essential boundary conditions. The domain integration is evaluated over More >

L. Codecasa¹, R. Specogna², F. Trevisan³

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 129-144, 2008, DOI:10.3970/cmes.2008.031.129

Abstract In the paper we introduce a methodology to construct discrete constitutive matrices relating magnetic fluxes with magneto motive forces (reluctance matrix) and electro motive forces with currents (conductance matrix) needed for discretizing eddy current problems over hexahedral primal grids by means of the Finite Integration Technique (FIT) and the Cell Method (CM). We prove that, unlike the mass matrices of Finite Elements, the proposed matrices ensure both the stability and the consistency of the discrete equations introduced in FIT and CM. More >

M. Steffen¹, P.C. Wallstedt², J.E. Guilkey^2,3, R.M. Kirby¹, M. Berzins¹

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.2, pp. 107-128, 2008, DOI:10.3970/cmes.2008.031.107

Abstract The Material Point Method (MPM) has shown itself to be a powerful tool in the simulation of large deformation problems, especially those involving complex geometries and contact where typical finite element type methods frequently fail. While these large complex problems lead to some impressive simulations and solutions, there has been a lack of basic analysis characterizing the errors present in the method, even on the simplest of problems. The large number of choices one has when implementing the method, such as the choice of basis functions and boundary treatments, further complicates this error analysis.\newline In More >

O. Buzzi¹, D. M. Pedroso², A. Giacomini¹

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.2, pp. 85-106, 2008, DOI:10.3970/cmes.2008.031.085

Abstract The material point method (MPM) is a numerical method for the solution of problems in continuum mechanics, including situations of large deformations. A generalization (GMPM) of this method was introduced by Bardenhagen and Kober (2004) in order to avoid some computational instabilities inherent to the original method (MPM). This generalization leads to a method more akin of the Petrov-Galerkin procedure. Although it is possible to find in the literature examples of the deduction and applications of the MPM/GMPM to specific problems, its detailed implementation is yet to be presented. Therefore, this paper attempts to describe… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.2, pp. 71-84, 2008, DOI:10.3970/cmes.2008.031.071

Abstract Iterative algorithms for solving a nonlinear system of algebraic equations of the type: F_i(x_j) = 0, i,j = 1,…,n date back to the seminal work of Issac Newton. Nowadays a Newton-like algorithm is still the most popular one due to its easy numerical implementation. However, this type of algorithm is sensitive to the initial guess of the solution and is expensive in the computations of the Jacobian matrix ∂ F_i/ ∂ x_j and its inverse at each iterative step. In a time-integration of a system of nonlinear Ordinary Differential Equations (ODEs) of the type B_ijx_j + F_i = 0… More >

T.Z. Desta¹, G. De Bruyne¹, J.-M. Aerts¹, M. Baelmans², D. Berckmans¹

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.2, pp. 61-70, 2008, DOI:10.3970/cmes.2008.031.061

Abstract This paper demonstrates the use of the concept of the local mean age of air (LMAA) to quantify ventilation effectiveness under bicycle rider's safety helmets. The specific objective is to study the effect of helmet openings on the resulting ventilation effectiveness. To quantify ventilation effectiveness using the concept of LMAA, dynamic tracer gas data are necessary. The data were generated using a Computational Fluid Dynamics (CFD) model. Two bicycle helmet designs were used and compared with respect to ventilation performance. The result showed that the helmet with more openings had better performance especially at the More >

A. Owatsiriwong¹, B. Phansri¹, K.H. Park^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.1, pp. 37-60, 2008, DOI:10.3970/cmes.2008.031.037

Abstract This study deals with the particular integral formulation for two (2D) and three (3D) dimensional elastoplastic analyses. The elastostatic equation is used for the complementary solution. The particular integrals for displacement, stress and traction rates are derived by introducing the concept of global shape function to approximate an initial stress rate term of the inhomogeneous equation. The Newton-Raphson algorithm for the plastic multiplier is used to solve the system equation. The developed program is integrated with the pre- and post-processor. The collapse analyses of the smooth flexible strip, square and circular footings are given by More >

Ömer Civalek ¹

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

Abstract In this paper, free vibration analysis of curvilinear or straight-sided quadrilateral membranes is presented. In the proposed approach, irregular physical domain is transformed into a rectangular domain by using geometric coordinate transformation. For demonstration of the accuracy and convergence of the method, some numerical examples are provided on membranes with different geometry such as skew, trapezoidal, sectorial, annular sectorial, and membranes with four curved edges. The results obtained by the DSC method are compared with those obtained by other numerical and analytical methods. More >

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

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.1, pp. 13-24, 2008, DOI:10.3970/cmes.2008.031.013

Abstract This paper presents large-scale finite element analyses of high frequency electromagnetic fields in commuter trains. The ADVENTURE_Magnetic is one of the main modules of the ADVENTURE system, which is an open source parallel finite element analyses system, and is able to solve eddy current and magnetostatic problems using the hierarchical domain decomposition method (HDDM) with an iterative linear algebraic solver. In this paper, we improve the module so as to solve a high frequency electromagnetic field of 500-1000 M[Hz]. A stationary Helmholtz equation for electromagnetic wave problems is solved taking an electric field as an More >