Ting Du¹, J. P. Xu¹, Y.Q. Liu², Z. B. Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.2, pp. 93-114, 2009, DOI:10.3970/cmes.2009.051.093

Abstract An eight-node hexahedral element with four-point quadrature in plane is developed using the assumed strain method, which can eliminate volumetric locking of incompressible material and absence of the portion of shear velocity strain related with hourglass mode to suppress hourglass mode and shear locking. In this approach, the radial return algorithm is adopted for more precise calculation of internal forces, stress and strain. In addition, a co-rotational coordinates system is established to make bending simulation much more effective, and the system is applicable to arbitrary 3D anisotropic yield criteria. A large elastic-plastic deformation of unconstrained thick plate bending example is… More >

Rubén Avila¹, Eduardo Ramos², S. N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.1, pp. 73-92, 2009, DOI:10.3970/cmes.2009.051.073

Abstract A Chebyshev Tau numerical algorithm is presented to solve the perturbation equations that result from the linear stability analysis of the convective motion of a fluid layer that appears when an unconfined solid melts in the presence of gravity. The system of equations that describe the phenomenon constitute an eigenvalue problem whose accurate solution requires a robust method. We solve the equations with our method and briefly describe examples of the results. In the limit where the liquid-solid interface recedes at zero velocity the Rayleigh-Bénard solution is recovered. We show that the critical Rayleigh number Ra_c and the critical wave… More >

G.S. Li¹, D. Yao², Y.Z. Wang³, H.Y. Jiang²

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

Abstract In this paper, an undisturbed soil-column infiltrating experiment is investigated, and a mathematical model describing multi-components solutes transport behaviors in the column is put forward by combing hydro-chemical analysis with advection dispersion mechanisms, which is a group of advection-dispersion-reaction partial differential equations. Since the model involving six reaction coefficients which can not be obtained directly, an optimal perturbation regularization algorithm of determining these parameters is performed, and numerical simulations under different conditions are carried out. Furthermore, the inversion algorithm is applied to solve the real inverse problem by utilizing the measured breakthrough data. The reconstruction data basically coincide with the… More >

Tae-Hyuk Ahn¹, Layne T. Watson^1,2, Yang Cao^1,1, Clifford A. Shaffer¹, William T. Baumann³

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.1, pp. 27-52, 2009, DOI:10.3970/cmes.2009.051.027

Abstract For biochemical systems, where some chemical species are represented by small numbers of molecules, discrete and stochastic approaches are more appropriate than continuous and deterministic approaches. The continuous deterministic approach using ordinary differential equations is adequate for understanding the average behavior of cells, while the discrete stochastic approach accurately captures noisy events in the growth-division cycle. Since the emergence of the stochastic simulation algorithm (SSA) by Gillespie, alternative algorithms have been developed whose goal is to improve the computational efficiency of the SSA. This paper explains and empirically compares the performance of some of these SSA alternatives on a realistic… More >

R. Zakhama^1,2,3, M.M. Abdalla², H. Smaoui^1,3, Z. Gürdal²

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.1, pp. 1-26, 2009, DOI:10.3970/cmes.2009.051.001

Abstract A multigrid accelerated cellular automata algorithm for two and three dimensional continuum topology optimization problems is presented. The topology optimization problem is regularized using the traditional SIMP approach. The analysis rules are derived from the principle of minimum total potential energy, and the design rules are derived based on continuous optimality criteria interpreted as local Kuhn-Tucker conditions. Three versions of the algorithm are implemented; a cellular automata based design algorithm, a baseline multigrid algorithm for analysis acceleration and a full multigrid integrated analysis and design algorithm. It is shown that the multigrid accelerated cellular automata scheme is a powerful tool… 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 2 and different arrangements regarding… 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 longer fundamental cell length is… More >

Zone-Ching Lin¹, Jia-Rong Ye²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.3, pp. 227-252, 2009, DOI:10.3970/cmes.2009.050.227

Abstract The paper develops a quasi-steady molecular statics model to analyze nanoscale cutting of copper materials by diamond cutters with different shapes. Cutting action, cutting force, equivalent strain and equivalent stress are discussed and compared. The quasi-steady molecular statics nanocutting model first assumes the trajectory of each atom of the copper workpiece being cut whenever the diamond cutter goes forward one step. It then uses the optimization search method to solve the force equilibrium equation of the Morse force in the X and Y directions when each atom moves a small distance, so as to find the new movement position of… More >

M.J. Galán Moreno, J.J. Sánchez Medina, L. Álvarez Álvarez E. Rubio Royo¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.3, pp. 191-226, 2009, DOI:10.3970/cmes.2009.050.191

Abstract Urban traffic is a key factor for the development of a city. There exist many different approaches facing traffic optimization. In our case we have focused on traffic lights optimization. We have designed and tested a new architecture to optimize traffic light cycle times. The purpose of this research is to demonstrate the good performance of our architecture in a congested scenario. We have simulated several congestion situations for a very large real world traffic network - "La Almozara" in Zaragoza, Spain. Our results seem encouraging in this extreme situation. As we increase the load in the network we get… More >

Nicolas Ali Libre^1,2, Arezoo Emdadi², Edward J. Kansa^3,4, Mohammad Shekarchi², Mohammad Rahimian²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 161-190, 2009, DOI:10.3970/cmes.2009.050.161

Abstract We present a wavelet based adaptive scheme and investigate the efficiency of this scheme for solving nearly singular potential PDEs over irregularly shaped domains. For a problem defined over Ω∈ℜ^d, the boundary of an irregularly shaped domain, Γ, is defined as a boundary curve that is a product of a Heaviside function along the normal direction and a piecewise continuous tangential curve. The link between the original wavelet based adaptive method presented in Libre, Emdadi, Kansa, Shekarchi, and Rahimian (2008, 2009) or LEKSR method and the generalized one is given through the use of simple Heaviside masking procedure. In addition… More >