Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mechanical robustness”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Game theoretic modeling and optimization of competition and collaboration in dual channel electronic waste supply chains

The rapid growth of electronic waste (e-waste) presents critical challenges for sustainable resource recovery and environmental protection. This study develops a dual-channel closed-loop supply chain (CLSC) model formulated as a hierarchical Stackelberg game, that integrates dynamic pricing and cost-sharing mechanisms to optimize both economic and environmental outcomes. The model explicitly captures strategic interactions between manufacturer-led and third-party recycling channels, accounting for consumer behavior, regulatory incentives, and market competition. Numerical simulations conducted (implemented over a four-iteration horizon using a commercial optimization solver) show that, relative to the baseline equilibrium, manufacturer profit increases from 11.6 thousand USD to 37.9 thousand USD (+226.8%), total recycled volume rises from 7,848 to 7,942 units (+1.2%), and collector profit nearly doubles under cost-sharing, enabling more equitable profit distribution. Furthermore, scenario-based simulations across Sub-Saharan Africa, high-income economies, and emerging Asian industrial countries reveal that infrastructure quality, policy intensity, and labor costs critically shape recycling efficiency and profit allocation. These findings demonstrate that subsidies alone are insufficient to ensure system efficiency. Instead, coordinated strategies that integrate internal incentive alignment with context-sensitive policy support are required. Overall, this study offers a robust framework for designing resilient, efficient, and regionally adaptable e-waste management systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Crystallization and Phase Transformations of Aluminum (Oxy)hydroxide Polymorphs in Caustic Aqueous Solution

Gibbsite, bayerite, and boehmite are important aluminum (oxy)hydroxide minerals in nature and have been widely deployed in various industrial applications. They are also major components in caustic nuclear wastes stored at various U.S. locations. Knowledge of their crystallization and phase transformation processes contributes to understanding their occurrence and could help optimize waste treatment processes. While it has been reported that partial conversion of bayerite and gibbsite to boehmite occurs in basic solutions at elevated temperatures, systematic studies of factors affecting the phase transformation as well as the underlying reaction mechanisms are non-existent, particularly in highly alkaline solutions. We explored the effects of sodium hydroxide concentrations (0.1~3 M), reaction temperature (60~100 ?) and aluminum concentrations (0.1~1 M) on the crystallization and transformation of these aluminum (oxy)hydroxides. Detailed structural and morphological characterization by X-ray diffraction (XRD), scanning electron microscopy (SEM), and nuclear magnetic resonance (NMR) spectrometry revealed that these processes depend largely on the reaction temperature and the Al/OH- ratio. When 1 = Al/OH- = 2.5, the reactions favor formation of high crystallinity precipitates, whereas at Al/OH- ratio ? 2.5 precipitation ceases unless the Al concentration is higher than 1 M. We identified pseudoboehmite, bayerite and gibbsite as intermediate phases to bayerite, gibbsite and boehmite, respectively, all of which transform via dissolution-reprecipitation. Gibbsite transforms to boehmite in both acidic and weak caustic environments at temperatures above 80 oC. However, a ‘bar-shaped’ gibbsite morphology dominates in highly caustic environments (3 M NaOH). The findings enable a robust basis for selection of various solid phases by tuning the reaction conditions.

aluminum (oxy)hydroxide polymorphs, crystallizatio↗
DOE OSTI2022

How a Formate Dehydrogenase Responds to Oxygen: Unexpected O 2 Insensitivity of an Enzyme Harboring Tungstopterin, Selenocysteine, and [4Fe–4S] Clusters

