H.K. Mebatsion^1,2, P. Verboven¹, P. T. Jancsók¹, Q.T. Ho¹, B.E. Verlinden³, B.M. Nicolaï^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.3, pp. 137-150, 2008, DOI:10.3970/cmes.2008.029.137

Abstract Transport processes of gas and moisture are among the most important physiological processes in plant tissue. Microscale transport models based on Navier-Stokes equations provide insight into such processes at the microscopic scale. Due to microscopic complexity, numerical solutions based on the finite element or finite volume methods are mandatory. Therefore, a 3D geometric model of the tissue is essential. In this article, a novel algorithm for geometric reconstruction of 2D slices of synchrotron tomographic images is presented. The boundaries of 2D cells on individual slices were digitized to establish a set of boundary coordinates and the slice index of individual… More >

H.M. Lee¹, H.S. Park^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 29-44, 2008, DOI:10.3970/cmes.2008.029.029

Abstract This paper presents a computational model to estimate deformed shapes of beam structures using 3D coordinate information from terrestrial laser scanning (TLS). The model is composed of five components: 1) formulation of polynomial shape function, 2) application of boundary condition, 3) inducement of compatibility condition, 4) application of the least square method and 5) evaluation of error vector and determination of reasonable polynomial shape function. In the proposed model, the optimal degree of polynomial function is selected based on the complexity of beam structures, instead of using a specific degree of polynomial function. The chosen polynomial function for estimation is… More >

G.F. Karlis¹, S.V. Tsinopoulos², D. Polyzos³, D.E. Beskos⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.3, pp. 189-208, 2008, DOI:10.3970/cmes.2008.026.189

Abstract A boundary element method, suitable for solving two and three dimensional gradient elastic fracture mechanics problems under static loading, is presented. A simple gradient elastic theory (a simplied version of Mindlin's Form-II general theory of gradient elasticity) is employed and the static gradient elastic fundamental solution is used to construct the boundary integral representation of the problem with the aid of a reciprocal integral identity. In addition to a boundary integral representation for the displacement, a boundary integral representation for its normal derivative is also necessary for the complete formulation of a well-posed problem. Surface quadratic line and quadrilateral boundary… More >

Y.Y. Shan¹, C. Shu^1,2, Z.L. Lu³

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 99-114, 2008, DOI:10.3970/cmes.2008.025.099

Abstract The local multiquadric-based differential quadrature (MQ-DQ) method proposed by [Shu, Ding, and Yeo (2003)] is a natural mesh-free approach for derivative approximation, which is easy to be implemented to solve problems with curved boundary. Previously, it has been well tested for the two-dimensional (2D) case. In this work, this mesh-free method was extended to simulate fluid flow problems with curved boundary in three-dimensional (3D) space. The main concern of this work is to numerically study the performance of the 3D local MQ-DQ method and demonstrate its capability and flexibility for simulation of 3D incompressible fluid flows with curved boundary. Fractional… More >

F.C. Araújo¹, L.J. Gray²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 103-122, 2008, DOI:10.3970/cmes.2008.024.103

Abstract In recent years, carbon nanotubes (CNTs) have been widely employed to build advanced composites. In this work, a Boundary Element Method (BEM) is applied to 3D representative volume elements (RVEs) to estimate mechanical properties of CNT-based composites. To model the thin-walled nanotubes, special integration procedures for calculating nearly-strongly-singular integrals have been developed. The generic BE substructuring algorithm allows modeling complex CNT-reinforced polymers, containing any number of nanotubes of any shape (straight or curved). The subregion-by-subregion strategy, based on Krylov solvers, makes the independent generation, assembly, and storage of the many parts of the complete BE model possible. Thus, significant memory… More >

E. Y. K. Ng¹, S. T. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.1, pp. 13-28, 2008, DOI:10.3970/cmes.2008.023.013

Abstract Microchannels have been recognized as a very effective chemical separation and heat transfer device. The electrical double layer (EDL) effect in a micro-scale flow is however anticipated to be critical. In this paper, Nernst-Planck model (NPM), is used to predict the ion concentration distribution as it is reported to be a more appropriate model for developing microchannel flow. The governing equations are discretised for developing rectangular microchannel flows in Cartesians coordinate. An additional body force source term that is relating to the electric potential, resulted from the EDL effect is introduced in the conventional z-axis momentum equation as a body… More >

Chun Yang¹, Dalin Tang², Chun Yuan³, Thomas S. Hatsukami⁴, Jie Zheng⁵, Pamela K. Woodard⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 233-246, 2007, DOI:10.3970/cmes.2007.019.233

Abstract It has been recognized that fluid-structure interactions (FSI) play an important role in cardiovascular disease initiation and development. However, in vivo MRI multi-component FSI models for human carotid atherosclerotic plaques with bifurcation and quantitative comparisons of FSI models with fluid-only or structure-only models are currently lacking in the literature. A 3D non-Newtonian multi-component FSI model based on in vivo/ex vivo MRI images for human atherosclerotic plaques was introduced to investigate flow and plaque stress/strain behaviors which may be related to plaque progression and rupture. Both artery wall and plaque components were assumed to be hyperelastic, isotropic, incompressible and homogeneous. Blood… More >

Chun-Ying. Wu¹, R. Al-Khoury², C. Kasbergen², Xue-Yan. Liu², A. Scarpas²

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.3, pp. 163-172, 2007, DOI:10.3970/cmes.2007.017.163

Abstract 3D spectral element method in the article is presented to predict the pavement layer modules using field measurement of Falling Weight Deflectometer (FWD). To improve the computational efficiency of the layer-condition assessment, one type of spectral element is proposed to develop the dynamic analysis of 3D multi-layered system subjected to an impulsive load. Each layer in structure is simulated as two-noded layer spectral element or one-noded spectral element in frequency domain. In order to identify the parameters of layered structures, a nonlinear optimization method called Powell hybrid algorithm is employed. The optimization procedure is performed in frequency domain and aims… More >

G.P. Nikishkov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 81-90, 2006, DOI:10.3970/cmes.2006.011.081

Abstract This paper presents the object oriented approach to programming the finite element method using the Java language. The developed finite element code consists of ten Java packages. Three main methods are related to generation of finite element models, solution of elastic and elastic-plastic boundary value problems, and visualization of models and solution results. Object-oriented model of the code is described. It is demonstrated that Java 1.5 new features are useful in development of the finite element code. Java 3D is used for visualization of models and results. More >

V.V. Mykhas’kiv¹, O.I. Stepanyuk²

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 45-64, 2005, DOI:10.3970/cmes.2005.010.045

Abstract The 3D elastostatic problem for an infinite remotely loaded matrix containing a finite number of arbitrarily located rigid disk-inclusions and plane cracks is solved by the boundary integral equation (BIE) method. Its boundary integral formulation is achieved by the superposition principle with the subsequent integral representations of superposition terms through surface integrals, which should satisfy the displacement linearity conditions in the inclusion domains and load-free conditions in the crack domains. The subtraction technique in the conjunction with mapping technique under taking into account the structure of the solution at the edges of inhomogeneities is applied for the regularization of BIE… More >