3,866 research outputs found
Error estimates for interpolation of rough data using the scattered shifts of a radial basis function
The error between appropriately smooth functions and their radial basis
function interpolants, as the interpolation points fill out a bounded domain in
R^d, is a well studied artifact. In all of these cases, the analysis takes
place in a natural function space dictated by the choice of radial basis
function -- the native space. The native space contains functions possessing a
certain amount of smoothness. This paper establishes error estimates when the
function being interpolated is conspicuously rough.Comment: 12 page
Extending the range of error estimates for radial approximation in Euclidean space and on spheres
We adapt Schaback's error doubling trick [R. Schaback. Improved error bounds
for scattered data interpolation by radial basis functions. Math. Comp.,
68(225):201--216, 1999.] to give error estimates for radial interpolation of
functions with smoothness lying (in some sense) between that of the usual
native space and the subspace with double the smoothness. We do this for both
bounded subsets of R^d and spheres. As a step on the way to our ultimate goal
we also show convergence of pseudoderivatives of the interpolation error.Comment: 10 page
Scattered Data Interpolation on Embedded Submanifolds with Restricted Positive Definite Kernels: Sobolev Error Estimates
In this paper we investigate the approximation properties of kernel
interpolants on manifolds. The kernels we consider will be obtained by the
restriction of positive definite kernels on , such as radial basis
functions (RBFs), to a smooth, compact embedded submanifold \M\subset \R^d.
For restricted kernels having finite smoothness, we provide a complete
characterization of the native space on \M. After this and some preliminary
setup, we present Sobolev-type error estimates for the interpolation problem.
Numerical results verifying the theory are also presented for a one-dimensional
curve embedded in and a two-dimensional torus
A High-Order Radial Basis Function (RBF) Leray Projection Method for the Solution of the Incompressible Unsteady Stokes Equations
A new projection method based on radial basis functions (RBFs) is presented
for discretizing the incompressible unsteady Stokes equations in irregular
geometries. The novelty of the method comes from the application of a new
technique for computing the Leray-Helmholtz projection of a vector field using
generalized interpolation with divergence-free and curl-free RBFs. Unlike
traditional projection methods, this new method enables matching both
tangential and normal components of divergence-free vector fields on the domain
boundary. This allows incompressibility of the velocity field to be enforced
without any time-splitting or pressure boundary conditions. Spatial derivatives
are approximated using collocation with global RBFs so that the method only
requires samples of the field at (possibly scattered) nodes over the domain.
Numerical results are presented demonstrating high-order convergence in both
space (between 5th and 6th order) and time (up to 4th order) for some model
problems in two dimensional irregular geometries.Comment: 34 pages, 8 figure
Surface Spline Approximation on SO(3)
The purpose of this article is to introduce a new class of kernels on SO(3)
for approximation and interpolation, and to estimate the approximation power of
the associated spaces. The kernels we consider arise as linear combinations of
Green's functions of certain differential operators on the rotation group. They
are conditionally positive definite and have a simple closed-form expression,
lending themselves to direct implementation via, e.g., interpolation, or
least-squares approximation. To gauge the approximation power of the underlying
spaces, we introduce an approximation scheme providing precise L_p error
estimates for linear schemes, namely with L_p approximation order conforming to
the L_p smoothness of the target function.Comment: 22 pages, to appear in Appl. Comput. Harmon. Ana
A Partition of Unity Method for Divergence-Free or Curl-Free Radial Basis Function Approximation
Divergence-free (div-free) and curl-free vector fields are pervasive in many areas of science and engineering, from fluid dynamics to electromagnetism. A common problem that arises in applications is that of constructing smooth approximants to these vector fields and/or their potentials based only on discrete samples. Additionally, it is often necessary that the vector approximants preserve the div-free or curl-free properties of the field to maintain certain physical constraints. Div/curl-free radial basis functions (RBFs) are a particularly good choice for this application as they are meshfree and analytically satisfy the div-free or curl-free property. However, this method can be computationally expensive due to its global nature. In this paper, we develop a technique for bypassing this issue that combines div/curl-free RBFs in a partition of unity framework, where one solves for local approximants over subsets of the global samples and then blends them together to form a div-free or curl-free global approximant. The method is applicable to div/curl-free vector fields in ℝ2 and tangential fields on two-dimensional surfaces, such as the sphere, and the curl-free method can be generalized to vector fields in ℝd. The method also produces an approximant for the scalar potential of the underlying sampled field. We present error estimates and demonstrate the effectiveness of the method on several test problems
Doctor of Philosophy
dissertationPlatelet aggregation, an important part of the development of blood clots, is a complex process involving both mechanical interaction between platelets and blood, and chemical transport on and o the surfaces of those platelets. Radial Basis Function (RBF) interpolation is a meshfree method for the interpolation of multidimensional scattered data, and therefore well-suited for the development of meshfree numerical methods. This dissertation explores the use of RBF interpolation for the simulation of both the chemistry and mechanics of platelet aggregation. We rst develop a parametric RBF representation for closed platelet surfaces represented by scattered nodes in both two and three dimensions. We compare this new RBF model to Fourier models in terms of computational cost and errors in shape representation. We then augment the Immersed Boundary (IB) method, a method for uid-structure interaction, with our RBF geometric model. We apply the resultant method to a simulation of platelet aggregation, and present comparisons against the traditional IB method. We next consider a two-dimensional problem where platelets are suspended in a stationary fluid, with chemical diusion in the fluid and chemical reaction-diusion on platelet surfaces. To tackle the latter, we propose a new method based on RBF-generated nite dierences (RBF-FD) for solving partial dierential equations (PDEs) on surfaces embedded in 2D domains. To robustly tackle the former, we remove a limitation of the Augmented Forcing method (AFM), a method for solving PDEs on domains containing curved objects, using RBF-based symmetric Hermite interpolation. Next, we extend our RBF-FD method to the numerical solution of PDEs on surfaces embedded in 3D domains, proposing a new method of stabilizing RBF-FD discretizations on surfaces. We perform convergence studies and present applications motivated by biology. We conclude with a summary of the thesis research and present an overview of future research directions, including spectrally-accurate projection methods, an extension of the Regularized Stokeslet method, RBF-FD for variable-coecient diusion, and boundary conditions for RBF-FD
- …