Predicting cave initiation and propagation

The most widely used method for predicting the onset of continuous caving is Laubscher's caving chart. A detailed examination of this method was undertaken which concluded that it had limitations which may impact on results, particularly when dealing with stronger rock masses that are outside current experience. These limitations relate to inadequate guidelines for adjustment factors to rock mass rating (RMR), concerns about the position on the chart of critical case history data, undocumented changes to the method and an inadequate number of data points to be confident of stability boundaries. A review was undertaken on the application and reliability of a numerical method of assessing cavability. The review highlighted a number of issues, which at this stage, make numerical continuum methods problematic for predicting cavability. This is in particular reference to sensitivity to input parameters that are difficult to determine accurately and mesh dependency. An extended version of the Mathews method for open stope design was developed as an alternative method of predicting the onset of continuous caving. A number of caving case histories were collected and analyzed and a caving boundary delineated statistically on the Mathews stability graph. The definition of the caving boundary was aided by the existence of a large and wide-ranging stability database from non-caving mines. A caving rate model was extrapolated from the extended Mathews stability graph but could only be partially validated due to a lack of reliable data.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Water table waves in an unconfined aquifer: Experiments and modeling

[1] Comprehensive measurements are presented of the piezometric head in an unconfined aquifer during steady, simple harmonic oscillations driven by a hydrostatic clear water reservoir through a vertical interface. The results are analyzed and used to test existing hydrostatic and nonhydrostatic, small-amplitude theories along with capillary fringe effects. As expected, the amplitude of the water table wave decays exponentially. However, the decay rates and phase lags indicate the influence of both vertical flow and capillary effects. The capillary effects are reconciled with observations of water table oscillations in a sand column with the same sand. The effects of vertical flows and the corresponding nonhydrostatic pressure are reasonably well described by small-amplitude theory for water table waves in finite depth aquifers. That includes the oscillation amplitudes being greater at the bottom than at the top and the phase lead of the bottom compared with the top. The main problems with respect to interpreting the measurements through existing theory relate to the complicated boundary condition at the interface between the driving head reservoir and the aquifer. That is, the small-amplitude, finite depth expansion solution, which matches a hydrostatic boundary condition between the bottom and the mean driving head level, is unrealistic with respect to the pressure variation above this level. Hence it cannot describe the finer details of the multiple mode behavior close to the driving head boundary. The mean water table height initially increases with distance from the forcing boundary but then decreases again, and its asymptotic value is considerably smaller than that previously predicted for finite depth aquifers without capillary effects. Just as the mean water table over-height is smaller than predicted by capillarity-free shallow aquifer models, so is the amplitude of the second harmonic. In fact, there is no indication of extra second harmonics ( in addition to that contained in the driving head) being generated at the interface or in the interior.

Selective entrainment of sediment graded by size and density under waves

The influence of near-bed sorting processes on heavy mineral content in suspension is discussed. Sediment concentrations above a rippled bed of mixed quartz and heavy mineral sand were measured under regular nonbreaking waves in the laboratory. Using the traditional gradient diffusion process, settling velocity would be expected to strongly affect sediment distribution. This was not observed during present trials. In fact, the vertical gradients of time-averaged suspension concentrations were found to be similar for the light and heavy minerals, despite their different settling velocities. This behavior implies a convective rather than diffusive distribution mechanism. Between the nonmoving bed and the lowest suspension sampling point, fight and heavy mineral concentration differs by two orders of magnitude. This discrimination against the heavy minerals in the pickup process is due largely to selective entrainment at the ripple face. Bed-form dynamics and the nature of quartz suspension profiles are found to be little affected by the trialed proportion of overall heavy minerals in the bed (3.8-22.1%).

Simultaneous FT-NIR and ESR analyses to yield propagation rate coefficients for polymerization of methyl methacrylate based monomers

A hyphenated instrumental approach has been used to obtain reliable values for the propagation rate coefficients as a function of conversion for polymerizations of methyl methacrylate (MMA) and a mixture of MMA and ethyleneglycol dimethacrylate (EGDMA) with a 1:1 concentration of double bonds, from near the onset of the Trommsdorf region into the glass region. ESR spectroscopy was used to measure the radical concentration while FT-NIR fibre-optic spectroscopy was employed to measure instantaneously the double-bond concentration within the temperature-controlled cavity of the ESR instrument during polymerization. The advantage of this approach to the measurement of the rate coefficient is that it is equally applicable to branching and linear polymerizations. For the polymerization of methyl methacrylate, the values of the rate coefficient at the lowest conversions at which reliable values could be obtained were in agreement with recently reported values obtained by the PLP-SEC method. For the lowest conversions, the values obtained were 403 1 mol(-1) s(-1) at 306 K for MMA and 5201 mol(-1) s(-1) at 310 K for a 1:1 mixture of MMA and EGDMA. (C) 2003 Society of Chemical Industry.