2,073 research outputs found
GRID2D/3D: A computer program for generating grid systems in complex-shaped two- and three-dimensional spatial domains. Part 1: Theory and method
An efficient computer program, called GRID2D/3D was developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2- and 3-D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation in which the distribution of grid points within the spatial domain is controlled by stretching functions. All single grid systems generated by GRID2D/3D can have grid lines that are continuous and differentiable everywhere up to the second-order. Also, grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For continuous composite grid systems, the grid lines are continuous and differentiable everywhere up to the second-order except at interfaces where different single grid systems meet. At interfaces where different single grid systems meet, the grid lines are only differentiable up to the first-order. For 2-D spatial domains, the boundary curves are described by using either cubic or tension spline interpolation. For 3-D spatial domains, the boundary surfaces are described by using either linear Coon's interpolation, bi-hyperbolic spline interpolation, or a new technique referred to as 3-D bi-directional Hermite interpolation. Since grid systems generated by algebraic methods can have grid lines that overlap one another, GRID2D/3D contains a graphics package for evaluating the grid systems generated. With the graphics package, the user can generate grid systems in an interactive manner with the grid generation part of GRID2D/3D. GRID2D/3D is written in FORTRAN 77 and can be run on any IBM PC, XT, or AT compatible computer. In order to use GRID2D/3D on workstations or mainframe computers, some minor modifications must be made in the graphics part of the program; no modifications are needed in the grid generation part of the program. This technical memorandum describes the theory and method used in GRID2D/3D
Grid generation for the solution of partial differential equations
A general survey of grid generators is presented with a concern for understanding why grids are necessary, how they are applied, and how they are generated. After an examination of the need for meshes, the overall applications setting is established with a categorization of the various connectivity patterns. This is split between structured grids and unstructured meshes. Altogether, the categorization establishes the foundation upon which grid generation techniques are developed. The two primary categories are algebraic techniques and partial differential equation techniques. These are each split into basic parts, and accordingly are individually examined in some detail. In the process, the interrelations between the various parts are accented. From the established background in the primary techniques, consideration is shifted to the topic of interactive grid generation and then to adaptive meshes. The setting for adaptivity is established with a suitable means to monitor severe solution behavior. Adaptive grids are considered first and are followed by adaptive triangular meshes. Then the consideration shifts to the temporal coupling between grid generators and PDE-solvers. To conclude, a reflection upon the discussion, herein, is given
Exponential Splines and Pseudo-Splines: Generation versus reproduction of exponential polynomials
Subdivision schemes are iterative methods for the design of smooth curves and
surfaces. Any linear subdivision scheme can be identified by a sequence of
Laurent polynomials, also called subdivision symbols, which describe the linear
rules determining successive refinements of coarse initial meshes. One
important property of subdivision schemes is their capability of exactly
reproducing in the limit specific types of functions from which the data is
sampled. Indeed, this property is linked to the approximation order of the
scheme and to its regularity. When the capability of reproducing polynomials is
required, it is possible to define a family of subdivision schemes that allows
to meet various demands for balancing approximation order, regularity and
support size. The members of this family are known in the literature with the
name of pseudo-splines. In case reproduction of exponential polynomials instead
of polynomials is requested, the resulting family turns out to be the
non-stationary counterpart of the one of pseudo-splines, that we here call the
family of exponential pseudo-splines. The goal of this work is to derive the
explicit expressions of the subdivision symbols of exponential pseudo-splines
and to study their symmetry properties as well as their convergence and
regularity.Comment: 25 page
Interpolation of equation-of-state data
Aims. We use Hermite splines to interpolate pressure and its derivatives
simultaneously, thereby preserving mathematical relations between the
derivatives. The method therefore guarantees that thermodynamic identities are
obeyed even between mesh points. In addition, our method enables an estimation
of the precision of the interpolation by comparing the Hermite-spline results
with those of frequent cubic (B-) spline interpolation.
Methods. We have interpolated pressure as a function of temperature and
density with quintic Hermite 2D-splines. The Hermite interpolation requires
knowledge of pressure and its first and second derivatives at every mesh point.
To obtain the partial derivatives at the mesh points, we used tabulated values
if given or else thermodynamic equalities, or, if not available, values
obtained by differentiating B-splines.
Results. The results were obtained with the grid of the SAHA-S
equation-of-state (EOS) tables. The maximum difference lies in the range
from to , and difference varies from to
. Specifically, for the points of a solar model, the maximum
differences are one order of magnitude smaller than the aforementioned values.
The poorest precision is found in the dissociation and ionization regions,
occurring at K. The best precision is achieved at
higher temperatures, K. To discuss the significance of the
interpolation errors we compare them with the corresponding difference between
two different equation-of-state formalisms, SAHA-S and OPAL 2005. We find that
the interpolation errors of the pressure are a few orders of magnitude less
than the differences from between the physical formalisms, which is
particularly true for the solar-model points.Comment: Accepted for publication in A&
- …