V. Giusti ¹, B. Montagnini ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.3, pp. 229-255, 2014, DOI:10.3970/cmes.2014.102.229

Abstract An application of a 3D Boundary Element Method (BEM), coupled with the Response Matrix (RM) technique, to solve the neutron diffusion and transport equations for multi-region domains is presented. The discussion is here limited to steady state problems, in which the neutrons have a wide energy spectrum, which leads to systems of several diffusion or transport equations. Moreover, the number of regions with different physical constants can be very large. The boundary integral equations concerning each region are solved via a polynomial moment expansion and, taking advantage of suitable recurrence formulas, the multi-fold integrals there… More >

E.F. Fontes Jr¹, J.A.F. Santiago¹, J.C.F. Telles¹

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.3, pp. 211-228, 2014, DOI:10.3970/cmes.2014.102.211

Abstract An iterative coupling procedure using different meshless methods is presented to solve linear elastic fracture mechanic (LEFM) problems. The domain of the problem is decomposed into two sub-domains, where each one is addressed using an appropriate meshless method. The method of fundamental solutions (MFS) based on the numerical Green’s function (NGF) procedure to generate the fundamental solution has been chosen for modeling embedded cracks in the elastic medium and the meshless local Petrov-Galerkin (MLPG) method has been chosen for modeling the remaining sub-domain. Each meshless method runs independently, coupled with an iterative update of interface More >

L. Rodríguez-Tembleque¹, M.H. Aliabadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.3, pp. 183-210, 2014, DOI:10.3970/cmes.2014.102.183

Abstract This work presents new contact constitutive laws for friction and wear modelling in fiber-reinforced plastics (FRP). These laws are incorporated to a numerical methodology which allows us to solve the contact problem taking into account the anisotropic tribological properties on the interfaces. This formulation uses the Boundary Element Method for computing the elastic influence coefficients. Furthermore, the formulation considers micromechanical models for FRP that also makes it possible to take into account the fiber orientation relative to the sliding direction, the fiber volume fraction, the aspect ratio of fibers, or the fiber arrangement. The proposed More >

Xiaomin Wang^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.2, pp. 169-182, 2014, DOI:10.3970/cmes.2014.102.169

Abstract In this paper, an efficient and robust wavelet Laplace inversion method of solving the fractional differential equations is proposed. Such an inverse function can be applied to any reasonable function categories and it is not necessary to know the properties of original function in advance. As an example, we have applied the proposed method to the solution of the Bagley–Torvik equations and Numerical examples are given to demonstrate the efficiency and accuracy of the proposed. More >

Z.W. Wu¹, J.K. Liu¹, Z.Q. Liu^1,2, Z.R. Lu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.2, pp. 149-168, 2014, DOI:10.3970/cmes.2014.102.149

Abstract This paper presents numerical analysis for the mooring system with nonlinear elastic mooring cables. The equation of motion for nonlinear elastic mooring cable is established by utilizing finite element method. A marine mooring system of floating rectangular box with nonlinear elastic cables is taken as an illustrative example. The dynamic analysis, static analysis, and uniformity analysis are carried out for the polyester mooring system and the results are compared with those of the steel wire and the chain mooring system. Results from the present study can provide valuable recommendations for the design and construction of More >

Jizeng Wang^1,2, Lei Zhang¹, Youhe Zhou¹

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.2, pp. 127-148, 2014, DOI:10.3970/cmes.2014.102.127

Abstract In this paper, we propose an efficient wavelet method for numerical solution of nonlinear integral equations with singular kernels. The proposed method is established based on a function approximation algorithm in terms of Coiflet scaling expansion and a special treatment of boundary extension. The adopted Coiflet bases in this algorithm allow each expansion coefficient being explicitly expressed by a single-point sampling of the function, which is crucially important for dealing with nonlinear terms in the equations. In addition, we use the technique of integration by parts to transform the original integral equations with non-smooth or More >

Jin Li^1,2, De-hao Yu^3,4

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.2, pp. 103-126, 2014, DOI:10.3970/cmes.2014.102.103

Abstract The generalized composite middle rectangle rule for the computation of Hilbert integral is discussed. The pointwise superconvergence phenomenon is presented, i.e., when the singular point coincides with some a priori known point, the convergence rate of the rectangle rule is higher than what is global possible. We proved that the superconvergence rate of the composite middle rectangle rule occurs at certain local coordinate of each subinterval and the corresponding superconvergence error estimate is obtained. By choosing the superconvergence point as the collocation points, a collocation scheme for solving the relevant Hilbert integral equation is presented More >

T. Ghisu¹, B. Arca², G. Pellizzaro², P. Duce²

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.1, pp. 83-102, 2014, DOI:10.3970/cmes.2014.102.083

Abstract Level-set approaches are efficient and versatile methods for solving interface tracking problems and have been used in recent years to describe wildland fire propagation. Being based on an Eulerian description of the spread problem, their numerical implementation offers improved computational agility and better portability to parallel computing environments with respect to vector-based simulators. The use of a continuous representation of the fire perimeter in place of the binary formulation used in Cellular Automata avoids the commonly observed distortion of the fire shape. This work presents an algorithm for fire-spread simulation based on a level-set formulation. More >

Weiwei Zhang¹, Yufeng Nie¹, Li Cai¹, Nan Qi²

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.1, pp. 55-82, 2014, DOI:10.3970/cmes.2014.102.055

Abstract In this paper, we derive an adaptive mesh generation method for discretizing the incompressible flow using node-based local grids. The flow problem is described by the Stokes equations which are solved by a stabilized low-order P1-P1 (linear velocity, linear pressure) mixed finite element method. The proposed node-based adaptive mesh generation method consists of four components: mesh size modification, a node placement procedure, a node-based local mesh generation strategy and an error estimation technique, which are combined so as to guarantee obtaining a conforming refined/coarsened mesh. The nodes are considered as particles with interaction forces, which… More >

G Devaraj¹, Shashi Narayan¹, Debasish Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.1, pp. 1-54, 2014, DOI:10.3970/cmes.2014.102.001

Abstract This work sets forth a 'hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (C^p-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization More >