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 271 records · Page 15

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗

Can Glacial Sea‐Level Drop‐Induced Gas Hydrate Dissociation Cause Submarine Landslides?

Abstract We conducted two‐dimensional numerical simulations to investigate the mechanisms underlying the strong spatiotemporal correlation observed between submarine landslides and gas hydrate dissociation due to glacial sea‐level drops. Our results suggest that potential plastic deformation or slip could occur at localized and small scales in the shallow‐water portion of the gas hydrate stability zone (GHSZ). This shallow‐water portion of the GHSZ typically lies within the area enclosed by three points: the BGHSZ–seafloor intersection, the seafloor at ∼600 m below sea level (mbsl), and the base of the GHSZ (BGHSZ) at ∼1,050 mbsl in low‐latitude regions. The deep BGHSZ (>1,050 mbsl) could not slip; therefore, the entire BGHSZ was not a complete slip surface. Glacial hydrate dissociation alone is unlikely to cause large‐scale submarine landslides. Observed deep‐water (much greater than 600 mbsl) turbidites containing geochemical evidence of glacial hydrate dissociation potentially formed from erosion or detachment in the GHSZ pinch‐out zone.

58 GEOSCIENCES↗

Deformation of Hydrous Phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and Stishovite [Si 1‐n H 4n O 2 ] Relevant to Anisotropy of the Earth's Mantle

Seismic anisotropy of the Earth's mantle has been mostly attributed to crystallographic preferred orientation (CPO) generated during subduction and convection of an anhydrous mantle. But some hydrous phases are also stable at mantle conditions. Here we present results from diamond-anvil cell deformation experiments at high pressure and temperature on hydrous phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and hydrous stishovite [Si 1-n H 4n O 2 ], transformed from the clay mineral kaolinite. They develop strong CPO during axial compression, suggesting that they likely contribute to seismic anisotropy and heterogeneity in the mantle. Comparing experimental results with viscoplastic polycrystal plasticity models suggest that phase Egg deforms dominantly by (001) slip, δ by (010) slip and stishovite by {100} slip which could be incorporated in future models of mantle geodynamics.

58 GEOSCIENCES↗

Enhancing fatigue life by ductile-transformable multicomponent B2 precipitates in a high-entropy alloy

Catastrophic accidents caused by fatigue failures often occur in engineering structures. Thus, a fundamental understanding of cyclic-deformation and fatigue-failure mechanisms is critical for the development of fatigue-resistant structural materials. Here we report a high-entropy alloy with enhanced fatigue life by ductile-transformable multicomponent B2 precipitates. Its cyclic-deformation mechanisms are revealed by real-time in-situ neutron diffraction, transmission-electron microscopy, crystal-plasticity modeling, and Monte-Carlo simulations. Multiple cyclic-deformation mechanisms, including dislocation slips, precipitation strengthening, deformation twinning, and reversible martensitic phase transformation, are observed in the studied high-entropy alloy. Its improved fatigue performance at low strain amplitudes, i.e., the high fatigue-crack-initiation resistance, is attributed to the high elasticity, plastic deformability, and martensitic transformation of the B2-strengthening phase. This study shows that fatigue-resistant alloys can be developed by incorporating strengthening ductile-transformable multicomponent intermetallic phases.

36 MATERIALS SCIENCE↗

Uncovering the phonon spectra and lattice dynamics of plastically deformable InSe van der Waals crystals

Stacking two-dimensional (2D) van der Waals (vdW) materials in a layered bulk structure provides an appealing platform for the emergence of exotic physical properties. As a vdW crystal with exceptional plasticity, InSe offers the opportunity to explore various effects arising from the coupling of its peculiar mechanical behaviors and other physical properties. Here, we employ neutron scattering techniques to investigate the correlations of plastic interlayer slip, lattice anharmonicity, and thermal transport in InSe crystals. Not only are the interlayer slip direction and magnitude well captured by shifts in the Bragg reflections, but we also observe a deviation from the expected Debye behaviour in the heat capacity and lattice thermal conductivity. Combining the experimental data with first-principles calculations, we tentatively attribute the observed evidence of strong phonon-phonon interactions to a combination of a large acoustic-optical frequency resonance and a nesting effect. These findings correlate the macroscopic plastic slip and the microscopic lattice dynamics, providing insights into the mechano-thermo coupling and modulation in 2D vdW materials.

36 MATERIALS SCIENCE↗

