David Stevens¹, Henry Power^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 111-146, 2010, DOI:10.3970/cmes.2010.055.111

Abstract Problems involving nonlinear time-dependent heat conduction in materials which have temperature-dependent thermal properties are solved with a novel meshless numerical solution technique using multiquadric radial basis functions (RBFs). Unlike traditional RBF collocation methods, the local Hermitian interpolation (LHI) method examined here can be scaled to arbitrarily large problems without numerical ill-conditioning or computational cost issues, due to the presence of small overlapping interpolation systems which grow in number but not in size as the global dataset grows. The flexibility of the full-domain multiquadric collocation method to directly interpolate arbitrary boundary conditions is maintained, via the More >

Ying Liu¹, Yufeng Nie², Weiwei Zhang², Lei Wang²

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.1, pp. 89-110, 2010, DOI:10.3970/cmes.2010.055.089

Abstract In the light of the ideas and treatment technologies about molecular dynamics simulation and bubble meshing, a new approach of node placement for the meshless method called node placement method by bubble simulation (NPBS method), is proposed. Nodes are seen as the centers of the bubbles which can be moved by their interacting forces. Through dynamic simulation, bubbles are placed into a near-optimal configuration, and the centers of bubbles will form a good-quality node distribution in the domain. This process doesn't need updating the mesh connection constantly, i.e., is totally meshfree. Some example results show… More >

Z. T. Ma¹, X. Zhang^1,2, P. Huang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.1, pp. 61-88, 2010, DOI:10.3970/cmes.2010.055.061

Abstract The Material Point Method (MPM) is more expensive in terms of storage than other methods, as MPM makes use of both mesh and particle data. Therefore, it is critical to develop an efficient MPM framework for engineering applications, such as impact and explosive simulations. This paper presents a new architecture for MPM computer code, developed using object-oriented design, which enables MPM analysis of a mass of grains, large deformation, high strain rates and complex material behavior. It is flexible, extendible, and easily modified for a variety of MPM analysis procedures. An MPM scheme combining contact More >

Z. Y. Yan^1,2, J. Zhang¹, W. Ye¹, T.X. Yu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.1, pp. 33-60, 2010, DOI:10.3970/cmes.2010.055.033

Abstract An 3-D precorrected-FFT accelerated BEM approach for the linear elastic analysis of porous solids with randomly distributed pores of arbitrary shape and size is described in this paper. Both the upper bound and the lower bound of elastic properties of solids with spherical pores are obtained using the developed fast BEM code. Effects of porosity and pore shape on the elastic properties are investigated. The performance of several theoretical models is evaluated by comparing the theoretical predictions with the numerical results. It is found that for porous solids with spherical pores, the performances of the More >

Surkay D. Akbarov^1,2,3, Mugan S. Guliev⁴, Ramazan Tekercioglu⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.1, pp. 1-32, 2010, DOI:10.3970/cmes.2010.055.001

Abstract Dispersion relations of axisymmetric longitudinal wave propagation in a finite pre-strained compound (bi-material) cylinder made from high elastic incompressible materials are investigated within the scope of a piecewise homogeneous body model utilizing three-dimensional linearized theory wave propagation in the initially stressed body. The materials of the inner and outer cylinders are assumed to be neo-Hookean. The numerical results regarding the influence of the initial strains in the inner and outer cylinders on the wave dispersion are presented and discussed. These results are obtained for the case where the material of the inner solid cylinder is More >

L. Dammak, E. Hadj-Taïeb

FDMP-Fluid Dynamics & Materials Processing, Vol.6, No.4, pp. 419-430, 2010, DOI:10.3970/fdmp.2010.006.419

Abstract This paper presents a short and focused analysis of the pressure development inside the fluid film related to a journal bearing (i.e. the pressure distribution in the the gap between the shaft, generally referred to as the "journal", and the bearing). The related flow is considered to be isotherm, laminar, steady and incompressible. The lubricant is assumed to be an isoviscous fluid. The Reynolds equation governing the lubricant pressure is derived from the coupled continuity and momentum balance equations written in the framework of the Stokes theory. The non linear system given by coupled equations More >

R. Mestiri¹, R.Hadaji¹, S. Ben Nasrallah¹

FDMP-Fluid Dynamics & Materials Processing, Vol.6, No.4, pp. 409-418, 2010, DOI:10.3970/fdmp.2010.006.409

Abstract The electroaerodynamic actuator or plasma actuator uses the characteristics of the non-thermal surface plasmas. These plasmas are created in atmospheric pressure by a DC electrical corona discharge at the surface of a dielectric material. The two electrodes are two conductive parallel wires. The applied voltage is of several kilovolts. The corona discharge creates a tangential electric wind that can modify the boundary layer flow properties. In this paper, we present the results found for two geometric configurations: the flat plate and the cylinder. In order to study the discharge specificity, we have found the current- More >

M. Chrigui^1,2, L. Schneider¹, A. Zghal², A. Sadiki¹, J. Janicka¹

FDMP-Fluid Dynamics & Materials Processing, Vol.6, No.4, pp. 399-408, 2010, DOI:10.3970/fdmp.2010.006.399

Abstract This paper deals with the numerical simulation of droplet dispersion and evaporation within an LPP (Lean Premix Prevaporized) burner. The Eulerian-Lagrangian approach was used for this purpose, and a fully two way-coupling was accounted for. For the phase transition, a non-equilibrium evaporation model was applied that differs strongly from the equilibrium one where there are high evaporation rates. The non-equilibrium conditions were fulfilled in the investigated configuration, as the droplets at the inlet had a mean diameter of 50mm. The numerical results of water droplet velocities, corresponding fluctuations, and diameters were compared with experimental data. More >

A. Radhouane¹, N. Mahjoub Said¹, H. Mhiri¹, G. Le Palec², P. Bournot²

FDMP-Fluid Dynamics & Materials Processing, Vol.6, No.4, pp. 385-398, 2010, DOI:10.3970/fdmp.2010.006.385

Abstract "Twin jets in crossflow" is a complex configuration that raises an increasing interest due to its presence in various common applications such as chimney stacks, film cooling, VSTOL aircrafts, etc... In the present paper, the twin jets were arranged inline with an oncoming crossflow;they were also inclined which resulted in similar elliptic cross sections of the nozzles' exits. The exploration of the flows in interaction was carried out numerically by means of the finite volume method together with the second order turbulent closure model, namely the Reynolds stress Model (RSM), and a non uniform grid… More >

S. Hamimid¹, A.Amroune¹

FDMP-Fluid Dynamics & Materials Processing, Vol.6, No.4, pp. 369-384, 2010, DOI:10.3970/fdmp.2010.006.369

Abstract The Navier-Stokes and energy equations are numerically solved to investigate two-dimensional convection (originating from the combined effect of buoyancy and surface tension forces) in a liquid metal subjected to transverse magnetic fields. In particular, a laterally heated horizontal cavity with aspect ratio (height/width) =1 and Pr=0.015 is considered (typically associated with the horizontal Bridgman crystal growth process and commonly used for benchmarking purposes). The effect of a uniform magnetic field with different magnitudes and orientations on the stability of the two distinct convective solution branches (with a single-cell or two-cell pattern) of the steady-state flows More >