Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “slip”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantifying the resistance to dislocation glide in single phase FeCrAl alloy

The resistance to dislocation glide associated with slip systems {110}<111> and {112}<111> in single phase FeCrAl alloy is measured via micromechanical testing in a scanning electronic microscopy at room temperature. Two important factors, the shape and orientation of a pillar, are discussed with respect to the glide resistance and stress-strain response. Maximizing Schmid factor of one specific slip system while minimizing others is recommended in order to diminish obvious dislocations-induced hardening during in-situ testing. Apparent Schmid factor analysis is conducted to select grains with preferred orientations. Two types of pillars with conventional cylindrical shape or dog-bone shape are tested under compression to estimate the resistance to dislocation glide and evaluate the effect of pillar shape on the compression stress-strain response. One dog-bone pillar is tested under tension to check the tension-compression isotropy of dislocation slip. We discover that the shape of a pillar to a smaller extent affects the measured resistance but strongly influences the stress strain response. Cylindrical pillars exhibit apparent hardening associated with early yielding due to stress or strain concentration at contact region, while dog-bone pillars show an obvious yielding and continuous shearing without hardening. The resistance is 220 MPa for slip system {110}<111> and 230 MPa for slip system {112}<111>. Finite element analysis is performed to account for the influence of pillar shape and contact condition on mechanical response.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2020

Analytical model for drag reduction on liquid-infused structured non-wetting surfaces

Liquid-infused structured non-wetting surfaces provide alternating no-slip and partial slip boundary conditions to the fluid flow, resulting in reduced friction at the interface. In this paper, an analytical model is developed for the evaluation of effective slip and, in turn, friction factor and drag reduction on liquid-infused structured non-wetting surfaces. By considering the entire range of anisotropy and heterogeneity of the surface structures as well as the full range of partial slip offered by the infusion liquid, the present model eliminates empirical fitting or correlations that are inherent in previous studies. Based on the effective slip length, drag reduction and skin friction coefficient values for Newtonian flow between two infinite parallel plates and flow in round tubes are presented. Extension of Moody charts for non-wetting surfaces and design maps of surface meso/micro/nano texturing for achieving desired drag reduction are presented for a broad range of engineering applications. The paper further presents independent validation of the model across experimental and computational data from the literature and brings together several previous studies in a unified manner.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High‐Temperature Mechanical Behavior of Single‐Crystal FeCrAl Alloy Under In Situ Micropillar Compression

FeCrAl cladding is one of the candidate materials for the near‐term accident‐tolerant fuel technologies under development. Research on high‐temperature mechanical behaviors of single‐crystal FeCrAl alloy is rather limited. Previous studies have reported the mechanical property of low‐index orientation in single‐crystal FeCrAl alloy at room temperature. However, the critical resolved shear stress to activate slip systems can be orientation and temperature dependent. Here, single‐crystal grains in a coarse‐grained FeCrAl alloy with different crystallographic orientations are selected to preferentially activate {110}<111> slip systems or {112}<111> slip systems. Micropillars are fabricated in the selected single‐crystal grains and tested at elevated temperatures in situ in a scanning electron microscope. The critical resolved shear stresses of {110}<111> slip systems and {112}<111> slip systems are determined at various temperatures. The critical resolved shear stress shows a temperature dependence and orientation independence. This study provides important insight for understanding the deformation mechanisms of FeCrAl alloys at elevated temperatures.

Sun, Tianyi↗
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Heterogeneous fatigue damage in a nickel-based single-crystal superalloy unraveled using correlative 3D X-ray technology

