Rate-and-State friction as a bulk visco-plastic flow law that includes generation, diffusion, and healing of distributed damage

2021 ◽  
Author(s):  
Casper Pranger ◽  
Patrick Sanan ◽  
Dave May ◽  
Alice Gabriel
Keyword(s):  
Numerical Models ◽  
Shear Bands ◽  
Theoretical Work ◽  
Small Scale ◽  
Theoretical Understanding ◽  
State Dependent ◽  
Tectonic Environment ◽  
Rate And State Friction ◽  
Original Laboratory ◽  
Flow Law

<p>The rate- and state-dependent friction (RSF) laws (Dieterich, 1979, JGR; Ruina, 1983, JGR-SE) have been widely successful in capturing the behavior of sliding surfaces in laboratory settings, as well as reproducing a range of natural fault slip phenomena in numerical models.</p><p>Studies of exhumed fault zones make it clear that faults are not two-dimensional features at the junction of two distinct bodies of rock, but instead evolve into complex damage zones that show clear signs of multi-scale fracturing, grain diminution, hydro-thermal effects and chemical and petrological changes. Many of these observed factors have been experimentally verified, and several studies have furthered our theoretical understanding of earthquakes and other seismic phenomena as volumetric, bulk-rock processes, including Sleep (1995, 1997), Lyakhovsky and Ben-Zion et al. (2011, J. Mech. Phys. Solids; 2014, PAGeoph; 2014,  J. Mech. Phys. Solids; 2016, GJI), Niemeijer, Chen, van den Ende et al. (2007, 2016, JGR-SE; 2018, Tectonophysics), Roubicek (2014, GJI), and Barbot (2019, Tectonophysics).</p><p>While the established numerical modeling approach of simulating faults as planar features undergoing friction can be a useful and powerful homogenization of small-scale volumetric processes, there are also cases where this practice falls short -- most notably when studying faults that grow and evolve in response to a changing tectonic environment. This is mainly due to the computational challenges associated with automating the construction of a fault-resolving conformal mesh.</p><p>Motivated by this issue, we formulate a generalization of RSF as a plastic or viscous flow law with generation, diffusion, and healing of damage that gives rise to mathematically and numerically well-behaved finite shear bands that closely mimic the behavior of the original laboratory-derived formulation. The proposed formulation includes the well-known RSF laws for an infinitely thin fault as a limit case as the damage diffusion length scale tends to zero. In contrast to previous theoretical work we focus only on a mathematical formalism that is used to generalize and regularize the existing RSF laws in order to retain close correspondence to existing experimental and numerical results. We will demonstrate the behavior of this new bulk RSF formulation with results of 1D and 2D numerical simulations, and hope to engage in a preliminary discussion of the physical implications.</p>

scholarly journals Chloride binding in concrete: recent investigations and recognised knowledge gaps: RILEM Robert L’Hermite Medal Paper 2021

2021 ◽  
Vol 54 (6) ◽  
Author(s):  
Klaartje De Weerdt
Keyword(s):  
Reinforced Concrete ◽  
Numerical Models ◽  
Concrete Structures ◽  
Theoretical Work ◽  
Reinforced Concrete Structures ◽  
Chloride Ingress ◽  
Chloride Binding ◽  
Theoretical Understanding ◽  
The Face ◽  
Composite Cements