The reversible two-electron interconversion of formate and CO 2 is catalyzed by both nonmetallo- and metallo-formate dehydrogenases (FDHs). The latter group comprises molybdenum- or tungsten-containing enzymes with the metal coordinated by two equivalents of a pyranopterin cofactor, a cysteinyl or selenocysteinyl (Sec) ligand supplied by the polypeptide, and a catalytically essential terminal sulfido ligand. In addition, these biocatalysts incorporate one or more [4Fe–4S] clusters for facilitating long-distance electron transfer. However, an interesting dichotomy arises when attempting to understand how the metallo-FDHs react with O 2 . Whereas existing scholarship portrays these enzymes as being unable to perform in air due to extreme O 2 lability of their metal centers, studies dating as far back as the 1930s emphasize that some of these systems exhibit formate oxidase (FOX) activity, coupling formate oxidation to O 2 reduction. Therefore, to reconcile these conflicting views, we explored context-dependent functional linkages between metallo-FDHs and their cognate electron acceptors within the same organism vis-à-vis catalysis under atmospheric O 2 . Here, we report the discovery and characterization of an O 2 -insensitive FDH2 from the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough (DvH) that ligates tungsten, Sec, and four [4Fe–4S] clusters. By advancing a robust expression platform for its recombinant production, we eliminate both the requirement of nitrate or azide during purification and reductive activation with thiols and/or formate prior to catalysis. Because the distinctive spectral signatures of formate-reduced DvH-FDH2 remain invariant under anaerobic and aerobic conditions, we benchmarked the enzyme activity in air, identifying CO 2 as the catalytic product. Full reaction progress curve analysis discloses a high catalytic efficiency when probed with a high-potential artificial electron acceptor. Furthermore, we show that DvH-FDH2 enables near-stoichiometric hydrogen peroxide production without superoxide release to achieve O 2 insensitivity. Notably, simultaneous electron transfer to cytochrome c and O 2 reveals that metal-based electron bifurcation is operational in this system. Taken together, our work proves the co-occurrence of redox bifurcated FDH and FOX activities within a metalloenzyme scaffold. These findings set the stage for uncovering previously unknown O 2 -insensitive flavin-based electron bifurcation mechanisms, as well as for developing authentic formate/air biofuel cells, engineering O 2 -stable FDHs and biohybrid metallocatalysts, and discerning formate bioenergetics of gut microbiota.

13C NMR↗
DOE OSTI2020

Almost medium-free measurement of the Hoyle state direct-decay component with a TPC

The structure of the Hoyle state, a highly α -clustered state at 7.65 MeV in 12 C, has long been the subject of debate. Understanding if the system comprises of three weakly interacting α particles in the 0s orbital, known as an α-condensate state, is possible by studying the decay branches of the Hoyle state. The direct decay of the Hoyle state into three α particles, rather than through the 8 Be ground state, can be identified by studying the energy partition of the three α particles arising from the decay. This paper provides details on the breakup mechanism of the Hoyle stating using a new experimental technique. Method: By using β-delayed charged-particle spectroscopy of 12 N using the Texas active target time-projection chamber, a high-sensitivity measurement of the direct 3α decay ratio can be performed without contributions from pileup events. A Bayesian approach to understanding the contribution of the direct components via a likelihood function shows that the direct component is <0.043% at the 95% confidence level. This value is in agreement with several other studies, and, here, we can demonstrate that a small nonsequential component with a decay fraction of about 10 –4 is most likely. Here, the measurement of the non-sequential component of the Hoyle state decay is performed in an almost medium-free reaction for the first time. The derived upper limit is in agreement with previous studies and demonstrates sensitivity to the absolute branching ratio. Further experimental studies would need to be combined with robust microscopic theoretical understanding of the decay dynamics to provide additional insight into the idea of the Hoyle state as an α condensate.

6 ≤ A ≤ 19↗
DOE OSTI2020

Characterizing photovoltaic backsheet adhesion degradation using the wedge and single cantilever beam tests, Part II: Accelerated tests

