Revisiting Event Horizon Finders

Event horizons are the defining physical features of black hole spacetimes, and are of considerable interest in studying black hole dynamics. Here, we reconsider three techniques to localise event horizons in numerical spacetimes: integrating geodesics, integrating a surface, and integrating a level-set of surfaces over a volume. We implement the first two techniques and find that straightforward integration of geodesics backward in time to be most robust. We find that the exponential rate of approach of a null surface towards the event horizon of a spinning black hole equals the surface gravity of the black hole. In head-on mergers we are able to track quasi-normal ringing of the merged black hole through seven oscillations, covering a dynamic range of about 10^5. Both at late times (when the final black hole has settled down) and at early times (before the merger), the apparent horizon is found to be an excellent approximation of the event horizon. In the head-on binary black hole merger, only {\em some} of the future null generators of the horizon are found to start from past null infinity; the others approach the event horizons of the individual black holes at times far before merger.Comment: 30 pages, 15 figures, revision

A comparison of 3D particle, fluid and hybrid simulations for negative streamers

In the high field region at the head of a discharge streamer, the electron energy distribution develops a long tail. In negative streamers, these electrons can run away and contribute to energetic processes such as terrestrial gamma-ray and electron flashes. Moreover, electron density fluctuations can accelerate streamer branching. To track energies and locations of single electrons in relevant regions, we have developed a 3D hybrid model that couples a particle model in the region of high fields and low electron densities with a fluid model in the rest of the domain. Here we validate our 3D hybrid model on a 3D (super-)particle model for negative streamers in overvolted gaps, and we show that it almost reaches the computational efficiency of a 3D fluid model. We also show that the extended fluid model approximates the particle and the hybrid model well until stochastic fluctuations become important, while the classical fluid model underestimates velocities and ionization densities. We compare density fluctuations and the onset of branching between the models, and we compare the front velocities with an analytical approximation

Hydroxyl vacancies in single-walled aluminosilicate and aluminogermanate nanotubes

We report the first theoretical study of hydroxyl vacancies in aluminosilicate and aluminogermanate single-walled metal-oxide nanotubes. The defects are modeled on both sides of the tube walls and lead to occupied and empty states in the band gap which are highly localized both in energy and in real space. We find different magnetization states depending on both the chemical composition and the specific side with respect to the tube cavity. The defect-induced perturbations to the pristine electronic structure are related to the electrostatic polarization across the tube walls and the ensuing change in Br{\o}nsted acid-base reactivity. Finally, the capacity to counterbalance local charge accumulations, a characteristic feature of these systems, is discussed in view of their potential application as insulating coatings for one-dimensional conducting nanodevices.Comment: manuscript: 4 pages, 4 figure

A computational study of the effect of windscreen shape and flow resistivity on turbulent wind noise reduction

This is the published version. Copyright 2011 Acoustical Society of AmericaIn this paper, numerical simulations are used to study the turbulentwind noise reduction effect of microphone windscreens with varying shapes and flow resistivities. Typical windscreen shapes consisting of circular, elliptical, and rectangular cylinders are investigated. A turbulent environment is generated by placing a solid circular cylinder upstream of the microphone. An immersed-boundary method with a fifth-order weighted essentially non-oscillatory scheme is implemented to enhance the simulation accuracy for high-Reynolds number flow around the solid cylinder as well as at the interface between the open air and the porous material comprising the windscreen. The Navier–Stokes equations for incompressible flow are solved in the open air. For the flow inside the porous material, a modified form of the Zwikker–Kosten equation is solved. The results show that, on average, the circular and horizontal ellipse windscreens have similar overall wind noise reduction performance, while the horizontal ellipse windscreen with medium flow resistivity provides the most effective wind noise reduction among all the considered cases. The vertical ellipse windscreen with high flow resistivity, in particular, increases the wind noise because of increased self-generation of turbulence