Adaptive computation of gravitational waves from black hole interactions

We construct a class of linear partial differential equations describing general perturbations of non-rotating black holes in 3D Cartesian coordinates. In contrast to the usual approach, a single equation treats all radiative $\ell -m$ modes simultaneously, allowing the study of wave perturbations of black holes with arbitrary 3D structure, as would be present when studying the full set of nonlinear Einstein equations describing a perturbed black hole. This class of equations forms an excellent testbed to explore the computational issues of simulating black spacetimes using a three dimensional adaptive mesh refinement code. Using this code, we present results from the first fully resolved 3D solution of the equations describing perturbed black holes. We discuss both fixed and adaptive mesh refinement, refinement criteria, and the computational savings provided by adaptive techniques in 3D for such model problems of distorted black holes.Comment: 16 Pages, RevTeX, 13 figure

Towards a Landau-Ginzburg-type Theory for Granular Fluids

In this paper we show how, under certain restrictions, the hydrodynamic equations for the freely evolving granular fluid fit within the framework of the time dependent Landau-Ginzburg (LG) models for critical and unstable fluids (e.g. spinodal decomposition). The granular fluid, which is usually modeled as a fluid of inelastic hard spheres (IHS), exhibits two instabilities: the spontaneous formation of vortices and of high density clusters. We suppress the clustering instability by imposing constraints on the system sizes, in order to illustrate how LG-equations can be derived for the order parameter, being the rate of deformation or shear rate tensor, which controls the formation of vortex patterns. From the shape of the energy functional we obtain the stationary patterns in the flow field. Quantitative predictions of this theory for the stationary states agree well with molecular dynamics simulations of a fluid of inelastic hard disks.Comment: 19 pages, LaTeX, 8 figure

Alpha, Betti and the Megaparsec Universe: on the Topology of the Cosmic Web

We study the topology of the Megaparsec Cosmic Web in terms of the scale-dependent Betti numbers, which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic. The richer information content of Betti numbers goes along the availability of fast algorithms to compute them. For continuous density fields, we determine the scale-dependence of Betti numbers by invoking the cosmologically familiar filtration of sublevel or superlevel sets defined by density thresholds. For the discrete galaxy distribution, however, the analysis is based on the alpha shapes of the particles. These simplicial complexes constitute an ordered sequence of nested subsets of the Delaunay tessellation, a filtration defined by the scale parameter, $\alpha$. As they are homotopy equivalent to the sublevel sets of the distance field, they are an excellent tool for assessing the topological structure of a discrete point distribution. In order to develop an intuitive understanding for the behavior of Betti numbers as a function of $\alpha$, and their relation to the morphological patterns in the Cosmic Web, we first study them within the context of simple heuristic Voronoi clustering models. Subsequently, we address the topology of structures emerging in the standard LCDM scenario and in cosmological scenarios with alternative dark energy content. The evolution and scale-dependence of the Betti numbers is shown to reflect the hierarchical evolution of the Cosmic Web and yields a promising measure of cosmological parameters. We also discuss the expected Betti numbers as a function of the density threshold for superlevel sets of a Gaussian random field.Comment: 42 pages, 14 figure

Drag reduction induced by superhydrophobic surfaces in turbulent pipe flow

The drag reduction induced by superhydrophobic surfaces is investigated in a turbulent pipe flow. Wetted superhydrophobic surfaces are shown to trap gas bubbles in their asperities. This stops the liquid from coming in direct contact with the wall in that location, allowing the flow to slip over the air bubbles. We consider a well-defined texture with streamwise grooves at the walls in which the gas is expected to be entrapped. This configuration is modeled with alternating no-slip and shear-free boundary conditions at the wall. With respect to the classical turbulent pipe flow, a substantial drag reduction is observed which strongly depends on the grooves’ dimension and on the solid fraction, i.e., the ratio between the solid wall surface and the total surface of the pipe’s circumference. The drag reduction is due to the mean slip velocity at the wall which increases the flow rate at a fixed pressure drop. The enforced boundary conditions also produce peculiar turbulent structures which on the contrary decrease the flow rate. The two concurrent effects provide an overall flow rate increase as demonstrated by means of the mean axial momentum balance. This equation provides the balance between the mean pressure gradient, the Reynolds stress, the mean flow rate, and the mean slip velocity contribution

Classical and quantum regimes of two-dimensional turbulence in trapped Bose-Einstein condensates

We investigate two-dimensional turbulence in finite-temperature trapped Bose-Einstein condensates within damped Gross-Pitaevskii theory. Turbulence is produced via circular motion of a Gaussian potential barrier stirring the condensate. We systematically explore a range of stirring parameters and identify three regimes, characterized by the injection of distinct quantum vortex structures into the condensate: (A) periodic vortex dipole injection, (B) irregular injection of a mixture of vortex dipoles and co-rotating vortex clusters, and (C) continuous injection of oblique solitons that decay into vortex dipoles. Spectral analysis of the kinetic energy associated with vortices reveals that regime (B) can intermittently exhibit a Kolmogorov $k^{-5/3}$ power law over almost a decade of length or wavenumber ($k$) scales. The kinetic energy spectrum of regime (C) exhibits a clear $k^{-3/2}$ power law associated with an inertial range for weak-wave turbulence, and a $k^{-7/2}$ power law for high wavenumbers. We thus identify distinct regimes of forcing for generating either two-dimensional quantum turbulence or classical weak-wave turbulence that may be realizable experimentally.Comment: 11 pages, 10 figures. Minor updates to text and figures 1, 2 and