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

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.3, pp. 263-284, 2008, DOI:10.3970/cmes.2008.038.263

Abstract We present a wavelet based adaptive scheme and investigate the efficiency of this scheme for solving nearly singular potential PDEs. Multiresolution wavelet analysis (MRWA) provides a firm mathematical foundation by projecting the solution of PDE onto a nested sequence of approximation spaces. The wavelet coefficients then were used as an estimation of the sensible regions for node adaptation. The proposed adaptation scheme requires negligible calculation time due to the existence of the fast Discrete Wavelet Transform (DWT). Certain aspects of the proposed adaptive scheme are discussed through numerical examples. It has been shown that the More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 65-76, 2008, DOI:10.3970/cmes.2008.030.065

Abstract Trefftz method (TM) is one of widely used meshless numerical methods in elliptic type boundary value problems, of which the approximate solution is expressed as a linear combination of T-complete bases, and the unknown coefficients are determined from boundary conditions by solving a linear equations system. However, the accuracy of TM is severely limited by its ill-conditioning. This paper is a continuation of the work of Liu (2007a). The collocation TM is modified and applied to the direct and inverse problems of biharmonic equation in a simply connected plane domain. Due to its well-conditioning of More >

Shu Li¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.1, pp. 37-56, 2008, DOI:10.3970/cmes.2008.030.037

Abstract In this paper, a method based on a combination of an optimization of directions of orthotropy, along with topology optimization, is applied to continuum orthotropic solids with the objective of minimizing their compliance. The spatial discretization algorithm is the so called Meshless Local Petrov-Galerkin (MLPG) "mixed collocation'' method for the design domain, and the material-orthotropy orientation angles and the nodal volume fractions are used as the design variables in material optimization and topology optimization, respectively. Filtering after each iteration diminishes the checkerboard effect in the topology optimization problem. The example results are provided to illustrate More >

Leevan Ling¹, Tomoya Takeuchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 45-54, 2008, DOI:10.3970/cmes.2008.029.045

Abstract The method of fundamental solutions is coupled with the boundary control technique to solve the Cauchy problems of the Laplace Equations. The main idea of the proposed method is to solve a sequence of direct problems instead of solving the inverse problem directly. In particular, we use a boundary control technique to obtain an approximation of the missing Dirichlet boundary data; the Tikhonov regularization technique and the L-curve method are employed to achieve such goal stably. Once the boundary data on the whole boundary are known, the numerical solution to the Cauchy problem can be More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.2, pp. 77-94, 2008, DOI:10.3970/cmes.2008.028.077

Abstract The method of fundamental solutions (MFS) is a truly meshless numerical method widely used in the elliptic type boundary value problems, of which the approximate solution is expressed as a linear combination of fundamental solutions and the unknown coefficients are determined from the boundary conditions by solving a linear equations system. However, the accuracy of MFS is severely limited by its ill-conditioning of the resulting linear equations system. This paper is motivated by the works of Chen, Wu, Lee and Chen (2007) and Liu (2007a). The first paper proved an equivalent relation of the Trefftz… More >

Chih-Ping Wu¹, Kuan-Hao Chiu, Yun-Ming Wang

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 163-186, 2008, DOI:10.3970/cmes.2008.027.163

Abstract A differential reproducing kernel particle (DRKP) method is proposed and developed for the analysis of simply supported, multilayered elastic and piezoelectric plates by following up the consistent concepts of reproducing kernel particle (RKP) method. Unlike the RKP method in which the shape functions for derivatives of the reproducing kernel (RK) approximants are obtained by directly taking the differentiation with respect to the shape functions of the RK approximants, we construct a set of differential reproducing conditions to determine the shape functions for the derivatives of RK approximants. On the basis of the extended Hellinger-Reissner principle, More >

Shu Li¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.1, pp. 61-74, 2008, DOI:10.3970/cmes.2008.026.061

Abstract The Meshless Local Petrov-Galerkin (MLPG) "mixed collocation'' method is applied to the problem of topology-optimization of elastic structures. In this paper, the topic of compliance minimization of elastic structures is pursued, and nodal design variables which represent nodal volume fractions at discretized nodes are adopted. A so-called nodal sensitivity filter is employed, to prevent the phenomenon of checkerboarding in numerical solutions to the topology-optimization problems. The example results presented in the paper demonstrate the suitability and versatility of the MLPG "mixed collocation'' method, in implementing structural topology-optimization. More >

Tienfuan Kerh^1,2, J.S. Lai¹, D. Gunaratnam², R. Saunders²

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.1, pp. 1-12, 2008, DOI:10.3970/cmes.2008.026.001

Abstract Taiwan frequently suffers from strong ground motion, and the current building code is essentially based on two seismic zones, A and B. The design value of horizontal acceleration for zone A is 0.33g, and the value for zone B is 0.23g. To check the suitability of these values, a series of actual earthquake records are considered for evaluating peak ground acceleration (PGA) for each of the zones by using neural network models. The input parameters are magnitude, epicenter distance, and focal depth for each of the checking stations, and the peak ground acceleration is calculated… More >

G. Kosec¹, B. Šarler¹

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.3, pp. 197-208, 2008, DOI:10.3970/cmes.2008.025.197

Abstract This paper explores the application of the mesh-free Local Radial Basis Function Collocation Method (LRBFCM) in solution of coupled heat transfer and fluid flow problems in Darcy porous media. The involved temperature, velocity and pressure fields are represented on overlapping sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective derivatives of the RBF's. The energy and momentum equations are solved through explicit time stepping. The pressure-velocity coupling is calculated iteratively, with pressure correction, predicted from the local continuity equation violation.… 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 >