A.E. Martínez-Castro¹, I.H. Faris¹, R. Gallego¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.3, pp. 177-206, 2012, DOI:10.3970/cmes.2012.087.177

Abstract In this paper, the identification of hidden defects inside a three-dimen -sional layer is set as an Identification Inverse Problem. This problem is solved by minimizing a cost functional which is linearized with respect to the volume defects, leading to a procedure that requires only computations at the host domain free of defects. The cost functional is stated as the misfit between experimental and computed displacements and spherical and/or ellipsoidal cavities are the defects to locate. The identification of these cavities is based on the measured displacements at a set of points due to time-harmonic point loads at an array… More >

Kentaro Doi^1,2, Ikumi Onishi¹, Satoyuki Kawano^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 113-136, 2011, DOI:10.3970/cmes.2011.077.113

Abstract Hydrogen technologies are currently one of the most actively researched topics. A lot of researches have tied to enhance their energy conversion efficiencies. In the present study, numerical analyses have been carried out focusing on hydrogen-storage carbon materials which are expected to realize high gravimetric and volumetric capacities. In particular, dissociative adsorption processes of H₂ molecules above graphene surfaces have been investigated by ab initio molecular dynamics. The present results indicate that a steric graphene surface plays an important role in enhancing the charge transfer which induces dissociation of H₂ and adsorption of H atoms on the surface. The dissociation… More >

Chein-Shan Liu¹, Chung-Lun Kuo²

CMES-Computer Modeling in Engineering & Sciences, Vol.76, No.2, pp. 109-132, 2011, DOI:10.3970/cmes.2011.076.109

Abstract The Tikhonov method is a famous technique for regularizing ill-posed systems. In this theory a regularization parameter α needs to be determined. Based-on an invariant-manifold defined in the space of (x,t) and from the Tikhonov minimization functional, we can derive an optimal vector driven system of nonlinear ordinary differential equations (ODEs). In the Optimal Vector Driven Algorithm (OVDA), the optimal regularization parameter αk is presented in the iterative solution of x, which means that a dynamical Tikhonov regularization method is involved in the solution of nonlinear ill-posed problem. The OVDA is an extension of the Landweber-Scherzer iterative algorithm. Numerical examples… More >

Cheng-Yu Ku^1,2,3,Weichung Yeih^1,2, Chein-Shan Liu⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.76, No.2, pp. 83-108, 2011, DOI:10.3970/cmes.2011.076.083

Abstract In this paper, a general dynamical method based on the construction of a scalar homotopy function to transform a vector function of Non-Linear Algebraic Equations (NAEs) into a time-dependent scalar function by introducing a fictitious time-like variable is proposed. With the introduction of a transformation matrix, the proposed general dynamical method can be transformed into several dynamical Newton-like methods including the Dynamical Newton Method (DNM), the Dynamical Jacobian-Inverse Free Method (DJIFM), and the Manifold-Based Exponentially Convergent Algorithm (MBECA). From the general dynamical method, we can also derive the conventional Newton method using a certain fictitious time-like function. The formulation presented… More >

Delfim Soares Jr.¹

CMES-Computer Modeling in Engineering & Sciences, Vol.76, No.1, pp. 61-82, 2011, DOI:10.3970/cmes.2011.076.061

Abstract This work proposes an iterative procedure to analyse pore-dynamic models discretized by time-domain Meshless Local Petrov-Galerkin formulations. By considering an iterative procedure based on a successive renew of variables, each phase of the coupled problem in focus can be treated separately, uncoupling the governing equations of the model. Thus, smaller and better conditioned systems of equations are obtained, rendering a more attractive methodology. A relaxation parameter is introduced here in order to improve the efficiency of the iterative procedure and an expression to compute optimal values for the relaxation parameter is discussed. Linear and nonlinear models are focused, highlighting that… More >

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 >

Wei-Xin Ren^1,2,3, Cong-Cong Su¹, Wang-Ji Yan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.70, No.1, pp. 67-92, 2010, DOI:10.3970/cmes.2010.070.067

Abstract A beam-arch segment assembly procedure is presented in this paper for the dynamic modelling and analysis of arch bridges. Such a beam-arch segment assembly is composed of different structural elements of arch bridges such as arch ribs (curved beams), suspenders, girders and floor beams. Based on the energy principle in structural dynamics, the stiffness matrix and mass matrix of such an assembly are formulated. The proposed procedure is then implemented to carry out the free vibration analysis of the Jian concrete filled tubular arch bridge. It is demonstrated that the proposed beam-arch segment assembly procedure is efficient with the advantages… 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 >

Yantao Zhang¹, Xiong Zhang^1,2, Yan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.2, pp. 143-166, 2010, DOI:10.3970/cmes.2010.069.143

Abstract Material point method(MPM) is a promising method in solving problems involving large deformations, especially explosion and penetration. In MPM, particles can move around the computing domain dynamically, which can result in load imbalance easily. In parallelizing MPM using OpenMP, data races will occur in the stage of grid node updating if we use loop-level parallelism for these loops. Huang et al. proposed a domain decomposition method to overcome data races [Huang, Zhang, Ma and Wang (2008)]. However, significant modifications of the original serial code are required. In this paper, we proposed a new alternated grid updating method to avoid data… More >