Yiorgos-Sokratis Smyrlis¹, Andreas Karageorghis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 535-550, 2003, DOI:10.3970/cmes.2003.004.535

Abstract In this study, we investigate the application of the Method of Fundamental Solutions for the solution of biharmonic Dirichlet problems on a disk. Modifications of the method for overcoming sources of inaccuracy are suggested. We also propose an efficient algorithm for the solution of the resulting systems which exploits the symmetries of the matrices involved. The techniques described in the paper are applied to standard test problems. More >

I. S. Raju¹, D. R. Phillips²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 141-160, 2003, DOI:10.3970/cmes.2003.004.141

Abstract An accurate and yet simple Meshless Local Petrov-Galerkin (MLPG) formulation for analyzing beam problems is presented. In the formulation, simple weight functions are chosen as test functions as in the conventional MLPG method. Linear test functions are also chosen, leading to a variation of the MLPG method that is computationally efficient compared to the conventional implementation. The MLPG method is evaluated by applying the formulation to a variety of patch tests, thin beam problems, and problems with load discontinuities. The formulation successfully reproduces exact solutions to machine accuracy when higher order power and spline functions More >

Jin Yeon Cho¹, Seung Jo Kim²

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.6, pp. 687-698, 2002, DOI:10.3970/cmes.2002.003.687

Abstract In this work, an explicit solution procedure for the recently developed discontinuous time integration method is proposed in order to reduce the computational cost while maintaining the desirable numerical characteristics of the discontinuous time integration method. In the present explicit solution procedure, a two-stage correction algorithm is devised to obtain the solution at the next time step without any matrix factorization. To observe the numerical characteristics of the proposed explicit solution procedure, stability and convergence analyses are performed. From the stability analysis, it is observed that the proposed algorithm gives a larger critical time step More >

Olivier Ricou¹, Michel Bercovier²

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.1, pp. 87-102, 2002, DOI:10.3970/cmes.2002.003.087

Abstract In this article, we present a method of reduction of the dimension of the Stokes equations by one in a quasi-cylindrical domain. It takes the special shape of the domain into account by the use of a projection onto a space of polynomials defined over the thickness. The polynomials are defined to fit as well as possible with the variables they approximate. Hence, this method restricted to the first polynomial, recovers the Hele-Shaw approximation.
The convergence of the approximate solution to the continuous one is shown. Under a regularity hypothesis, we also obtain error estimates. More >

Weimin Han, Xueping Meng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.1, pp. 65-76, 2002, DOI:10.3970/cmes.2002.003.065

Abstract Interests in meshfree (or meshless) methods have grown rapidly in the recent years in solving boundary value problems arising in mechanics, especially in dealing with difficult problems involving large deformation, moving discontinuities, etc. Rigorous error estimates of a meshfree method, the reproducing kernel particle method, for smooth solutions have been theoretically derived and experimentally tested in Han, Meng (2001). In this paper, we provide an error analysis of the meshfree method for solving problems with singular solutions. The results are presented in the context of one-dimensional problems. The error estimates are of optimal order and More >

Shuyao Long¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.1, pp. 53-64, 2002, DOI:10.3970/cmes.2002.003.053

Abstract Meshless methods have been extensively popularized in literature in recent years, due to their flexibility in solving boundary value problems. The meshless local Petrov-Galerkin(MLPG) method for solving the bending problem of the thin plate is presented and discussed in the present paper. The method uses the moving least-squares approximation to interpolate the solution variables, and employs a local symmetric weak form. The present method is a truly meshless one as it does not need a mesh, either for the purpose of interpolation of the solution or for the integration of the energy. All integrals can More >

António Tadeu, Julieta António¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.4, pp. 477-496, 2001, DOI:10.3970/cmes.2001.002.477

Abstract This paper presents analytical solutions, together with explicit expressions, for the steady state response of homogeneous three-dimensional layered acoustic and elastic formations subjected to a spatially sinusoidal harmonic line load. These formulas are theoretically interesting in themselves and they are also useful as benchmark solutions for numerical applications. In particular, they are very important in formulating three-dimensional elastodynamic problems in layered fluid and solid formations using integral transform methods and/or boundary elements, avoiding the discretization of the solid-fluid interfaces. The proposed Green's functions will allow the solution to be obtained for high frequencies, for which More >

F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 365-372, 2001, DOI:10.3970/cmes.2001.002.365

Abstract The aim of this paper is to present a methodology for solving the plane elasticity problem using the Cell Method. It is shown that with the use of a parabolic interpolation in a vectorial problem, a convergence rate of 3.5 is obtained. Such a convergence rate compares with, or is even better than, the one obtained with FEM with the same interpolation – depending on the integration technique used by the FEM application. The accuracy of the solution is also comparable or better. More >

M.Ganesh¹, D. Sheen²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 183-194, 2001, DOI:10.3970/cmes.2001.002.183

Abstract A parallel numerical algorithm is introduced and analyzed for solving inhomogeneous initial-boundary value parabolic problems. The scheme is based on the method recently introduced in Sheen, Sloan, and Thomée (2000) for homogeneous problems. We give a method based on a suitable choice of multiple parameters. Our scheme allows one to compute solutions in a wide range of time. Instead of using a standard time-marching method, which is not easily parallelizable, we take the Laplace transform in time of the parabolic problems. The resulting elliptic problems can be solved in parallel. Solutions are then computed by More >

S.T. Lie¹, Guoyou Yu, Z. Zhao²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 171-182, 2001, DOI:10.3970/cmes.2001.002.171

Abstract The BEM/FEM coupling procedure is applied to 2-D time domain structural-acoustic interaction problems. The acoustic domain for fluid or air is modeled by BEM scheme that is suitable for both finite and infinite domains, while the structure is modeled by FEM scheme. The input impact, which can be either plane waves or non-plane waves, can either be forces acting directly on the structural-acoustic system or be explosion sources. The far field or near field explosion sources which are difficult to be simulated by finite element modeling, can be simulated exactly by boundary element modeling as More >