603 research outputs found
B\'ezier curves that are close to elastica
We study the problem of identifying those cubic B\'ezier curves that are
close in the L2 norm to planar elastic curves. The problem arises in design
situations where the manufacturing process produces elastic curves; these are
difficult to work with in a digital environment. We seek a sub-class of special
B\'ezier curves as a proxy. We identify an easily computable quantity, which we
call the lambda-residual, that accurately predicts a small L2 distance. We then
identify geometric criteria on the control polygon that guarantee that a
B\'ezier curve has lambda-residual below 0.4, which effectively implies that
the curve is within 1 percent of its arc-length to an elastic curve in the L2
norm. Finally we give two projection algorithms that take an input B\'ezier
curve and adjust its length and shape, whilst keeping the end-points and
end-tangent angles fixed, until it is close to an elastic curve.Comment: 13 pages, 15 figure
Recommended from our members
A rational cubic spline with tension
A rational cubic spline curve is described which has tension control parameters for manipulating the shape of the curve. The spline is presented in both interpolatory and rational B-spline forms, and the behaviour of the resulting representations is analysed with respect to variation of the control parameters
Fast Isogeometric Boundary Element Method based on Independent Field Approximation
An isogeometric boundary element method for problems in elasticity is
presented, which is based on an independent approximation for the geometry,
traction and displacement field. This enables a flexible choice of refinement
strategies, permits an efficient evaluation of geometry related information, a
mixed collocation scheme which deals with discontinuous tractions along
non-smooth boundaries and a significant reduction of the right hand side of the
system of equations for common boundary conditions. All these benefits are
achieved without any loss of accuracy compared to conventional isogeometric
formulations. The system matrices are approximated by means of hierarchical
matrices to reduce the computational complexity for large scale analysis. For
the required geometrical bisection of the domain, a strategy for the evaluation
of bounding boxes containing the supports of NURBS basis functions is
presented. The versatility and accuracy of the proposed methodology is
demonstrated by convergence studies showing optimal rates and real world
examples in two and three dimensions.Comment: 32 pages, 27 figure
Discontinuities in numerical radiative transfer
Observations and magnetohydrodynamic simulations of solar and stellar
atmospheres reveal an intermittent behavior or steep gradients in physical
parameters, such as magnetic field, temperature, and bulk velocities. The
numerical solution of the stationary radiative transfer equation is
particularly challenging in such situations, because standard numerical methods
may perform very inefficiently in the absence of local smoothness. However, a
rigorous investigation of the numerical treatment of the radiative transfer
equation in discontinuous media is still lacking. The aim of this work is to
expose the limitations of standard convergence analyses for this problem and to
identify the relevant issues. Moreover, specific numerical tests are performed.
These show that discontinuities in the atmospheric physical parameters
effectively induce first-order discontinuities in the radiative transfer
equation, reducing the accuracy of the solution and thwarting high-order
convergence. In addition, a survey of the existing numerical schemes for
discontinuous ordinary differential systems and interpolation techniques for
discontinuous discrete data is given, evaluating their applicability to the
radiative transfer problem
A discrete methodology for controlling the sign of curvature and torsion for NURBS
This paper develops a discrete methodology for approximating the so-called convex domain of a NURBS curve, namely the domain in the ambient space, where a user-specified control point is free to move so that the curvature and torsion retains its sign along the NURBS parametric domain of definition. The methodology provides a monotonic sequence of convex polyhedra, converging from the interior to the convex domain. If the latter is non-empty, a simple algorithm is proposed, that yields a sequence of polytopes converging uniformly to the restriction of the convex domain to any user-specified bounding box. The algorithm is illustrated for a pair of planar and a spatial Bézier configuration
- …