Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic fields”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Enhanced relaxed physical factorization preconditioner for coupled poromechanics

The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFF T ), where C is a symmetric positive definite matrix, FF T a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Disentangling Coexisting Structural Order Through Phase Lock-In Analysis of Atomic-Resolution STEM Data

As a real-space technique, atomic-resolution STEM imaging contains both amplitude and geometric phase information about structural order in materials, with the latter encoding important information about local variations and heterogeneities present in crystalline lattices. Such phase information can be extracted using geometric phase analysis (GPA), a method which has generally focused on spatially mapping elastic strain. Here we demonstrate an alternative phase demodulation technique and its application to reveal complex structural phenomena in correlated quantum materials. As with other methods of image phase analysis, the phase lock-in approach can be implemented to extract detailed information about structural order and disorder, including dislocations and compound defects in crystals. Extending the application of this phase analysis to Fourier components that encode periodic modulations of the crystalline lattice, such as superlattice or secondary frequency peaks, we extract the behavior of multiple distinct order parameters within the same image, yielding insights into not only the crystalline heterogeneity but also subtle emergent order parameters such as antipolar displacements. When applied to atomic-resolution images spanning large (~0.5 × 0.5 μm 2 ) fields of view, this approach enables vivid visualizations of the spatial interplay between various structural orders in novel materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2022

Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators

Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. This white paper discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Geometric compatibility measure m' for twin transmission: A predictor or descriptor?

In this work, the geometric compatibility factor m' is critically analyzed to assess whether it can be used to interpret/predict twin transmission (TT) across grain boundaries (GBs). This geometric measure is widely used to relate the likelihood of TT to the misalignment of both the shear and plane-normal directions within a twin set (i.e., incoming and outgoing twin). Here, using a large set of electron back scattering diffraction (EBSD) data, a detailed statistical analysis of twin-GB interactions is performed for {${1\bar{01}}2$} tensile twins in hexagonal close-packed (HCP) metals Mg, Zr, and Ti at different strain levels. In addition, a full-field crystal plasticity model is employed to quantify the role of local stresses and the applicability of m' as a criterion for the TT process. This combined study addresses the following three main questions: (i) What is the fidelity of m' in describing experimentally observed TTs? (ii) Can m' be used as a metric to predict/anticipate TT? (iii) Does m' naturally capture local stress effects? As a descriptor, m' cannot rationalize ~25% of TT events observed in Mg or more than 50% of TT events in Zr and Ti. As a predictor, the m'-measure does not predict TT events in over ~50% of twin-GB interactions analyzed. Further, the applicability of m' to describe and predict TT events decreases with an increase in elastic anisotropy, plastic anisotropy, and macroscopic strain levels. Finally, the twinning simulations reveal that m' does not capture the key effects of local stresses on variant selection upon twin transmission. The local stress induced by the twinning shear transformation plays a dominant role in driving the TT process compared to the geometric alignment of the constituting twins, i.e., m'.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimizing energetic light dark matter searches in dark matter and neutrino experiments

Neutrino and dark matter experiments with large-volume ($ \underset{\sim }{>} $,2 1 ton) detectors can provide excellent sensitivity to signals induced by energetic light dark matter coming from the present universe. Taking boosted dark matter as a concrete example of energetic light dark matter, we scrutinize two representative search channels, electron scattering and proton scattering including deep inelastic scattering processes, in the context of elastic and inelastic boosted dark matter, in a completely detector-independent manner. In this work, a dark gauge boson is adopted as the particle to mediate the interactions between the Standard Model particles and boosted dark matter. We find that the signal sensitivity of the two channels highly depends on the (mass-)parameter region to probe, so search strategies and channels should be designed sensibly especially at the earlier stage of experiments. In particular, the contribution from the boosted-dark-matter-initiated deep inelastic scattering can be subleading (important) compared to the quasi-elastic proton scattering, if the mass of the mediator is below (above) $ \mathcal{O} $(GeV). We demonstrate how to practically perform searches and relevant analyses, employing example detectors such as DarkSide-20k, DUNE, Hyper-Kamiokande, and DeepCore, with their respective detector specifications taken into consideration. For other potential detectors we provide a summary table, collecting relevant information, from which similar studies can be fulfilled readily.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Mechanical response and microstructure evolution from multiaxial stress relaxation in textured zircaloy nuclear cladding

Modern energy markets provide motivation for nuclear reactors to increase operational flexibility and load following capacity of existing reactor designs. The operational power changes result in fuel pellet expansion and contraction resulting in pellet-cladding mechanical interaction (PCMI) and high cladding stresses that are strongly related to stressed corrosion cracking (SCC) which can lead to fuel rod failure. This work aims to quantify the stress fields in cladding with a focus on multiaxial stress relaxation response from an imposed strain and the resulting microstructure changes in textured zircaloy claddings. Stress relaxation tests resulting in multiaxial stress fields are performed to measure the material response of Zircaloy-4 nuclear cladding. Multiple experiments using samples from the same piece of cladding are performed so samples with the same processing history may be removed at strategic points throughout mechanical testing for microstructure characterization. The microstructure evolution is analyzed with mechanical response to evaluate the static recovery impact on crystallographic texture and grain morphology. Experiments are performed at Idaho National Laboratory using an experimental device capable of applying internal pressure and independent axial load to commercial nuclear cladding at elevated temperatures with stress- or strain-controlled experimental capability. X-ray diffraction and electron backscatter diffraction are used for microstructure analysis. Rheological models are used to visualize the multiaxial stress-strain relationships and isolate elastic, plastic, and viscoplastic material properties for analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Using Deep Learning Techniques to Search for the MiniBooNE Low Energy Excess in MicroBooNE with > 3$\sigma$ Sensitivity

