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 >

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 >

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 >

C. Gato and Y. Shie¹

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.3, pp. 269-292, 2008, DOI:10.3970/cmes.2008.033.269

Abstract Numerical simulations of large deformation dynamic fracture in thin shell structures using 3-D meshfree method is presented. Due to the smoothness of the meshfree shape functions, they are well suited to simulate large deformation of thin shell structures while avoiding ill-conditioning as well as stiffening in numerical computations. Dynamic fracture is modeled by simple criterion, i.e. removing connectivity between adjacent nodes once a fracture criterion is met. The main advantage of such 3-D meshfree continuum approach is its simplicity in both formulation and implementation as compared to shell theory approach, or degenerated continuum approach. Moreover, it is believed that the… More >

K. Han ¹, Y. T. Feng ¹, D. R. J. Owen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.2, pp. 87-100, 2007, DOI:10.3970/cmes.2007.018.087

Abstract Numerical procedures are introduced for simulations of irregular particle transport in turbulent flows using the coupled lattice Boltzmann method (LBM) and the discrete element method (DEM). The fluid field is solved by the extended LBM with the incorporation of the Smagorinsky turbulence approach, while particle interaction is modeled by the DEM. The hydrodynamic interactions between fluid and particles are realised through an immersed boundary condition, which gives rise to a coupled solution strategy to model the fluid-particle system under consideration. Main computational aspects comprise the lattice Boltzmann formulation for the solution of fluid flows; the incorporation of the large eddy… More >

Charles E. Anderson, Jr.¹, Matthew S. Burkins², James D. Walker¹, William A. Gooch²

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 91-104, 2005, DOI:10.3970/cmes.2005.008.091

Abstract A modification of Wilkins computational ceramics model is used to simulate experiments of the impact of the APM2 bullet against boron carbide/aluminum targets. Flash radiography provides time-resolved penetration histories. The simulation results are compared to the experimental data; generally, agreement is very good, including capturing dwell and then the onset of penetration. Crater width and debris diameter are also reproduced by the simulations reasonably well. A critical discussion of deficiencies of this computational engineering model is provided. More >

Joseph G. Lawrence, Arunan Nadarajah¹

FDMP-Fluid Dynamics & Materials Processing, Vol.3, No.3, pp. 247-254, 2007, DOI:10.3970/fdmp.2007.003.247

Abstract A simplified 2D pseudo steady state model was developed for an atmospheric chemical vapor deposition (CVD) process on glass. This is used to study the feasibility of converting a continuous coating process to one with discrete glass plates with a gap between them. A preliminary estimate employing mass transfer correlations suggested that there would be significant concentration variations due to the gap between the plates. More detailed studies were done by solving the model numerically employing a finite difference scheme with a vorticity-stream function formulation, and employing the commercial computational fluid dynamics program FIDAP which employs a finite element scheme.… More >

Alamusi¹, Y.L. Liu¹, N. Hu^1,2

CMC-Computers, Materials & Continua, Vol.20, No.2, pp. 101-118, 2010, DOI:10.3970/cmc.2010.020.101

Abstract In this work, we propose a 3 dimensional (3D) numerical model to predict the piezoresistivity behaviors of a nanocomposite material made from an insulating polymer filled by carbon nanotubes (CNTs). This material is very hopeful for its application in highly sensitive strain sensor by measuring its piezoresistivity, i.e., the ratio of resistance change versus applied strain. In this numerical approach, a 3D resistor network model is firstly proposed to predict the electrical conductivity of the nanocomposite with a large amount of randomly dispersed CNTs under the zero strain state. By focusing on the fact that the piezoresistivity of the nanocomposite… More >