Nives Brajčić Kurbaša¹, Blaž Gotovac¹, Vedrana Kozulić¹

CMES-Computer Modeling in Engineering & Sciences, Vol.111, No.6, pp. 493-530, 2016, DOI:10.3970/cmes.2016.111.493

Abstract This paper presents exponential Atomic Basis Functions (ABF), which are called Eup(x,ω). These functions are infinitely differentiable finite functions that unlike algebraic up(x) basis functions, have an unspecified parameter - frequency w. Numerical experiments show that this class of atomic functions has good approximation properties, especially in the case of large gradients (Gibbs phenomenon). In this work, for the first time, the properties of exponential ABF are thoroughly investigated and the expression for calculating the value of the basis function at an arbitrary point of the domain is given in a form suitable for implementation in numerical analysis. Application of… More >

Guoqing Yang¹, Baotong Li^2,3, Yang Wang², Jun Hong²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.4, pp. 251-279, 2014, DOI:10.3970/cmes.2014.103.251

Abstract Mechanical behaviors arising at the contact interface largely depend on its surface topographies, particularly when it comes to rough surfaces. A numerical simulation based on an appropriate characterization of rough surfaces especially in terms of three dimensional can be of great significance when it comes to capturing the deformation patterns of micro-scale contacts. In this paper, a simple and practical scheme is developed to generate 3D rough surfaces and to analyze and evaluate the contact characteristics. Firstly amplitude and spatial statistical characterizations of asperities are introduced to avert from the redundancy of topography data caused by traditional measuring methods. A… More >

Xuefei He¹, Kian-Meng Lim^1,2,3, Siak-Piang Lim^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.1, pp. 21-48, 2008, DOI:10.3970/cmes.2008.035.021

Abstract The boundary element method (BEM) is known to have the advantage of reducing the dimension of problem by discretizing only the boundary of the domain. But it becomes less attractive for solving Poisson-type equations, due to the need to evaluate the domain integral which is computationally expensive. In this paper, we present the extension of a recently developed fast algorithm for Laplace equation, based on fast Fourier transform on multipoles (FFTM), to solve large scale 3D Poisson-type equations. We combined the Laplace solver with two fast methods for handling the domain integral based on fast Fourier transform (FFT). The first… More >

Fabian M. E. Duddeck¹

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.1, pp. 1-14, 2006, DOI:10.3970/cmes.2006.016.001

Abstract In general, the use of Boundary Element Methods (BEM) is restricted to physical cases for which a fundamental solution can be obtained. For simple differential operators (e.g. isotropic elasticity) these special solutions are known in their explicit form. Hence, the realization of the BEM is straight forward. For more complicated problems (e.g. anisotropic materials), we can only construct the fundamental solution numerically. This is normally done before the actual problem is tackled; the values of the fundamental solutions are stored in a table and all values needed later are interpolated from these entries. The drawbacks of this approach lie in… More >

B. Yang², V. K. Tewary³

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 165-178, 2006, DOI:10.3970/cmes.2006.015.165

Abstract Green's function (GF) modeling of defects may take effect only if the GF as well as its various integrals over a line, a surface and/or a volume can be efficiently evaluated. The GF is needed in modeling a point defect, while integrals are needed in modeling line, surface and volumetric defects. In a matrix of multilayered, generally anisotropic and linearly elastic and piezoelectric materials, the GF has been derived by applying 2D Fourier transforms and the Stroh formalism. Its use involves another two dimensions of integration in the Fourier inverse transform. A semi-analytical scheme has been developed previously for efficient… More >

Chun-Ying Wu^1,3, Xue-Yan Liu², A. Scarpas², Xiu-Run Ge³

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.2, pp. 149-158, 2006, DOI:10.3970/cmes.2006.012.149

Abstract For the spectral analysis of the three-dimensional multi-layered pavement, 3D layer spectral element method is presented to solve the problems of bounded layer system subjected to a transient load pulse. In spectral element, each layer is treated as one spectral element. The wave propagation inside each layer element is achieved by the superposition of the incident wave and the reflection wave. Fast Fourier transformation is used to transform FWD datum from time domain to frequency domain. The accuracy and efficiency of 3D layer spectral element approach were verified by analyzing the Falling weight deflectometer(FWD) testing model with the spectral methods… More >

J.C. Michel¹, H. Moulinec, P. Suquet

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 79-88, 2000, DOI:10.3970/cmes.2000.001.239

Abstract An iterative numerical method based on Fast Fourier Transforms has been proposed by \cite{MOU98} to investigate the effective properties of periodic composites. This iterative method is based on the exact expression of the Green function for a linear elastic, homogeneous reference material. When dealing with linear phases, the number of iterations required to reach convergence is proportional to the contrast between the phases properties, and convergence is therefore not ensured in the case of composites with infinite contrast (those containing voids or rigid inclusions or highly nonlinear materials). It is proposed in this study to overcome this difficulty by using… More >

A. Abbassi^1,2, N. Kechiche¹, H. Ben Aissia¹

FDMP-Fluid Dynamics & Materials Processing, Vol.10, No.3, pp. 319-341, 2014, DOI:10.3970/fdmp.2014.010.319

Abstract Jet transition towards a turbulent state is an interesting topic requiring a detailed analysis of the process leading to the onset and amplification of small flow disturbances. Here we examine experimentally the transition process for an isothermal laminar round free jet at low values of the Reynolds number. Close to the inlet nozzle, the turbulence intensity is assumed to be small enough so that the initial shear layer can be considered laminar and the velocity profile uniform. Experimental data are obtained using a Laser Doppler Anemometry (LDA) technique at various longitudinal and transversal coordinates, (x,y). Spectral analysis of the instantaneous… More >

D.T. Fullwood¹, S.R. Kalidindi², B.L. Adams¹, S. Ahmadi¹

CMC-Computers, Materials & Continua, Vol.9, No.1, pp. 25-40, 2009, DOI:10.3970/cmc.2009.009.025

Abstract Localization relations arise naturally in the formulation of multi-scale models. They facilitate statistical analysis of local phenomena that may contribute to failure related properties. The computational burden of dealing with such relations is high and recent work has focused on spectral methods to provide more efficient models. Issues with the inherent integrations in the framework have led to a tendency towards calibration-based approaches. In this paper a discrete Fourier transform framework is introduced, leading to an extremely efficient basis for the localization relations. Previous issues with the Green's function integrals are resolved, and the method is validated against finite element… More >

B. Yang², V. K. Tewary³

CMC-Computers, Materials & Continua, Vol.8, No.1, pp. 23-32, 2008, DOI:10.3970/cmc.2008.008.023

Abstract A three-dimensional Green's function for a material system consisting of anisotropic and linearly elastic planar multilayers with interfacial membrane and flexural rigidities has been derived. The Stroh formalism and two-dimensional Fourier transforms are applied to derive the general solution for each homogeneous layer. The Green's function for the multilayers is then solved by imposing the surface boundary condition, the interfacial displacement continuity condition, and the interfacial traction discontinuity condition. The last condition is given by the membrane and bending equilibrium equations of the interphases modeled as Kirchhoff plates. Numerical results that demonstrate the validity and efficiency of the formulation are… More >