R.D. Phillips¹, M.S. Hossain¹, L.T. Watson^1,2, R.H. Wynne³, Naren Ramakrishnan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.2, pp. 175-197, 2014, DOI:10.3970/cmes.2014.097.175

Abstract Clusters, typically mined by modeling locality of attribute spaces, are often evaluated for their ability to demonstrate ‘enrichment’ of categorical features. A cluster enrichment procedure evaluates the membership of a cluster for significant representation in predefined categories of interest. While classical enrichment procedures assume a hard clustering definition, this paper introduces a new statistical test that computes enrichments for soft clusters. Application of the new test to several scientific datasets is given. More >

Yuping Zhu^1,2, Tao Shi¹, Yuanbing Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.6, pp. 501-517, 2013, DOI:10.3970/cmes.2013.095.501

Abstract Single crystal ferromagnetic shape memory alloy is a kind of new intelligent materials, it shows obvious anisotropy. Micromechanics theory has been used to analyze the whole mechanical behaviors of this material. However, Eshelby’s tensor of this material which plays an important role has still not solved efficiently. Based on the existing micromechanics constitutive model, this paper analyzes the numerical calculation formula of Eshelby’s tensor of anisotropic ferromagnetic shape memory alloy. Adopting the way of Gauss integral, the optimal Gaussian integral points for different inclusion shapes and the corresponding numerical solution of Eshelby’s tensor are obtained.Furthermore, the influence of inclusion shapes… More >

G. Kosec^1,2, P. Zinterhof³

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.5, pp. 377-396, 2013, DOI:10.3970/cmes.2013.091.377

Abstract This paper deals with the implementation of the local meshless numerical method (LMM) on general purpose graphics processing units (GPU) in solving partial differential equations (PDE). The local meshless solution procedure is formulated in a way suitable for parallel execution and has been implemented on multiple GPUs. The implementation is tested on a solution of diffusion equation in a 2D domain. Different setups of the meshless approach regarding the selection of basis functions are tested on an interval up to 2.5 million of computational points. It is shown that monomials are a good selection of the basis when working with… More >

A.R.Ervilha^1,2, F.A.Sousa¹, J.M.C.Pereira¹, J.C.F.Pereira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.89, No.5, pp. 361-387, 2012, DOI:10.3970/cmes.2012.089.361

Abstract Faster than real time stochastic fire spread predictions are reported using a Non-Intrusive Spectral Projection (NISP) method based on Polynomial Chaos expansion and Graphic Processing Units (GPUs). The fireLib BEHAVE model together with a raster surface fire growth algorithm was implemented using the Compute Unified Device Architecture (CUDA) programming language. The uncertainty generated by the four random variables considered (wind speed, wind direction, fuel moisture, and fuel load) is quantified in the stochastic solution. Stochastic simulation of an idealized vegetation fire in a realistic complex terrain is obtained with speed-ups as high as 176 when compared to Central Processing Unit… More >

W.W. Zhang¹, Y.F. Nie¹, Y.Q. Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.6, pp. 507-530, 2012, DOI:10.3970/cmes.2012.088.507

Abstract A new algorithm for anisotropic triangular meshes generation in two dimension is presented. The inputs to the algorithm are the boundary geometry information and a metric tensor that specifies the desired element size and shape. The initial nodes are firstly distributed according to the above mentioned geometrical information, after bubble simulation, the optimized nodes set that meets the requirements of the metric tensor is obtained quickly. Then taking full advantage of the nodes set and the adjacency lists information provided by the process of node placement, a handful of non-satellite nodes are removed from the adjacency lists of the nodes… More >

J.A. García^1,2, J.M. Peña¹, S. McHugh², A. Jérusalem^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.5, pp. 411-432, 2012, DOI:10.3970/cmes.2012.087.411

Abstract Neuronal growth is a complex process involving many intra- and extracellular mechanisms which are collaborating conjointly to participate to the development of the nervous system. More particularly, the early neocortical development involves the creation of a multilayered structure constituted by neuronal growth (driven by axonal or dendritic guidance cues) as well as cell migration. The underlying mechanisms of such structural lamination not only implies important biochemical changes at the intracellular level through axonal microtubule (de)polymerization and growth cone advance, but also through the directly dependent stress/stretch coupling mechanisms driving them. Efforts have recently focused on modeling approaches aimed at accounting… More >

Ali Bekir Yildiz¹

CMES-Computer Modeling in Engineering & Sciences, Vol.75, No.1, pp. 33-42, 2011, DOI:10.3970/cmes.2011.075.033

Abstract The generalized method for determination of time constants of linear circuits is introduced. Nodal and mesh analysis, conventional methods whose applications are simpler than the state-space formulation, are used in obtaining the system equations. The approach is based on the use of the relationship between transfer functions and system equations of linear circuits, obtained by the conventional methods. The examples of active and passive circuits are given to illustrate the method. More >

Peter J. Ireman, Anders Klarbring¹, Niclas Strömberg

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.2, pp. 125-158, 2009, DOI:10.3970/cmes.2009.052.125

Abstract In this paper finite element approaches for fretting fatigue are proposed on the basis of a non-local model of continuum damage coupled to friction and wear. The model is formulated in the frame-work of a standard material. In a previous paper this was done in the spirit of Maugin, where an extra entropy flux is introduced in the second law in order to include the gradient of the internal variable in a proper manner. In this paper we follow instead the ideas of Frémond and others, where this extra entropy flux is no longer needed, but instead new non-classical balance… More >

Z. C. Li², T. T. Lu³, H. T. Huang⁴, A. H.-D. Cheng⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.1, pp. 39-82, 2009, DOI:10.3970/cmes.2009.052.039

Abstract For Laplace's equation and other homogeneous elliptic equations, when the particular and fundamental solutions can be found, we may choose their linear combination as the admissible functions, and obtain the expansion coefficients by satisfying the boundary conditions only. This is known as the Trefftz method (TM) (or boundary approximation methods). Since the TM is a meshless method, it has drawn great attention of researchers in recent years, and Inter. Workshops of TM and MFS (i.e., the method of fundamental solutions). A number of efficient algorithms, such the collocation algorithms, Lagrange multiplier methods, etc., have been developed in computation. However, there… More >

Paola Gardano¹, Peter Dabnichki^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.2, pp. 87-98, 2006, DOI:10.3970/cmes.2006.015.087

Abstract The work presents a numerical simulation of hydrodynamic forces generated in front crawl swimming. The three dimensional Laplace's equation is used for the analysis of the flow around a moving body in an infinite domain and considers the effect of the added mass and the acceleration on the hydrodynamic forces (Drag and Lift) generated by the interaction between the flow and the body at different geometric configurations of the arm -- variable elbow angle. Boundary Element Method (BEM) was used to obtain the solution of the three dimensional equation numerically. The aim of the work was two-fold:
1) to… More >