No Arabic abstract
We investigate the nature of friction in granular layers by means of numerical simulation focusing on the critical slip distance, over which the system relaxes to a new stationary state. Analyzing a transient process in which the sliding velocity is instantaneously changed, we find that the critical slip distance is proportional to the sliding velocity. We thus define the relaxation time, which is independent of the sliding velocity. It is found that the relaxation time is proportional to the layer thickness and inversely proportional to the square root of the pressure. An evolution law for the relaxation process is proposed, which does not contain any length constants describing the surface geometry but the relaxation time of the bulk granular matter. As a result, the critical slip distance is scaled with a typical length scale of a system. It is proportional to the layer thickness in an instantaneous velocity change experiment, whereas it is scaled with the total slip distance in a spring-block system on granular layers.
The way granular materials response to an applied shear stress is of the utmost relevance to both human activities and natural environment. One of the their most intriguing and less understood behavior, is the stick-instability, whose most dramatic manifestation are earthquakes, ultimately governed by the dynamics of rocks and debris jammed within the fault gauge. Many of the features of earthquakes, i.e. intermittency, broad times and energy scale involved, are mimicked by a very simple experimental set-up, where small beads of glass under load are slowly sheared by an elastic medium. Analyzing data from long lasting experiments, we identify a critical dynamical regime, that can be related to known theoretical models used for crackling-noise phenomena. In particular, we focus on the average shape of the slip velocity, observing a breakdown of scaling: while small slips show a self-similar shape, large does not, in a way that suggests the presence of subtle inertial effects within the granular system. In order to characterise the crossover between the two regimes, we investigate the frictional response of the system, which we trat as a stochastic quantity. Computing different averages, we evidence a weakening effect, whose Stribeck threshold velocity can be related to the aforementioned breaking of scaling.
YBaCuO nanowires were reproducibly fabricated down to widths of 50 nm. A Au/Ti cap layer on YBCO yielded high electrical performance up to temperatures above 80 K in single nanowires. Critical current density of tens of MA/cm2 at T = 4.2 K and of 10 MA/cm2 at 77 K were achieved that survive in high magnetic fields. Phase-slip processes were tuned by choosing the size of the nanochannels and the intensity of the applied external magnetic field. Data indicate that YBCO nanowires are rather attractive system for the fabrication of efficient sensors, supporting the notion of futuristic THz devices.
The coupled mechanics of fluid-filled granular media controls the behavior of many natural systems such as saturated soils, fault gouge, and landslides. The grain motion and the fluid pressure influence each other: It is well established that when the fluid pressure rises, the shear resistance of fluid-filled granular systems decreases, and as a result catastrophic events such as soil liquefaction, earthquakes, and accelerating landslides may be triggered. Alternatively, when the pore pressure drops, the shear resistance of these systems increases. Despite the great importance of the coupled mechanics of grains-fluid systems, the basic physics that controls this coupling is far from understood. We developed a new multi-scaled model based on the discrete element method, coupled with a continuum model of fluid pressure, to explore this dynamical system. The model was shown recently to capture essential feedbacks between porosity changes arising from rearrangement of grains, and local pressure variations due to changing pore configurations. We report here new results from numerical experiments of a continuously shearing layer of circular two-dimensional grains, trapped between two parallel rough boundaries. The experiments use a fixed confining stress on the boundary walls, and a constant velocity applied to one of the boundaries, as if this system was the interior of a sliding geological fault filled with fault gouge. In addition, we control the layer permeability and the drainage boundary conditions. This paper presents modeling results showing that the localization of shear (into a narrow shear band within the shearing layer) is strongly affected by the presence of fluids. While in dry granular layers there is no preferred position for the onset of localization, drained systems tend to localize shear on their boundary. We propose a scaling argument to describe the pressure deviations in a shear band, and use that to predict the allowable positions of shear localizations as a function of the fault and gouge properties.
Data-driven machine-learning for predicting instantaneous and future fault-slip in laboratory experiments has recently progressed markedly due to large training data sets. In Earth however, earthquake interevent times range from 10s-100s of years and geophysical data typically exist for only a portion of an earthquake cycle. Sparse data presents a serious challenge to training machine learning models. Here we describe a transfer learning approach using numerical simulations to train a convolutional encoder-decoder that predicts fault-slip behavior in laboratory experiments. The model learns a mapping between acoustic emission histories and fault-slip from numerical simulations, and generalizes to produce accurate results using laboratory data. Notably slip-predictions markedly improve using the simulation-data trained-model and training the latent space using a portion of a single laboratory earthquake-cycle. The transfer learning results elucidate the potential of using models trained on numerical simulations and fine-tuned with small geophysical data sets for potential applications to faults in Earth.
Earthquakes cause lasting changes in static equilibrium, resulting in global deformation fields that can be observed. Consequently, deformation measurements such as those provided by satellite based InSAR monitoring can be used to infer an earthquakes faulting mechanism. This inverse problem requires a numerical forward model that is both accurate and fast, as typical inverse procedures require many evaluations. The Weakly-enforced Slip Method (WSM) was developed to meet these needs, but it was not before applied in an inverse problem setting. Consequently, it was unknown what effect particular properties of the WSM, notably its inherent continuity, have on the inversion process. Here we show that the WSM is able to accurately recover slip distributions in a Bayesian-inference setting, provided that data points in the vicinity of the fault are removed. In a representative scenario, an element size of 2 km was found to be sufficiently fine to generate a posterior probability distribution that is close to the theoretical optimum. For rupturing faults a masking zone of 20 km sufficed to avoid numerical disturbances that would otherwise be induced by the discretization error. These results demonstrate that the WSM is a viable forward method for earthquake inversion problems. While our synthesized scenario is basic for reasons of validation, our results are expected to generalize to the wider gamut of scenarios that finite element methods are able to capture. This has the potential to bring modeling flexibility to a field that if often forced to impose model restrictions in a concession to computability.