Photovoltaic (PV) backsheets provide critical moisture, mechanical, and electrical insulation to the backside of PV modules, but their continued functionality depends upon their ability to remain well adhered over years of harsh environmental exposure. A study of adhesive strength was conducted on several PV backsheet types exposed to indoor accelerated weathering. Two adhesion tests – the wedge test and single cantilever beam test – were used to measure adhesion energy in four backsheets: two with fluoropolymer-based outer (airside) layers – polyvinyl fluoride (PVF), polyvinylidene fluoride (PVDF) – and two with non-fluoropolymer-based outer layers–polyamide (PA) and polyethylene terephthalate (PET). The effects of changes in temperature, ultraviolet (UV) irradiance, relative humidity, and a periodic water spray were studied. The PET-, PVF- and PVDF-based backsheets all showed instances of debonding in an adhesive layer, which is commonly reported. Furthermore, for the PA-, PET-, and PVDF-based backsheets, adhesion energy dropped fairly rapidly with exposure, reaching very low levels by 4000 h. The PVF-based backsheet was relatively robust to weathering. Pull-off of an outer chalking layer in the PA- and PET-based backsheets suggested a sensitivity to UV. Adhesion in these two backsheets was also highly sensitive to changes in moisture level. Changes in UV irradiance and temperature did not show a statistically significant effect on adhesion loss for the exposure levels used. The present work is part of a two-part adhesion study on both field-weathered and indoor-exposed backsheets, and forms a basis for understanding adhesion degradation across a variety of backsheet types and degradation factors.

14 SOLAR ENERGY↗
DOE OSTI2020

A robust and multifunctional calcium coordination polymer as a selective fluorescent sensor for acetone and iron (+3) and as a tunable proton conductor

A multifunctional calcium coordination polymer (Ca-CP) formulated as [(CH 3 ) 2 NH 2 ][Ca(Me 2 tcpbH)(H 2 O)] (1) has been synthesized and characterized. Compound 1 features a robust and novel three-dimensional (3D) network with cage-like cavities. The high thermal and moisture stability of the title compound makes it a good candidate for possible environment or energy related applications. Fluorescence (FL) studies demonstrate that 1 shows emission with a maximum at 350 nm under 270 nm excitation. The title compound acts as a selective and sensitive fluorescent sensor for the detection of acetone and Fe 3+ ions at low concentrations. The FL intensity quenching percentage and the concentration ratio of Fe 3+ /Fe 2+ showed a linear relation with a Fe 3+ concentration as low as 50 μM in Fe 2+ aqueous solutions, offering an easy way to probe the extent of Fe 2+ oxidation. Moreover, due to the relatively high hydration character of Ca 2+ , the terminal water molecules and undeprotonated carboxylate groups form a strong hydrogen-bonded network within the structure, facilitating its tunable proton conduction under various humidity and temperature conditions. The proton conducting mechanism has been investigated using PXRD, IR, and EA methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Dynamic subcanopy leaf traits drive resistance of net primary production across a disturbance severity gradient

Across the globe, the forest carbon sink is increasingly vulnerable to an expanding array of low- to moderate-severity disturbances. However, some forest ecosystems exhibit functional resistance (i.e., the capacity of ecosystems to continue functioning as usual) following disturbances such as extreme weather events and insect or fungal pathogen outbreaks. Unlike severe disturbances (e.g., stand-replacing wildfires), moderate severity disturbances do not always result in near-term declines in forest production because of the potential for compensatory growth, including enhanced subcanopy production. Community-wide shifts in subcanopy plant functional traits, prompted by disturbance-driven environmental change, may play a key mechanistic role in resisting declines in net primary production (NPP) up to thresholds of canopy loss. However, the temporal dynamics of these shifts, as well as the upper limits of disturbance for which subcanopy production can compensate, remain poorly characterized. In this study, we leverage a 4-year dataset from an experimental forest disturbance in northern Michigan to assess subcanopy community trait shifts as well as their utility in predicting ecosystem NPP resistance across a wide range of implemented disturbance severities. Through mechanical girdling of stems, we achieved a gradient of severity from 0% (i.e., control) to 45, 65, and 85% targeted gross canopy defoliation, replicated across four landscape ecosystems broadly representative of the Upper Great Lakes ecoregion. We found that three of four examined subcanopy community weighted mean (CWM) traits including leaf photosynthetic rate (p = 0.04), stomatal conductance (p = 0.07), and the red edge normalized difference vegetation index (p < 0.0001) shifted rapidly following disturbance but before widespread changes in subcanopy light environment triggered by canopy tree mortality. Surprisingly, stimulated subcanopy production fully compensated for upper canopy losses across our gradient of experimental severities, achieving complete resistance (i.e., no significant interannual differences from control) of whole ecosystem NPP even in the 85% disturbance treatment. Additionally, we identified a probable mechanistic switch from nutrient-driven to light-driven trait shifts as disturbance progressed. Our findings suggest that remotely sensed traits such as the red edge normalized difference vegetation index (reNDVI) could be particularly sensitive and robust predictors of production response to disturbance, even across compositionally diverse forests. The potential of leaf spectral indices to predict post-disturbance functional resistance is promising given the capabilities of airborne to satellite remote sensing. We conclude that dynamic functional trait shifts following disturbance can be used to predict production response across a wide range of disturbance severities.

