Kazuhiko Suga^1,2, Takahiko Ito¹

CMES-Computer Modeling in Engineering & Sciences, Vol.75, No.2, pp. 141-172, 2011, DOI:10.3970/cmes.2011.075.141

Abstract The multiple-relaxation-time lattice Boltzmann method for micro-scale flows (MRT µ-flow LBM) is extensively evaluated in this study. Following the study of Chai, Shi, Guo and Lu (2010), the diffusive bounce-back wall boundary condition and the collision matrix are modeled. To determine the model parameters, the first-order, 1.5-order and second-order slip-flow models are discussed. Since the mean free path of gas molecules is considered to be influenced by the wall in micro flow systems, the effects of a correction function after Stops (1970) are also evaluated. As the increase of the Knudsen number (Kn), it is necessary to introduce the regularization… More >

L.Q. Moreira¹, F.P. Mariano², A. Silveira-Neto¹

CMES-Computer Modeling in Engineering & Sciences, Vol.75, No.2, pp. 113-140, 2011, DOI:10.3970/cmes.2011.075.113

Abstract Turbulence in fluid flow is one of the most challenging problems in classical physics. It is a very important research problem because of its numerous implications, such as industrial applications that involve processes using mixtures of components, heat transfer and lubrication and injection of fuel into the combustion chambers and propulsion systems of airplanes. Turbulence in flow presents characteristics that are fully nonlinear and that occur at high Reynolds numbers. Because of the nonlinear nature of turbulent flow, an increase in the Reynolds number implies an increase in the Kolmogorov wave numbers, and the flow spectrum becomes larger in both… More >

V.V. Zozulya¹

CMES-Computer Modeling in Engineering & Sciences, Vol.70, No.3, pp. 253-350, 2010, DOI:10.3970/cmes.2010.070.253

Abstract In this article the divergent integrals, which arise when the boundary integral equation (BIE) methods are used for solution of the 3-D elastostatic problems is considered. The same approach for weakly singular, singular and hypersingular integral regularization is developed. The approach is based on theory of distribution and Green's theorems. This approach is applied for regularization of the divergent integrals over convex polygonal boundary elements (BE) in the case of piecewise constant approximation and over rectangular and triangular BE for piecewise linear approximation. The divergent integrals are transformed into the regular contour integrals that can be easily calculated analytically. Proposed… More >

Xianqiao Wang¹, James D. Lee²

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.3, pp. 199-222, 2010, DOI:10.3970/cmes.2010.069.199

Abstract This paper presents an atom-based continuum (ABC) method aiming at a seamless transition from the atomistic to the continuum description of multi-element crystalline solids (which has more than one kind of atom in the unit cell). Contrary to many concurrent multiscale approaches, ABC method is naturally suitable for the analysis of multi-element crystals within a finite element (FE) framework. Taking both efficiency and accuracy into account, we adopt a cluster-based summation rule for atomic force calculations in the FE formulations. Single-crystals MgO, BaTiO₃ and Cu under mechanical loading are modeled and simulated. With a coarse-grained mesh, ABC method is shown… More >

A.P.S.Selvadurai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.3, pp. 259-278, 2009, DOI:10.3970/cmes.2009.052.259

Abstract The paper deals with a computational approach for modelling the fragmentation of ice sheets during their impact with stationary structures. The modelling takes into consideration the intact continuum behaviour of the ice as a rate-sensitive elastoplastic material. During impact, the ice sheet can undergo fragmentation, which is controlled by a brittle strength criterion based on the current stress state. The fragmentation allows the generation of discrete elements of the ice sheet, the movements of which are governed by the equations of motion. The contact between individual fragments is governed by a Coulomb criterion. The individual fragments can themselves undergo further… More >

A. Ferrari¹, L. Laloui^1,2, Ch. Bonnard^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.3, pp. 217-236, 2009, DOI:10.3970/cmes.2009.052.217

Abstract A coupled hydro-mechanical formulation is presented for the analysis of landslide motion during crisis episodes. The mathematical formulation is used to model a natural slope affected by a multiple slip surface failure mechanism, in which pore water pressure evolution was identified as the main cause for movement accelerations. An elasto-plastic constitutive model is adopted for the behaviour of slip surfaces. Material parameters are obtained by combining the available laboratory tests and the back analysis of some crisis episodes. After being calibrated and validated, the model is applied to improve the understanding of the physical processes involved and to predict the… More >

Anil Misra¹, Shiping Huang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.2, pp. 197-216, 2009, DOI:10.3970/cmes.2009.052.197

Abstract The behavior of contact between solid bodies with rough surfaces under combined normal and shear loading remains a problem of interest in many areas of engineering. In this paper, we have utilized a micromechanical methodology to derive an expression of stress-displacement relationship applicable to combined normal and shear loading conditions. The micromechanical methodology considers the mechanics of asperity contacts and the interface roughness in terms of asperity height and asperity contact orientation distribution. A numerical procedure is implemented to evaluate the derived expressions under complex and mixed loading conditions using an incremental approach. We find that the proposed numerical procedure… More >

E. Mesquita¹, J.Labaki ¹ and L.O.S.Ferreira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.2, pp. 143-168, 2009, DOI:10.3970/cmes.2009.051.143

Abstract There is a growing trend towards solving problems of computational mechanics by parallelization strategies. The traditional approach is to implement the parallelization procedures on CPUs based on the MPI or OpenMP paradigms. Recent efforts have been made to implement computational tasks on general-purpose programmable graphics hardware (GPGPU). The GPU is specially well-suited to address problems that can be formulated in form of data-parallel computations with high arithmetic intensity. This work addresses the implementation of the Longman's integration method on graphics hardware. A serial implementation of Longman's method was rewritten under the SIMD (Single Input Multiple Data) parallel programming paradigm. The… More >

G. Veble^1,2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.3, pp. 243-260, 2008, DOI:10.3970/cmes.2008.036.243

Abstract A simple cost function is proposed that depends on the inviscid pressure distribution around an airfoil and that, when minimized, results in airfoils that promote laminar flow. Additional constraints specify the design point of the airfoil. The method allows for straightforward inclusion of multiple design points. The resulting airfoils are quantitatively similar to those already successfully used in practice. More >

G.J.A. Loeven^1,2, H. Bijl³

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.3, pp. 193-212, 2008, DOI:10.3970/cmes.2008.036.193

Abstract In this paper a Two-Step approach is presented for uncertainty quantification for expensive problems with multiple uncertain parameters. Both steps are performed using the Probabilistic Collocation method. The first step consists of a sensitivity analysis to identify the most important parameters of the problem. The sensitivity derivatives are obtained using a first or second order Probabilistic Collocation approximation. For the most important parameters the probability distribution functions are propagated using the Probabilistic Collocation method using higher order approximations. The Two-Step approach is demonstrated for flow around a NACA0012 airfoil with eight uncertain parameters in the free stream conditions and geometry.… More >