AbstractA theoretical understanding of chloride binding is urgently needed if we are to use the new low-CO2 composite cements in reinforced concrete structures exposed to chloride-loaded environments. The empirical models and “deemed-to-satisfy” rules currently in use will not help in the face of the wide variety of new SCMs now being proposed. Engineers need generic chloride ingress models that can handle the variations in binder chemistry and exposure conditions. This paper reviews a selection of recent investigations by a team of NTNU researchers and international collaborators on various Portland composite cements using a combination of equilibrium chloride-binding experiments and thermodynamic modelling. One of our main findings is that “leaching” simulated by varying the pH and the calcium concentration has a dominant effect on chloride binding. However, challenges regarding the quantification and characterization of the chloride binding phases have yet to be overcome. To advance in this field we need additional crucial thermodynamic data for chloride-binding hydrates formed by low-CO2 cements containing SCMs, and we need to overcome some experimental challenges. A major break-through would come from understanding the interactions between chlorides and C(-A)-S-H. Part of the answer might be found at the atomic level. Finally, mechanistic numerical models are needed to upscale the findings from chloride binding to chloride ingress models. We conclude by showing the significance of such theoretical work on chloride binding for engineers tasked with the construction and maintenance of the reinforced concrete structures that are so vital a part of modern infrastructure.

scholarly journals The Effectiveness of Ensemble-Neural Network Techniques to Predict Peak Uplift Resistance of Buried Pipes in Reinforced Sand

2021 ◽  
Vol 11 (3) ◽  
pp. 908
Author(s):  
Jie Zeng ◽  
Panagiotis G. Asteris ◽  
Anna P. Mamou ◽  
Ahmed Salih Mohammed ◽  
Emmanuil A. Golias ◽  
...  
Keyword(s):  
Neural Network ◽  
Numerical Models ◽  
Correlation Coefficients ◽  
Network Models ◽  
Small Scale ◽  
Offshore Platforms ◽  
Neural Network Models ◽  
Buried Pipes ◽  
Ensemble Techniques ◽  
Uplift Resistance

Buried pipes are extensively used for oil transportation from offshore platforms. Under unfavorable loading combinations, the pipe’s uplift resistance may be exceeded, which may result in excessive deformations and significant disruptions. This paper presents findings from a series of small-scale tests performed on pipes buried in geogrid-reinforced sands, with the measured peak uplift resistance being used to calibrate advanced numerical models employing neural networks. Multilayer perceptron (MLP) and Radial Basis Function (RBF) primary structure types have been used to train two neural network models, which were then further developed using bagging and boosting ensemble techniques. Correlation coefficients in excess of 0.954 between the measured and predicted peak uplift resistance have been achieved. The results show that the design of pipelines can be significantly improved using the proposed novel, reliable and robust soft computing models.

scholarly journals Observational aspects of prominence oscillations

2007 ◽  
Vol 3 (S247) ◽  
pp. 152-157 ◽  
Author(s):  
Oddbjørn Engvold
Keyword(s):  
Magnetic Structure ◽  
Numerical Models ◽  
Small Scale ◽  
Solar Prominences ◽  
Scale Structures ◽  
Oscillations And Waves ◽  
Available Information ◽  
Small Scale Structures

AbstractSeismology has become a powerful tool in studies of the magnetic structure of solar prominences and filaments. Reversely, analytical and numerical models are guided by available information about the spatial and thermodynamical structure of these enigmatic structures. The present invited paper reviews recent observational results on oscillations and waves as well as details about small-scale structures and dynamics of prominences and filaments.

Coupled Thermo-Hydro-Geochemical Models of Engineered Barrier Systems: The Febex Project

2000 ◽  
Vol 663 ◽  
Author(s):  
J. Samper ◽  
R. Juncosa ◽  
V. Navarro ◽  
J. Delgado ◽  
L. Montenegro ◽  
...  
Keyword(s):  
Large Scale ◽  
Reference Model ◽  
Numerical Models ◽  
Full Scale ◽  
Small Scale ◽  
Model Parameters ◽  
In Situ Test ◽  
Wide Range ◽  
Engineered Barrier

