The lattice solid model to simulate the physics of rocks and earthquakes: Incorporation of friction

The particle-based lattice solid model developed to study the physics of rocks and the nonlinear dynamics of earthquakes is refined by incorporating intrinsic friction between particles. The model provides a means for studying the causes of seismic wave attenuation, as well as frictional heat generation, fault zone evolution, and localisation phenomena. A modified velocity-Verlat scheme that allows friction to be precisely modelled is developed. This is a difficult computational problem given that a discontinuity must be accurately simulated by the numerical approach (i.e., the transition from static to dynamical frictional behaviour). This is achieved using a half time step integration scheme. At each half time step, a nonlinear system is solved to compute the static frictional forces and states of touching particle-pairs. Improved efficiency is achieved by adaptively adjusting the time step increment, depending on the particle velocities in the system. The total energy is calculated and verified to remain constant to a high precision during simulations. Numerical experiments show that the model can be applied to the study of earthquake dynamics, the stick-slip instability, heat generation, and fault zone evolution. Such experiments may lead to a conclusive resolution of the heat flow paradox and improved understanding of earthquake precursory phenomena and dynamics. (C) 1999 Academic Press.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

Discrete element modelling of the influence of lifters on power draw of tumbling mills

Crushing and grinding are the most energy intensive part of the mineral recovery process. A major part of rock size reduction occurs in tumbling mills. Empirical models for the power draw of tumbling mills do not consider the effect of lifters. Discrete element modelling was used to investigate the effect of lifter condition on the power draw of tumbling mill. Results obtained with PFC3D code show that lifter condition will have a significant influence on the power draw and on the mode of energy consumption in the mill. Relatively high lifters will consume less power than low lifters, under otherwise identical conditions. The fraction of the power that will be consumed as friction will increase as the height of the lifters decreases. This will result in less power being used for high intensity comminution caused by the impacts. The fraction of the power that will be used to overcome frictional resistance is determined by the material's coefficient of friction. Based on the modelled results, it appears that the effective coefficient of friction for in situ mill is close to 0.1. (C) 2003 Elsevier Science Ltd. All rights reserved.

Dynamic rupture in a 3-D particle-based simulation of a rough planar fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.

Transducer for direct measurement of skin friction in hypervelocity impulse facilities

An acceleration compensated transducer was developed to enable the direct measurement of skin friction in hypervelocity impulse facilities. The gauge incorporated a measurement and acceleration element that employed direct shear of a piezoelectric ceramic. The design integrated techniques to maximize rise time and shear response while minimizing the affects of acceleration, pressure, heat transfer, and electrical interference. The arrangement resulted in a transducer natural frequency near 40 kHz. The transducer was calibrated for shear and acceleration in separate bench tests and was calibrated for pressure within an impulse facility. Uncertainty analyses identified only small experimental errors in the shear and acceleration calibration techniques. Although significant errors were revealed in the method of pressure calibration, total skin-friction measurement errors as low as +/-7-12% were established. The transducer was successfully utilized in a shock tunnel, and sample measurements are presented for flow conditions that simulate a flight Mach number near 8.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Friction factors and rheological behavior of Australian honey in a straight pipe

The pumping characteristics of four Australian honey samples were investigated in a straight pipe. Six flow rates (100-500 kg h(-1)) were studied at three temperatures (35-50degreesC). The pressure loss increased with an increase in the length of the pipe, as the low rate was increased and as the temperature was reduced. In the 25.4 mm-pipe, the Reynolds number ranged from 0.2-32.0 and are substantially less than the critica value (2040-2180) for laminar condition in the system. The relationship between the wall shear stress and shear rate approximated power-law behaviour, and the power-law index was not significantly (p>0.05) different from 1.0. The honey samples exhibited Newtonian behaviour at all the temperatures and this was confirmed by rheometric studies using Couette geometry. A friction chart was generated independent of temperature and the type of honey. An equation was developed to predict the pressure loss of the honey in a typical pipeline at any temperature once the viscosity-temperature relationship had been established.

Control of hypersonic turbulent skin friction by boundary-layer combustion of hydrogen