This thesis describes an analysis developed for the MicroBooNE experiment to investigate an anomalous excess of electron-like events observed in the MiniBooNE detector. The hypothesis investigated here is that the MiniBooNE anomaly represents appearance of electron neutrinos. Using an amalgam of novel Deep Learning and standard algorithmic techniques this analysis reconstructs and identifies a highly pure sample of charged current quasi-elastic muon neutrino and electron neutrino interactions. This thesis describes the steps in the analysis chain and provides data-to-simulation comparisons for each step that establish confidence in the final prediction. When interpreted in the context of a $\nu e$ appearance like model, this analysis predicts a 3.2$\sigma$ sensitivity to exclude a standard model fluctuation which would appear as a MiniBooNE like anomaly using $7\times10^{20}$ protons on target of MicroBooNE Data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Distinguishing Charged Lepton Flavor Violation Scenarios with Inelastic 𝜇 → 𝑒 Conversion

The Mu2e and COMET experiments are expected to improve existing limits on charged lepton flavor violation (CLFV) by roughly 4 orders of magnitude. 𝜇 → 𝑒 conversion experiments are typically optimized for electrons produced without nuclear excitation, as this maximizes the electron energy and minimizes backgrounds from the free decay of the muon. Here we argue that Mu2e and COMET will be able to extract additional constraints on CLFV from inelastic 𝜇 → 𝑒 conversion, given the 27 Al target they have chosen and backgrounds they anticipate. We describe CLFV scenarios in which inelastic CLFV can induce measurable distortions in the near-endpoint spectrum of conversion electrons, including cases where certain contributing operators cannot be probed in elastic 𝜇 → 𝑒 conversion. We extend the nonrelativistic EFT treatment of elastic 𝜇 → 𝑒 conversion to include the new nuclear operators needed for the inelastic process, evaluate the associated nuclear response functions, and describe several new-physics scenarios where the inelastic process can provide additional information on CLFV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Unavoidable multilevel biaxial symmetry breaking in chiral hybrid liquid crystals

Chiral nematic or cholesteric liquid crystals (LCs) are fluid mesophases with long-ranged orientational order featuring a quasilayered periodicity imparted by a helical director configuration but lacking long-range positional order. Doping molecular cholesteric LCs with strongly anisotropic uniaxial colloidal particles adds another level of complexity because of the interplay between weak surface-anchoring boundary conditions and bulk-based elastic distortions near the particle-LC interface. Using cylindrical colloidal disks and rods with different geometric shapes and surface conditions, we demonstrate that these colloidal inclusions generically exhibit biaxial orientational probability distributions which may impart anomalously strong local biaxiality onto the hybrid cholesteric LC structure. Unlike nonchiral hybrid molecular-colloidal LCs, where biaxial order emerges only at critical colloid volume fractions exceeding some uniaxial-biaxial transition value, the orientational probability of the colloidal inclusions immersed in chiral nematic hosts is unambiguously biaxial even at infinite dilution. We demonstrate that the colloids induce local biaxial perturbations within the molecular orientational order of the LC host medium which strongly enhances the weak but native biaxial order of chiral nematic LC induced by the chiral symmetry breaking of the director field. With the help of analytical modeling and computer simulations based on the Landau–de Gennes free energy of the host LC around the colloids, we rationalize the observed multilevel biaxial order and conclude that it is not only unavoidable but also strongly enhanced compared to both achiral hybrid LCs and purely molecular cholesteric LCs.

biaxial liquid crystals↗
DOE OSTI2020

Unified model of nucleon elastic form factors and implications for neutrino-oscillation experiments

Precise knowledge of the nucleon’s axial-current form factors is crucial for modeling GeV-scale neutrino-nucleus interactions. Unfortunately, the axial form factor remains insufficiently constrained to meet the precision requirements of upcoming long-baseline neutrino-oscillation experiments. This work studies the nucleon’s axial and vector form factors using the light-front approach to build a quark-diquark model of the nucleon with an explicit pion cloud. The light-front wave functions in both the quark and pion-baryon Fock spaces are first calibrated to existing experimental information on the nucleon’s electromagnetic form factors and then used to predict the axial form factor. The resulting squared charge radius of the axial pseudovector form factor is predicted to be $r$$^{2}_{A}$ = 0.29 ± 0.03 fm 2 , where the small error accounts for the model’s parametric uncertainty. We use our form factor results to explore the (quasi)elastic scattering of neutrinos by (nuclei)nucleons, with the result that the widely implemented dipole ansatz is an inadequate approximation of the full form factor for modeling both processes. Here, the approximation leads to a 5%–10% overestimation of the total cross section, depending on the (anti)neutrino energy. We project overestimations of similar size in the flux-averaged cross sections for the upcoming DUNE long-baseline neutrino-oscillation experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Non-Equilibrium relaxation, aging scaling, and critical depinning dynamics of skyrmions in disordered magnetic films. Final Report

