1,229 research outputs found
On the motion of hairy black holes in Einstein-Maxwell-dilaton theories
Starting from the static, spherically symmetric black hole solutions in
massless Einstein-Maxwell-dilaton (EMD) theories, we build a "skeleton" action,
that is, we phenomenologically replace black holes by an appropriate effective
point particle action, which is well suited to the formal treatment of the
many-body problem in EMD theories. We find that, depending crucially on the
value of their scalar cosmological environment, black holes can undergo steep
"scalarization" transitions, inducing large deviations to the general
relativistic two-body dynamics, as shown, for example, when computing the first
post-Keplerian Lagrangian of EMD theories
A Path-integral for the Master Constraint of Loop Quantum Gravity
In the present paper, we start from the canonical theory of loop quantum
gravity and the master constraint programme. The physical inner product is
expressed by using the group averaging technique for a single self-adjoint
master constraint operator. By the standard technique of skeletonization and
the coherent state path-integral, we derive a path-integral formula from the
group averaging for the master constraint operator. Our derivation in the
present paper suggests there exists a direct link connecting the canonical Loop
quantum gravity with a path-integral quantization or a spin-foam model of
General Relativity.Comment: 19 page
A fast semi-direct least squares algorithm for hierarchically block separable matrices
We present a fast algorithm for linear least squares problems governed by
hierarchically block separable (HBS) matrices. Such matrices are generally
dense but data-sparse and can describe many important operators including those
derived from asymptotically smooth radial kernels that are not too oscillatory.
The algorithm is based on a recursive skeletonization procedure that exposes
this sparsity and solves the dense least squares problem as a larger,
equality-constrained, sparse one. It relies on a sparse QR factorization
coupled with iterative weighted least squares methods. In essence, our scheme
consists of a direct component, comprised of matrix compression and
factorization, followed by an iterative component to enforce certain equality
constraints. At most two iterations are typically required for problems that
are not too ill-conditioned. For an HBS matrix with
having bounded off-diagonal block rank, the algorithm has optimal complexity. If the rank increases with the spatial dimension as is
common for operators that are singular at the origin, then this becomes
in 1D, in 2D, and
in 3D. We illustrate the performance of the method on
both over- and underdetermined systems in a variety of settings, with an
emphasis on radial basis function approximation and efficient updating and
downdating.Comment: 24 pages, 8 figures, 6 tables; to appear in SIAM J. Matrix Anal. App
A Phase Field Model for Continuous Clustering on Vector Fields
A new method for the simplification of flow fields is presented. It is based on continuous clustering. A well-known physical clustering model, the Cahn Hilliard model, which describes phase separation, is modified to reflect the properties of the data to be visualized. Clusters are defined implicitly as connected components of the positivity set of a density function. An evolution equation for this function is obtained as a suitable gradient flow of an underlying anisotropic energy functional. Here, time serves as the scale parameter. The evolution is characterized by a successive coarsening of patterns-the actual clustering-during which the underlying simulation data specifies preferable pattern boundaries. We introduce specific physical quantities in the simulation to control the shape, orientation and distribution of the clusters as a function of the underlying flow field. In addition, the model is expanded, involving elastic effects. In the early stages of the evolution shear layer type representation of the flow field can thereby be generated, whereas, for later stages, the distribution of clusters can be influenced. Furthermore, we incorporate upwind ideas to give the clusters an oriented drop-shaped appearance. Here, we discuss the applicability of this new type of approach mainly for flow fields, where the cluster energy penalizes cross streamline boundaries. However, the method also carries provisions for other fields as well. The clusters can be displayed directly as a flow texture. Alternatively, the clusters can be visualized by iconic representations, which are positioned by using a skeletonization algorithm.
Hierarchical interpolative factorization for elliptic operators: differential equations
This paper introduces the hierarchical interpolative factorization for
elliptic partial differential equations (HIF-DE) in two (2D) and three
dimensions (3D). This factorization takes the form of an approximate
generalized LU/LDL decomposition that facilitates the efficient inversion of
the discretized operator. HIF-DE is based on the multifrontal method but uses
skeletonization on the separator fronts to sparsify the dense frontal matrices
and thus reduce the cost. We conjecture that this strategy yields linear
complexity in 2D and quasilinear complexity in 3D. Estimated linear complexity
in 3D can be achieved by skeletonizing the compressed fronts themselves, which
amounts geometrically to a recursive dimensional reduction scheme. Numerical
experiments support our claims and further demonstrate the performance of our
algorithm as a fast direct solver and preconditioner. MATLAB codes are freely
available.Comment: 37 pages, 13 figures, 12 tables; to appear, Comm. Pure Appl. Math.
arXiv admin note: substantial text overlap with arXiv:1307.266
Hierarchical interpolative factorization for elliptic operators: integral equations
This paper introduces the hierarchical interpolative factorization for
integral equations (HIF-IE) associated with elliptic problems in two and three
dimensions. This factorization takes the form of an approximate generalized LU
decomposition that permits the efficient application of the discretized
operator and its inverse. HIF-IE is based on the recursive skeletonization
algorithm but incorporates a novel combination of two key features: (1) a
matrix factorization framework for sparsifying structured dense matrices and
(2) a recursive dimensional reduction strategy to decrease the cost. Thus,
higher-dimensional problems are effectively mapped to one dimension, and we
conjecture that constructing, applying, and inverting the factorization all
have linear or quasilinear complexity. Numerical experiments support this claim
and further demonstrate the performance of our algorithm as a generalized fast
multipole method, direct solver, and preconditioner. HIF-IE is compatible with
geometric adaptivity and can handle both boundary and volume problems. MATLAB
codes are freely available.Comment: 39 pages, 14 figures, 13 tables; to appear, Comm. Pure Appl. Mat
Combined 3D thinning and greedy algorithm to approximate realistic particles with corrected mechanical properties
The shape of irregular particles has significant influence on micro- and
macro-scopic behavior of granular systems. This paper presents a combined 3D
thinning and greedy set-covering algorithm to approximate realistic particles
with a clump of overlapping spheres for discrete element method (DEM)
simulations. First, the particle medial surface (or surface skeleton), from
which all candidate (maximal inscribed) spheres can be generated, is computed
by the topological 3D thinning. Then, the clump generation procedure is
converted into a greedy set-covering (SCP) problem.
To correct the mass distribution due to highly overlapped spheres inside the
clump, linear programming (LP) is used to adjust the density of each component
sphere, such that the aggregate properties mass, center of mass and inertia
tensor are identical or close enough to the prototypical particle. In order to
find the optimal approximation accuracy (volume coverage: ratio of clump's
volume to the original particle's volume), particle flow of 3 different shapes
in a rotating drum are conducted. It was observed that the dynamic angle of
repose starts to converge for all particle shapes at 85% volume coverage
(spheres per clump < 30), which implies the possible optimal resolution to
capture the mechanical behavior of the system.Comment: 34 pages, 13 figure
- …