ABSTRACTFEBEX (Full-scale Engineered Barrier EXperiment) is a demonstration and research project dealing with the bentonite engineered barrier designed for sealing and containment of waste in a high level radioactive waste repository (HLWR). It includes two main experiments: an situ full-scale test performed at Grimsel (GTS) and a mock-up test operating since February 1997 at CIEMAT facilities in Madrid (Spain) [1,2,3]. One of the objectives of FEBEX is the development and testing of conceptual and numerical models for the thermal, hydrodynamic, and geochemical (THG) processes expected to take place in engineered clay barriers. A significant improvement in coupled THG modeling of the clay barrier has been achieved both in terms of a better understanding of THG processes and more sophisticated THG computer codes. The ability of these models to reproduce the observed THG patterns in a wide range of THG conditions enhances the confidence in their prediction capabilities. Numerical THG models of heating and hydration experiments performed on small-scale lab cells provide excellent results for temperatures, water inflow and final water content in the cells [3]. Calculated concentrations at the end of the experiments reproduce most of the patterns of measured data. In general, the fit of concentrations of dissolved species is better than that of exchanged cations. These models were later used to simulate the evolution of the large-scale experiments (in situ and mock-up). Some thermo-hydrodynamic hypotheses and bentonite parameters were slightly revised during TH calibration of the mock-up test. The results of the reference model reproduce simultaneously the observed water inflows and bentonite temperatures and relative humidities. Although the model is highly sensitive to one-at-a-time variations in model parameters, the possibility of parameter combinations leading to similar fits cannot be precluded. The TH model of the “in situ” test is based on the same bentonite TH parameters and assumptions as for the “mock-up” test. Granite parameters were slightly modified during the calibration process in order to reproduce the observed thermal and hydrodynamic evolution. The reference model captures properly relative humidities and temperatures in the bentonite [3]. It also reproduces the observed spatial distribution of water pressures and temperatures in the granite. Once calibrated the TH aspects of the model, predictions of the THG evolution of both tests were performed. Data from the dismantling of the in situ test, which is planned for the summer of 2001, will provide a unique opportunity to test and validate current THG models of the EBS.

scholarly journals Numerical models of slab migration in continental collision zones

Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 293-306 ◽  
Author(s):  
V. Magni ◽  
J. van Hunen ◽  
F. Funiciello ◽  
C. Faccenna
Keyword(s):  
Subduction Zone ◽  
Plate Tectonics ◽  
Numerical Models ◽  
Continental Collision ◽  
Force Balance ◽  
Small Scale ◽  
Dip Angle ◽  
External Forces ◽  
Stress Dependent ◽  
Subduction System

Abstract. Continental collision is an intrinsic feature of plate tectonics. The closure of an oceanic basin leads to the onset of subduction of buoyant continental material, which slows down and eventually stops the subduction process. In natural cases, evidence of advancing margins has been recognized in continental collision zones such as India-Eurasia and Arabia-Eurasia. We perform a parametric study of the geometrical and rheological influence on subduction dynamics during the subduction of continental lithosphere. In our 2-D numerical models of a free subduction system with temperature and stress-dependent rheology, the trench and the overriding plate move self-consistently as a function of the dynamics of the system (i.e. no external forces are imposed). This setup enables to study how continental subduction influences the trench migration. We found that in all models the slab starts to advance once the continent enters the subduction zone and continues to migrate until few million years after the ultimate slab detachment. Our results support the idea that the advancing mode is favoured and, in part, provided by the intrinsic force balance of continental collision. We suggest that the advance is first induced by the locking of the subduction zone and the subsequent steepening of the slab, and next by the sinking of the deepest oceanic part of the slab, during stretching and break-off of the slab. These processes are responsible for the migration of the subduction zone by triggering small-scale convection cells in the mantle that, in turn, drag the plates. The amount of advance ranges from 40 to 220 km and depends on the dip angle of the slab before the onset of collision.

scholarly journals The South Georgia Wave Experiment: A Means for Improved Analysis of Gravity Waves and Low-Level Wind Impacts Generated from Mountainous Islands

