Chia-Cheng Tsai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 151-162, 2008, DOI:10.3970/cmes.2008.027.151

Abstract In this paper we develop analytical particular solutions for the polyharmonic and the products of Helmholtz-type partial differential operators with Chebyshev polynomials at right-hand side. Our solutions can be written explicitly in terms of either monomial or Chebyshev bases. By using these formulas, we can obtain the approximate particular solution when the right-hand side has been represented by a truncated series of Chebyshev polynomials. These formulas are further implemented to solve inhomogeneous partial differential equations (PDEs) in which the homogeneous solutions are complementarily solved by the method of fundamental solutions (MFS). Numerical experiments, which include eighth order PDEs and three-dimensional… 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. This formulation does not require… 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 compare favourably with those obtained… 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 >

P. Orsini¹, H. Power^1,2, H. Morvan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 187-208, 2008, DOI:10.3970/cmes.2008.023.187

Abstract In this work, we present a modified Control Volume (CV) method that uses a Radial Basis Function (RBF) interpolation to improve the prediction of the flux accuracy at the faces of the CV. The method proposed differs from classical CV methods in the way that the flux at the cell surfaces is computed. A local RBF interpolation of the field variable is performed at the centres of the cell being integrated and its neighbours. This interpolation is then used to reconstruct the solution and its gradient in the integration points which support the flux computation. In addition, it is required… More >

E.J. Sellountos¹, A. Sequeira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 127-148, 2008, DOI:10.3970/cmes.2008.023.127

Abstract In this work a hybrid velocity-vorticity scheme for the solution of the 2D Navier-Stokes equations is presented. The multi-region Local Boundary Integral Equation (LBIE) combined with Radial Basis Functions (RBF) interpolation is used for the solution of the kinematics and the multi-region BEM for the solution of the transport kinetics. The final system of equations is in band form for both methods. The issue of RBF discontinuities is resolved by constructing the RBF matrix locally in every region. The kinematics integral equation is used in three different forms, for coupling the velocity field on the boundary, on interior points and… More >

P. H. Wen¹, Y. C. Hon²

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.3, pp. 177-192, 2007, DOI:10.3970/cmes.2007.021.177

Abstract In this paper, we perform a geometrically nonlinear analysis of Reissner-Mindlin plate by using a meshless collocation method. The use of the smooth radial basis functions (RBFs) gives an advantage to evaluate higher order derivatives of the solution at no cost on extra-interpolation. The computational cost is low and requires neither the connectivity of mesh in the domain/boundary nor integrations of fundamental/particular solutions. The coupled nonlinear terms in the equilibrium equations for both the plane stress and plate bending problems are treated as body forces. Two load increment schemes are developed to solve the nonlinear differential equations. Numerical verifications are… More >

S.Y. Wang^1,2, K.M. Lim^2,3, B.C. Khoo^2,3, M.Y. Wang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.1, pp. 1-40, 2007, DOI:10.3970/cmes.2007.021.001

Abstract The level set method is a numerical technique for simulating moving interfaces. In this paper, an unconditionally BIBO (Bounded-Input-Bounded-Output) time-stable consistent meshfree level set method is proposed and applied as a more effective approach to simultaneous shape and topology optimization. In the present level set method, the meshfree infinitely smooth inverse multiquadric Radial Basis Functions (RBFs) are employed to discretize the implicit level set function. A high level of smoothness of the level set function and accuracy of the solution to the Hamilton-Jacobi partial differential equation (PDE) can be achieved. The resulting dynamic system of coupled Ordinary Differential Equations (ODEs)… More >

N. Mai-Duy¹, L. Mai-Cao², T. Tran-Cong³

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.1, pp. 59-78, 2007, DOI:10.3970/cmes.2007.018.059

Abstract In this paper, an indirect/integrated radial-basis-function network (IRBFN) method is further developed to solve transient partial differential equations (PDEs) governing fluid flow problems. Spatial derivatives are discretized using one- and two-dimensional IRBFN interpolation schemes, whereas temporal derivatives are approximated using a method of lines and a finite-difference technique. In the case of moving interface problems, the IRBFN method is combined with the level set method to capture the evolution of the interface. The accuracy of the method is investigated by considering several benchmark test problems, including the classical lid-driven cavity flow. Very accurate results are achieved using relatively low numbers… More >

S. N. Atluri¹, Z. D. Han¹, S. Shen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 507-518, 2003, DOI:10.3970/cmes.2003.004.507

Abstract The general Meshless Local Petrov-Galerkin (MLPG) type weak-forms of the displacement & traction boundary integral equations are presented, for solids undergoing small deformations. These MLPG weak forms provide the most general basis for the numerical solution of the non-hyper-singular displacement and traction BIEs [given in Han, and Atluri (2003)], which are simply derived by using the gradients of the displacements of the fundamental solutions [Okada, Rajiyah, and Atluri (1989a,b)]. By employing the various types of test functions, in the MLPG-type weak-forms of the non-hyper-singular dBIE and tBIE over the local sub-boundary surfaces, several types of MLPG/BIEs are formulated, while also… More >