SEARCH · Engineering Papers
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.
In this report, we summarize the progress made in FY23 for experiments fielded at Triangle Universities Nuclear Laboratory (TUNL) of Duke University on the neutron scattering measurements of 238 U. The preliminary results on the elastic and inelastic reaction channels will be presented in addition to details of experiment and data analysis.
We fit standard power-law models for ground acceleration and particle velocity using scaled ground motion observations from underground nuclear explosions at Aqueduct Mesa on the Nevada National Security Site. The models are then validated using small-scale high-explosive (HE) test data from nearby Rainier Mesa. The comparison between model predictions and observations is good and can be made more favorable if the assumed yield for the HE tests are doubled. Further validation is made with a numerical experiment that additionally shows a transition in ground motion attenuation to elastic propagation (1 = r). An extended model is provided to incorporate the long-range transition in attenuation, which results in increased ground motions at farther range. We then make predictions for the field experiment, Physical Experiment One (PE1), to take place in P-tunnel on Aqueduct Mesa, and suggest that the models could be used in planning various operations for PE1.
We report the low-temperature magnetic properties of Nd 3+ pyrochlore Nd 2 ScNbO 7 . Susceptibility and magnetization show an easy-axis moment, and heat capacity reveals a phase transition to long-range order at T N =371(2) mK with a fully recovered ΔS=Rln(2), 53% of it recovered for T>T N . Elastic neutron scattering shows a long-range all-in all-out magnetic order with low-Q diffuse elastic scattering. Inelastic neutron scattering shows a low-energy flat band, indicating a magnetic Hamiltonian similar to Nd 2 Zr 2 O 7 . Nuclear hyperfine excitations measured by ultra-high-resolution neutron backscattering indicate a distribution of static electronic moments below T N , which may be due to B-site disorder influencing Nd crystal electric fields. Analysis of heat-capacity data shows an unexpected T-linear or T 3/2 term which is inconsistent with conventional magnon quasiparticles, but is consistent with fractionalized spinons or gapless local spin excitations. We use legacy data to show similar behavior in Nd 2 Zr 2 O 7 . Comparing local static moments also reveals a suppression of the nuclear Schottky anomaly in temperature, evidencing a fraction of Nd sites with nearly zero static moment, consistent with exchange-disorder-induced random singlet formation. Finally, taken together, these measurements suggest an unusual fluctuating magnetic ground state which mimics a spin liquid, but may not actually be one.
Abstract The calibration of solid constitutive models with full‐field experimental data is a long‐standing challenge, especially in materials that undergo large deformations. In this paper, we propose a physics‐informed deep‐learning framework for the discovery of hyperelastic constitutive model parameterizations given full‐field surface displacement data and global force‐displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic material models suitable for different material classes is considered including the Neo–Hookean, Gent, and Blatz–Ko constitutive models as exemplars for general non‐linear elastic behaviour, elastomer behaviour with finite strain lock‐up, and compressible foam behaviour, respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full‐field experimental data are utilized to calibrate constitutive models under finite deformations.
Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).
While the contour method for residual stress assessment has developed rapidly, no published study documents its interlaboratory reproducibility. Objective: Here we report an initial reproducibility experiment focused on contour method data analysis and residual stress calculation. The experiment uses surface topography data from a physical process simulation of elastic-plastic beam bending. The simulation provides surface topography, for input to the contour method data analysis, as well as a known residual stress field with 130 MPa peak magnitude. Additionally, to increase realism, noise and specific artifacts are added to the topography data. A group of participants received the topography data (without the known residual stress), independently analyzed the data, and submitted results as a two-dimensional residual stress field. Analysis of submissions provides a group average residual stress field and the spatial distribution of reproducibility standard deviation. The group average residual stress agrees with the known stress in magnitude and spatial trend. The reproducibility standard deviation ranges from 2 to 54 MPa over the measurement plane, with an average of 5.4 MPa. Reproducibility standard deviation is smaller in the cross-section interior (≤ 5 MPa), modest near local extrema in the stress field (5 to 10 MPa), and larger near the cross-section boundaries (10 to 30 MPa). Overall, the largest values of reproducibility standard deviation (up to 54 MPa) occur in limited areas where artifacts had been added to the topography data; while some participants identified and removed these artifacts, some did not, leading to systematic differences that elevated the standard deviation.
A bridgmanite/periclase aggregate is investigated due to its prevalence in the Earth’s lower mantle and its importance for understanding geodynamic processes. In dual-phase polycrystalline aggregates, both strength contrast between phases and single crystal elastic and plastic anisotropy are known to influence the development of intragrain misorientation, and the evolution of crystallographic texture. In this study, full-field crystal plasticity simulations are performed with a finite element solution method, and applied to an aggregate of a mechanically strong orthorhombic phase with perovskite structure (bridgmanite, MgSiO 3 ), and a relatively weaker cubic phase (periclase, MgO). The relative strengths of the phases are parameterized to elucidate the effect that strength contrast has on texture evolution and single-phase simulations are performed for comparison. Altogether, results indicate that the relative strength between the two phases influences the development of plasticity, and the overall texture evolution. Results are discussed in light of trends related to the evolution of plasticity and misorientation in the aggregate, and their dependence on both the strength contrast between phases as well as the spatial distribution of grains and phases.
Abstract YbPtBi is one of the heavy-fermion systems with largest Sommerfeld coefficient γ and is thus classified as a ‘super’-heavy fermion material. In this work, we resolve the long-debated question about the hierarchy of relevant energy scales, such as crystal-electric field (CEF) levels, Kondo and magnetic ordering temperature, in YbPtBi. Through measurements of the a.c. elastocaloric effect and generic symmetry arguments, we identify an elastic level splitting that is unambiguously associated with the symmetry-allowed splitting of a quartet CEF level. This quartet, which we identify to be the first excited state at Δ/ k B ≈ 1.6 K above the doublet ground state at ambient pressure, is well below the proposed Kondo temperature T K ≈ 10 K. Consequently, this analysis of the energy scheme can provide support models that predict that the heavy electron mass is a result of an enhanced degeneracy of the CEF ground state, i.e., a quasi-sextet in YbPtBi. At the same time, our study shows the potential of the a.c. elastocaloric effect to control and quantify strain-induced changes of the CEF schemes, opening a different route to disentangle the CEF energy scales from other relevant energy scales in correlated quantum materials.
Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.
Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.
In this article, we describe the seismoacoustic wavefield recorded outdoors but inside the facility fence of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (Tennessee). HFIR is a research nuclear reactor that generates neutrons for scattering, irradiation research, and isotope production. This reactor operates at a nominal power of 85 MW, with a full-power period between 24 and 26 days. This study uses data from a single seismoacoustic station that operated for 60 days and sampled a full operating reactor cycle, that is, full-power operation and end-of-cycle outage. The analysis presented here is based on identifying signals that characterize the steady, that is, full-power operation and end-of-cycle outage, and transitional, that is, start-up and shutdown, states of the reactor. We found that the overall seismoacoustic energy closely follows the main power cycle of the reactor and identified spectral regions excited by specific reactor operational conditions. In particular, we identified a tonal noise sequence with a fundamental frequency around 21.4 Hz and multiple harmonics that emerge as the reactor reaches 90% of nominal power in both seismic and acoustic channels. We also utilized temperature measurements from the monitoring system of the reactor to suggest links between the operation of reactor’s subsystems and seismoacoustic signals. We demonstrate that seismoacoustic monitoring of an industrial facility can identify and track some industrial processes and detect events related to operations that involve energy transport.
Neutrino oscillations are so far the only experimental observation beyond the standard model since its development four decades ago. The remaining fundamental questions that can be answered by long-baseline neutrino oscillations are (1) the CP violation (whether neutrinos and antineutrinos behave the same way in oscillation), and (2) the mass hierarchy ($m_3>m_{1,2}$ or $m_{1,2}>m_3$), and (3) the octant of $\theta_{23}$ ($\theta_{23}>45^\circ, =45^\circ$ or $<45^\circ$). The PI's group at UCI focuses on long-baseline neutrino experiments DUNE and NOvA, key components of DOE's Intensity Frontier program. NOvA is the current major U.S.-based long-baseline neutrino experiment aiming to solve the mass hierarchy and $\theta_{23}$ octant. It is taking $\nu$ and $\bar{\nu}$ data from the NuMI beam at Fermilab. DUNE is the next-generation flagship neutrino experiment in the U.S. which is designed to decisively determine neutrino CP violation, mass hierarchy and $\theta_{23}$ octant. DUNE is Fermilab's first priority, and has been established as an international collaboration. During the period covered by the report, the PI's group has been exceptionally productive, with leading contributions to DUNE argon-purity monitoring detectors, NOvA/DUNE deep-learning reconstruction, and NOvA oscillation and near detector (ND) analyses. The group also is actively involved in DUNE/NOvA detector operation, data production, calibration, software upgrade and maintenance, and ProtoDUNE cold electronic QA/QC. Specifically, we have made accomplishments in the following four research projects: 1) Purity monitor R\&D and data analysis for DUNE, 2) Deep-learning-based energy reconstruction at NOvA and DUNE, 3) Oscillation analysis at NOvA, and 4) $\nu$-e elastic scattering analysis at NOvA.
Neutrinos continue to provide a testing ground for the structure of the standard model of particle physics as well as hints towards the physics beyond the standard model. Neutrinos of energies spanning over several orders of magnitude, originating in many terrestrial and astrophysical processes, have been detected via various decay and interaction mechanisms. At MeV scales, there has been one elusive process, until a few years ago, known as coherent elastic neutrino-nucleus scattering (CEvNS) that was theoretically predicted over five decades ago but was never observed experimentally. The recent experimental observation of the CEvNS process by the COHERENT collaboration at a stopped pion neutrino source has inspired physicists across many subfields. This has vital implications for nuclear physics, high-energy physics, astrophysics, and beyond. CEvNS, being a low-energy process, provides a natural window to study light, weakly-coupled, new physics in the neutrino sector. In this review, we intend to provide the current status of low energy neutrino scattering physics and its implications for the standard and beyond the standard model physics. We discuss the general formalism of calculating the tree-level CEvNS cross section and present estimated theoretical uncertainties on the CEvNS cross section stemming from different sources. We also discuss the inelastic scattering of tens of MeV neutrinos that have implications for supernova detection in future neutrino experiments. We discuss how the CEvNS experiments can be used as a testing ground for the Standard Model (SM) weak physics as well as in searching for the Beyond the Standard Model (BSM) physics signals. Any deviation from the SM predicted event rate either with a change in the total event rate or with a change in the shape of the recoil spectrum, could indicate new contributions to the interaction cross-section. The SM implications include the study of weak nuclear form factor and weak mixing angle. The BSM studies include non-standard interactions, neutrino electromagnetic properties, and sterile neutrino searches. In conclusion, stopped pion facilities are also a copious source of neutral and changed mesons that allow study of several dark sector physics scenarios such as vector portal models, leptophobic dark matter as well as axion-like particle searches.
Here, we present the analysis and results of the first dataset collected with the MARS neutron detector deployed at the Oak Ridge National Laboratory Spallation Neutron Source (SNS) for the purpose of monitoring and characterizing the beam-related neutron (BRN) background for the COHERENT collaboration. MARS was positioned next to the COH-CsI coherent elastic neutrino-nucleus scattering detector in the SNS basement corridor. This is the basement location of closest proximity to the SNS target and thus, of highest neutrino flux, but it is also well shielded from the BRN flux by infill concrete and gravel. Furthermore, these data show the detector registered roughly one BRN per day. Using MARS' measured detection efficiency, the incoming BRN flux is estimated to be 1.20 ± 0.56 neutrons/m 2 /MWh for neutron energies above ~3.5 MeV and up to a few tens of MeV. We compare our results with previous BRN measurements in the SNS basement corridor reported by other neutron detectors.
High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.
As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.
The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.
In the era of big data, materials science workflows need to handle large-scale data distribution, storage, and computation. Any of these areas can become a performance bottleneck. We present a framework for analyzing internal material structures (e.g., cracks) to mitigate these bottlenecks. We demonstrate the effectiveness of our framework for a workflow performing synchrotron X-ray computed tomography reconstruction and segmentation of a silica-based structure. Our framework provides a cloud-based, cutting-edge solution to challenges such as growing intermediate and output data and heavy resource demands during image reconstruction and segmentation. Specifically, our framework efficiently manages data storage, scaling up compute resources on the cloud. The multi-layer software structure of our framework includes three layers. A top layer uses Jupyter notebooks and serves as the user interface. A middle layer uses Ansible for resource deployment and managing the execution environment. A low layer is dedicated to resource management and provides resource management and job scheduling on heterogeneous nodes (i.e., GPU and CPU). At the core of this layer, Kubernetes supports resource management, and Dask enables large-scale job scheduling for heterogeneous resources. The broader impact of our work is four-fold: through our framework, we hide the complexity of the cloud’s software stack to the user who otherwise is required to have expertise in cloud technologies; we manage job scheduling efficiently and in a scalable manner; we enable resource elasticity and workflow orchestration at a large scale; and we facilitate moving the study of nonporous structures, which has wide applications in engineering and scientific fields, to the cloud. While we demonstrate the capability of our framework for a specific materials science application, it can be adapted for other applications and domains because of its modular, multi-layer architecture.