Qiang Li¹, Jie Ouyang¹, Guorong Wu¹, Xiaoyang Xu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.82, No.3&4, pp. 215-232, 2011, DOI:10.32604/cmes.2011.082.215

Abstract The governing equations of two-phase flows including gas and polymer melt are presented, which are solved using finite volume and domain extension methods with SIMPLEC technology. The melt filling and primary gas penetration in gas-assisted injection molding (GAIM) process are simulated, where the Cross-viscosity model is employed to describe the melt rheological behavior, and the CLSVOF(coupled Level Set and Volume of fluid) method is employed to capture the moving interfaces. In order to test and verify the coupled methods, melt filling in a rectangular plate with an insert is simulated, and the numerical results are in good agreement with those… More >

M.A. Argenta¹, A.P. Gebert², E.S. Filho³, B.A. Felizari⁴, M.B. Hecke⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.79, No.3&4, pp. 159-182, 2011, DOI:10.3970/cmes.2011.079.159

Abstract Various methods in the literature proposesequations to calculate the stiffness as a function of density of bone tissue such as apparent density and ash density among others[Helgason, Perilli, Schileo, Taddei, Brynjolfsson and Viceconti, 2008]. Other ones present a value of an equivalent elasticity modulus, obtained by statistical adjustments of curves generated through mechanical compression tests over various specimens[Chevalier, Pahr, Allmer, Charlebois and Zysset, 2007; Cuppone, Seedhom, Berry and Ostell, 2004]. Bone tissue is a material withdifferent behaviors according to the scale of observation. It has a complex composite hierarchical structure, which is responsible for assign optimal mechanical properties. Its characteristics,… More >

S. Abbasbandy^1,2, V. Sladek³, A. Shirzadi¹, J. Sladek³

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.1, pp. 15-38, 2011, DOI:10.3970/cmes.2011.071.015

Abstract This paper deals with the development of a new method for solution of initial-boundary value problems governed by a couple of nonlinear diffusion equations occurring in the theory of self-organization in non-equilibrium systems. The time dependence is treated by linear interpolation using the finite difference method and the semi-discrete partial differential equations are considered in a weak sense by using the local integral equation method with approximating 2-d spatial variations of the field variables by the Moving Least Squares. The evaluation techniques are discussed and the applicability of the presented method is demonstrated on two illustrative examples with exact solutions… More >

D. Kakogiannis¹, D. Van Hemelrijck¹, J. Wastiels¹, S. Palanivelu², W. Van Paepegem², J. Vantomme³, A. Kotzakolios⁴, V. Kostopoulos⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.1, pp. 75-94, 2010, DOI:10.3970/cmes.2010.065.075

Abstract In the present study the estimation of the blast wave by two types of finite element methods is investigated: Eulerian multi-material modeling and pure Lagrangian. The main goal is to compare and study their ability to predict the clearing effect during blast. Element shape and improvements on the codes are also considered. For the Lagrangian finite element models the load is applied by using an empirical method, deriving from databases, for the time-spatial distribution of the pressure profiles. In the ideal case of the above method the blast load is applied as an equivalent triangular pulse to represent the decay… More >

Y. Inoue¹, R. Kobayashi¹, S. Ogata¹, T. Gotoh¹

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 137-162, 2010, DOI:10.3970/cmes.2010.063.137

Abstract Vibration of a single graphene and a pair of graphenes at micro meter scale induced by air flow is numerically simulated and examined by using a hybrid computational method starting from a microscopic level of description for the graphene. In order to bridge a huge gap in spatial and time scales in their motions, the carbon atoms of the graphene are represented by a small number of coarse grained particles, the fluid motion is described by the lattice Boltzmann equation and the momentum exchange at the boundary is treated by the time averaged immersed boundary method. It is found that… More >

F.P. Mariano¹, L.Q. Moreira¹, A. Silveira-Neto¹, C.B. da Silva², J.C.F. Pereira²

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.2, pp. 181-216, 2010, DOI:10.3970/cmes.2010.059.181

Abstract A new numerical methodology combining Fourier pseudo-spectral and immersed boundary methods - IMERSPEC - is developed for fluid flow problems governed by the incompressible Navier-Stokes equations. The numerical algorithm consists in a classical Fourier pseudo-spectral methodology using the collocation method where wall boundary conditions are modelled by using an immersed boundary method (IBM). The performance of that new methodology is exemplified in two-dimensional numerical simulations of Green-Taylor decaying vortex, lid-driven cavity and flow over a square cylinder. The convergence rate, the accuracy, the influence of the Reynolds number and the external domain size are analyzed. This new method combines some… More >

Y. T. Gu¹, P. Zhuang², F. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.3, pp. 303-334, 2010, DOI:10.3970/cmes.2010.056.303

Abstract Recently, the numerical modelling and simulation for anomalous subdiffusion equation (ASDE), which is a type of fractional partial differential equation(FPDE) and has been found with widely applications in modern engineering and sciences, are attracting more and more attentions. The current dominant numerical method for modelling ASDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings. This paper aims to develop an implicit meshless approach based on the radial basis functions (RBF) for numerical simulation of the non-linear ASDE. The discrete system of equations is obtained by using the meshless shape functions… 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 >

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 >

Q.W. Ma^1,2, J.T. Zhou¹

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 277-304, 2009, DOI:10.3970/cmes.2009.043.277

Abstract Following our previous work, the Meshless Local Petrov-Galerin me -thod based on Rankine source solution (MLPG_R) will be extended in this paper to deal with breaking waves. For this purpose, the governing equation for pressure is improved and a new technique called Mixed Particle Number Density and Auxiliary Function Method (MPAM) is suggested for identifying the free surface particles. Due to complexity of the problem, two dimensional (2D) breaking waves are only concerned here. Various cases are investigated and some numerical results are compared with experimental data available in literature to show the newly developed method is robust. More >