Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

Cyclic extrusion of a lava dome based on a stick-slip mechanism

Lava dome eruptions are sometimes characterised by large periodic fluctuations in extrusion rate over periods of hours that may be accompanied by Vulcanian explosions and pyroclastic flows. We consider a simple system of nonlinear equations describing a 1D flow of lava extrusion through a deep elastic dyke feeding a shallower cylindrical conduit in order to simulate this short-period cyclicity. Stick-slip conditions depending on a critical shear stress are assumed at the wall boundary of the cylindrical conduit. By analogy with the behaviour of industrial polymers in a plastic extruder, the elastic dyke acts like a barrel and the shallower cylindrical portion of the conduit as a die for the flow of magma acting as a polymer. When we applied the model to the Soufrière Hills Volcano, Montserrat, for which the key parameters have been evaluated from previous studies, cyclic extrusions with periods from 3 to 30 h were readily simulated, matching observations. The model also reproduces the reduced period of cycles observed when a major unloading event occurs due to lava dome collapse.

On stick-slip homoclinic chaos and bifurcations in a mechanical system with dry friction

In this paper we consider a self-excited mechanical system by dry friction in order to study the bifurcational behavior of the arisen vibrations. The oscillating system consists of a mass block-belt-system which is self-excited by static and Coulomb friction. We analyze the system behavior numerically through bifurcation diagrams, phase portraits, frequency spectra and Poincare maps, which show the existence of nonhomoclinic and homoclinic chaos and a route to homoclinic chaos. The homoclinic chaos is also analyzed analytically via the Melnikov prediction method. The system dynamic is characterized by the existence of two potential wells in the phase plane which exhibit rich bifurcational and chaotic behavior.

Stick-slip chaos in a non-ideal self-excited system

In engineering practical systems the excitation source is generally dependent on the system dynamic structure. In this paper we analyze a self-excited oscillating system due to dry friction which interacts with an energy source of limited power supply (non ideal problem). The mechanical system consists of an oscillating system sliding on a moving belt driven by a limited power supply. In the oscillating system considered here, dry friction acts as an excitation mechanism for stick-slip oscillations. The stick-slip chaotic oscillations are investigated because the knowledge of their dynamic characteristics is an important step in system design and control. Many engineering systems present stick-slip chaotic oscillations such as machine tools, oil well drillstrings, car brakes and others.

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.

ACcESS: Australia's contribution to the iSERVO Institute's development

The new Australian Computational Earth Systems Simulator research facility provides a virtual laboratory for studying the solid earth and its complex system behavior. The facility's capabilities complement those developed by overseas groups, thereby creating the infrastructure for an international computational solid earth research virtual observatory.

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.

Correlazioni tra proprietá reologiche, struttura e processabilitá di blend di LLDPE/LDPE

Polymer blends constitute a valuable way to produce relatively low cost new materials. A still open question concerns the miscibility of polyethylene blends. Deviations from the log-additivity rule of the newtonian viscosity are often taken as a signature of immiscibility of the two components. The aim of this thesis is to characterize the rheological behavior in shear and elongation of five series of LLDPE/LDPE blends whose parent polymers have been chosen with different viscosity and SCB content and length. Synergistic effects have been measured for both zero shear viscosity and melt strength. Both SCB length and viscosity ratio between the components have been found to be key parameters for the miscibility of the pure polymers. In particular the miscibility increases with increasing SCB length and with decreasing the LDPE molecular weight and viscosity. This rheological behavior has significant effects on the processability window of these blends when the uni or biaxial elongational flows are involved. The film blowing is one of the processes for which the synergistic effects above mentioned can be crucial. Small scale experiments of film blowing performed for one of the series of blends has demonstrated that the positive deviation of the melt strength enlarges the processability window. In particular, the bubble stability was found to improve or disappear when the melt strength of the samples increased. The blending of LDPE and LLDPE can even reduce undesired melt flow instability phenomena widening, as a consequence, the processability window in extrusion. One of the series of blends has been characterized by means of capillary rheometry in order to allow a careful morphological analysis of the surface of the extruded polymer jets by means of Scanning Electron Microscopy (SEM) with the aim to detect the very early stages of the small scale melt instabilty at low shear rates (sharksin) and to follow its subsequent evolution as long as the shear rate was increased. With this experimental procedure it was possible to evaluate the shear rate ranges corresponding to different flow regions: smooth extrudate surface (absence of instability), sharkskin (small scale instability produced at the capillary exit), stick-slip transition (instability involving the whole capillary wall) and gross melt fracture (i.e. a large scale "upstream" instability originating from the entrance region of the capillary). A quantitative map was finally worked out using which an assessment of the flow type for a given shear rate and blend composition can be predicted.

Slip complexity in dynamic models of earthquake faults.

We summarize recent evidence that models of earthquake faults with dynamically unstable friction laws but no externally imposed heterogeneities can exhibit slip complexity. Two models are described here. The first is a one-dimensional model with velocity-weakening stick-slip friction; the second is a two-dimensional elastodynamic model with slip-weakening friction. Both exhibit small-event complexity and chaotic sequences of large characteristic events. The large events in both models are composed of Heaton pulses. We argue that the key ingredients of these models are reasonably accurate representations of the properties of real faults.

Comparison of friction models applied to a control valve

Eight different models to represent the effect of friction in control valves are presented: four models based on physical principles and four empirical ones. The physical models, both static and dynamic, have the same structure. The models are implemented in Simulink/Matlab (R) and compared, using different friction coefficients and input signals. Three of the models were able to reproduce the stick-slip phenomenon and passed all the tests, which were applied following ISA standards. (C) 2008 Elsevier Ltd. All rights reserved.

On the frictional contacts in multibody system dynamics

"Series title: Computational methods in applied sciences, ISSN1871-3033, vol. 42"

Velocity fluctuations of fracturelike disruptions of associating polymer solutions

Velocity has been measured as a function of time for propagating crack tips as water is injected into solutions of end-capped associating polymers in a rectanguar Hele-Shaw cell. Measurements were performed for flows with different values of cell gap, channel width, polymer molecular weight, and polymer concentration. The condition for the onset of fracturelike behavior is well described by a Deborah number which uses the shear-thinning shear rate of the polymer solution as a characteristic frequency for network relaxation. At low molecular weight, the onset of fracturelike pattern evolution is accompanied by an abrupt jump in tip velocity, followed by a lower and approximately constant acceleration. At high molecular weight, the transition to fracturelike behavior involves passing through a regime that may be understood in terms of stick-slip dynamics. The crack-tip wanders from side to side and fluctuates (in both speed and velocity along the channel) with a characteristic frequency which depends linearly on the invading fluid injection rate.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.