No description available.

No description available.

No description available.

No description available.

No description available.

No description available.

No description available.

This thesis deals with the study of rock falls using a mathematical model, codified for computer use, entitled GeoFall. GeoFall, which was developed by the author, allows predictions to be made of rock fall trajectories, run out distances and kinetic energies for a rock of any arbitrary shape. Its main purpose is to assist in the design of remedial works. The mathematical model is based on rigid body mechanics, and analyses a fall in 2D space using a new theory of impact dynamics developed by Brach (1991). The main features and algorithms of the program are presented in this thesis. The performance of GeoFall was evaluated by comparing actual rock fall events described in several published papers with the output created by GeoFall. Also the output from GeoFall has been compared with the output from other rock fall simulation programs used to simulate the documented rockfalls. A new rock slope inventory system entitled the Rock Fall Risk Assessment System (RFRAS) has been developed by the author to determine the rock fall risk at specific rock fall sites. It consists of three phases of inspection, the slope survey, and the preliminary and detailed rating phases. The detailed rating phase uses 13 parameters that when assessed, evaluated and totalled, numerically differentiates slopes from the least to the most hazardous producing an overall rating in the range 21-1926. It not only allows the relative risk of rockfall between slopes to be assessed but it also categorises the rock fall risk and the potential number of future rockfalls. It has been tested on 18 slopes at ten locations in County Durham. The final part of the thesis details a new laboratory based procedure that can be used to determine the coefficients of restitution for any type of rock material. The normal coefficient of restitution has been determined for seven different types of rock, and the tangential coefficient of restitution has been determined for a local sandstone. Some tentative correlations between the normal coefficient of restitution and the rocks physical properties, such as its Unconfined Compressive Strength (UCS) have been presented.

No description available.

Understanding the hydraulic properties of rock fractures is an issue of great importance in fields such as petroleum engineering, groundwater hydrology, and underground waste isolation. Traditionally, models of fluid flow through rock fractures have been based on the cubic law, which asserts that the local value of the fracture transmissivity is equal to /i^/12, where h is the fracture aperture. The local cubic law is mathematically equivalent to assuming that flow in the fracture is governed by the Reynolds lubrication equation. However, this equation is only applicable if the aperture does not change too abruptly, and flow-rates are suitably low. Recent previous investigations showed that the Reynolds equation may over-predict the transmissivity of a fracture by as much as 100%. Other analyses, both theoretical and computational, have estimated that the cubic law will also break down if the Reynolds numbers reach some critical value, variously estimated to be between 1-10. The implication of these results has been that the full three-dimensional, nonlinear Navier-Stokes equations are actually needed to accurately simulate fluid flow in a rock fracture. This conjecture, however, has never been verified, and its verification is the main objective of this research. A surface profilometer was used to measure fracture profiles every 20 nm over the surface of a replica of a fracture in a red Permian sandstone of size 4 cm^. These surface data were used as input to two finite element codes that solve the Navier-Stokes equations and the Reynolds equation, respectively. Numerical simulations of flow through these measured aperture fields were then carried out at different values of the mean aperture, which corresponds to different values of the relative roughness (which is defined as the ratio of the roughness to the mean aperture). At low Reynolds numbers, the Navier-Stokes simulations yielded transmissivities that were 10-100% lower (depending on the relative roughness) than those predicted by the Reynolds simulations, in close agreement with the range of discrepancy between the Reynolds equation and experimental flow measurements reported by previous investigators. At Reynolds numbers << 1, the computed transmissivity is constant. Appreciable deviations from Darcy's law began to be observed when the Reynolds number (defined using the mean aperture as the length scale) exceeded unity. In the regime Re > 20, the computed transmissivities could be fit very well to a Forchheimer-type equation, in which the additional pressure drop varies quadratically with the Reynolds number. The initial deviations from linearity, for Reynolds number around 1, are consistent with the "weak inertia" model of Mei and Auriault (1991). Experimental flow measurements were conducted, on the same fracture replica for which the surface measurements were taken and on which the simulations were conducted. The experiments showed that the measured transmissivities were in close agreement to the values predicted by the Navier-Stokes equations, whereas the Reynolds equation over-estimated fracture transmissivity by as much as 30%. Appreciable deviations from Darcy's law began to be observed when the Reynolds number exceeded unity. When the Reynolds number exceeds about 15, measured transmissivities could be described by the Forchheimer equation, which is consistent with the numerical solutions of the Navier-Stokes equations.