Yiming Chen, Mingxu Yi, Chen Chen, Chunxiao Yu

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.6, pp. 639-654, 2012, DOI:10.3970/cmes.2012.083.639

Abstract In this paper, Bernstein polynomials method is proposed for the numerical solution of a class of space-time fractional convection-diffusion equation with variable coefficients. This method combines the definition of fractional derivatives with some properties of Bernstein polynomials and are dispersed the coefficients efficaciously. The main characteristic behind this method is that the original problem is translated into a Sylvester equation. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples show that the method is effective. More >

Ping Hu¹, Yang Xia¹, Limin Tang²

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.2, pp. 103-136, 2011, DOI:10.3970/cmes.2011.073.103

Abstract In this paper, an assumed displacement quasi-conforming finite element method with truncated polynomial expansions of in-domain displacements and derived expansions of strains is introduced. Based on the method a four-node quadrilateral flat shell element with complete quadratic polynomials for membrane and bending displacement fields is developed. Numerical tests are carried out for validation of the present element. The results show that the present element preserves all the advantages of the quasi-conforming i.e., explicit stiffness matrix, convenient post processing and free from membrane locking and shear locking. The tests also prove that the present element gives excellent results, especially for the… More >

J.L. Sereno¹, J.C.F. Pereira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.1, pp. 23-44, 2011, DOI:10.3970/cmes.2011.073.023

Abstract The evolution of a two-dimensional vortex pair in ground effect was studied under the influence of random initial inputs comprising vortex strength (circulation) or initial vortex position. The paper addresses the questions of how do variations and uncertainties of initial conditions translate to the variability of vortex pair evolution. The stochastic solutions were obtained recurring to the Polynomial Chaos Expansion method of random processes applied to the Navier-Stokes equations for a laminar flow. The method quantifies the extent, dependence and propagation of uncertainty through the model system and, in particular, a methodology for the calculation of the vortices trajectory variability,… More >

Sunilkumar N¹, D Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.3, pp. 249-284, 2010, DOI:10.3970/cmes.2010.066.249

Abstract We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and 1D NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of… More >

Sunilkumar N¹, D Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.2, pp. 107-154, 2010, DOI:10.3970/cmes.2010.065.107

Abstract The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries. Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials. There is thus a case for combining these advantages in a so-called hybrid scheme or a 'smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C^p(p ≥ 1) continuity. One such recent attempt, a… More >

Jun-Sheng Duan¹, Ai-Ping Guo²

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.2, pp. 139-150, 2010, DOI:10.3970/cmes.2010.060.139

Abstract Adomian polynomials are constituted of reduced polynomials and derivatives of nonlinear operator. The reduced polynomials are independent of the form of the nonlinear operator. A recursive algorithm of the reduced polynomials is discovered and its symbolic implementation by the software Mathematica is given. As a result, a new and convenient algorithm for the Adomian polynomials is obtained. More >

Chia-Cheng Tsai¹, Chein-Shan Liu², Wei-Chung Yeih³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.2, pp. 131-152, 2010, DOI:10.3970/cmes.2010.056.131

Abstract The fictitious time integration method (FTIM) previously developed by Liu and Atluri (2008a) is combined with the method of fundamental solutions and the Chebyshev polynomials to solve Poisson-type nonlinear PDEs. The method of fundamental solutions with Chebyshev polynomials (MFS-CP) is an exponentially-convergent meshless numerical method which is able to solving nonhomogeneous partial differential equations if the fundamental solution and the analytical particular solutions of the considered operator are known. In this study, the MFS-CP is extended to solve Poisson-type nonlinear PDEs by using the FTIM. In the solution procedure, the FTIM is introduced to convert a Poisson-type nonlinear PDE into… More >

M. A. Kelmanson¹ and M. C. Tenwick¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 191-210, 2010, DOI:10.3970/cmes.2010.055.191

Abstract A method is presented for improving the accuracy of the widely used Gauss-Legendre Nyström method for determining approximate solutions of Fredholm integral equations of the second kind on finite intervals. The authors' recent continuous-kernel approach is generalised in order to accommodate kernels that are either singular or of limited continuous differentiability at a finite number of points within the interval of integration. This is achieved by developing a Gauss-Jacobi Nyström method that moreover includes a mean-value estimate of the truncation error of the Hermite interpolation on which the quadrature rule is based, making it particularly accurate at low orders. A… More >

Hua Li¹, Shantanu S. Mulay¹, Simon See²

CMES-Computer Modeling in Engineering & Sciences, Vol.54, No.2, pp. 147-200, 2009, DOI:10.3970/cmes.2009.054.147

Abstract In this paper, the stability characteristics of a novel strong-form meshless method, called the random differential quadrature (RDQ), are studied using the location of zeros or roots of its characteristic polynomials with respect to unit circle in complex plane by discretizing the domain with the uniform or random field nodes. This is achieved by carrying out the RDQ method stability analysis for the 1st-order wave, transient heat conduction and transverse beam deflection equations using both the analytical and numerical approaches. The RDQ method extends the applicability of the differential quadrature (DQ) method over irregular domain, discretized by randomly or uniformly… More >

Chia-Cheng Tsai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.45, No.3, pp. 249-272, 2009, DOI:10.3970/cmes.2009.045.249

Abstract Analytical particular solutions of Chebyshev polynomials are obtained for problems of Reissner plates under arbitrary loadings, which are governed by three coupled second-ordered partial differential equation (PDEs). Our solutions can be written explicitly in terms of monomials. By using these formulas, we can obtain the approximate particular solution when the arbitrary loadings have been represented by a truncated series of Chebyshev polynomials. In the derivations of particular solutions, the three coupled second-ordered PDE are first transformed into a single six-ordered PDE through the Hörmander operator decomposition technique. Then the particular solutions of this six-ordered PDE can be found in the… More >