Irradiation-enhanced torsional buckling capacity of carbon nanotube bundles

Molecular dynamics simulations are used to understand the torsional buckling of pristine and irradiated carbon nanotube (CNT) bundles. Irradiation-induced inter-tube defects are shown to significantly increase the critical buckling torque and critical buckling angle, while slightly increasing the torsional stiffness. In contrast, intra-tube defects are found to degrade the torsional properties. Such competing interactions cause irradiation enhancement to occur in large bundles where significant inter-tube bonding can occur. However, the irradiation enhancement effect becomes weak for very large bundles in which enhanced inter-tube interactions already exist in unirradiated bundles. In pristine CNT bundles of all sizes under torsional loading, CNTs can slip via the weakly interacting van der Waals force, whereas in the irradiated bundles, the inter-tube defects prevent slipping. The study further shows that the formation of one-dimensional carbon chain defects contributes to enhanced friction under slipping.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chiral propulsion: The method of effective boundary conditions

We propose to apply an “effective boundary condition” method to the problem of chiral propulsion. For the case of a rotating helix moving through a fluid at a low Reynolds number, the method amounts to replacing the original helix (in the limit of small pitch) by a cylinder, but with a special kind of partial slip boundary conditions replacing the non-slip boundary conditions on the original helix. These boundary conditions are constructed to reproduce far-field velocities of the original problem and are defined by a few parameters (slipping lengths) that can be extracted from a problem in planar rather than cylindrical geometry. We derive the chiral propulsion coefficients for spirals, helicoids, helically modulated cylinders and some of their generalizations using the introduced method. In the case of spirals, we compare our results with the ones derived by Lighthill and find a very good agreement. Here, the proposed method is general and can be applied to any helical shape in the limit of a small pitch. Furthermore, we have established that for a broad class of helical surfaces the dependence of the chiral propulsion on the helical angle θ is universal, x ~ cos θ sin 2θ with the maximal propulsion achieved at the universal angle θ m = tan –1 (1/√2) ≈ 35.26°.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shape effects on the local dynamics of suspensions of spheroidal particles

The effect of shape on the dynamics of suspensions of non-spherical heavy particles is examined by fully resolved numerical simulations of oblate and prolate spheroids, as well as spheres, for a density ratio of ten, volume fractions ranging from 0.5% to 5%, and Reynolds numbers between 20 and 30. The dynamics is determined both by the interactions of the particles with the fluid as well as by collisions, with the number and importance of collisions increasing with volume fractions. A single isolated oblate or prolate spheroid falling under gravity generally falls broadside on, for the governing parameters examined here, and at low-volume fractions, the majority of particles in a suspension fall that way. At higher-volume fractions, the orientation is more random. The slip velocity decreases as the volume fraction increases for all shapes, as expected, but the effect of the shape is much less than seen for a single particle. Furthermore, this seems to be due to two effects. For all volume fractions, the most deformed particles cluster more than spheres and less deformed particles, which increases their slip velocity. As the concentration increases, the increased particle interactions also causes more particles to fall short side-on, which reduces the frontal area and the resulting drag, increasing the slip velocity. This second effect is, of course, absent for spherical particles.

42 ENGINEERING↗

Spherical and sessile droplet dynamics by fluctuating hydrodynamics

We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.

Bell, John B↗

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗

High-Fidelity CFD Modeling of Cryogenic Hydrogen Isotope Extrusion for Fusion Reactor Pellet Fueling

This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.

ANSYS-Polyflow↗

Hypothesis tests on Rayleigh wave radiation pattern shapes: a theoretical assessment of idealized source screening

SUMMARY Shallow seismic sources excite Rayleigh wave ground motion with azimuthally dependent radiation patterns. We place binary hypothesis tests on theoretical models of such radiation patterns to screen cylindrically symmetric sources (like explosions) from non-symmetric sources (like non-vertical dip-slip or non-VDS faults). These models for data include sources with several unknown parameters, contaminated by Gaussian noise and embedded in a layered half-space. The generalized maximum likelihood ratio tests that we derive from these data models produce screening statistics and decision rules that depend on measured, noisy ground motion at discrete sensor locations. We explicitly quantify how the screening power of these statistics increase with the size of any dip-slip and strike-slip components of the source, relative to noise (faulting signal strength) and how they vary with network geometry. As applications of our theory, we apply these tests to (1) find optimal sensor locations that maximize the probability of screening non-circular radiation patterns and (2) invert for the largest non-VDS faulting signal that could be mistakenly attributed to an explosion with damage, at a particular attribution probability. Finally, we quantify how certain errors that are sourced by opening cracks increase screening rate errors. While such theoretical solutions are ideal and require future validation, they remain important in underground explosion monitoring scenarios because they provide fundamental physical limits on the discrimination power of tests that screen explosive from non-VDS faulting sources.

