Microlensing in dark matter haloes

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.     

Stochastic computational modelling of highly heterogeneous poroelastic media with long‐range correlations

The compaction of highly heterogeneous poroelastic reservoirs with the geology characterized by long‐range correlations displaying fractal character is investigated within the framework of the stochastic computational modelling. The influence of reservoir heterogeneity upon the magnitude of the stresses induced in the porous matrix during fluid withdrawal and rock consolidation is analysed by performing ensemble averages over realizations of a log‐normally distributed stationary random hydraulic conductivity field. Considering the statistical distribution of this parameter characterized by a coefficient of variation governing the magnitude of heterogeneity and a correlation function which decays with a power‐law scaling behaviour we show that the combination of these two effects result in an increase in the magnitude of effective stresses of the rock during reservoir depletion. Further, within the framework of a perturbation analysis we show that the randomness in the hydraulic conductivity gives rise to non‐linear corrections in the upscaled poroelastic equations. These corrections are illustrated by a self‐consistent recursive hierarchy of solutions of the stochastic poroelastic equations parametrized by a scale parameter representing the fluctuating log‐conductivity standard deviation. A classical example of land subsidence caused by fluid extraction of a weak reservoir is numerically simulated by performing Monte Carlo simulations in conjunction with finite elements discretizations of the poroelastic equations associated with an ensemble of geologies. Numerical results illustrate the effects of the spatial variability and fractal character of the permeability distribution upon the evolution of the Mohr–Coulomb function of the rock. Copyright © 2004 John Wiley & Sons, Ltd.