2018 ◽  
Vol 99 (5) ◽  
pp. 1027-1040 ◽  
Author(s):  
D. R. Jackson ◽  
A. Gadian ◽  
N. P. Hindley ◽  
L. Hoffmann ◽  
J. Hughes ◽  
...  
Keyword(s):  
High Resolution ◽  
Gravity Waves ◽  
Numerical Models ◽  
Small Island ◽  
Small Scale ◽  
The South ◽  
South Georgia ◽  
Low Level ◽  
Wave Experiment ◽  
Realistic Representation

AbstractGravity waves (GWs) play an important role in many atmospheric processes. However, the observation-based understanding of GWs is limited, and representing them in numerical models is difficult. Recent studies show that small islands can be intense sources of GWs, with climatologically significant effects on the atmospheric circulation. South Georgia, in the South Atlantic, is a notable source of such “small island” waves. GWs are usually too small scale to be resolved by current models, so their effects are represented approximately using resolved model fields (parameterization). However, the small-island waves are not well represented by such parameterizations, and the explicit representation of GWs in very-high-resolution models is still in its infancy. Steep islands such as South Georgia are also known to generate low-level wakes, affecting the flow hundreds of kilometers downwind. These wakes are also poorly represented in models.We present results from the South Georgia Wave Experiment (SG-WEX) for 5 July 2015. Analysis of GWs from satellite observations is augmented by radiosonde observations made from South Georgia. Simulations were also made using high-resolution configurations of the Met Office Unified Model (UM). Comparison with observations indicates that the UM performs well for this case, with realistic representation of GW patterns and low-level wakes. Examination of a longer simulation period suggests that the wakes generally are well represented by the model. The realism of these simulations suggests they can be used to develop parameterizations for use at coarser model resolutions.

Hydrostatic Collapse Pressure and Radial Collapse Force Comparisons for Ultra-Deepwater Pipelines

2014 ◽  
Author(s):  
Marco A. P. Rosas ◽  
Ana Paula F. Souza ◽  
Marcos V. Rodrigues ◽  
Danilo Machado L. da Silva
Keyword(s):  
Numerical Models ◽  
Research Work ◽  
Compressive Force ◽  
Quadratic Function ◽  
Radial Force ◽  
Small Scale ◽  
Collapse Pressure ◽  
First Results ◽  
Specimen Test ◽  
The Relationship

In this paper the behavior and the relationship between hydrostatic collapse pressure and diametrically opposed radial compressive force for pipelines were analyzed. This study presents an introduction of a research work aimed to assess the pipeline collapse pressure based on the radial collapse force. Initially the hydrostatic collapse pressure is analyzed, for pipes with different diameter to wall thickness ratio (D/t) and ovalities, using classical assessment (DNV method) and numerical models (FE). Then, the compressive radial force is also analyzed using numerical models validated by a small-scale ring specimen test. After that, the relationship between hydrostatic collapse pressure and compressive radial force is discussed. These first results show that the radial force is a quadratic function of the collapse pressure.

scholarly journals Turbulent vertical diffusivity in the sub-tropical stratosphere

2007 ◽  
Vol 7 (3) ◽  
pp. 6603-6629 ◽  
Author(s):  
I. Pisso ◽  
B. Legras
Keyword(s):  
Transport Model ◽  
Numerical Models ◽  
Small Scale ◽  
Chemistry Transport Model ◽  
Weather Forecasts ◽  
Chemical Tracers ◽  
Vertical Diffusivity ◽  
Large Ensembles ◽  
Local Equivalent ◽  
Using Data

Abstract. Vertical (cross-isentropic) mixing is produced by small-scale turbulent processes which are still poorly understood and parametrized in numerical models. In this work we provide estimates of local equivalent diffusion in the lower stratosphere by comparing balloon borne high-resolution measurements of chemical tracers with reconstructed mixing ratio from large ensembles of random Lagrangian backward trajectories using European Center for Medium-range Weather Forecasts analysed winds and a chemistry-transport model (REPROBUS). We have investigated cases in subtropical latitudes using data from HIBISCUS campaign. Upper bound on the vertical diffusivity is found to be of the order of 0.5 m2 s−1 in the subtropical region, which is larger than the estimates at higher latitudes. The relation between diffusion and dispersion is studied by estimating Lyapunov exponents and studying their variation according to the presence of active dynamical structures.

