P. C. Wallstedt¹, J. E. Guilkey¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 223-232, 2007, DOI:10.3970/cmes.2007.019.223

Abstract The standard velocity projection scheme for the Material Point Method (MPM) and a typical form of the GIMP Method are examined. It is demonstrated that the fidelity of information transfer from a particle representation to the computational grid is strongly dependent on particle density and location. In addition, use of non-uniform grids and even non-uniform particle sizes are shown to introduce error. An enhancement to the projection operation is developed which makes use of already available velocity gradient information. This enhancement facilitates exact projection of linear functions and reduces the dependence of projection accuracy on particle location and density for… More >

D.L. Young¹, K.H. Chen², J.T. Chen³, J.H. Kao⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 197-222, 2007, DOI:10.3970/cmes.2007.019.197

Abstract A boundary-type method for solving the Laplace problems using the modified method of fundamental solutions (MMFS) is proposed. The present method (MMFS) implements the singular fundamental solutions to evaluate the solutions, and it can locate the source points on the real boundary as contrasted to the conventional MFS, where a fictitious boundary is needed to avoid the singularity of diagonal term of influence matrices. The diagonal term of influence matrices for arbitrary domain can be novelly determined by relating the MFS with the indirect BEM and are also solved for circular domain analytically by using separable kernels and circulants. The… More >

S. Chen^1,2, Y. Liu^1,3, B.C. Khoo⁴, X.J. Fan⁵, J.T. Fan⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 181-196, 2007, DOI:10.3970/cmes.2007.019.181

Abstract The mesoscopic simulation for fluids flow in a square microchannel is investigated using dissipative particle dynamics. The velocity distribution for single fluid in a square channel is compared with the solutions of CFD solver, which is found to be in good agreement with each other. The no-slip boundary condition could be well held for the repulsive coefficient ranged from 9.68 to 18.0. For the same range of repulsive coefficient, various wettabilities could be obtained by changing the repulsive coefficient for binary immiscible fluids, in which the immiscible fluids are achieved by increasing the repulsive force between species. The typical motion… More >

A. Arockiarajan¹, A. Menzel²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 163-178, 2007, DOI:10.3970/cmes.2007.019.163

Abstract To study rate-dependent properties of piezoelectric materials a micro-mechanically motivated model is applied in this work. The developed framework is embedded into a coupled three-dimensional finite element setting, whereby each element is assumed to represent one grain and, moreover, possesses a random initialisation of the underlying polarisation direction. Furthermore, an energy-based criterion is used for the initiation of the onset of domain switching and the subsequent propagation of domain wall motion during the switching process is modelled via a linear kinetics theory. The interaction between individual grains is thereby incorporated by means of a probabilistic approach -- a purely phenomenologically… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 145-162, 2007, DOI:10.3970/cmes.2007.019.145

Abstract A new method is developed to solve the Dirichlet problems for the two-dimensional Laplace equation in the doubly-connected domains, namely the meshless regularized integral equations method (MRIEM), which consists of three portions: Fourier series expansion, the Fredholm integral equations, and linear equations to determine the unknown boundary conditions onartificial circles. The boundary integral equations on artificial circles are singular-free and the kernels are degenerate. When boundary-type methods are inefficient to treat the problems with complicated domains, the new method can be applicable for such problems. The new method by using the Fourier series and the Fourier coefficients can be adopted… More >

Y.F. Nie¹, S. Chang¹, X.K. Fan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 135-144, 2007, DOI:10.3970/cmes.2007.019.135

Abstract A new parallel mechanism for Node-based Seamless Finite Element Method was proposed in this paper, which possessed the following three prominent points: realizing the workload balance for the parallel processes naturally, achieving synchronization of all the schedules under complex parallel environment, and filling up the gap between pre-processing and main processing. To support the scheme, three specific solutions of the parallel mechanism were proposed in this paper, all of which achieved the highly efficient parallel seamless connection between the FEM mesh generation process and structure analysis process. Two of the three schemes, i.e. dynamic scheme and self-adaptive scheme, enabled both… More >

T.C. Theodosiou¹, D.A. Saravanos²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 121-134, 2007, DOI:10.3970/cmes.2007.019.121

Abstract In this paper a new method is introduced for carbon nanotube modeling combining features of Molecular Mechanics and Finite Element Analysis. Repetitive atomic cells are treated as finite elements, whose internal energy is determined by the semi-empirical Brenner molecular potential model; internal forces and linearized stiffness matrices are formulated analytically in order to gain in speed and accuracy, and the resultant discrete system is formulated and solved using the Newton-Raphson method. The presented method is validated through comparisons to numerical and experimental results provided by other researchers. The bending and shearing of CNTs is also simulated. More >

Shunde Yin¹, Leo Rothenburg¹, Maurice B. Dusseault¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 111-120, 2007, DOI:10.3970/cmes.2007.019.111

Abstract Subsidence problem is of great importance in petroleum engineering and environmental engineering. In this paper, we firstly apply a hybrid Displacement Discontinuity-FEM modeling to this classic problem: the evaluation of subsidence over a compacting oil reservoir. We use displacement discontinuity method to account for the reservoir surrounding area, and finite element methods in the fully coupled simulation of the reservoir itself. This approach greatly reduces the number of degrees of freedom compared to an analyzing fully coupled problem using only a finite element or finite difference discretization. More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 99-110, 2007, DOI:10.3970/cmes.2007.019.099

Abstract A new meshless regularized integral equation method (MRIEM) is developed to solve the interior and exterior Dirichlet problems for the two-dimensional Laplace equation, which consists of three parts: Fourier series expansion, the second kind Fredholm integral equation and an analytically regularized solution of the unknown boundary condition on an artificial circle. We find that the new method is powerful even for the problem with complex boundary shape and with random noise disturbing the boundary data. More >

Amit Shaw¹, D Roy²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 69-98, 2007, DOI:10.3970/cmes.2007.019.069

Abstract A novel discretization strategy and derivative reproduction based on reproducing kernel (RK) particle approximations of functions are proposed. The proposed scheme is in the form of an RK interpolation that offers significant numerical advantages over a recent version of the strategy by Chen et al. (2003), wherein the authors added a set of primitive functions to the reproducing kernel (enrichment) functions. It was also required that the support size of the primitive function be less than the smallest distance between two successive grid points. Since the primitive function was required to vary from 0 to 1 within half of this… More >