11,196 research outputs found
A Simple Method for Computing Singular or Nearly Singular Integrals on Closed Surfaces
We present a simple, accurate method for computing singular or nearly
singular integrals on a smooth, closed surface, such as layer potentials for
harmonic functions evaluated at points on or near the surface. The integral is
computed with a regularized kernel and corrections are added for regularization
and discretization, which are found from analysis near the singular point. The
surface integrals are computed from a new quadrature rule using surface points
which project onto grid points in coordinate planes. The method does not
require coordinate charts on the surface or special treatment of the
singularity other than the corrections. The accuracy is about , where
is the spacing in the background grid, uniformly with respect to the point
of evaluation, on or near the surface. Improved accuracy is obtained for points
on the surface. The treecode of Duan and Krasny for Ewald summation is used to
perform sums. Numerical examples are presented with a variety of surfaces.Comment: to appear in Commun. Comput. Phy
Integration over curves and surfaces defined by the closest point mapping
We propose a new formulation for integrating over smooth curves and surfaces
that are described by their closest point mappings. Our method is designed for
curves and surfaces that are not defined by any explicit parameterization and
is intended to be used in combination with level set techniques. However,
contrary to the common practice with level set methods, the volume integrals
derived from our formulation coincide exactly with the surface or line
integrals that one wish to compute. We study various aspects of this
formulation and provide a geometric interpretation of this formulation in terms
of the singular values of the Jacobian matrix of the closest point mapping.
Additionally, we extend the formulation - initially derived to integrate over
manifolds of codimension one - to include integration along curves in three
dimensions. Some numerical examples using very simple discretizations are
presented to demonstrate the efficacy of the formulation.Comment: Revised the pape
A High-Order Kernel Method for Diffusion and Reaction-Diffusion Equations on Surfaces
In this paper we present a high-order kernel method for numerically solving
diffusion and reaction-diffusion partial differential equations (PDEs) on
smooth, closed surfaces embedded in . For two-dimensional
surfaces embedded in , these types of problems have received
growing interest in biology, chemistry, and computer graphics to model such
things as diffusion of chemicals on biological cells or membranes, pattern
formations in biology, nonlinear chemical oscillators in excitable media, and
texture mappings. Our kernel method is based on radial basis functions (RBFs)
and uses a semi-discrete approach (or the method-of-lines) in which the surface
derivative operators that appear in the PDEs are approximated using
collocation. The method only requires nodes at "scattered" locations on the
surface and the corresponding normal vectors to the surface. Additionally, it
does not rely on any surface-based metrics and avoids any intrinsic coordinate
systems, and thus does not suffer from any coordinate distortions or
singularities. We provide error estimates for the kernel-based approximate
surface derivative operators and numerically study the accuracy and stability
of the method. Applications to different non-linear systems of PDEs that arise
in biology and chemistry are also presented
Twist operators in higher dimensions
We study twist operators in higher dimensional CFT's. In particular, we
express their conformal dimension in terms of the energy density for the CFT in
a particular thermal ensemble. We construct an expansion of the conformal
dimension in power series around n=1, with n being replica parameter. We show
that the coefficients in this expansion are determined by higher point
correlations of the energy-momentum tensor. In particular, the first and second
terms, i.e. the first and second derivatives of the scaling dimension, have a
simple universal form. We test these results using holography and free field
theory computations, finding agreement in both cases. We also consider the
`operator product expansion' of spherical twist operators and finally, we
examine the behaviour of correlators of twist operators with other operators in
the limit n ->1.Comment: 44 pages, 2 figure
A high-order Nystrom discretization scheme for boundary integral equations defined on rotationally symmetric surfaces
A scheme for rapidly and accurately computing solutions to boundary integral
equations (BIEs) on rotationally symmetric surfaces in R^3 is presented. The
scheme uses the Fourier transform to reduce the original BIE defined on a
surface to a sequence of BIEs defined on a generating curve for the surface. It
can handle loads that are not necessarily rotationally symmetric. Nystrom
discretization is used to discretize the BIEs on the generating curve. The
quadrature is a high-order Gaussian rule that is modified near the diagonal to
retain high-order accuracy for singular kernels. The reduction in
dimensionality, along with the use of high-order accurate quadratures, leads to
small linear systems that can be inverted directly via, e.g., Gaussian
elimination. This makes the scheme particularly fast in environments involving
multiple right hand sides. It is demonstrated that for BIEs associated with the
Laplace and Helmholtz equations, the kernel in the reduced equations can be
evaluated very rapidly by exploiting recursion relations for Legendre
functions. Numerical examples illustrate the performance of the scheme; in
particular, it is demonstrated that for a BIE associated with Laplace's
equation on a surface discretized using 320,800 points, the set-up phase of the
algorithm takes 1 minute on a standard laptop, and then solves can be executed
in 0.5 seconds.Comment: arXiv admin note: substantial text overlap with
arXiv:1012.56301002.200
A boundary element regularised Stokeslet method applied to cilia and flagella-driven flow
A boundary element implementation of the regularised Stokeslet method of
Cortez is applied to cilia and flagella-driven flows in biology.
Previously-published approaches implicitly combine the force discretisation and
the numerical quadrature used to evaluate boundary integrals. By contrast, a
boundary element method can be implemented by discretising the force using
basis functions, and calculating integrals using accurate numerical or analytic
integration. This substantially weakens the coupling of the mesh size for the
force and the regularisation parameter, and greatly reduces the number of
degrees of freedom required. When modelling a cilium or flagellum as a
one-dimensional filament, the regularisation parameter can be considered a
proxy for the body radius, as opposed to being a parameter used to minimise
numerical errors. Modelling a patch of cilia, it is found that: (1) For a fixed
number of cilia, reducing cilia spacing reduces transport. (2) For fixed patch
dimension, increasing cilia number increases the transport, up to a plateau at
cilia. Modelling a choanoflagellate cell it is found that the
presence of a lorica structure significantly affects transport and flow outside
the lorica, but does not significantly alter the force experienced by the
flagellum.Comment: 20 pages, 7 figures, postprin
Comments on the Sign and Other Aspects of Semiclassical Casimir Energies
The Casimir energy of a massless scalar field is semiclassically given by
contributions due to classical periodic rays. The required subtractions in the
spectral density are determined explicitly. The so defined semiclassical
Casimir energy coincides with that obtained using zeta function regularization
in the cases studied. Poles in the analytic continuation of zeta function
regularization are related to non-universal subtractions in the spectral
density. The sign of the Casimir energy of a scalar field on a smooth manifold
is estimated by the sign of the contribution due to the shortest periodic rays
only. Demanding continuity of the Casimir energy under small deformations of
the manifold, the method is extended to integrable systems. The Casimir energy
of a massless scalar field on a manifold with boundaries includes contributions
due to periodic rays that lie entirely within the boundaries. These
contributions in general depend on the boundary conditions. Although the
Casimir energy due to a massless scalar field may be sensitive to the physical
dimensions of manifolds with boundary, its sign can in favorable cases be
inferred without explicit calculation of the Casimir energy.Comment: 39 pages, no figures, references added, some correction
- …