58 GEOSCIENCES↗

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗

Androgynous { 10 1 ¯ 2 } twin in zinc

Under ambient conditions, Zn is a hexagonal metal with a large c / a ratio of 1.856. Plastic deformation is predominantly accommodated by basal 〈 a 〉 slip and compression twins on the { 10 1 ¯ 2 } planes. Increasing hydrostatic pressure drastically reduces the c / a ratio of Zn and, when a critical threshold of c / a = 3 at about 10 GPa is crossed, the { 10 1 ¯ 2 } twins are predicted to change from compressive to tensile in nature. What happens at the transition point, when c / a = 3 , remains unknown. Here, we strain-cycle a textured polycrystalline sample of pure Zn at uniform hydrostatic pressures ranging between 2 and 17 GPa, over which the c / a ratio crosses the c / a = 3 compressive-tensile transition for { 10 1 ¯ 2 } twins. During deformation, the state of the sample is monitored through x-ray diffraction to extract texture and internal strain evolution. By comparing the experimental results with the predictions of an elastoviscoplastic polycrystal simulation, we confirm the androgynous nature of { 10 1 ¯ 2 } twin response at low and high pressures. When c / a = 3 , polycrystalline Zn does not display any evidence of twinning and its plastic behavior is controlled by mostly basal and pyramidal 〈 c + a 〉 slip activity, with a very small contribution of prismatic 〈 a 〉 slip. Evidence for the activity of other { 10 1 ¯ n } twinning modes, which have been suggested for Zn under high pressure, are not observed. Published by the American Physical Society 2024

Merkel, Sébastien (ORCID:000000032767581X)↗

Helmholtz vibrations in bowed strings

We report for almost 160 years, it has been known that Helmholtz oscillations, unique to vibrating strings in bowed instruments (violin, cello, etc.), have two distinct regimes: “slip” and “stick.” During the slip regime, the force at the bow-string interaction is attributed to friction between the sliding bow hair and the vibrating string, with a friction coefficient that decreases with increasing relative velocity. Yet the hair-string interaction during the stick regime is less understood. We propose that the interaction force during the stick regime is proportional to the product of the longitudinal acoustic impedance of the bow hair to the relative bow-string velocity. We validate this hypothesis by solving the string's differential equation of motion, including an enhanced formulation to avoid parasitic high-frequency oscillations. This physical model enables us to analyze, in real time, the characteristics of the Helmholtz oscillations, including the string shape, excitation of harmonics, Schelleng ripples, and string energy, showing that the bowed string gains energy during the stick regime and loses energy during the slip regime.

36 MATERIALS SCIENCE↗

Seismic fault weakening via CO 2 pressurization enhanced by mechanical deformation of dolomite fault gouges

Carbon dioxide emissions from dolomite decarbonation play an essential role in the weakening of carbonate faults by lowering the effective normal stress, which is thermally activated at temperatures above 600–700 °C. However, the mechanochemical effect of low-crystalline ultrafine fault gouge on the decarbonation and slip behavior of dolomite-bearing faults remains unclear. In this study, we obtained a series of artificial dolomite fault gouges with systematically varying particle sizes and dolomite crystallinities using a high-energy ball mill. The laboratory-scale pulverization of dolomite yielded MgO at temperatures below 50 °C, indicating that mechanical decarbonation without significant heating occurred due to the collapse of the crystalline structure, as revealed by X-ray diffraction and solid-state nuclear magnetic resonance results. Furthermore, the onset temperature of thermal decarbonation decreased to ~400 °C. Numerical modeling reproduced this two-stage decarbonation, where the pore pressure increased due to low-temperature thermal decarbonation, leading to slip weakening on the fault plane even at 400–500 °C; i.e., 200–300 °C lower than previously reported temperatures. Thus, the presence of small amounts of low-crystalline dolomite in a fault plane may lead to a severely reduced shear strength due to thermal decomposition at ~400 °C with a small slip weakening distance.

58 GEOSCIENCES↗