464 research outputs found

    Local bisection refinement for nn-simplicial grids generated by reflection

    Get PDF
    A simple local bisection refinement algorithm for the adaptive refinement of nn-simplicial grids is presented. The algorithm requires that the vertices of each simplex be ordered in a special way relative to those in neighboring simplices. It is proven that certain regular simplicial grids on [0,1]n[0,1]^n have this property, and the more general grids to which this method is applicable are discussed. The edges to be bisected are determined by an ordering of the simplex vertices, without local or global computation or communication. Further, the number of congruency classes in a locally refined grid turns out to be bounded above by nn, independent of the level of refinement. Simplicial grids of higher dimension are frequently used to approximate solution manifolds of parametrized equations, for instance, as in [W. C. Rheinboldt, Numer. Math., 53 (1988), pp. 165–180] and [E. Allgower and K. Georg, Utilitas Math., 16 (1979), pp. 123–129]. They are also used for the determination of fixed points of functions from Rn{\bf R}^n to Rn{\bf R}^n , as described in [M. J. Todd, Lecture Notes in Economic and Mathematical Systems, 124, Springer-Verlag, Berlin, 1976]. In two and three dimensions, such grids of triangles, respectively, tetrahedrons, are used for the computation of finite element solutions of partial differential equations, for example, as in [O. Axelsson and V. A. Barker, Finite Element Solution of Boundary Value Problems, Academic Press, Orlando, 1984], [R. E. Bank and B. D. Welfert, SIAM J. Numer. Anal., 28 (1991), pp. 591–623], [W. F. Mitchell, SIAM J. Sci. Statist. Comput., 13 (1992), pp. 146–147], and [M. C. Rivara, J. Comput. Appl. Math., 36 (1991), pp. 79–89]. The new method is applicable to any triangular grid and may possibly be applied to many tetrahedral grids using additional closure refinement to avoid incompatibilities

    Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

    Get PDF
    One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. In the past, limited computer resources have motivated the development of efficient numerical schemes in computational fluid dynamics (CFD) utilizing structured meshes. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Unstructured grids on the other hand offer an alternative to modeling complex geometries. Unstructured meshes have irregular connectivity and usually contain combinations of triangles, quadrilaterals, tetrahedra, and hexahedra. The generation and use of unstructured grids poses new challenges in CFD. The purpose of this note is to present recent developments in the unstructured grid generation and flow solution technology

    Master of Science

    Get PDF
    thesisWe present a procedure for generating a coarse, high-quality, tetrahedral mesh whose exterior surface encloses and approximates a given triangle mesh. A tetrahedral mesh is useful for computing perturbation of the triangle mesh based on continuum mechanics: perturbation such as plastic flow, fracture, and elastic deformation. The computer graphics community has long used this physics-based simulation to produce animations of objects exhibiting such physical phenomena. Interactive animation applications such as industrial design, medical training, and computer entertainment require meshes that are particularly efficient and robust, and our meshing procedure targets these properties. We begin with a BCC background lattice and sculpt an initial mesh from it whose tetrahedra occupy some of the volume bounded by the triangle mesh. We then refine this initial mesh with an iterative optimization procedure that simultaneously minimizes the distance from the triangle mesh to the surface of the tetrahedral mesh and maximizes the numerical quality of the tetrahedra. Our procedure provides a trade-off among the mesh's quality, resolution, and degree of approximation of the triangle mesh

    Parallel software tool for decomposing and meshing of 3d structures

    Get PDF
    An algorithm for automatic parallel generation of three-dimensional unstructured computational meshes based on geometrical domain decomposition is proposed in this paper. Software package build upon proposed algorithm is described. Several practical examples of mesh generation on multiprocessor computational systems are given. It is shown that developed parallel algorithm enables us to reduce mesh generation time significantly (dozens of times). Moreover, it easily produces meshes with number of elements of order 5 · 107, construction of those on a single CPU is problematic. Questions of time consumption, efficiency of computations and quality of generated meshes are also considered

    Activities of the Research Institute for Advanced Computer Science

    Get PDF
    The Research Institute for Advanced Computer Science (RIACS) was established by the Universities Space Research Association (USRA) at the NASA Ames Research Center (ARC) on June 6, 1983. RIACS is privately operated by USRA, a consortium of universities with research programs in the aerospace sciences, under contract with NASA. The primary mission of RIACS is to provide research and expertise in computer science and scientific computing to support the scientific missions of NASA ARC. The research carried out at RIACS must change its emphasis from year to year in response to NASA ARC's changing needs and technological opportunities. Research at RIACS is currently being done in the following areas: (1) parallel computing; (2) advanced methods for scientific computing; (3) high performance networks; and (4) learning systems. RIACS technical reports are usually preprints of manuscripts that have been submitted to research journals or conference proceedings. A list of these reports for the period January 1, 1994 through December 31, 1994 is in the Reports and Abstracts section of this report

    Foundations of space-time finite element methods: polytopes, interpolation, and integration

    Full text link
    The main purpose of this article is to facilitate the implementation of space-time finite element methods in four-dimensional space. In order to develop a finite element method in this setting, it is necessary to create a numerical foundation, or equivalently a numerical infrastructure. This foundation should include a collection of suitable elements (usually hypercubes, simplices, or closely related polytopes), numerical interpolation procedures (usually orthonormal polynomial bases), and numerical integration procedures (usually quadrature rules). It is well known that each of these areas has yet to be fully explored, and in the present article, we attempt to directly address this issue. We begin by developing a concrete, sequential procedure for constructing generic four-dimensional elements (4-polytopes). Thereafter, we review the key numerical properties of several canonical elements: the tesseract, tetrahedral prism, and pentatope. Here, we provide explicit expressions for orthonormal polynomial bases on these elements. Next, we construct fully symmetric quadrature rules with positive weights that are capable of exactly integrating high-degree polynomials, e.g. up to degree 17 on the tesseract. Finally, the quadrature rules are successfully tested using a set of canonical numerical experiments on polynomial and transcendental functions.Comment: 34 pages, 18 figure

    Trefftz Difference Schemes on Irregular Stencils

    Full text link
    The recently developed Flexible Local Approximation MEthod (FLAME) produces accurate difference schemes by replacing the usual Taylor expansion with Trefftz functions -- local solutions of the underlying differential equation. This paper advances and casts in a general form a significant modification of FLAME proposed recently by Pinheiro & Webb: a least-squares fit instead of the exact match of the approximate solution at the stencil nodes. As a consequence of that, FLAME schemes can now be generated on irregular stencils with the number of nodes substantially greater than the number of approximating functions. The accuracy of the method is preserved but its robustness is improved. For demonstration, the paper presents a number of numerical examples in 2D and 3D: electrostatic (magnetostatic) particle interactions, scattering of electromagnetic (acoustic) waves, and wave propagation in a photonic crystal. The examples explore the role of the grid and stencil size, of the number of approximating functions, and of the irregularity of the stencils.Comment: 28 pages, 12 figures; to be published in J Comp Phy

    Diamond-based models for scientific visualization

    Get PDF
    Hierarchical spatial decompositions are a basic modeling tool in a variety of application domains including scientific visualization, finite element analysis and shape modeling and analysis. A popular class of such approaches is based on the regular simplex bisection operator, which bisects simplices (e.g. line segments, triangles, tetrahedra) along the midpoint of a predetermined edge. Regular simplex bisection produces adaptive simplicial meshes of high geometric quality, while simplifying the extraction of crack-free, or conforming, approximations to the original dataset. Efficient multiresolution representations for such models have been achieved in 2D and 3D by clustering sets of simplices sharing the same bisection edge into structures called diamonds. In this thesis, we introduce several diamond-based approaches for scientific visualization. We first formalize the notion of diamonds in arbitrary dimensions in terms of two related simplicial decompositions of hypercubes. This enables us to enumerate the vertices, simplices, parents and children of a diamond. In particular, we identify the number of simplices involved in conforming updates to be factorial in the dimension and group these into a linear number of subclusters of simplices that are generated simultaneously. The latter form the basis for a compact pointerless representation for conforming meshes generated by regular simplex bisection and for efficiently navigating the topological connectivity of these meshes. Secondly, we introduce the supercube as a high-level primitive on such nested meshes based on the atomic units within the underlying triangulation grid. We propose the use of supercubes to associate information with coherent subsets of the full hierarchy and demonstrate the effectiveness of such a representation for modeling multiresolution terrain and volumetric datasets. Next, we introduce Isodiamond Hierarchies, a general framework for spatial access structures on a hierarchy of diamonds that exploits the implicit hierarchical and geometric relationships of the diamond model. We use an isodiamond hierarchy to encode irregular updates to a multiresolution isosurface or interval volume in terms of regular updates to diamonds. Finally, we consider nested hypercubic meshes, such as quadtrees, octrees and their higher dimensional analogues, through the lens of diamond hierarchies. This allows us to determine the relationships involved in generating balanced hypercubic meshes and to propose a compact pointerless representation of such meshes. We also provide a local diamond-based triangulation algorithm to generate high-quality conforming simplicial meshes
    • …
    corecore