H. Nguyen-Van¹, N. Mai-Duy², T. Tran-Cong³

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 209-222, 2008, DOI:10.3970/cmes.2008.023.209

Abstract This paper reports a study of linear elastic analysis of two-dimensional piezoelectric structures using a smoothed four-node piezoelectric element. The element is built by incorporating the strain smoothing method of mesh-free conforming nodal integration into the standard four-node quadrilateral piezoelectric finite element. The approximations of mechanical strains and electric potential fields are normalized using a constant smoothing function. This allows the field gradients to be directly computed from shape functions. No mapping or coordinate transformation is necessary so that the element can be used in arbitrary shapes. Through several examples, the simplicity, efficiency and reliability More >

P. Orsini¹, H. Power^1,2, H. Morvan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 187-208, 2008, DOI:10.3970/cmes.2008.023.187

Abstract In this work, we present a modified Control Volume (CV) method that uses a Radial Basis Function (RBF) interpolation to improve the prediction of the flux accuracy at the faces of the CV. The method proposed differs from classical CV methods in the way that the flux at the cell surfaces is computed. A local RBF interpolation of the field variable is performed at the centres of the cell being integrated and its neighbours. This interpolation is then used to reconstruct the solution and its gradient in the integration points which support the flux computation. More >

Shueei-Muh Lin¹, Sen-Yung Lee², Yu-Sheng Lin³

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 175-186, 2008, DOI:10.3970/cmes.2008.023.175

Abstract The blade of a horizontal-axis wind power turbine is modeled as a rotating beam with pre-cone angles and setting angles. Based on the Bernoulli-Euler beam theory, without considering the axial extension deformation and the Coriolis forces effect, the governing differential equations for the bending vibration of the beam are derived. It is pointed out that if the geometric and the material properties of the beam are in polynomial forms, then the exact solution for the system can be obtained. Based on the frequency relations as revealed, without tedious numerical analysis, one can reach many general More >

Liviu Marin¹, Henry Power¹, Richard W. Bowtell², Clemente Cobos Sanchez², Adib A. Becker¹, Paul Glover²，Arthur Jones¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 149-174, 2008, DOI:10.3970/cmes.2008.023.149

Abstract We investigate the reconstruction of a divergence-free surface current distribution from knowledge of the magnetic flux density in a prescribed region of interest in the framework of static electromagnetism. This inverse problem is motivated by the design of gradient coils for use in magnetic resonance imaging (MRI) and is formulated using its corresponding integral representation according to potential theory. A novel boundary element method (BEM) which employs linear interpolation on quadratic surfaces and also satisfies the continuity equation for the current density, i.e. a divergence-free BEM, is presented. Since the discretised BEM system is ill-posed More >

E.J. Sellountos¹, A. Sequeira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 127-148, 2008, DOI:10.3970/cmes.2008.023.127

Abstract In this work a hybrid velocity-vorticity scheme for the solution of the 2D Navier-Stokes equations is presented. The multi-region Local Boundary Integral Equation (LBIE) combined with Radial Basis Functions (RBF) interpolation is used for the solution of the kinematics and the multi-region BEM for the solution of the transport kinetics. The final system of equations is in band form for both methods. The issue of RBF discontinuities is resolved by constructing the RBF matrix locally in every region. The kinematics integral equation is used in three different forms, for coupling the velocity field on the More >

C. Xu¹, X. Jia², R. A. Williams², E. H. Stitt³, M. Nijemeisland³, S. El-Bachir⁴, A. J. Sederman⁴, L. F. Gladden⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 117-126, 2008, DOI:10.3970/cmes.2008.023.117

Abstract Existing theories and computer models for packed columns are either incapable of handling complex pellet shapes or based on over-simplified packing geometry. A digital packing algorithm, namely DigiPac, has recently been developed to fill the gap. It is capable of packing of particles of any shapes and sizes in a container of arbitrary geometry, and is a first step towards a practical computational tool for reliable predictions of packed column properties based on the actual pellet shapes. DigiPac can operate in two modes: a Monte Carlo mode in which particles undergo directional diffusive motions; and… More >

Nagalinga Rajan, Soumyendu Raha¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 91-116, 2008, DOI:10.3970/cmes.2008.023.091

Abstract The article describes a numerical method for time domain integration of noisy dynamical systems originating from engineering applications. The models are second order stochastic differential equations (SDE). The stochastic process forcing the dynamics is treated mainly as multiplicative noise involving a Wiener Process in the Itô sense. The developed numerical integration method is a drift implicit strong order 2.0 method. The method has user-selectable numerical dissipation properties that can be useful in dealing with both multiplicative noise and stiffness in a computationally efficient way. A generalized analysis of the method including the multiplicative noise is More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.023.075

Abstract In the MLPG_R (Meshless Local Petrove-Galerkin based on Rankine source solution) method, one needs a meshless interpolation scheme for an unknown function to discretise the governing equation. The MLS (moving least square) method has been used for this purpose so far. The MLS method requires inverse of matrix or solution of a linear algebraic system and so is quite time-consuming. In this paper, a new scheme, called simplified finite difference interpolation (SFDI), is devised. This scheme is generally as accurate as the MLS method but does not need matrix inverse and consume less CPU time More >

A.P.S. Selvadurai¹, H. Ghiabi²

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

Abstract This paper deals with the problem of the consolidation of a composite consisting of alternate layers of soft clay and a granular material. A series of experiments were conducted on components to develop the constitutive models that can be implemented in a computational approach. The constitutive response of the soft clay is represented by a poro-elasto-plastic Cam clay-based model and the granular medium by an elasto-plastic model with a Drucker-Prager type failure criterion and a non-associated flow rule. The computational poro-elasto-plastic model is used to calibrate the experimental results derived from the one-dimensional tests and More >

L. Parussini¹, V. Pediroda²

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

Abstract In this paper different Polynomial Chaos methods coupled to Fictitious Domain approach have been applied to one- and two- dimensional elliptic problems with multi uncertain variables in order to compare the accuracy and convergence of the methodologies. Both intrusive and non-intrusive methods have been considered, with particular attention to their employment for quantification of geometric uncertainties. A Fictitious Domain approach with Least-Squares Spectral Element approximation has been employed for the analysis of differential problems with uncertain boundary domains. Its main advantage lies in the fact that only a Cartesian mesh, that represents the enclosure, needs More >