Instability of Extremal Relativistic Charged Spheres

With the question, ``Can relativistic charged spheres form extremal black holes?" in mind, we investigate the properties of such spheres from a classical point of view. The investigation is carried out numerically by integrating the Oppenheimer-Volkov equation for relativistic charged fluid spheres and finding interior Reissner-Nordstr\"om solutions for these objects. We consider both constant density and adiabatic equations of state, as well as several possible charge distributions, and examine stability by both a normal mode and an energy analysis. In all cases, the stability limit for these spheres lies between the extremal ($Q = M$) limit and the black hole limit ($R = R_+$). That is, we find that charged spheres undergo gravitational collapse before they reach $Q = M$, suggesting that extremal Reissner-Nordtr\"om black holes produced by collapse are ruled out. A general proof of this statement would support a strong form of the cosmic censorship hypothesis, excluding not only stable naked singularities, but stable extremal black holes. The numerical results also indicate that although the interior mass-energy $m(R)$ obeys the usual $m/R < 4/9$ stability limit for the Schwarzschild interior solution, the gravitational mass $M$ does not. Indeed, the stability limit approaches $R_+$ as $Q \to M$. In the Appendix we also argue that Hawking radiation will not lead to an extremal Reissner-Nordstr\"om black hole. All our results are consistent with the third law of black hole dynamics, as currently understood

A Phase Space Approach to Gravitational Enropy

We examine the definition S = ln Omega as a candidate "gravitational entropy" function. We calculate its behavior for gravitationl and density perturbations in closed, open and flat cosmologies and find that in all cases it increases monotonically. Using the formalism to calculate the gravitational entropy produced during inflation gives the canonical answer. We compare the behavior of S with the behavior of the square of the Weyl tensor. Applying the formalism to black holes has proven more problematical.Comment: Talk delivered at South African Relativistic Cosmology Symposium, Feb 1999. Some new results over Rothman and Anninos 97. To appear in GRG, 17 page

A Two-Threshold Model for Scaling Laws of Non-Interacting Snow Avalanches

The sizes of snow slab failure that trigger snow avalanches are power-law distributed. Such a power-law probability distribution function has also been proposed to characterize different landslide types. In order to understand this scaling for gravity driven systems, we introduce a two-threshold 2-d cellular automaton, in which failure occurs irreversibly. Taking snow slab avalanches as a model system, we find that the sizes of the largest avalanches just preceeding the lattice system breakdown are power law distributed. By tuning the maximum value of the ratio of the two failure thresholds our model reproduces the range of power law exponents observed for land-, rock- or snow avalanches. We suggest this control parameter represents the material cohesion anisotropy.Comment: accepted PR

Two-dimensional hydrodynamic lattice-gas simulations of binary immiscible and ternary amphiphilic fluid flow through porous media

