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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

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↗
DOE OSTI2022

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↗
DOE OSTI2024

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)↗
DOE OSTI2019

Cycloidal paths in physics as superpositions of translational and rotational motions

Cycloidal paths are ubiquitous in physics. Here, we show that representative cycloidal paths in physics can be described as superpositions of translations and rotations of a point through space. Using this unifying principle, the parametric equations of the path of a point on a rolling disk are derived for rolling without slipping, rolling with frictionless slipping, and when kinetic solid-on-solid friction is present during rolling with slipping. In a similar way, the parametric equations versus time for the orbit with respect to a star of a moon in a circular orbit about a planet that is in a circular orbit about a star are derived, where the orbits are coplanar. The parametric equations versus time for the path of the magnetization vector during undamped electron-spin resonance are found using the same principle, which show that cycloidal paths can occur under specified conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

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↗
DOE OSTI2019

Mixed deformation styles observed on a shallow subduction thrust, Hikurangi margin, New Zealand

Geophysical observations show spatial and temporal variations in fault slip style on shallow subduction thrust faults, but geological signatures and underlying deformation processes remain poorly understood. International Ocean Discovery Program (IODP) Expeditions 372 and 375 investigated New Zealand’s Hikurangi margin in a region that has experienced both tsunami earthquakes and repeated slow-slip events. We report direct observations from cores that sampled the active Papaku splay fault at 304 m below the seafloor. This fault roots into the plate interface and comprises an 18-m-thick main fault underlain by ~30 m of less intensely deformed footwall and an ~10-m-thick subsidiary fault above undeformed footwall. Fault zone structures include breccias, folds, and asymmetric clasts within transposed and/or dismembered, relatively homogeneous, silty hemipelagic sediments. The data demonstrate that the fault has experienced both ductile and brittle deformation. As a result, this structural variation indicates that a range of local slip speeds can occur along shallow faults, and they are controlled by temporal, potentially far-field, changes in strain rate or effective stress.

58 GEOSCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2018

In situ neutron diffraction study on tensile deformation behavior of carbon-strengthened CoCrFeMnNi high-entropy alloys at room and elevated temperatures

In this paper, carbon is doped into a CoCrFeMnNi high-entropy alloy as an interstitial atom, improving the single phase solid solution alloy with a good combination of strength and ductility at room temperature by introducing deformation twins. In situ neutron diffraction (ND) is applied to investigate the carbon-doped CoCrFeMnNi deformation mechanism and micromechanical behaviors during uniaxial tension at room and elevated temperatures. With in situ results accompanied with the microstructure and texture measurement, it is found that the plastic deformation is dominated by dislocation slip at an early stage at both temperatures. However, at high strain level, deformation is mediated simultaneously by deformation twins and microbands at room temperature, while it is governed solely by microbands at elevated temperature of 573 K. The evolution of lattice strain, peak intensity, and peak width from in situ ND elucidates the micromechanical behaviors regarding the role of slips and twins. Finally, the texture represented by orientation distribution function indicates that the initial specimen possesses a relatively strong {112} <110> texture component, and the room-temperature tension deformed texture comprises of slip-induced fiber texture and twinning-induced {115} <552> texture component.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Multi-Objective Optimization of a Wrought Magnesium Alloy for High Strength and Ductility

