Tao Jiang^1,2, Yuan-Sheng Tang¹, Jin-Lian Ren^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.4, pp. 249-297, 2014, DOI:10.3970/cmes.2014.101.249

Abstract In this work, a corrected three-dimensional smoothed particle hydrodynamics (CSPH-3D) method is proposed to simulate the polymer free surface flows in the filling process based on the eXtended Pom-Pom (XPP) model, and some complex deformation phenomena are also numerically predicted. The proposed CSPH-3D method is mainly motivated by a coupled concept that an extended kernel-gradient-corrected SPH (KGC-SPH) method is used in the interior of fluid flow and the traditional SPH (TSPH) method is used near the boundary domain. The present 3D particle method has higher accuracy and better stability than the TSPH-3D method. Meanwhile, a density diffusive term is introduced… More >

N. Mitsume¹, S. Yoshimura¹, K. Murotani¹, T. Yamada¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.4, pp. 229-247, 2014, DOI:10.3970/cmes.2014.101.229

Abstract The MPS-FE method, which adopts the Finite Element (FE) method for structure computation and the Moving Particle Simulation (MPS) method for fluid computation involving free surfaces, was developed to solve fluid-structure interaction problems with free surfaces. The conventional MPS-FE method, in which MPS wall boundary particles and finite elements are overlapped in order to exchange information at a fluid-structure interface, is not versatile and reduces the advantages of the software modularity. In this study, we developed a nonoverlapping approach in which the interface in the fluid computation corresponds to the interface in the structure computation through an MPS polygon wall… More >

W. Kaewjuea¹, T. Senjuntichai², R.K.N.D. Rajapakse³

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.3, pp. 207-228, 2014, DOI:10.3970/cmes.2014.101.207

Abstract Dynamic response of an infinite cylindrical borehole in a poroelastic medium with an excavation disturbed zone is investigated in this paper. The borehole is subjected to axisymmetric time-harmonic loads and fluid sources applied to its surface, which is either fully permeable or impermeable. The governing equations based on Biot’s poroelastodynamics theory are solved by using two scalar potentials and two vector potentials. The general solutions are then derived through the application of Fourier integral transform with respect to the vertical coordinate. An exact stiffness matrix scheme is established from the derived general solutions to include the excavation disturbed zone. Boundary… More >

Yuan Li¹, Jianming Zhang^1,2, Guizhong Xie¹, Xingshuai Zheng¹, Shuaiping Guo¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.3, pp. 187-206, 2014, DOI:10.3970/cmes.2014.101.187

Abstract In this paper, a time-domain boundary element method formulation for the analysis of three-dimensional elastodynamic problems with arbitrary, non-null initial conditions is presented. The formulation is based on the convolution quadrature method, by which the numerical stability is improved significantly. In order to take into account the non-null initial conditions in this formulation, a general method is developed to replace the initial conditions by equivalent pseudo-forces based on the pseudo-force method. The original governing equation is transformed into a new one subjected to null initial conditions. In the numerical examples, longitudinal vibrations of a free beam and a cantilevered beam… More >

W.M. Abd- Elhameed^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.3, pp. 159-185, 2014, DOI:10.3970/cmes.2014.101.159

Abstract This paper is concerned with developing two new algorithms for direct solutions of linear and nonlinear sixth-order two point boundary value problems. These algorithms are based on the application of the two spectral methods namely, collocation and Petrov-Galerkin methods. The suggested algorithms are completely new and they depend on introducing a novel operational matrix of derivatives which is expressed in terms of the well-known harmonic numbers. The basic idea for the suggested algorithms rely on reducing the linear or nonlinear sixth-order boundary value problem governed by its boundary conditions to a system of linear or nonlinear algebraic equations which can… More >

Shantaram Parab¹, Shreya Srivastava², Pijush Samui³, A. Ramachandra Murthy⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.2, pp. 139-158, 2014, DOI:10.3970/cmes.2014.101.139

Abstract This paper studies the applicability of Gaussian Process Regression (GPR) and Least Squares Support Vector Machines (LSSVM) to predict fracture parameters and failure load (P_max) of high strength and ultra-high strength concrete beams. Fracture characteristics include fracture energy (G_F), critical stress intensity factor (K_IC) and critical crack tip opening displacement (CTOD_C) Mathematical models have been developed in the form of relation between several input variables such as beam dimensions, water cement ratio, compressive strength, split tensile strength, notch depth, modulus of elasticity and output fracture parameters. Four GPR and four LSSVM models have been developed using MATLAB software for training… More >

M. Mohammadi¹, R. Mokhtari^2,3, R. Schaback⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.2, pp. 113-138, 2014, DOI:10.3970/cmes.2014.101.113

Abstract In this paper, the two-dimensional (2D) Brusselator reaction-diffusion system is simulated numerically by the method of lines. The proposed method is implemented as a meshless method based on spatial trial functions in the reproducing kernel Hilbert spaces. For efficiency and stability reasons, we use the Newton basis introduced recently by Müller and Schaback. The method is shown to work in all interesting situations described by Hopf bifurcations and Turing patterns. More >

Lifeng Wang¹, Yunpeng Ma^1,2, Yongqiang Yang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.2, pp. 97-111, 2014, DOI:10.3970/cmes.2014.101.097

Abstract In this paper, a numerical method based on the Legendre polynomials is presented for a class of variable order fractional differential equation. We adopt the Coimbra variable order fractional operator, which can be viewed as a Caputo-type definition. Three different kinds of operational matrixes with Legendre polynomials are derived. A truncated the Legendre polynomials series together with the products of several dependent matrixes are utilized to reduce the variable order fractional differential equation to a system of algebraic equations. The solution of this system gives the approximation solution for the truncated limited n. An error analysis technique is also given.… More >

S.A. Khuri¹, A. Sayfy¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.2, pp. 81-96, 2014, DOI:10.3970/cmes.2014.101.081

Abstract In this article, we introduce a numerical scheme to solve a class of nonlinear two-point BVPs on a semi-infinite domain that arise in engineering applications and the physical sciences. The strategy is based on replacing the boundary condition at infinity by an asymptotic boundary condition (ABC) specified over a finite interval that approaches the given value at infinity. Then, the problem complimented with the resulting ABC is solved using a fourth order spline collocation approach constructed over uniform meshes on the truncated domain. A number of test examples are considered to confirm the accuracy, efficient treatment of the boundary condition… More >

ShiXu Guo¹, Shili Chen¹, Xinjing Huang¹, Yu Zhang¹, Shijiu Jin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.1, pp. 59-80, 2014, DOI:10.3970/cmes.2014.101.059

Abstract Spherical leak detectors can detect very tiny leakage in pipelines and have low risk of blockage. In this paper the passing ability of the detector in the vertical segment of a pipe was studied using CFD simulations and experiments. The Reynolds number for the sphere exceeds 10⁴ at the economical velocity range for oil pipelines, and there were few researches related to the hydrodynamic force on the sphere by the pipe flow at high Reynolds number. For sphere with different sizes and density, and at different flow rates, more than 100 3-D steady numerical simulations were carried out. The simulation… More >