V.G. Yakhno¹, H.Çerdik Yaslan²

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.3&4, pp. 295-310, 2011, DOI:10.3970/cmes.2011.081.295

Abstract The time-dependent differential equations of elasticity for 3D quasicrystals are considered in the paper. These equations are written in the form of a vector partial differential equation of the second order with symmetric matrix coefficients. The Green's function is defined for this vector partial differential equation. A new method of the numerical computation of values of the Green's function is proposed. This method is based on the Fourier transformation and some matrix computations. Computational experiments confirm the robustness of our method for the computation of the time-dependent Green's function in icosahedral quasicrystals. More >

C.B. Du^1,2, S.Y Jiang², W. Qin³, Y.M. Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.3&4, pp. 229-248, 2011, DOI:10.3970/cmes.2011.081.229

Abstract A numerical analysis of concrete composites at the mesoscale based on three-dimensional (3D) reconstruction technology of X-ray computed tomography (CT) images is presented in this paper. For X-ray CT images of concrete, morphology processing was used to recover complete image information, including borders, and the median filtering method was applied to eliminate potential impurities in the images. The final X-ray CT images obtained after processing for a concrete section were composed of three-value pixels that indicated aggregate particles, mortar matrix and air voids, and the 3D structure of the concrete specimen was reconstructed using the volume data method. The mapping… More >

Hui-Fen Guo^1,2, Bin-Gang Xu^1,3, Sheng-Yan Li¹, Chong-Wen Yu²

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.2, pp. 87-112, 2011, DOI:10.3970/cmes.2011.080.087

Abstract Three-dimensional transient simulation is presented for swirling recirculating flow in a nozzle with a slotted-tube (different grooves) and the effect of the groove number is also investigated. The numerical results on the streamline angles are validated by experimental visualization using the surface oil flow technology. In the downstream center of the injectors, the vortex breakdown experiences a transition from bubble- to spiral- breakdown as time is increased. For all cases under study, as the sizes of two recirculation zones near the injector upstream wall and the step retain almost constant, the spiral breakdown shows a periodic development. The more the… More >

Kun Luan¹, Baozhong Sun¹, Bohong Gu^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.79, No.1, pp. 31-62, 2011, DOI:10.3970/cmes.2011.079.031

Abstract This paper reports finite element multi-scale simulations of ballistic impact damage of three-dimensional angle-interlock woven composite (3DAWC) penetrated under a hemispherical rigid projectile. A multi-scale geometrical model of the 3DAWC was established for the numerical simulation. The multi-scale geometrical model of the 3DAWC consists two parts: one is the microstructure model and another is the continuum model. The microstructure model has the same microstructure with that of the 3DAWC composite panel, including the fiber tows' diameter, fiber tow configuration and fiber volume fraction. The continuum model has the same mechanical properties with the 3DAWC. The commercial-available finite element software package… More >

Y.C. Shiah¹, C. L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 95-108, 2011, DOI:10.3970/cmes.2011.078.095

Abstract By differentiating the Green function of Ting and Lee (1997) for 3D general anisotropic elastotatics in a spherical coordinate system as an intermediate step, and then using the chain rule, derivatives of up to the second order of this fundamental solution are obtained in exact, explicit, algebraic forms. No tensors of order higher than two are present in these derivatives, thereby allowing these quantities to be numerically evaluated quite expeditiously. These derivatives are required for the computation of the internal point displacements and stresses via Somigliana's identity in BEM analysis. Some examples are presented to demonstrate their successful implementation to… More >

E. Schnack¹, W. Weber², Y. Zhu³

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

Abstract Three dimensional fracture mechanics was done by several groups in the past. One topic for these three dimensional fracture mechanics is to consider re-entrant corners or wedges for isotropic material. An algorithm was developed in the past to compute the dominant eigenvalues for those problems with high accuracy. Based on Kondratiev's Lemma for elliptic boundary value problems it is started with the asymptotic for the displacement and stress distribution around these three dimensional corners. By considering the mixed boundary value problem, the field quantities in the vicinity of corner points are computed by using a special finite element formulation, which… More >

Chun Yang^1,2, Dalin Tang², Satya Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.1, pp. 53-78, 2011, DOI:10.3970/cmes.2011.072.053

Abstract Previously, we introduced a computational procedure based on three-dimensional meshless generalized finite difference (MGFD) method and serial magnetic resonance imaging (MRI) data to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. Structure-only models were used in our previous report. In this paper, fluid-stricture interaction (FSI) was added to improve on prediction accuracy. One participating patient was scanned three times (T1, T2, and T3, at intervals of about 18 months) to obtain plaque progression data. Blood flow was assumed to laminar, Newtonian, viscous and incompressible. The Navier-Stokes equations with arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing… More >

G.A. Gravvanis¹, K.M. Giannoutakis²

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.4, pp. 305-330, 2011, DOI:10.3970/cmes.2011.071.305

Abstract In this paper we present parallel explicit approximate inverse matrix techniques for solving sparse linear systems on shared memory systems, which are derived using the finite element method for biharmonic equations in three space variables. Our approach for solving such equations is by considering the biharmonic equation as a coupled equation approach (pair of Poisson equation), using a FE approximation scheme, yielding an inner-outer iteration method. Additionally, parallel approximate inverse matrix algorithms are introduced for the efficient solution of sparse linear systems, based on an anti-diagonal computational approach that eliminates the data dependencies. Parallel explicit preconditioned conjugate gradient-type schemes in… More >

Sunilkumar N¹, D Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.3, pp. 249-284, 2010, DOI:10.3970/cmes.2010.066.249

Abstract We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and 1D NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of… More >

P. Bettini¹, M. Midrio², R. Specogna²

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.2, pp. 117-134, 2010, DOI:10.3970/cmes.2010.066.117

Abstract In this paper we propose a geometric formulation to solve 3D electromagnetic wave problems in unbounded regions in the frequency domain. An absorbing boundary condition (ABC) is introduced to limit the size of the computational domain by means of anisotropic Perfectly Matched Layers (PML) absorbing media in the outer layers of an unstructured mesh. The numerical results of 3D benchmark problems are presented and the effect of the PML parameters and scaling functions on PML effectiveness are discussed. More >