An optimization technique is coupled with crystal plasticity based finite element (CPFE) computations to aid the microstructural design of a wrought magnesium alloy for improved strength and ductility. The initial microstructure consists of a collection of sub-micron sized grains containing deformation twins. The variables used in the simulations are crystallographic texture, and twin spacing within the grains. It is assumed that plastic deformation occurs mainly by dislocation slip on two sets of slip systems classified as hard and soft modes. The hard modes are those slip systems that are inclined to the twin planes and the soft mode consists of dislocation glide along the twin plane. The CPFE code calculates the stress-strain response of the microstructure as a function of the microstructural parameters and the length-scale of the features. A failure criterion based on a critical shear strain and a critical hydrostatic stress is used to define ductility. The optimization is based on the sequential generation of an initial population defined by the texture and twin spacing variables. The CPFE code and the optimizer are coupled in parallel so that new generations are created and analyzed dynamically. In each successive generation, microstructures that satisfy at least 90% of the mean strength and mean ductility in the current generation are retained. Multiple generation runs based on the above procedure are carried out in order to obtain maximum strength-ductility combinations. Finally, the implications of the computations for the design of a wrought magnesium alloy are discussed. Research sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U. S. Department of Energy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of the Earthquake Rupture on Ground-Motion Variability of the 24 August 2016 M w 6.2 Amatrice, Italy, Earthquake

Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.

3D Ground Motions Simulations↗
DOE OSTI2024

Rupture Model of the 5 April 2024 Tewksbury, New Jersey, Earthquake Based on Regional Lg -Wave Data

On 5 April 2024, an earthquake of magnitude 4.8 occurred in Tewksbury, New Jersey. It was the largest instrumentally recorded event since 1900 in New Jersey and southern New York. Millions of people around New York City, ~65 km east–northeast of Tewksbury, felt the shaking from the mainshock, but the epicentral area experienced no known significant property damages. We determine the focal mechanism, which is oblique faulting, and retrieve the Lg–wave relative source time functions (RSTFs) from the stations at regional distances to understand rupture processes and ground motions. Our fault–slip models well explain azimuthal variations of the RSTFs. The models show the rupture propagating toward the east–northeast (~50° to 60°), not along the fault strike. The slip distribution on the nodal plane striking north and dipping to the east shows a slip area of 1.1 km radius with the rupture propagating down–dip. The down–dip rupture may account for the observed lack of strong shaking in the epicentral area.

58 GEOSCIENCES↗
DOE OSTI2013

Statistical Mechanics Modeling of Mesoscale Deformation in Metals

The research under this project focused on a theoretical and computational modeling of dislocation dynamics of mesoscale deformation of metal single crystals. Specifically, the work aimed to implement a continuum statistical theory of dislocations to understand strain hardening and cell structure formation under monotonic loading. These aspects of crystal deformation are manifestations of the evolution of the underlying dislocation system under mechanical loading. The project had three research tasks: 1) Investigating the statistical characteristics of dislocation systems in deformed crystals. 2) Formulating kinetic equations of dislocations and coupling these kinetics equations and crystal mechanics. 3) Computational solution of coupled crystal mechanics and dislocation kinetics. Comparison of dislocation dynamics predictions with experimental results in the area of statistical properties of dislocations and their field was also a part of the proposed effort. In the first research task, the dislocation dynamics simulation method was used to investigate the spatial, orientation, velocity, and temporal statistics of dynamical dislocation systems, and on the use of the results from this investigation to complete the kinetic description of dislocations. The second task focused on completing the formulation of a kinetic theory of dislocations that respects the discrete nature of crystallographic slip and the physics of dislocation motion and dislocation interaction in the crystal. Part of this effort also targeted the theoretical basis for establishing the connection between discrete and continuum representation of dislocations and the analysis of discrete dislocation simulation results within the continuum framework. This part of the research enables the enrichment of the kinetic description with information representing the discrete dislocation systems behavior. The third task focused on the development of physics-inspired numerical methods of solution of the coupled dislocation kinetics and crystal mechanics framework. To a large extent, this task has also been successfully started. We have developed a custom finite-element approach with mesh points being a subset of the underlying crystal structure. When used to predict the evolution of the dislocation system, the planar motion of dislocations is naturally captured for all slip systems, thus minimizing numerical errors and providing simple ways to investigate cross slip and dislocation reactions. Preliminary results in this direction show that we are closer than ever in building a predictive framework for dislocation dynamics and mesoscale plasticity based on the first principles of dislocation dynamics. The rest of the report gives and overview of the research performed under this project and highlights the key results and open questions left for future investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2017