scholarly journals The role of Edge-Driven Convection in the generation of volcanism I: a 2D systematic study

2020 ◽  
Author(s):  
Antonio Manjón-Cabeza Córdoba ◽  
Maxim Ballmer
Keyword(s):  
Numerical Models ◽  
High Volume ◽  
Oceanic Lithosphere ◽  
Companion Paper ◽  
Small Scale ◽  
Continental Lithosphere ◽  
Transient Phenomenon ◽  
Mantle Viscosity ◽  
Wide Range

Abstract. The origin of intraplate volcanism is not explained by the plate tectonic theory, and several models have been put forward for explanation. One of these models involves Edge-Driven Convection (EDC), in which cold and thick continental lithosphere is juxtaposed to warm and thin oceanic lithosphere to trigger convective instability. To test whether EDC can produce long-lived high-volume magmatism, we run numerical models of EDC for a wide range of mantle properties and edge (i.e., the oceanic-continental transition) geometries. We find that the most important parameters that govern EDC are the rheological paramaters mantle viscosity η0 and activation energy Ea. However, even the maximum melting volumes found in our models are insufficient to account for island-building volcanism on old seafloor, such as at the Canary Islands and Cape Verde. Also, beneath old seafloor, localized EDC-related melting commonly transitions into widespread melting due to small-scale sublithospheric convection, inconsistent with the distribution of volcanism at these volcanic chains. In turn, EDC is a good candidate to sustain the formation of small seamounts on young seafloor, as it is a highly transient phenomenon that occurs in all our models soon after initiation. In a companion paper, we investigate the implications of interaction of EDC with mantle-plume activity.

scholarly journals Numerical Simulation of Deep-Sea Sediment Transport Induced by a Dredge Experiment in the Northeastern Pacific Ocean

2021 ◽  
Vol 8 ◽  
Author(s):  
Kaveh Purkiani ◽  
Benjamin Gillard ◽  
André Paul ◽  
Matthias Haeckel ◽  
Sabine Haalboom ◽  
...  
Keyword(s):  
Sediment Transport ◽  
Pacific Ocean ◽  
Suspended Sediment ◽  
Deep Sea ◽  
Numerical Models ◽  
Sediment Concentration ◽  
Ocean Model ◽  
Small Scale ◽  
Plume Dispersion ◽  
Deep Sea Sediment

Predictability of the dispersion of sediment plumes induced by potential deep-sea mining activities is still very limited due to operational limitations on in-situ observations required for a thorough validation and calibration of numerical models. Here we report on a plume dispersion experiment carried out in the German license area for the exploration of polymetallic nodules in the northeastern tropical Pacific Ocean in 4,200 m water depth. The dispersion of a sediment plume induced by a small-scale dredge experiment in April 2019 was investigated numerically by employing a sediment transport module coupled to a high-resolution hydrodynamic regional ocean model. Various aspects including sediment characteristics and ocean hydrodynamics were examined to obtain the best statistical agreement between sensor-based observations and model results. Results show that the model is capable of reproducing suspended sediment concentration and redeposition patterns observed during the dredge experiment. Due to a strong southward current during the dredging, the model predicts no sediment deposition and plume dispersion north of the dredging tracks. The sediment redeposition thickness reaches up to 9 mm directly next to the dredging tracks and 0.07 mm in about 320 m away from the dredging center. The model results suggest that seabed topography and variable sediment release heights above the seafloor cause significant changes especially for the low sedimentation pattern in the far-field area. Near-bottom mixing is expected to strongly influence vertical transport of suspended sediment.

Sign in / Sign up

Export Citation Format

Share Document