A. Colombo¹, U. Galvanetto²

CMES-Computer Modeling in Engineering & Sciences, Vol.53, No.3, pp. 235-254, 2009, DOI:10.3970/cmes.2009.053.235

Abstract The paper addresses the problem of computing the stable manifolds of equilibria and limit cycles of saddle type in piecewise smooth dynamical systems. All singular points that are generically present along one-dimensional or two-dimensional manifolds are classified and such a classification is then used to define a method for the numerical computation of the stable manifolds. Finally the proposed method is applied to the case of a stick-slip oscillator. More >

Yi W. Wan¹, Huan J. Keh²

CMES-Computer Modeling in Engineering & Sciences, Vol.53, No.1, pp. 73-94, 2009, DOI:10.3970/cmes.2009.053.073

Abstract The problem of the rotation of a rigid particle of revolution about its axis in a viscous fluid is studied theoretically in the steady limit of low Reynolds number. The fluid is allowed to slip at the surface of the particle. A singularity method based on the principle of distribution of a set of spherical singularities along the axis of revolution within a prolate particle or on the fundamental plane within an oblate particle is used to find the general solution for the fluid velocity field that satisfies the boundary condition at infinity. The slip… More >

Yu Miao¹, Qiao Wang¹, Bihai Liao^1,2, Junjie Zheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.53, No.1, pp. 1-22, 2009, DOI:10.3970/cmes.2009.053.001

Abstract As a truly meshless method, the Hybrid Boundary Node method (Hybrid BNM) does not require a `boundary element mesh', either for the purpose of interpolation of the solution variables or for the integration of `energy'. This paper presents a further development of the Hybrid BNM to the 2D elastodynamics. Based on the radial basis function (RBF) and the Hybrid BNM, it presents an inherently meshless, boundary-only technique, which named dual hybrid boundary node method (DHBNM), for solving 2D elastodynamics. In this study, the RBFs are employed to approximate the inhomogeneous terms via dual reciprocity method More >

A.R. da Silva¹, A. Silveira-Neto^2,3, D.A. Rade^2,4, R.Francis⁴, E.A. Santos⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.3, pp. 285-304, 2009, DOI:10.3970/cmes.2009.050.285

Abstract In the context of computational fluid dynamics a numerical investigation of incompressible flow around fixed pairs of rigid circular cylinders was carried out. The two-dimensional filtered Navier-Stokes equations with the Smagorinsky sub-grid scale model were solved using a Cartesian non-uniform grid. The immersed Boundary Method with the Virtual Physical Model was used in order to model the presence of two circular cylinders embedded in the flow. The fractional time step method was used to couple pressure and velocity fields. The simulations were carried out for Reynolds number equal to 72,000 for pitch ratio equal to More >

Chun-Lang Yeh ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.3, pp. 253-284, 2009, DOI:10.3970/cmes.2009.050.253

Abstract This paper investigates the instability of a liquid thread by molecular dynamics (MD) simulation. The influences of liquid thread radius, fundamental cell length, and temperature are discussed. Snapshots of molecules, number of liquid particles formed, and density field are analyzed. Two linear stability criteria, namely Rayleigh's stability criterion and Kim's stability criterion, are accessed for their validity in molecular scale. It is found that a liquid thread is more unstable and produces more liquid particles in the fundamental cell when it is thinner or at a higher temperature. In addition, a liquid thread with a More >

D. Soares Jr.¹, J. Sladek², V. Sladek²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 115-140, 2009, DOI:10.3970/cmes.2009.050.115

Abstract In recent years the idea of using mesh-free or mesh-less methods for numerical solution of partial differential equations has received much attention, due to their potential advantage in eliminating the costly effort of mesh generation and re-meshing. A variety of meshless methods has been proposed so far. Many of them are derived from a weak-form formulation on global domain or a set of local subdomains. In the global formulation background cells are required for the integration of the weak form. In methods based on local weak-form formulation no cells are required and therefore they are… More >

Tai-Ming Chang¹, Chien-Chou Weng¹, Mei-Jiau Huang^1,2, Chun-KaiLiu², Chih-Kuang Yu²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.1, pp. 47-66, 2009, DOI:10.3970/cmes.2009.050.047

Abstract We employ the non-equilibrium molecular dynamics (NEMD) simulation to calculate the in-plane thermal conductivity of silicon thin films of thickness 2.2nm and 11nm. To eliminate the finite-size effect, samples of various lengths are simulated and an extrapolation technique is applied. To perform the quantum correction which is necessary as the MD simulation temperature is lower than Debye temperature, the confined phonon spectra are obtained in advance via the EMD simulations. The investigation shows the thermal conductivities corrected based on the bulk and thin-film phonon densities of states are very close and they agree excellently with More >

Z.H.Yao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.1, pp. 21-46, 2009, DOI:10.3970/cmes.2009.050.021

Abstract The traditional time domain boundary integral equation (TDBIE) of elastodynamics is formulated based on the time dependent fundamental solution and the reciprocal theorem of elastodynamics. The time dependent fundamental solution of the elastodynamics is the response of the infinite elastic medium under a unit concentrate impulsive force subjected at a point and at an instant, including not only the pressure wave and shear wave, but also the Laplace wave with speed between that of P and S waves. In this paper, a new TDBIE is derived directly from the initial boundary value problem of the… More >

Haiping Ye¹, Yongsheng Ding²

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.3, pp. 255-268, 2009, DOI:10.3970/cmes.2009.049.255

Abstract A fractional-order model for the immunological and therapeutic control of HIV is studied qualitatively. The equilibria are found and their local stability are investigated. Also the global stability of the infection-free equilibrium is established. The optimal efficacy level of anti-retroviral therapy needed to eradicate HIV from the body of an HIV-infected individual is obtained. More >

J.M.C.Pereira¹, N.A.R.Maia¹, J.C.F.Pereira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.2, pp. 113-142, 2009, DOI:10.3970/cmes.2009.049.113

Abstract We present a 2D incompressible Navier-Stokes numerical simulation of a virtual model of an elliptic, or flat plate, foil in hovering flight configuration. Computations obtained with a general purpose solver were validated against reference data on forward flapping flight, normal or dragonfly hovering. The moving mesh technique allows airfoil translation and angular mesh movement accompaining the airfoil stroke motion. Close to the ground the mesh deforms to occupy the narrow computational domain formed between the airfoil and the ground. Computations have been carried out for some parameters, including the distances h between the foil center and… More >