Shvab-Zeldovich coupling of flow variables has been used to extend Van Driest's theory of turbulent boundary-layer skin friction to include injection and combustion of hydrogen in the boundary layer. The resulting theory is used to make predictions of skin friction and heat transfer that are found to be consistent with experimental and numerical results. Using the theory to extrapolate to larger downstream distances at the same experimental conditions, it is found that the reduction in skin-friction drag with hydrogen mixing and combustion is three times that with mixing alone. In application to flow on a flat plate at mainstream velocities of 2, 4, and 6 knits, and Reynolds numbers from 3 X 10(6) to 1 x 10(8), injection and combustion of hydrogen yielded values of skin-friction drag that were less than one-half of the no-injection skin-friction drag, together with a net reduction in heat transfer when the combustion heat release in air was less than the stagnation enthalpy. The mass efficiency of hydrogen injection, as measured by effective specific impulse values, was approximately 2000 s.

A unified friction description and its application to the simulation of frictional instability using the finite element method

This article first summarizes some available experimental results on the frictional behaviour of contact interfaces, and briefly recalls typical frictional experiments and relationships, which are applicable for rock mechanics, and then a unified description is obtained to describe the entire frictional behaviour. It is formulated based on the experimental results and applied with a stick and slip decomposition algorithm to describe the stick-slip instability phenomena, which can describe the effects observed in rock experiments without using the so-called state variable, thus avoiding related numerical difficulties. This has been implemented to our finite element code, which uses the node-to-point contact element strategy proposed by the authors to handle the frictional contact between multiple finite-deformation bodies with stick and finite frictional slip, and applied here to simulate the frictional behaviour of rocks to show its usefulness and efficiency.

Analytical solution of forced convection in a duct of rectangular cross-section saturated by a porous medium

A theoretical analysis is presented to investigate fully developed (both thermally and hydrodynamically) forced convection in a duct of rectangular cross-section filled with a hyper-porous medium. The Darcy-Brinkman model for flow through porous media was adopted in the present analysis. A Fourier series type solution is applied to obtain the exact velocity and temperature distribution within the duct. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford [1], is treated. Values of the Nusselt number and the friction factor as a function of the aspect ratio, the Darcy number, and the viscosity ratio are reported.

Effect of sodium dodecylbenzene sulfonate on the motion of three-phase contact lines on the Wilhelmy plate surface

The combined approach of the molecular-kinetic and hydrodynamic theories for description of the motion of three-phase gas-liquid-solid contact lines has been examined using the Wilhelmy plate method. The whole dynamic meniscus has been divided into molecular, hydrodynamic, and static-like regions. The Young-Laplace equation and the molecular-kinetic and hydrodynamic dewetting theories have been applied to describe the meniscus profiles and contact angle. The dissipative forces accompanying the dynamic dewetting have also been investigated. The experiments with a Wilhelmy plate made from an acrylic polymer sheet were carried out using a computerized apparatus for contact angle analysis (OCA 20, DataPhysics, Germany). The extrapolated dynamic contact angle versus velocity of the three-phase contact line for Milli-Q water and 5 x 10(-4) M SDBS solution was experimentally obtained and compared with the combined MHD models with low and moderate Reynolds numbers. The models predict similar results for the extrapolated contact angle. SDBS decreases the equilibrium contact angle and increases the molecular jumping length but does not affect the molecular frequency significantly. The hydrodynamic deformation of the meniscus, viscous dissipation, and friction were also influenced by the SDBS surfactant. (c) 2005 Elsevier Inc. All rights reserved.

Effect of hardness on three very different forms of wear

Different abrasive wear tests have been applied to materials with hardnesses ranging from 80 HV (aluminium) to 1700 HV (tungsten carbide). The tests were: dry sand rubber wheel (DSRbrW); a similar test using a steel wheel (DSStlW); a new combined impact-abrasion test (FIA). The DSRbrW results were as expected, giving generally decreasing wear with increasing hardness. White cast irons and tool steels containing coarse, hard carbide particles performed better than more homogeneous materials of comparable hardness. When normalized to load and distance, the DSStlW results for the homogeneous materials were similar to the DSRbrW results. The multi-phase materials performed poorly in the DSStlW test, with volume loss for high-speed steel (880 HV) higher than that of aluminium. Within this group, wear increased with increasing hardness. These unexpected results are explained in terms of (a) differential friction coefficients of wheel and specimen, (b) increased fracture of sand, and (c) introduction of microfracture wear mechanisms. The FIA combined impact-abrasion results lacked clear correlations with hardness. The span of relative wear rates was similar to that reported for materials in ball mills. White cast irons at maximum hardness performed fairly poorly and showed evidence of microfracture. (C) 1997 Elsevier Science S.A.