39,947 research outputs found
Low-l CMB Analysis and Inpainting
Reconstruction of the CMB in the Galactic plane is extremely difficult due to
the dominant foreground emissions such as Dust, Free-Free or Synchrotron. For
cosmological studies, the standard approach consists in masking this area where
the reconstruction is not good enough. This leads to difficulties for the
statistical analysis of the CMB map, especially at very large scales (to study
for e.g., the low quadrupole, ISW, axis of evil, etc). We investigate in this
paper how well some inpainting techniques can recover the low- spherical
harmonic coefficients. We introduce three new inpainting techniques based on
three different kinds of priors: sparsity, energy and isotropy, and we compare
them. We show that two of them, sparsity and energy priors, can lead to
extremely high quality reconstruction, within 1% of the cosmic variance for a
mask with Fsky larger than 80%.Comment: Submitte
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
A Study of Different Modeling Choices For Simulating Platelets Within the Immersed Boundary Method
The Immersed Boundary (IB) method is a widely-used numerical methodology for
the simulation of fluid-structure interaction problems. The IB method utilizes
an Eulerian discretization for the fluid equations of motion while maintaining
a Lagrangian representation of structural objects. Operators are defined for
transmitting information (forces and velocities) between these two
representations. Most IB simulations represent their structures with
piecewise-linear approximations and utilize Hookean spring models to
approximate structural forces. Our specific motivation is the modeling of
platelets in hemodynamic flows. In this paper, we study two alternative
representations - radial basis functions (RBFs) and Fourier-based
(trigonometric polynomials and spherical harmonics) representations - for the
modeling of platelets in two and three dimensions within the IB framework, and
compare our results with the traditional piecewise-linear approximation
methodology. For different representative shapes, we examine the geometric
modeling errors (position and normal vectors), force computation errors, and
computational cost and provide an engineering trade-off strategy for when and
why one might select to employ these different representations.Comment: 33 pages, 17 figures, Accepted (in press) by APNU
Solving Einstein's Equations With Dual Coordinate Frames
A method is introduced for solving Einstein's equations using two distinct
coordinate systems. The coordinate basis vectors associated with one system are
used to project out components of the metric and other fields, in analogy with
the way fields are projected onto an orthonormal tetrad basis. These field
components are then determined as functions of a second independent coordinate
system. The transformation to the second coordinate system can be thought of as
a mapping from the original ``inertial'' coordinate system to the computational
domain. This dual-coordinate method is used to perform stable numerical
evolutions of a black-hole spacetime using the generalized harmonic form of
Einstein's equations in coordinates that rotate with respect to the inertial
frame at infinity; such evolutions are found to be generically unstable using a
single rotating coordinate frame. The dual-coordinate method is also used here
to evolve binary black-hole spacetimes for several orbits. The great
flexibility of this method allows comoving coordinates to be adjusted with a
feedback control system that keeps the excision boundaries of the holes within
their respective apparent horizons.Comment: Updated to agree with published versio
- …