resistance↗
DOE OSTI2022

Contech to Accelerate Cleantech: Seeding Emerging Innovation Programs for Construction Productivity and Energy Efficiency Integration; Preprint

Investments in U.S.-based start-ups that focus on advanced building construction technologies to increase construction productivity (contech) surged to approximately $3.1 billion in 2018 (as per Crunchbase data). More recently, emerging programs by government funding agencies, philanthropic foundations, and venture capitalists have been instrumental in supporting contech start-ups for innovations that increase productivity of energy efficiency integration and accelerate clean energy technologies (cleantech) for the buildings sector. These programs include R&D support and funding mechanisms for contech and cleantech. Traditionally, contech and cleantech are considered as two different innovation ecosystems. To enhance and scale up energy efficiency in buildings, creative programs that bring together contech and cleantech are critical. This paper provides a landscape assessment of the "contech-for-cleantech" innovation ecosystem in the U.S. and its impact in accelerating technology readiness and the development pipeline of "contech-for-cleantech". Technologies highlighted are robotics for retrofits, prefabrication of energy-efficient products, and advanced manufacturing of low-carbon net-zero buildings construction. Programs discussed include those led by (1) government funding agencies: American-Made Challenges with prizes like E-ROBOT for retrofits with robotics, (2) philanthropic foundations: Wells Fargo Innovation Incubator (IN2) that includes focus on energy efficiency and prefabrication, and (3) venture capitalists: Shadow Ventures Green Building Accelerator program which provides funding support to start-ups with ambitious plans for decarbonizing the built environment. This paper will also expand upon robust processes and criteria involved in judging and down-selection of start-ups through vetting and feedback from national lab researchers and industry experts in both cleantech and contech.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Regioselective On-Surface Synthesis of [3]Triangulene Graphene Nanoribbons

The integration of low-energy states into bottom-up engineered graphene nanoribbons (GNRs) is a robust strategy for realizing materials with tailored electronic band structure for nanoelectronics. Low-energy zero-modes (ZMs) can be introduced into nanographenes (NGs) by creating an imbalance between the two sublattices of graphene. This phenomenon is exemplified by the family of [n]triangulenes (n ϵ $\mathbb{N}$). Here, we demonstrate the synthesis of [3]triangulene-GNRs, a regioregular one-dimensional (1D) chain of [3]triangulenes linked by five-membered rings. Hybridization between ZMs on adjacent [3]triangulenes leads to the emergence of a narrow band gap, E g,exp ~ 0.7 eV, and topological end states that are experimentally verified using scanning tunneling spectroscopy. Tight-binding and first-principles density functional theory calculations within the local density approximation corroborate our experimental observations. Our synthetic design takes advantage of a selective on-surface head-to-tail coupling of monomer building blocks enabling the regioselective synthesis of [3]triangulene-GNRs. Detailed ab initio theory provides insights into the mechanism of on-surface radical polymerization, revealing the pivotal role of Au-C bond formation/breakage in driving selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Nonlinear nanoelectrodynamics of a Weyl metal