2017 GTO Project review Laboratory Evaluation of EGS Shear Stimulation.

The objectives and purpose of this research has been to produce laboratory-based experimental and numerical analyses to provide a physics-based understanding of shear stimulation phenomena (hydroshearing) and its evolution during stimulation. Water was flowed along fractures in hot and stressed fractured rock, to promote slip. The controlled laboratory experiments provide a high resolution/high quality data resource for evaluation of analysis methods developed by DOE to assess EGS “behavior” during this stimulation process. Segments of the experimental program will provide data sets for model input parameters, i.e., material properties, and other segments of the experimental program will represent small scale physical models of an EGS system, which may be modeled. The coupled lab/analysis project has been a study of the response of a fracture in hot, water-saturated fractured rock to shear stress experiencing fluid flow. Under this condition, the fracture experiences a combination of potential pore pressure changes and fracture surface cooling, resulting in slip along the fracture. The laboratory work provides a means to assess the role of “hydroshearing” on permeability enhancement in reservoir stimulation. Using the laboratory experiments and results to define boundary and input/output conditions of pore pressure, thermal stress, fracture shear deformation and fluid flow, and models were developed and simulations completed by the University of Oklahoma team. The analysis methods are ones used on field scale problems. The sophisticated numerical models developed contain parameters present in the field. The analysis results provide insight into the role of fracture slip on permeability enhancement-“hydroshear” is to be obtained. The work will provide valuable input data to evaluate stimulation models, thus helping design effective EGS.

58 GEOSCIENCES↗
DOE OSTI2019

Nanoindentation of Micrograins in Polycrystals under Multiaxial Stress: Control of Abrupt & Stochastic Plastic Events

