A. Grimaldi¹, G. Pascazio², M. Napolitano³

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.1, pp. 45-56, 2006, DOI:10.3970/cmes.2006.014.045

Abstract This paper provides a numerical method for solving two- and three-dimensional unsteady incompressible flows. The vorticity-velocity formulation of the Navier--Stokes equations is considered, employing the vorticity transport equation and a second-order Poisson equation for the velocity. Second-order-accurate centred finite differences on a staggered grid are used for the space discretization. The vorticity equation is discretized in time using a fully implicit three-level scheme. At each physical time level, a dual-time stepping technique is used to solve the coupled system of non linear algebraic equations by various efficient relaxation schemes. Steady flows are computed by dropping More >

A. N. Guz, V. A. Dekret¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 165-170, 2006, DOI:10.3970/cmes.2006.013.165

Abstract Stability problem of composite material reinforced by two parallel short fibers is solved. The problem is formulated with application of equations of linearized three-dimensional theory of stability. The composite is modeled as piecewise-homogeneous medium. The influence of geometrical and mechanical parameters of the material on critical strain is investigated. More >

Sunil Punjabi¹, K. Muralidhar², P. K. Panigrahi²

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.2, pp. 95-106, 2006, DOI:10.3970/fdmp.2006.002.095

Abstract Convection in a differentially heated two-layer system consisting of air and water was studied experimentally, using laser-interferometry. The cavity used for flow visualization was square in cross-section and rectangular in-plan having dimensions of 447 × 32 × 32 mm³. Experiments performed over different layer thicknesses of water filled in a square cross-section cavity, the rest being air, are reported in the present work. The following temperature differences for each layer height were imposed across the hot and the cold walls of the superposed fluid layers: (i) ΔT=10K and (ii)ΔT =18 K. The present study was aimed… More >

B.B. Karki¹, D. Bhattarai¹, L. Stixrude²

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 107-118, 2006, DOI:10.3970/cmc.2006.003.107

Abstract Computer modeling of liquid phase poses tremendous challenge: It requires a relatively large simulation size, long simulation time and accurate interatomic interaction and as such, it produces massive amounts of data. Recent advances in hardware and software have made it possible to accurately simulate the liquid phase. This paper reports the details of methodology used in the context of liquid simulations and subsequent analysis of the output data. For illustration purpose, we consider the results for the liquid phases of two geophysically relevant materials, namely MgO and MgSiO₃. The simulations are performed using the parallel first-principles More >

J. Ma¹, H. Lu¹, B. Wang¹, S. Roy¹, R. Hornung², A. Wissink², R. Komanduri^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 135-152, 2005, DOI:10.3970/cmes.2005.008.135

Abstract In the simulation of a wide range of mechanics problems including impact/contact/penetration and fracture, the material point method (MPM), Sulsky, Zhou and Shreyer (1995), demonstrated its computational capabilities. To resolve alternating stress sign and instability problems associated with conventional MPM, Bardenhagen and Kober (2004) introduced recently the generalized interpolation material point (GIMP) method and implemented for one-dimensional simulations. In this paper we have extended GIMP to 2D and applied to simulate simple tension and indentation problems. For simulations spanning multiple length scales, based on the continuum mechanics approach, we present a parallel GIMP computational method… More >

J. Bielak¹, O. Ghattas², E.-J. Kim³

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.2, pp. 99-112, 2005, DOI:10.3970/cmes.2005.010.099

Abstract We present a parallel octree-based finite element method for large-scale earthquake ground motion simulation in realistic basins. The octree representation combines the low memory per node and good cache performance of finite difference methods with the spatial adaptivity to local seismic wavelengths characteristic of unstructured finite element methods. Several tests are provided to verify the numerical performance of the method against Green's function solutions for homogeneous and piecewise homogeneous media, both with and without anelastic attenuation. A comparison is also provided against a finite difference code and an unstructured tetrahedral finite element code for a More >

Marcello Lappa^{1, 2}

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.1, pp. 1-20, 2005, DOI:10.3970/fdmp.2005.001.001

Abstract Possible natural transport mechanisms in cubical and shallow cavities with different heating conditions (from below or from the side) are investigated by means of numerical solution of the non-linear model equations and multiprocessor computations. Attention is focused on a variety of three-dimensional steady effects that can arise in such configurations in the case of low-Pr liquids (silicon melt) even for relatively small values of the temperature gradient due to localized boundary effects and/or true instabilities of the flow. Such aspects are still poorly known or completely ignored owing to the fact that most of the More >

T. Furumura¹, L. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.2, pp. 153-168, 2004, DOI:10.3970/cmes.2004.006.153

Abstract Recent developments of the Earth Simulator, a high-performance parallel computer, has made it possible to realize realistic 3D simulations of seismic wave propagations on a regional scale including higher frequencies. Paralleling this development, the deployment of dense networks of strong ground motion instruments in Japan (K-NET and KiK-net) has now made it possible to directly visualize regional seismic wave propagation during large earthquakes. Our group has developed an efficient parallel finite difference method (FDM) code for modeling the seismic wavefield and a 3D visualization technique, both suitable for implementation on the Earth Simulator. Large-scale 3D… More >

A. Rama Mohan Rao¹, T.V.S.R. Appa Rao², B. Dattaguru³

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.3, pp. 213-234, 2004, DOI:10.3970/cmes.2004.005.213

Abstract This paper presents an algorithm for automatic partitioning of unstructured meshes for parallel finite element computations employing float-encoded genetic algorithms (FEGA). The problem of mesh partitioning is represented in such a way that the number of variables considered in the genome (chromosome) construction is constant irrespective of the size of the problem. In order to accelerate the computational process, several acceleration techniques like constraining the search space, local improvement after initial global partitioning have been attempted. Finally, micro float-encoded genetic algorithms have been developed to accelerate the computational process. More >

John L. Volakis¹, Kubilay Sertel¹, Erik Jørgensen², Rick W. Kindt¹

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.5, pp. 463-476, 2004, DOI:10.3970/cmes.2004.005.463

Abstract In this paper we address several topics relating to the development and implementation of volume integral and hybrid finite element methods for electromagnetic modeling. Comparisons of volume integral equation formulations with the finite element-boundary integral method are given in terms of accuracy and computing resources. We also discuss preconditioning and parallelization of the multilevel fast multipole method, and propose higher-order basis functions for curvilinear quadrilaterals and volumetric basis functions for curvilinear hexahedra. The latter have the desirable property of vanishing divergence within the element but non-zero curl. In addition, a new domain decomposition is introduced More >