Chiral Weyl fermions with linear energy-momentum dispersion in the bulk accompanied by Fermi-arc states on the surfaces prompt a host of enticing optical effects. While new Weyl semimetal materials keep emerging, the available optical probes are limited. In particular, isolating bulk and surface electrodynamics in Weyl conductors remains a challenge. We devised an approach to the problem based on near-field photocurrent imaging at the nanoscale and applied this technique to a prototypical Weyl semimetal TaIrTe 4 . As a first step, we visualized nano-photocurrent patterns in real space and demonstrated their connection to bulk nonlinear conductivity tensors through extensive modeling augmented with density functional theory calculations. Notably, our nanoscale probe gives access to not only the in-plane but also the out-of-plane electric fields so that it is feasible to interrogate all allowed nonlinear tensors including those that remained dormant in conventional far-field optics. Surface- and bulk-related nonlinear contributions are distinguished through their “symmetry fingerprints” in the photocurrent maps. Robust photocurrents also appear at mirror-symmetry breaking edges of TaIrTe 4 single crystals that we assign to nonlinear conductivity tensors forbidden in the bulk. Here, nano-photocurrent spectroscopy at the boundary reveals a strong resonance structure absent in the interior of the sample, providing evidence for elusive surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Understanding the Breakdown Behavior of Ultrawide-Bandgap Boron Nitride Power Diodes Using Device Modeling

Herein, a device study using technology computer-aided design simulation to theoretically analyze the electrical performance of ultrawide-bandgap boron nitride (BN)-based vertical junction devices is performed, including h-BN Schottky diode, h-BN pn diode, and h-BN/AlN pn diode; this is also the first demonstration of the BN power devices in simulation. The material properties of BN are defined with recently reported data, and the physical mechanisms of the device performance are systematically investigated. Additionally, the h-BN junctions in this simulation shows excellent performance, especially for breakdown behaviors. Schottky diode shows a turn-on voltage of 0.6 V for Pt Schottky contact and breakdown voltages over 450 V for 5 μm, 6 × 10 15 cm –3 p-type-doped drift layer; The h-BN pn diode shows a turn-on voltage of 6 V and breakdown voltages over 3 kV with a critical electric field of 13.6 MV cm –1 for 2.5 μm, 2 × 10 16 cm –3 p-type-doped drift layer. The h-BN/AlN heterojunction pn diode shows a turn-on voltage of 5.8 V and breakdown voltage over 2 kV for 2.5 μm, 2 × 10 16 cm –3 n-type-doped AlN drift layer. Herein, an understanding of the device principles of vertical BN junctions is provided, which can serve as a reference for the future development of robust BN power electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The essential character of the neutrino mechanism of core-collapse supernova explosions

Abstract Calibrating with detailed 2D core-collapse supernova (CCSN) simulations, we derive a simple CCSN explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behaviour of the one hundred 2D state-of-the-art models we performed to gauge its usefulness. The goal is to predict, without resort to detailed simulation, the explodability of a given massive star. We find that the simple maximum fractional ram pressure jump discriminant we define works well ∼90 per cent of the time and we speculate on the origin of the few false positives and false negatives we witness. The maximum ram pressure jump generally occurs at the time of accretion of the silicon/oxygen interface, but not always. Our results depend upon the fidelity with which the current implementation of our code F ornax adheres to Nature and issues concerning the neutrino–matter interaction, the nuclear equation of state, the possible effects of neutrino oscillations, grid resolution, the possible role of rotation and magnetic fields, and the accuracy of the numerical algorithms employed remain to be resolved. Nevertheless, the explodability condition we obtain is simple to implement, shows promise that it might be further generalized while still employing data from only the unstable Chandrasekhar progenitors, and is a more credible and robust simple explosion predictor than can currently be found in the literature.

Astronomy & Astrophysics↗
DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The role of international collaboration in the United States geologic disposal research program