Deforming objects in the world of the small has proven a much more complex task than just bending a common-sized object like a spoon. Over the last decade, it has become clear that the application of stress on micron-sized objects, such as micropillars, leads to abrupt and stochastic plastic strain jumps that do not satisfy Weibull statistics (as expected for uncorrelated defect yielding) but instead seem to follow highly wide, power-law distributions, pointing towards cooperative dislocation plasticity. Are power-law distributed plastic events the rule or the exception in crystals? If the rule, can we design protocols to control and suppress them in order to achieve deterministic and smooth microforming processes? We propose to combine experiments and theory to unravel the mysteries of abrupt microplastic deformation in a natural dislocation environment: micron-sized grains (\micrograins") of polycrystals. We propose a novel experimental design that utilizes nanoindentation and 4-point bending techniques on the substrate of polycrystalline thin films of two typical FCC metals (Al and Cu), and make concrete theoretical conjectures, based on intuition built from previous theoretical and experimental efforts in micropillar compression and nanoindentation. Based on a wide set of observations over the last decade, we put forward two distinct hypotheses that we propose to investigate in both theory and experiments: Our first conjecture is that crystal plasticity in the central region of micron-sized grains has micropillar-like stochastic features: power-law plastic strain jumps at stresses above the macroscale yield strength, displaying strong size and rate effects. We will demonstrate this similarity in spherical nanoindentation experiments in the central-region of multi-orientation grains of Al and Cu thin films. Al and Cu are selected for their disparity in stacking fault energy and consequently, their different propensity for cross-slip relaxation; however, our protocol may be applied to any crystalline material of technological interest. We will test our hypothesis by performing two-dimensional extended discrete dislocation plasticity (2.5D-DDP) simulations, inspired by comparable successful studies of micropillar plasticity, and extended to model nanoindentation with two slip systems and on a surface of a bicrystal with microhard grain-boundaries (GB), ie. not allowing dislocation absorption (116). Our model will include effective three dimensional dislocation mechanisms such as dislocation multiplication and double-cross-slip assisted glide. Second, we hypothesize that the application of uniform in-plane stress on our thin films below the yield point can smoothen the micrograin plastic response during nanoindentation. On a technological level, we propose that indentation-induced localization (sink-in) may be diminished by the application of concurrent in-plane stress that activates slow, viscoplastic relaxations near GBs. We aim to demonstrate our hypothesis through a 4-point bending apparatus that introduces a homogeneous in-plane stress state before indentation. We will perform simulations of a coarse-grained mesoscale theory of plasticity in a system with microhard GBs, benchmarked through 2.5D-DDP simulations, which will capture the effects of slow viscoplastic relaxations. Our nanoindentation abrupt event statistics studies, either in the micrograin interior or near GBs, will be correlated to direct observations of induced surface dislocation microstructures. For we will use TEM-based methods (PACOM/TKD) in combination with EBSD. Indents placed in the center of grains will be analyzed based on the orientation of the grain, while those near to GBs will be analyzed based on the orientation of each grain, the distance from the GB and the 5 degrees of freedom of the GB.

36 MATERIALS SCIENCE↗
DOE OSTI2020

NH 4 OH Looping with Membrane CO 2 Absorber and Distributed Stripper for Enhanced Algae Growth

The University of Kentucky Center for Applied Energy (UK CAER) has devised a unique, integrated CO2 capture and utilization technology. CO2 from coal-fired power generation flue gas is first captured at half the operating cost of a typical aqueous CO2 capture system (CCS), distributed in an aqueous stream and then fixed by algae in bioreactors where the algae production is increased by 50% over that with a typical intermittent nutrient feeding system. Lower CCS operating cost is achieved by eliminating the flue gas pretreatment step for cooling and SO2 removal, eliminating steam extraction from the power generation steam cycle for solvent regeneration, and eliminating CO2 compression. Higher algae production is achieved by continuous, just-in-time nutrient feed to the bioreactors directly from a distributed solvent regenerator, which maintains the bioreactor pH for optimum growth. The process starts with a uniquely configured membrane absorber, where the flue gas is indirectly contacted with an ammonium hydroxide (NH4OH) solvent. Dissolved NH3 is attractive for both CO2 capture and as an algae nutrient. For CO2 capture it is inexpensive, has a low regeneration energy, is thermally- and oxidatively-stable and has a viscosity near that of water, which makes is easy to transport. Numerous studies have shown that the scrubbing capacity of NH3 is approximately 0.9-1.2 kg of CO2/kg of NH3, with a CO2 removal efficiency of ~99% and half the solvent regeneration energy than that of 30 wt% MEA[1, 2, 3]. NH3 is attractive as an algae nutrient due to its low cost. The rich NH4OH solvent is pumped to a set of distributed regenerators which are co-located with the algae bioreactors. Solvent pumping, transport and distribution reduces the balance of plant (BOP) cost compared to a typical aqueous CCS related to the flue gas duct and boost fan required to transport the flue gas. The energy required for the distributed solvent regeneration is supplied by solar-thermal panels eliminating the need for steam extraction from the power generation steam cycle. After solvent regeneration, the product stream contains both the CO2 captured from the flue gas and volatized NH3 from the solvent. This product stream is fed directly to the bioreactors, eliminating the need for compression of the CO2 stream. The relative amounts of CO2 and NH3 in the product stream are adjusted and controlled by a controlling the regeneration conditions (pressure and temperature). The continuous feed of the right ratio of nutrients overcomes the typical inhibition of algae growth resulting from frequent pH swings in the bioreactor due to unbalanced (intermittent) feeding systems for CO2 and N. Also, because the regenerators will operate at pressure and be located in close proximity to the bioreactors, there is no worry about pressure drop when sparging the gas into the algae. Sparging produces small bubbles which is beneficial for mass transfer efficiency. One known challenge when using an NH4OH solvent is high NH3 emission. Hydrophobic membranes are used for CO2 capture using an aqueous NH3 solution[4, 5] without the direct contact between flue gas and aqueous solution. Additionally, UK CAER CO2 capture and utilization process manages NH3 slip in three extra measures. First, NH3 slip is minimized by working with minimal species partial pressure, which is proportional to the concentration in the liquid. Hence, lowering the capture solvent concentration will lower the NH3 partial pressure. Second, UK CAER’s previous work has demonstrated that the addition of Zn2+ into NH3 solutions to chelate the NH3 can reduce NH3 volatility. Third, the configuration of the membrane CO2 absorber utilizes condensed water from the flue gas to continually wash the gas-side of the membrane to reduce fouling and recapture NH3 slip. Additional details about the UK CAER unique, integrated CO2 capture and utilization technology will be presented along with technology development plans. Diao, N., Q. Li, and Z. Fang. 2004. Heat transfer in ground heat exchangers with groundwater advection. International Journal of Thermal Sciences. 43: 1203-1211, He, Q., M. Chen, L. Meng, K. Liu, and W. Pan. 2004. Study on Carbon Dioxide Removal from Flue Gas by Absorption of Aqueous Ammonia. Western Kentucky University. Yeh, A.C., and H. Bai. 1999. Comparison of ammonia and monoethanolamine solvents to reduce CO2 greenhouse gas emissions. The Science of the Total Environment. 228: 121-133, Villeneuve, K., D. Roizard, J.C. Remigy, M. Iacono, and S. Rode. 2018. CO2 capture by aqueous ammonia with hollow fiber membrane contactors: Gas phase reactions and performance stability. Separation and Purification Technology, 199: 189-197, Toro Molina, C., and C. Bouallou. 2016. Carbon dioxide absorption by ammonia intensified with membrane contactors. Clean Techn Environ Policy 18, 2133–2146 (2016)

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Data Challenges for Structural Health Monitoring of Electrical Machines [Slides]

Induction motors (ID) consist of a stator and a rotor. Induction motor slip decreases with decreasing load. Induction motor slip decreases with decreasing load. This talk will focus on removing closely coupled power grid frequencies in an effort to observe the true dynamic response signature of induction motors. Using the time difference between the zeros, we found the frequency as it changes over time for each of the 3 phases. The frequencies of the 3 phases were averaged to reduce the effects of sidebands on the instantaneous frequency. The frequencies of the 3 phases were averaged to reduce the effects of sidebands on the instantaneous frequency. The Kalman filter rapidly synchronizes with the expected output. The sidebands were recreated using spectral subtraction. The frequency domain shows the reduction in magnitude of the power grid frequency the Kalman output was subtracted out. Kalman filtering also reveals important data about the sidebands in healthy motors. The Continuous Wavelet Transform is an alternative to a Bandpass IIR filter used in the analysis of the acceleration data. Analysis of the acceleration data reveals 4 sidebands around the slip frequency.

42 ENGINEERING↗
DOE OSTI2022

Characterization and Modeling of Deformation Induced Damage in Titanium Alloys

Over the lifetime of their DOE support, and specifically during the past funding cycle, the PIs have primarily focused on ambient temperature studies of the mechanisms of heterogeneous deformation and shear transfer at GBs. The PIs have developed and applied a broad array of experimental approaches complemented by advanced computational methods and simulations to enable new understanding of GB shear transfer, and a new approach for modeling this behavior at the mechanistic level has been established. The PIs have examined how slip behavior differs on the surface and interior of polycrystals using a variety of approaches. Some of the associated data has been used to generate 3-D models of experimentally-characterized microstructures that both support interpretation of deformation mechanisms as well as reveal where modeling needs to be improved to more closely agree with what is experimentally observed. The PIs progress in constitutive modeling of crystal plasticity has introduced explicit GB resistance to shear that depends on the mutual geometry of active slip systems. This past work places the PIs in a good position to take on the challenges associated with mechanistically understanding the interplay between GBS and slip interactions.

42 ENGINEERING↗
DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