13 research outputs found
Analysis of and workarounds for element reversal for a finite element-based algorithm for warping triangular and tetrahedral meshes
We consider an algorithm called FEMWARP for warping triangular and
tetrahedral finite element meshes that computes the warping using the finite
element method itself. The algorithm takes as input a two- or three-dimensional
domain defined by a boundary mesh (segments in one dimension or triangles in
two dimensions) that has a volume mesh (triangles in two dimensions or
tetrahedra in three dimensions) in its interior. It also takes as input a
prescribed movement of the boundary mesh. It computes as output updated
positions of the vertices of the volume mesh. The first step of the algorithm
is to determine from the initial mesh a set of local weights for each interior
vertex that describes each interior vertex in terms of the positions of its
neighbors. These weights are computed using a finite element stiffness matrix.
After a boundary transformation is applied, a linear system of equations based
upon the weights is solved to determine the final positions of the interior
vertices. The FEMWARP algorithm has been considered in the previous literature
(e.g., in a 2001 paper by Baker). FEMWARP has been succesful in computing
deformed meshes for certain applications. However, sometimes FEMWARP reverses
elements; this is our main concern in this paper. We analyze the causes for
this undesirable behavior and propose several techniques to make the method
more robust against reversals. The most successful of the proposed methods
includes combining FEMWARP with an optimization-based untangler.Comment: Revision of earlier version of paper. Submitted for publication in
BIT Numerical Mathematics on 27 April 2010. Accepted for publication on 7
September 2010. Published online on 9 October 2010. The final publication is
available at http://www.springerlink.co
The Sudbury Neutrino Observatory
The Sudbury Neutrino Observatory is a second generation water Cherenkov
detector designed to determine whether the currently observed solar neutrino
deficit is a result of neutrino oscillations. The detector is unique in its use
of D2O as a detection medium, permitting it to make a solar model-independent
test of the neutrino oscillation hypothesis by comparison of the charged- and
neutral-current interaction rates. In this paper the physical properties,
construction, and preliminary operation of the Sudbury Neutrino Observatory are
described. Data and predicted operating parameters are provided whenever
possible.Comment: 58 pages, 12 figures, submitted to Nucl. Inst. Meth. Uses elsart and
epsf style files. For additional information about SNO see
http://www.sno.phy.queensu.ca . This version has some new reference