The behaviour of two dimensional binary and ternary amphiphilic fluids under flow conditions is investigated using a hydrodynamic lattice gas model. After the validation of the model in simple cases (Poiseuille flow, Darcy's law for single component fluids), attention is focussed on the properties of binary immiscible fluids in porous media. An extension of Darcy's law which explicitly admits a viscous coupling between the fluids is verified, and evidence of capillary effects are described. The influence of a third component, namely surfactant, is studied in the same context. Invasion simulations have also been performed. The effect of the applied force on the invasion process is reported. As the forcing level increases, the invasion process becomes faster and the residual oil saturation decreases. The introduction of surfactant in the invading phase during imbibition produces new phenomena, including emulsification and micellisation. At very low fluid forcing levels, this leads to the production of a low-resistance gel, which then slows down the progress of the invading fluid. At long times (beyond the water percolation threshold), the concentration of remaining oil within the porous medium is lowered by the action of surfactant, thus enhancing oil recovery. On the other hand, the introduction of surfactant in the invading phase during drainage simulations slows down the invasion process -- the invading fluid takes a more tortuous path to invade the porous medium -- and reduces the oil recovery (the residual oil saturation increases).Comment: 48 pages, 26 figures. Phys. Rev. E (in press

Simulating Three-Dimensional Hydrodynamics on a Cellular-Automata Machine

We demonstrate how three-dimensional fluid flow simulations can be carried out on the Cellular Automata Machine 8 (CAM-8), a special-purpose computer for cellular-automata computations. The principal algorithmic innovation is the use of a lattice-gas model with a 16-bit collision operator that is specially adapted to the machine architecture. It is shown how the collision rules can be optimized to obtain a low viscosity of the fluid. Predictions of the viscosity based on a Boltzmann approximation agree well with measurements of the viscosity made on CAM-8. Several test simulations of flows in simple geometries -- channels, pipes, and a cubic array of spheres -- are carried out. Measurements of average flux in these geometries compare well with theoretical predictions.Comment: 19 pages, REVTeX and epsf macros require

Leukemia-related chromosomal loss detected in hematopoietic progenitor cells of benzene-exposed workers.

Benzene exposure causes acute myeloid leukemia and hematotoxicity, shown as suppression of mature blood and myeloid progenitor cell numbers. As the leukemia-related aneuploidies monosomy 7 and trisomy 8 previously had been detected in the mature peripheral blood cells of exposed workers, we hypothesized that benzene could cause leukemia through the induction of these aneuploidies in hematopoietic stem and progenitor cells. We measured loss and gain of chromosomes 7 and 8 by fluorescence in situ hybridization in interphase colony-forming unit-granulocyte-macrophage (CFU-GM) cells cultured from otherwise healthy benzene-exposed (n=28) and unexposed (n=14) workers. CFU-GM monosomy 7 and 8 levels (but not trisomy) were significantly increased in subjects exposed to benzene overall, compared with levels in the control subjects (P=0.0055 and P=0.0034, respectively). Levels of monosomy 7 and 8 were significantly increased in subjects exposed to &lt;10 p.p.m. (20%, P=0.0419 and 28%, P=0.0056, respectively) and ≥ 10 p.p.m. (48%, P=0.0045 and 32%, 0.0354) benzene, compared with controls, and significant exposure-response trends were detected (P(trend)=0.0033 and 0.0057). These data show that monosomies 7 and 8 are produced in a dose-dependent manner in the blood progenitor cells of workers exposed to benzene, and may be mechanistically relevant biomarkers of early effect for benzene and other leukemogens

Fisher's arrow of `time' in cosmological coherent phase space

Fisher's arrow of `time' in a cosmological phase space defined as in quantum optics (i.e., whose points are coherent states) is introduced as follows. Assuming that the phase space evolution of the universe starts from an initial squeezed cosmological state towards a final thermal one, a Fokker-Planck equation for the time-dependent, cosmological Q phase space probability distribution can be written down. Next, using some recent results in the literature, we derive an information arrow of time for the Fisher phase space cosmological entropy based on the Q function. We also mention the application of Fisher's arrow of time to stochastic inflation modelsComment: 10 pages, LaTex, Honorable Mention at GRF-199

Computer simulations of domain growth and phase separation in two-dimensional binary immiscible fluids using dissipative particle dynamics

We investigate the dynamical behavior of binary fluid systems in two dimensions using dissipative particle dynamics. We find that following a symmetric quench the domain size R(t) grows with time t according to two distinct algebraic laws R(t) = t^n: at early times n = 1/2, while for later times n = 2/3. Following an asymmetric quench we observe only n = 1/2, and if momentum conservation is violated we see n = 1/3 at early times. Bubble simulations confirm the existence of a finite surface tension and the validity of Laplace's law. Our results are compared with similar simulations which have been performed previously using molecular dynamics, lattice-gas and lattice-Boltzmann automata, and Langevin dynamics. We conclude that dissipative particle dynamics is a promising method for simulating fluid properties in such systems.Comment: RevTeX; 22 pages, 5 low-resolution figures. For full-resolution figures, connect to http://www.tcm.phy.cam.ac.uk/~ken21/tension/tension.htm

Fluctuations of elastic interfaces in fluids: Theory and simulation

We study the dynamics of elastic interfaces-membranes-immersed in thermally excited fluids. The work contains three components: the development of a numerical method, a purely theoretical approach, and numerical simulation. In developing a numerical method, we first discuss the dynamical coupling between the interface and the surrounding fluids. An argument is then presented that generalizes the single-relaxation time lattice-Boltzmann method for the simulation of hydrodynamic interfaces to include the elastic properties of the boundary. The implementation of the new method is outlined and it is tested by simulating the static behavior of spherical bubbles and the dynamics of bending waves. By means of the fluctuation-dissipation theorem we recover analytically the equilibrium frequency power spectrum of thermally fluctuating membranes and the correlation function of the excitations. Also, the non-equilibrium scaling properties of the membrane roughening are deduced, leading us to formulate a scaling law describing the interface growth, W^2(L,T)=L^3 g[t/L^(5/2)], where W, L and T are the width of the interface, the linear size of the system and the temperature respectively, and g is a scaling function. Finally, the phenomenology of thermally fluctuating membranes is simulated and the frequency power spectrum is recovered, confirming the decay of the correlation function of the fluctuations. As a further numerical study of fluctuating elastic interfaces, the non-equilibrium regime is reproduced by initializing the system as an interface immersed in thermally pre-excited fluids.Comment: 15 pages, 11 figure