Technical applications of type-II superconductors, especially high-Tc compounds, in external magnetic fields require an effective flux pinning mechanism to reduce Ohmic losses due to flux creep and flow. Naturally abundant point defects, dislocations, and twin boundaries as well as artificially introduced disorder may serve as effective pinning centers for magnetic vortices. The properties of interacting flux lines subject to strong thermal fluctuations and point-like or extended pinning centers have therefore been a major research focus in condensed matter physics, and considerable progress has been achieved in the characterization and mathematical description of the equilibrium phases and non-equilibrium steady states of vortex matter. Detailed studies of the rich phase diagrams that reveal novel glassy states as well as the associated complex dynamical phenomena have in turn enriched the theoretical understanding of condensed matter systems, elucidating specifically the kinetics of elastic objects driven through disordered media and the important role of topological defects such as vortices and dislocations in various materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Extended isogeometric analysis of multi-material and multi-physics problems using hierarchical B-splines

Here, this paper presents an immersed, isogeometric finite element framework to predict the response of multi-material, multi-physics problems with complex geometries using locally refined discretizations. To circumvent the need to generate conformal meshes, this work uses an extended finite element method (XFEM) to discretize the governing equations on non-conforming, embedding meshes. A flexible approach to create truncated hierarchical B-splines discretizations is presented. This approach enables the refinement of each state variable field individually to meet field-specific accuracy requirements. To obtain an immersed geometry representation that is consistent across all hierarchically refined B-spline discretizations, the geometry is immersed into a single mesh, the XFEM background mesh, which is constructed from the union of all hierarchical B-spline meshes. An extraction operator is introduced to represent the truncated hierarchical B-spline bases in terms of Lagrange shape functions on the XFEM background mesh without loss of accuracy. The truncated hierarchical B-spline bases are enriched using a generalized Heaviside enrichment strategy to accommodate small geometric features and multi-material problems. The governing equations are augmented by a formulation of the face-oriented ghost stabilization enhanced for locally refined B-spline bases. We present examples for two- and three-dimensional linear elastic and thermo-elastic problems. The numerical results validate the accuracy of our framework. The results also demonstrate the applicability of the proposed framework to large, geometrically complex problems.

42 ENGINEERING↗
DOE OSTI2023

Externally Heated Diamond ANvil Cell Experimentation (EH-DANCE) for studying materials and processes under extreme conditions

Externally heated diamond anvil cells provide a stable and uniform thermal environment, making them a versatile device to simultaneously generate high-pressure and high-temperature conditions in various fields of research, such as condensed matter physics, materials science, chemistry, and geosciences. The present study features the Externally Heated Diamond ANvil Cell Experimentation (EH-DANCE) system, a versatile configuration consisting of a diamond anvil cell with a customized microheater for stable resistive heating, bidirectional pressure control facilitated by compression and decompression membranes, and a water-cooled enclosure suitable for vacuum and controlled atmospheres. This integrated system excels with its precise control of both pressure and temperature for mineral and materials science research under extreme conditions. Here, we showcase the capabilities of the system through its successful application in the investigation of the melting temperature and thermal equation of state of high-pressure ice-VII at temperatures up to 1400 K. The system was also used to measure the elastic properties of solid ice-VII and liquid H 2 O using Brillouin scattering and Raman spectra of carbonates using Raman spectroscopy, highlighting the potential of the EH-DANCE system in high-pressure research.

Instruments & Instrumentation↗
DOE OSTI2026

Structural characterization of neutron irradiated hexagonal boron-10 nitride-15 single crystals

The negatively charged boron vacancy ($V$$^{–}_{⁠B}$⁠) in hexagonal boron nitride (hBN) is a promising quantum defect that can be used to sense pressure, temperature, and magnetic field with high spatial resolution. hBN enriched with the boron-10 and nitrogen-15 isotopes, denoted h 10 B 15 N, has good contrast and coherence for quantum sensing because nitrogen-15 has nuclear spin (1/2), reducing hyperfine interactions. Boron vacancies can be generated by neutron irradiation, which causes the transmutation of the boron-10 isotope to lithium-7. In this study, the ancillary structural, compositional, and mechanical properties of h 10 B 15 N crystals that have been subjected to neutron irradiation fluences from 1.4 × 10 16 to 8.4 × 10 17 n/cm 2 were thoroughly characterized. Besides creating $V$$^{–}_{⁠B}$⁠, the process also induces other defects that generate strain in the crystal lattice. In turn, the mechanical properties of these crystals change drastically. Investigated here are the visual changes, lattice integrity, composition, crystal strain, and elastic constants (C 33 and C 66 ) to assess how these characteristics change as a function of neutron fluence.

36 MATERIALS SCIENCE↗