4,855 research outputs found
There are 174 Subdivisions of the Hexahedron into Tetrahedra
This article answers an important theoretical question: How many different
subdivisions of the hexahedron into tetrahedra are there? It is well known that
the cube has five subdivisions into 6 tetrahedra and one subdivision into 5
tetrahedra. However, all hexahedra are not cubes and moving the vertex
positions increases the number of subdivisions. Recent hexahedral dominant
meshing methods try to take these configurations into account for combining
tetrahedra into hexahedra, but fail to enumerate them all: they use only a set
of 10 subdivisions among the 174 we found in this article.
The enumeration of these 174 subdivisions of the hexahedron into tetrahedra
is our combinatorial result. Each of the 174 subdivisions has between 5 and 15
tetrahedra and is actually a class of 2 to 48 equivalent instances which are
identical up to vertex relabeling. We further show that exactly 171 of these
subdivisions have a geometrical realization, i.e. there exist coordinates of
the eight hexahedron vertices in a three-dimensional space such that the
geometrical tetrahedral mesh is valid. We exhibit the tetrahedral meshes for
these configurations and show in particular subdivisions of hexahedra with 15
tetrahedra that have a strictly positive Jacobian
On the dimension of spline spaces on planar T-meshes
We analyze the space of bivariate functions that are piecewise polynomial of
bi-degree \textless{}= (m, m') and of smoothness r along the interior edges of
a planar T-mesh. We give new combinatorial lower and upper bounds for the
dimension of this space by exploiting homological techniques. We relate this
dimension to the weight of the maximal interior segments of the T-mesh, defined
for an ordering of these maximal interior segments. We show that the lower and
upper bounds coincide, for high enough degrees or for hierarchical T-meshes
which are enough regular. We give a rule of subdivision to construct
hierarchical T-meshes for which these lower and upper bounds coincide. Finally,
we illustrate these results by analyzing spline spaces of small degrees and
smoothness
Finding Hexahedrizations for Small Quadrangulations of the Sphere
This paper tackles the challenging problem of constrained hexahedral meshing.
An algorithm is introduced to build combinatorial hexahedral meshes whose
boundary facets exactly match a given quadrangulation of the topological
sphere. This algorithm is the first practical solution to the problem. It is
able to compute small hexahedral meshes of quadrangulations for which the
previously known best solutions could only be built by hand or contained
thousands of hexahedra. These challenging quadrangulations include the
boundaries of transition templates that are critical for the success of general
hexahedral meshing algorithms.
The algorithm proposed in this paper is dedicated to building combinatorial
hexahedral meshes of small quadrangulations and ignores the geometrical
problem. The key idea of the method is to exploit the equivalence between quad
flips in the boundary and the insertion of hexahedra glued to this boundary.
The tree of all sequences of flipping operations is explored, searching for a
path that transforms the input quadrangulation Q into a new quadrangulation for
which a hexahedral mesh is known. When a small hexahedral mesh exists, a
sequence transforming Q into the boundary of a cube is found; otherwise, a set
of pre-computed hexahedral meshes is used.
A novel approach to deal with the large number of problem symmetries is
proposed. Combined with an efficient backtracking search, it allows small
shellable hexahedral meshes to be found for all even quadrangulations with up
to 20 quadrangles. All 54,943 such quadrangulations were meshed using no more
than 72 hexahedra. This algorithm is also used to find a construction to fill
arbitrary domains, thereby proving that any ball-shaped domain bounded by n
quadrangles can be meshed with no more than 78 n hexahedra. This very
significantly lowers the previous upper bound of 5396 n.Comment: Accepted for SIGGRAPH 201
Low-order continuous finite element spaces on hybrid non-conforming hexahedral-tetrahedral meshes
This article deals with solving partial differential equations with the
finite element method on hybrid non-conforming hexahedral-tetrahedral meshes.
By non-conforming, we mean that a quadrangular face of a hexahedron can be
connected to two triangular faces of tetrahedra. We introduce a set of
low-order continuous (C0) finite element spaces defined on these meshes. They
are built from standard tri-linear and quadratic Lagrange finite elements with
an extra set of constraints at non-conforming hexahedra-tetrahedra junctions to
recover continuity. We consider both the continuity of the geometry and the
continuity of the function basis as follows: the continuity of the geometry is
achieved by using quadratic mappings for tetrahedra connected to tri-affine
hexahedra and the continuity of interpolating functions is enforced in a
similar manner by using quadratic Lagrange basis on tetrahedra with constraints
at non-conforming junctions to match tri-linear hexahedra. The so-defined
function spaces are validated numerically on simple Poisson and linear
elasticity problems for which an analytical solution is known. We observe that
using a hybrid mesh with the proposed function spaces results in an accuracy
significantly better than when using linear tetrahedra and slightly worse than
when solely using tri-linear hexahedra. As a consequence, the proposed function
spaces may be a promising alternative for complex geometries that are out of
reach of existing full hexahedral meshing methods
Subdivision Directional Fields
We present a novel linear subdivision scheme for face-based tangent
directional fields on triangle meshes. Our subdivision scheme is based on a
novel coordinate-free representation of directional fields as halfedge-based
scalar quantities, bridging the finite-element representation with discrete
exterior calculus. By commuting with differential operators, our subdivision is
structure-preserving: it reproduces curl-free fields precisely, and reproduces
divergence-free fields in the weak sense. Moreover, our subdivision scheme
directly extends to directional fields with several vectors per face by working
on the branched covering space. Finally, we demonstrate how our scheme can be
applied to directional-field design, advection, and robust earth mover's
distance computation, for efficient and robust computation
Cut Size Statistics of Graph Bisection Heuristics
We investigate the statistical properties of cut sizes generated by heuristic
algorithms which solve approximately the graph bisection problem. On an
ensemble of sparse random graphs, we find empirically that the distribution of
the cut sizes found by ``local'' algorithms becomes peaked as the number of
vertices in the graphs becomes large. Evidence is given that this distribution
tends towards a Gaussian whose mean and variance scales linearly with the
number of vertices of the graphs. Given the distribution of cut sizes
associated with each heuristic, we provide a ranking procedure which takes into
account both the quality of the solutions and the speed of the algorithms. This
procedure is demonstrated for a selection of local graph bisection heuristics.Comment: 17 pages, 5 figures, submitted to SIAM Journal on Optimization also
available at http://ipnweb.in2p3.fr/~martin
- …