Abstract. More than a decade ago, the United States disposal program discontinued all research activities focused on the unsaturated fractured tuff formation at Yucca Mountain as the geologic disposal site for spent fuel and high-level radioactive waste. A new research and development (R&D) program was initiated to provide a sound technical basis for alternative disposal options across clay, crystalline, and salt rocks. The goals of this broad program were (and still are) to (1) increase confidence in the robustness of generic disposal concepts, (2) develop the science and engineering tools needed to support disposal concept implementation, and (3) conduct R&D on the direct disposal of existing dual-purpose (storage and transportation) canisters. Recognizing the benefits of international collaboration toward the common goal of safely and efficiently managing the back end of the nuclear fuel cycle, the program emphasized international cooperation as an effective strategy for sharing information and knowledge. In a multi-laboratory effort coordinated by Lawrence Berkeley National Laboratory, the United States Department of Energy (DOE) program established formal and informal cooperation partnerships with several international initiatives and institutions and developed a number of collaborative R&D activities in important research areas, such as engineered barrier integrity, near-field perturbations, radionuclide transport, performance assessment, and methods for characterization and monitoring of engineered and natural barriers. This presentation gives an overview of these R&D activities, with a specific focus on activities that improve our current understanding of the coupled thermal–hydrological–mechanical–chemical (THMC) processes occurring in engineered and natural barriers. We start with a brief review of selected international collaboration initiatives and then describe a few specific collaboration projects. We focus specifically on such studies that use experimental data sets provided by international research cooperation for joint modeling work to increase confidence in performance-relevant predictions of coupled processes. Overall, the focus on international collaboration has allowed deep engagement of US researchers with the international waste management R&D community in terms of best practices, new scientific advances, state-of-the-art simulation tools, new monitoring and performance confirmation approaches, and lessons learned. The joint R&D with international researchers, worldwide sharing of knowledge and experience, and access to relevant data and experiments from a variety of host rocks have helped our researchers to significantly improve their understanding of the current technical basis for disposal in a range of potential host rock environments. International collaboration also provides ample opportunity for training and educating junior staff that are well suited to move the United States disposal research program forward into the next decades, a promising avenue for developing a next-generation workforce of disposal scientists.

Birkholzer, Jens T.↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2020

SO(3)-invariance of informed-graph-based deep neural network for anisotropic elastoplastic materials

This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.

42 ENGINEERING↗
DOE OSTI2022

Method to account for natural fracture induced elastic anisotropy in geomechanical characterization of shale gas reservoirs

Shale has been usually recognized as a transverse isotropic (TI) medium in conventional geomechanical log interpretation due to its laminated nature. However, when natural fractures exist in the shale rock, additional elastic anisotropy is introduced, converting laminated Shale to an orthorhombic (OB) medium. Previous studies illustrate that neglecting the natural fracture induced anisotropy in shale geomechanical log interpretation could lead to inaccurate evaluations of elastic moduli and in-situ stresses. In this paper, a new method is developed to account for the natural fracture induced anisotropy in geomechanical log interpretation based upon the TI acoustic model developed by the author and a characterization technique of elastic wave anisotropy (Sayers, 1991). The new OB model incorporates the four acoustic log data inputs and five modeling constraints in a nonlinear optimization algorithm to solve for the nine independent stiffness coefficients of an OB rock, and further to solve for the geomechanical properties and in-situ stress profiles in an OB formation. The new method was validated with a Marcellus Gas Shale field case. Both the new OB model and the conventional TI model were applied to interpret the minimum horizontal stress profile for the same formation. By comparing the results, the OB model is more robust than the TI from two aspects. First, the average stress magnitude predicted by the OB model is closer to the one measured by the Diagnostic Fracture Injection Test (DFIT). Second, the OB model predicts a more obvious stress barrier between the lower Marcellus and upper Onondaga Limestone than the TI model does. Finally, the predicted stress barrier is consistent with the observation of the microseismic events of a horizontal well drilled and completed nearby, which reveals that no hydraulic fracture propagates downward through the bottom boundary of Marcellus Shale into the underlying Onondaga Limestone.

03 NATURAL GAS↗
DOE OSTI2022

Layered complex networks as fluctuation amplifiers

Abstract In complex networked systems theory, an important question is how to evaluate the system robustness to external perturbations. With this task in mind, I investigate the propagation of noise in multi-layer networked systems. I find that, for a two layer network, noise originally injected in one layer can be strongly amplified in the other layer, depending on how well-connected are the complex networks in each layer and on how much the eigenmodes of their Laplacian matrices overlap. These results allow to predict potentially harmful conditions for the system and its sub-networks, where the level of fluctuations is important, and how to avoid them. The analytical results are illustrated numerically on various synthetic networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