Nickel-based single-crystal (Ni-SX) superalloys under cyclic stress are susceptible to cracking at stress-concentration sites, eventually leading to low-cycle fatigue (LCF) failure. LCF cracks typically originate from intrinsic defects (e.g., voids and carbides) within solidified dendrites. However, systematic quantitative experimental analyses of defect-mediated local damage remain limited. To thoroughly understand the microscopic origins and evolution of LCF damage, correlated 3D mapping of dendrites across various regions is essential. Here, in this study, macroscale micro-computed tomography (μ-CT) was initially used to capture internal interdendritic secondary cracks within bulk DD413 superalloy after LCF testing at 760 °C. Subsequently, a multimodal methodology combining synchrotron 3D microdiffraction (3D-μXRD), high-resolution μ-CT, and electron microscopy was established. This approach allowed precise localization of internal damage zones near interdendritic secondary cracks and detailed mapping of the 3D correlated distributions of dendrites, defects, and residual stress/strain fields within these zones at submicron spatial resolution. Finally, the same approach was applied to specimens subjected to interrupted loading at approximately 40 % of the fatigue life to uncover the early damage states of dendrites. The dendrite cores (DCs) and interdendritic regions (IDs) exhibit microscale heterogeneous mechanical responses: nearly defect-free DCs accumulate local irreversible slip along specific slip systems to generate slip bands, while the IDs containing various defects accommodate local microplasticity through the activation of multiple slip systems around these defects. The local tensile stress near defects in the IDs exceeds that in the DC slip band regions by more than threefold, leading to the generation of local damage zones within the IDs. Chain-like defect distributions facilitate the interconnection of these local zones into a continuous damage region, further elevating the overall tensile stress in the IDs. Additionally, geometrically necessary dislocations alone are insufficient as indicators of LCF damage; both the internal stress state and its magnitude must be considered. These experimental results provide critical data and insights for the development of multi-physics fatigue models.

Localized deformation↗
DOE OSTI2020

Dislocation density-based plasticity model from massive discrete dislocation dynamics database

In this work, we present a dislocation density-based strain hardening model for single crystal copper through a systematic coarse-graining analysis of more than 200 discrete dislocation dynamics (DDD) simulations of plastic deformation under uniaxial tension. The proposed constitutive model has two components: a generalized Taylor relation connecting resolved shear stresses to dislocation densities on individual slip systems, and a generalized Kocks-Mecking model for dislocation multiplication. The DDD data strongly suggests a logarithmic dependence of flow stress on the plastic shear strain rate on each slip system, and, equivalently, an exponential dependence of the plastic shear strain rate on the resolved shear stress. Hence the proposed generalized Taylor relation subsumes the Orowan relation for plastic flow. The DDD data also calls for a correction to the Kocks-Mecking model of dislocation multiplication to account for the increase of dislocation density on slip systems with negligible plastic shear strain rate. This is accomplished by allowing the multiplication rate on each slip system to include contributions from the plastic strain rates of the two coplanar slip systems. The resulting constitutive model successfully captures the strain hardening rate dependence on the loading orientation as predicted by the DDD simulations, which is also consistent with existing experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

Subduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well-studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present-day megathrust at depths of ~5–>25 km, and we shear these at their peak in situ pressure-temperature (P-T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity-strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity-weakening behavior was observed at 350°C, with velocity-strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a-b) increases with slip velocity at temperatures up to ~200°C and remains constant or decreases with slip velocity at higher temperatures. Altogether, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

58 GEOSCIENCES↗
DOE OSTI2023

Foreshock properties illuminate nucleation processes of slow and fast laboratory earthquakes

Understanding the connection between seismic activity and the earthquake nucleation process is a fundamental goal in earthquake seismology with important implications for earthquake early warning systems and forecasting. We use high-resolution acoustic emission (AE) waveform measurements from laboratory stick-slip experiments that span a spectrum of slow to fast slip rates to probe spatiotemporal properties of laboratory foreshocks and nucleation processes. We measure waveform similarity and pairwise differential travel-times (DTT) between AEs throughout the seismic cycle. AEs broadcasted prior to slow labquakes have small DTT and high waveform similarity relative to fast labquakes. We show that during slow stick-slip, the fault never fully locks, and waveform similarity and pairwise differential travel times do not evolve throughout the seismic cycle. In contrast, fast laboratory earthquakes are preceded by a rapid increase in waveform similarity late in the seismic cycle and a reduction in differential travel times, indicating that AEs begin to coalesce as the fault slip velocity increases leading up to failure. These observations point to key differences in the nucleation process of slow and fast labquakes and suggest that the spatiotemporal evolution of laboratory foreshocks is linked to fault slip velocity.

58 GEOSCIENCES↗
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARY Observational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

, RSQSim↗
DOE OSTI2024

Development of prechamber enabled mixing-controlled combustion strategy for ultra-low methane emissions from lean burn natural gas engines

