6,976 research outputs found
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
A priori error for unilateral contact problems with Lagrange multiplier and IsoGeometric Analysis
In this paper, we consider unilateral contact problem without friction
between a rigid body and deformable one in the framework of isogeometric
analysis. We present the theoretical analysis of the mixed problem using an
active-set strategy and for a primal space of NURBS of degree and for
a dual space of B-Spline. A inf-sup stability is proved to ensure a good
property of the method. An optimal a priori error estimate is demonstrated
without assumption on the unknown contact set. Several numerical examples in
two- and three-dimensional and in small and large deformation demonstrate the
accuracy of the proposed method
Einstein equations in the null quasi-spherical gauge III: numerical algorithms
We describe numerical techniques used in the construction of our 4th order
evolution for the full Einstein equations, and assess the accuracy of
representative solutions. The code is based on a null gauge with a
quasi-spherical radial coordinate, and simulates the interaction of a single
black hole with gravitational radiation. Techniques used include spherical
harmonic representations, convolution spline interpolation and filtering, and
an RK4 "method of lines" evolution. For sample initial data of "intermediate"
size (gravitational field with 19% of the black hole mass), the code is
accurate to 1 part in 10^5, until null time z=55 when the coordinate condition
breaks down.Comment: Latex, 38 pages, 29 figures (360Kb compressed
Local interpolation schemes for landmark-based image registration: a comparison
In this paper we focus, from a mathematical point of view, on properties and
performances of some local interpolation schemes for landmark-based image
registration. Precisely, we consider modified Shepard's interpolants,
Wendland's functions, and Lobachevsky splines. They are quite unlike each
other, but all of them are compactly supported and enjoy interesting
theoretical and computational properties. In particular, we point out some
unusual forms of the considered functions. Finally, detailed numerical
comparisons are given, considering also Gaussians and thin plate splines, which
are really globally supported but widely used in applications
Orthonormal bases of regular wavelets in spaces of homogeneous type
Adapting the recently developed randomized dyadic structures, we introduce
the notion of spline function in geometrically doubling quasi-metric spaces.
Such functions have interpolation and reproducing properties as the linear
splines in Euclidean spaces. They also have H\"older regularity. This is used
to build an orthonormal basis of H\"older-continuous wavelets with exponential
decay in any space of homogeneous type. As in the classical theory, wavelet
bases provide a universal Calder\'on reproducing formula to study and develop
function space theory and singular integrals. We discuss the examples of
spaces, BMO and apply this to a proof of the T(1) theorem. As no extra
condition {(like 'reverse doubling', 'small boundary' of balls, etc.)} on the
space of homogeneous type is required, our results extend a long line of works
on the subject.Comment: We have made improvements to section 2 following the referees
suggestions. In particular, it now contains full proof of formerly Theorem
2.7 instead of sending back to earlier works, which makes the construction of
splines self-contained. One reference adde
- …