E. Y. K. Ng¹, S. T. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.1, pp. 13-28, 2008, DOI:10.3970/cmes.2008.023.013

Abstract Microchannels have been recognized as a very effective chemical separation and heat transfer device. The electrical double layer (EDL) effect in a micro-scale flow is however anticipated to be critical. In this paper, Nernst-Planck model (NPM), is used to predict the ion concentration distribution as it is reported to be a more appropriate model for developing microchannel flow. The governing equations are discretised for developing rectangular microchannel flows in Cartesians coordinate. An additional body force source term that is relating to the electric potential, resulted from the EDL effect is introduced in the conventional z-axis momentum More >

Xiang Jiawei¹, Chen Xuefeng², Yang Lianfa³, He Zhengjia⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.1, pp. 1-12, 2008, DOI:10.3970/cmes.2008.023.001

Abstract A class of flat shell elements is constructed by using the scaling functions of two-dimensional tensor product B-spline wavelet on the interval (BSWI). Unlike the process of direct wavelets adding in the wavelet Galerkin method, the element displacement field represented by the coefficients of wavelets expansions was transformed from wavelet space into physical space via the constructed two-dimensional transformation matrix. Then, the BSWI flat shell element is constructed by the assembly of BSWI plane elastomechanics and Mindlin plate elements. Because of the good character of BSWI scaling functions, the BSWI flat shell element combine the More >

Ya Yan Lu¹, Jianxin Zhu²

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 235-248, 2007, DOI:10.3970/cmes.2007.022.235

Abstract The perfectly matched layer (PML) is a widely used technique for truncating unbounded domains in numerical simulations of wave propagation problems. In this paper, the PML technique is used with a standard one-way model to solve a benchmark problem for underwater acoustics modeling. Accurate solutions are obtained with a PML layer with a thickness of only a quarter of the wavelength. The effect of a PML is analyzed in a perturbation analysis for waveguides. More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², P. Solek³

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 217-234, 2007, DOI:10.3970/cmes.2007.022.217

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for the solution of boundary value problems for coupled thermo-electro-mechanical fields. Transient dynamic governing equations are considered here. To eliminate the time-dependence in these equations, the Laplace-transform technique is applied. Material properties of piezoelectric materials are influenced by a thermal field. It is leading to an induced nonhomogeneity and the governing equations are more complicated than in a homogeneous counterpart. Two-dimensional analyzed domain is subdivided into small circular subdomains surrounding nodes randomly spread over the whole domain. A unit step function is used as More >

Cha’o-Kuang Chen¹, Shing-Cheng Chang¹, Szu-Yu Sun¹

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 203-216, 2007, DOI:10.3970/cmes.2007.022.203

Abstract In this study, the channel flow is discussed by the LBM simulations. In the cases of channel with obstacles inside, the square pillars play the role of causing interruption within the fluid field, and hence change the direction of fluid flow. The recirculation region is formed behind the obstacles and influences the fluid passed through not only in the velocity field but also in the temperature field. Therefore, heat transfer is enhanced in local region.
The field synergy principle is applied in the research to demonstrate that the increased interruption within the fluid increases the More >

S. Namilae¹, U. Ch,ra², A Srinivasan³, N. Ch,ra⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 189-202, 2007, DOI:10.3970/cmes.2007.022.189

Abstract Molecular dynamics (MD) simulations play an important predictive role in understanding the behavior of nanoscale systems. In this paper, parallel MD simulations are used to understand the mechanical behavior of interfaces in CNT based composites. We present an algorithm for parallel implementation of MD simulations of carbon nanotube (CNT) based systems using reactive bond order potentials. We then use that algorithm to model the CNT-polymer interfaces with various levels of interaction as (a) described only by long range Van Der Waals interactions (b) chemically bonded with fixed matrix and (c) chemically bonded with matrix explicitly More >

J.N. Johnson¹, J.M. Owen²

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 165-188, 2007, DOI:10.3970/cmes.2007.022.165

Abstract In this paper, we propose a Meshless Local Petrov-Galerkin method for studying the diffusion of a magnetic field within a non-magnetic (μ = μ₀) conducting medium with non-homogeneous and anisotropic electrical resistivity. We derive a local weak form for the magnetic diffusion equation and discuss the effects of different trial/test functions and nodal spacings on its solution. We then demonstrate that the method produces convergent results for several relevant one-dimensional test problems for which solutions are known. This method has the potential to be combined with other mesh-free methods such as Smoothed Particle Hydrodynamics (SPH) to More >

R. Criado¹, J.E. Ortiz¹, V. Mantič¹, L.J. Gray^1,2, F. París¹

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.2, pp. 151-164, 2007, DOI:10.3970/cmes.2007.022.151

Abstract A numerical implementation of the Somigliana identity in displacements for the solution of 3D elastic problems in exponentially graded isotropic solids is presented. An expression for the fundamental solution in displacements, U_jl, was deduced by Martin et al. (Proc. R. Soc. Lond. A, 458, pp. 1931--1947, 2002). This expression was recently corrected and implemented in a Galerkin indirect 3D BEM code by Criado et al. (Int. J. Numer. Meth. Engng., 2008). Starting from this expression of U_jl, a new expression for the fundamental solution in tractions T_jl has been deduced in the present work. These quite complex expressions of More >

D. Moens¹, M. De Munck, D. Vandepitte

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.2, pp. 129-150, 2007, DOI:10.3970/cmes.2007.022.129

Abstract Recently, an interval finite element methodology has been developed to calculate envelope frequency response functions of uncertain structures with interval parameters. The methodology is based on a hybrid interval implementation of the modal superposition principle. This hybrid procedure consists of a preliminary optimization step, followed by an interval arithmetic procedure. The final envelope frequency response functions have been proved to give a very good approximation of the actual response range of the interval problem. Initially, this method was developed for undamped structures. Based on the theoretical principles of this approach, this paper introduces a new… More >

A.N. Guz¹, V.V. Zozulya²

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.2, pp. 119-128, 2007, DOI:10.3970/cmes.2007.022.119

Abstract The frictional contact interaction of the edges of a finite plane crack is studied for the case of normal incidence of a harmonic SH-shear wave which produces antiplane deformation. The forces of contact interaction and displacement discontinuity are analyzed. Influence of the wave frequency on the stress intensity factor for different coefficients of friction is studied here. More >