This numerical study explores the optimization of Prechamber Enabled Mixing-Controlled Combustion (PC-MCC) using natural gas in heavy-duty engines, aiming to enhance combustion efficiency to minimize methane slip and NOx emissions. The approach involves a prechamber ignition system, distinct from conventional spark ignition (SI) systems, to initiate combustion of direct injected natural gas. By leveraging the robust ignition characteristics of the prechamber, the PC-MCC method demonstrates significant potential in achieving efficient combustion akin to diesel engines but with lower greenhouse gas emissions. The research evaluates the effects of various geometric and operational parameters on the combustion process and emissions, including prechamber volume, nozzle diameter, direct injector (DI) geometry, and engine operating strategies. Computational Fluid Dynamics (CFD) simulations are utilized, focusing on a heavy-duty, single-cylinder engine modeled after the Caterpillar C9.3B engine. Key findings indicate that a prechamber volume of 3 cc, coupled with a nozzle diameter of 2.75 mm for two prechamber holes, strikes an optimal balance between combustion efficiency and emissions reduction. This configuration ensures robust combustion across a range of operating conditions while maintaining methane slip within targeted limits. Further investigation into DI geometry shows the significance of the injector umbrella angle and nozzle diameter in shaping the fuel-air mixing and combustion dynamics. An umbrella angle of 130° and a nozzle diameter of 300 microns are identified as optimal, promoting rapid and efficient combustion with minimized methane and NOx emissions. The study also investigates the impact of injection timing and pressure, highlighting their roles in controlling combustion timing and influencing emissions levels. Advanced injection timing is found to be crucial in achieving the desired low methane slip, whereas retarded injection timing assists to reduce NOx emissions while having a slight increase in methane emissions. Operating strategies incorporating various levels of Exhaust Gas Recirculation (EGR) are assessed for their effectiveness in further reducing emissions. The research demonstrates that a judicious combination of internal hot EGR and careful calibration of DI pressure and SOI timing can achieve significant reductions in NOx emissions while keeping methane slip under control. Specifically, an internal EGR level of 15%–25%, combined with DI pressures of 200–300 bar and injection timings at or after top dead center, is recommended. These findings contribute valuable insights into the development of advanced combustion techniques for natural gas engines, offering a viable pathway to reduce methane slip without compromising engine efficiency or performance. The PC-MCC system presents a promising solution for the future of heavy-duty natural gas engine technology to reduce methane emissions.

Nsaif, Osama↗
DOE OSTI2023

Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake

An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.

58 GEOSCIENCES↗
DOE OSTI2024

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr –1 . The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

58 GEOSCIENCES↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Investigation of crossed-twin structure formation in magnesium and magnesium alloys

In this work, the effect of alloying addition on the propensity for twin-twin interactions to transform into crossed-twin structures in magnesium alloys is investigated. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) framework combined with a discrete twin model and dislocation density-based hardening law for slip strengths is used to calculate the micromechanical fields in the crystalline matrix around the interacting twins. AZ31 and MgLi alloys are selected along with pure Mg to study the influence of plastic anisotropy in connection with alloying elements. These alloys were selected since their plastic anisotropy measure, which is defined as the ratio between the critical resolved shear stress for pyramidal $\langle$c +a$\rangle$ and basal $\langle$a$\rangle$ slip modes, spanned a wide range. To quantify the role of twin thicknesses, we probe a range of impinging twin thicknesses while fixing the recipient twin thickness. The analysis reveals that: (i) the local driving stress for crossed twin structure formation generated from the interaction of the two twins is lower in a low plastically anisotropic material, like a MgLi alloy, than a high plastically anisotropic material like pure Mg and (ii) the critical impinging twin thickness needed to form the crossed twin structure in pure Mg, AZ31 and MgLi alloys is ~0.5, ~0.75 and ~1.5 times the recipient twin thickness. We propose a relationship between the tendency for crossed twin structure formation and the experimentally observed higher ductility in MgLi alloys compared to pure Mg. Finally, one key implication of the findings is that crossed- twin structure formation can be hindered and the ductility of magnesium alloy thereby improved by properly choosing alloying elements that lower the slip strength for pyramidal $\langle$c + a$\rangle$ slip.

36 MATERIALS SCIENCE↗