SEARCH · Engineering Papers
Results for “model types”
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.
A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282
Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.
Site heterogeneity and broad surface-binding isotherms in modern catalysis: Building intuition beyond the Sabatier principle
Learning the science of heterogeneous catalysis and electrocatalysis always starts with the simple case of a flat, uniform surface with an ideal adsorbate. It has of course been recognized for a century that real catalysts are more complicated. For the increasingly complex catalysts of the 21st century, this Perspective argues that surface heterogeneity and non-ideal binding isotherms are central features, and their implications need to be incorporated in current thinking. A variety of systems are described herein where catalyst complexity leads to broad, non-Langmuirian surface isotherms for the binding of hydrogen atoms – and this occurs even for ideal, flat Pt(111) surfaces. Modern catalysis employs nanoscale materials whose surfaces have substantial step, edge, corner, impurity, and other defect sites, and they increasingly have both metallic and non-metallic elements M n X m , including metal oxides, chalcogenides, pnictides, carbides, doped carbons, etc. The surfaces of such catalysts are often not crystal facets of the bulk phase underneath, and they typically have a variety of potential active sites. Catalytic surfaces in operando are often non-stoichiometric, amorphous, dynamic, and impure, and often vary from one part of the surface to another. Understanding of the issues that arise at such nanoscale, multi-element catalysts is just beginning to emerge. Yet these catalysts are widely discussed using Brønsted/Bell-Evans-Polanyi (BEP) relations, volcano plots, Tafel slopes, the Butler-Volmer equation, and other linear free energy relations (LFERs), which all depend on the implicit assumption that the active sites are “similar” and that surface adsorption is close to ideal. These assumptions underly the ubiquitous intuition based on the Sabatier Principle, that the fastest catalysis will occur when key intermediates have free energies of adsorption that are not too strong nor too weak. Current catalysis research often aims to minimize the complexity of non-ideal isotherms through experimental and computational design (e.g., the use of single crystal surfaces), and these studies are the foundation of the field. In contrast, this Perspective argues that the heterogeneity of binding sites and binding energies is an inherent strength of these catalysts. Here, this diversity makes many nanoscale catalysts inherently a high-throughput screen wrapped in a tiny package. Only by making the heterogeneity part of the foundation of catalysis models, sorting the types of active sites and dissecting non-ideal binding isotherms, will modern catalysis learn to harness the inherent diversity of real catalysts. Controlling and exploiting diversity rather than avoiding it will help to optimize complex modern catalysts and catalytic conditions.
Accelerated radiation transport modeling techniques for pencil beam computed tomography using gamma rays
Monte Carlo radiation transport modeling studies were performed for a compact, and high-resolution gamma-ray computed tomography system designed for imaging irradiated nuclear fuel. The system comprises a 60 Co source – chosen for its highly penetrating 1173 keV and 1332 keV gamma rays – a pair of high-aspect-ratio pencil beam collimators, and an inorganic scintillator detector. Here two acceleration methods are proposed to rapidly model a transmission type gamma-ray tomography system. The first, a variance reduction technique, is based on performing Monte Carlo simulations with a monodirectionally-biased source, sampled from a characteristic sub-volume of the full source volume. The second acceleration method is based on the deterministic calculations using the Beer–Lambert law and detector response characteristics. Comparison of simulations using acceleration approaches with analog simulations of the fully isotropic, full-volume equivalent, show that the Monte Carlo variance reduction technique gives quantitatively accurate predictions for large collimator aspect ratios while the deterministic calculations are semi-quantitative but converge close to the correct result as the collimator aspect ratio increases. As such, these techniques can be used to reduce the computational cost in generating simulated radiographs and tomographs by several orders of magnitude. Experimental validation efforts are currently underway and will be demonstrated in future work.
Numerical validation of scaling laws for stratified turbulence
Recent theoretical progress using multiscale asymptotic analysis has revealed various possible regimes of stratified turbulence. Notably, buoyancy transport can either be dominated by advection or diffusion, depending on the effective Péclet number of the flow. Two types of asymptotic models have been proposed, which yield measurably different predictions for the characteristic vertical velocity and length scale of the turbulent eddies in both diffusive and non-diffusive regimes. The first, termed a ‘single-scale model’, is designed to describe flow structures having large horizontal and small vertical scales, while the second, termed a ‘multiscale model’, additionally incorporates flow features with small horizontal scales, and reduces to the single-scale model in their absence. By comparing predicted vertical velocity scaling laws with direct numerical simulation data, we show that the multiscale model correctly captures the properties of strongly stratified turbulence within regions dominated by small-scale isotropic motions, whose volume fraction decreases as the stratification increases. Meanwhile its single-scale reduction accurately describes the more orderly, layer-like, quiescent flow outside those regions.
Effects of Temperature on Flow Fouling of Smooth and Nonwetting Surfaces
Mineral fouling (scaling) of heat transfer surfaces is a pervasive problem in heat exchangers, chemical reactors, and other equipment in energy, environmental, and process industrial applications. Many of the applications involve dynamic flow of fluids with the impurity minerals and, furthermore, operate at elevated temperature. Strategies to mitigate fouling under these conditions are of much value in industrial applications. This paper presents a comparative study of temperature-dependent mineral fouling deposition on smooth surfaces and nonwetting superhydrophobic and lubricant-infused surfaces under dynamic flow conditions. The surfaces are represented in a unified manner using the viscosity ratio of the infused material within the porous asperities on a surface to that of the flowing fluid such that the spectrum of surfaces from superhydrophobic to smooth is captured by the range of viscosity ratios from 0 to ∞. Using a forced convection experimental setup, deposition of calcium sulfate on the surfaces is quantified in terms of asymptotic fouling resistance over a range of temperature, Reynolds number, and mineral foulant supersaturation. Through a systematic set of accelerated fouling experiments, an empirical relationship for the asymptotic fouling resistance is developed in terms of Reynolds number, foulant concentration, temperature, and surface type. The empirical model is validated with a comprehensive set of experimental data from this study as well as from the literature. Optimum nonwetting surface designs for minimizing fouling resistance compared to conventional smooth surfaces are developed as a function of temperature. The results of the study offer insight into the temperature-dependent fouling of surfaces under flow conditions and a rational design of fouling-resistant nonwetting surfaces that can be readily translated to practice.
Slater transition methods for core-level electron binding energies
Methods for computing core-level ionization energies using self-consistent field (SCF) calculations are evaluated and benchmarked. These include a “full core hole” (or “ΔSCF”) approach that fully accounts for orbital relaxation upon ionization, but also methods based on Slater’s transition concept in which the binding energy is estimated from an orbital energy level that is obtained from a fractional-occupancy SCF calculation. A generalization that uses two different fractional-occupancy SCF calculations is also considered. The best of the Slater-type methods afford mean errors of 0.3–0.4 eV with respect to experiment for a dataset of K-shell ionization energies, a level of accuracy that is competitive with more expensive many-body techniques. An empirical shifting procedure with one adjustable parameter reduces the average error below 0.2 eV. Here, this shifted Slater transition method is a simple and practical way to compute core-level binding energies using only initial-state Kohn–Sham eigenvalues. It requires no more computational effort than ΔSCF and may be especially useful for simulating transient x-ray experiments where core-level spectroscopy is used to probe an excited electronic state, for which the ΔSCF approach requires a tedious state-by-state calculation of the spectrum. As an example, we use Slater-type methods to model x-ray emission spectroscopy.
Constraining the Synthesis of the Lightest 𝑝 Nucleus 74 Se
We provide the first experimental cross section of the 73 As(p, γ) 74 Se reaction to constrain one of the main destruction mechanisms of the p nucleus 74 Se in explosive stellar environments. The measurement was done using a radioactive 73 As beam at incident energies of 3.1 and 3.7 MeV/nucleon. Along with the total cross section measurement, statistical properties of the 74 Se compound nucleus were extracted, constraining the reaction cross section across the entire Gamow window of the γ process. The impact of the experimentally constrained reaction rate on 74 Se production in Type II supernovae is investigated through Monte Carlo one-zone network simulations. The results indicate that the overproduction of 74 Se by Type II supernova models cannot be resolved by nuclear physics and point towards the description of the astrophysical conditions of the γ process.
Operator dynamics in Floquet many-body systems
We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025
Search for heavy right-handed Majorana neutrinos in the decay of top quarks produced in proton-proton collisions at s = 13 TeV with the ATLAS detector
A search for heavy right-handed Majorana neutrinos is performed with the ATLAS detector at the CERN Large Hadron Collider, using the 140 fb -1 of proton–proton collision data at $\sqrt{s}$ =13 TeV collected during Run 2. This search targets $t\overline{t}$ production, in which both top quarks decay into a bottom quark and a W boson, where one of the W bosons decays hadronically and the other decays into an electron or muon and a heavy neutral lepton. The heavy neutral lepton is identified through a decay into an electron or muon and another W boson, resulting in a pair of same-charge same-flavor leptons in the final state. This paper presents the first search for heavy neutral leptons in the mass range of 15–75 GeV using $t\overline{t}$ events. No significant excess is observed over the background expectation, and upper limits are placed on the signal cross sections. Assuming a benchmark scenario of the phenomenological type-I seesaw model, these cross section limits are then translated into upper limits on the mixing parameters of the heavy Majorana neutrino with Standard Model neutrinos.
Unified origin of inflation, baryon asymmetry, and neutrino mass
In this work, we present a unified theoretical framework that simultaneously addresses some of the most intriguing puzzles in particle physics and cosmology, namely the origins of neutrino mass, baryon asymmetry, and cosmic inflation. In our model, inflation is driven by a combination of the Standard Model Higgs, the type II seesaw Higgs responsible for neutrino mass generation, and the unified symmetry-breaking Higgs field. During inflation, nonzero values of the latter field ensure the absence of the monopole problem. The baryon asymmetry is generated through the Affleck-Dine mechanism, facilitated by the nonzero angular motion in the phase of a complex scalar field, which is part of the inflaton. We find that the successful parameter region for generating baryon asymmetry through a renormalizable term in the scalar potential requires a rather heavy type II seesaw triplet, with a mass well beyond the TeV scale. Inflationary observables, in particular, the spectral index is in excellent agree with experimental observation, whereas tensor-to scalar ratio is expected to be probed by the future LiteBIRD and CMB-S4 missions. Published by the American Physical Society 2024
Reactor scale simulations of ALD and ALE: Ideal and non-ideal self-limited processes in a cylindrical and a 300 mm wafer cross-flow reactor
Here, we have developed a simulation tool to model self-limited processes such as atomic layer deposition (ALD) and atomic layer etching inside reactors of arbitrary geometry as well the output of in situ quartz crystal microbalance and mass spectrometry. We have applied this model to two standard types of cross-flow reactors: a cylindrical tube reactor and a model 300 mm wafer reactor, and explored both ideal and nonideal self-limited kinetics. The model results are in agreement with experimental results and analytic expressions obtained using a simple plug-flow model for the cylindrical tube reactor. We also extended the simulations to consider two nonideal self-limited processes: soft-saturating processes characterized by a slow reaction pathway and processes where surface by-products can compete with the precursor for the same pool of adsorption sites. Our results show that it is possible to have a self-limited process with saturated yet inhomogeneous growth profiles due to the competition of reactor by-products. This is in agreement with experimental observations for titanium dioxide ALD from titanium tetraisopropoxide and titanium tetrachloride precursors.
Search for doubly charged Higgs boson production in multi-lepton final states using 139 fb$^{-1}$ of proton–proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector
A search for pair production of doubly charged Higgs bosons ($H$ ±± ), each decaying into a pair of prompt, isolated, and highly energetic leptons with the same electric charge, is presented. The search uses a proton–proton collision data sample at a centre-of-mass energy of 13 TeV corresponding to an integrated luminosity of 139 fb -1 recorded by the ATLAS detector during Run 2 of the Large Hadron Collider (LHC). This analysis focuses on same-charge leptonic decays, $H$ ±± → $\ell$ ± $\ell$ '± where $\ell$,$\ell$ ' = $e, μ, τ$, in two-, three-, and four-lepton channels, but only considers final states which include electrons or muons. No evidence of a signal is observed. Corresponding upper limits on the production cross-section of a doubly charged Higgs boson are derived, as a function of its mass $m$($H$ ±± ), at 95% confidence level. Assuming that the branching ratios to each of the possible leptonic final states are equal, $\mathcal{B}$($H$ ±± → $e$ ± $e$ ± ) = $\mathcal{B}$($H$ ±± → $e$ ± $μ$ ± ) = $\mathcal{B}$($H$ ±± → $μ$ ± $μ$ ± ) = $\mathcal{B}$($H$ ±± → $e$ ± $τ$ ± ) = $\mathcal{B}$($H$ ±± → $μ$ ± $τ$ ± ) = $\mathcal{B}$($H$ ±± → $τ$ ± $τ$ ± ) = 1/6, the observed (expected) lower limit on the mass of a doubly charged Higgs boson is 1080 GeV (1065 GeV) within the left right symmetric type-II seesaw model, which is the strongest limit to date produced by the ATLAS Collaboration. Additionally, this paper provides the first direct test of the Zee– Babu neutrino mass model at the LHC, yielding an observed (expected) lower limit of $m$($H$ ±± ) = 900 GeV (880 GeV).
BayFlux v1.0
This is an open source software library for performing Bayesian genome scale 13C Metabolic Flux Analysis (Bayesian 13C MFA). This is a new type of metabolic modeling that quantifies flux in a probabilistic manner at the genome scale.
In vitro toxicity assessment of uranium particulates on different human lung epithelial cell models
Inhalation of uranium aerosols produced via human activities such as mining can pose a threat to human respiratory systems. Uranium oxide particulates emit short-range alpha particles that elicit DNA and direct damage, beyond associated physiochemical heavy-metal toxicity, to internal epithelial tissues. The availability of reliable in vitro models to study radiation exposure can greatly enhance our ability to understand and combat the biological impacts of exposure. However, the toxicological effects of alpha emissions and/or the oxidation states of uranium particulates vary across different human lung epithelial cell models and have not been systematically compared. We have endeavored to address this limitation by comparing impacts in three different human lung cell models: primary human bronchial and tracheal epithelial cells, primary human small airway epithelial cells, and human adenocarcinoma alveolar basal epithelial cells. Other studies have mainly investigated the toxicity of depleted uranium. Here, we compared the exposure of uranium oxide particulates (U 3 O 8 and UO 3 ) of different enrichment states on the chosen cell systems. Each cell model was exposed to 0.1, 1, 10, 50, 100, and 500 µg/mL of depleted U 3 O 8 , highly-enriched U 3 O 8 , and natural UO 3 particulates for 24 hours in submerged monolayer cultures. We compared viability and superoxide dismutase activity results across cell lines and uranium enrichment/ oxidative states. The results showed that 1) the oxide state of the particulates affected cell viability, implying that uranium’s different oxidation states contribute to different toxicological responses, and 2) each cell model reacts differently when exposed to uranium oxides, which may provide insights into the mechanistic processes associated with the exposure of radiological particulates on different biological systems. For instance, increased uranium enrichment corresponds to increased toxicity for the primary cells, but not for the immortalized cells. Our study shows that a holistic approach that incorporates similarities between model systems and types of radionuclides is required to truly develop empirical solutions for radiation exposure.
Review of Landscaping Equipment Vibration Risk Exposure, LLNL Site 200 (Various Locations), Whole Body Vibration Seated and Standing Mowers
WBV (Whole-Body Vibration) for operators of seated and standing lawn mowers was conducted during landscape activities from April 2019 through August 2019 at various locations around LLNL Site 200. Mowing is generally performed by each landscaper for 2-3 days every other week with a maximum duration of 3-5 hours each day. On mowing days, other tasks are performed for the remainder of the day that do not expose the worker to WBV such as irrigation, blowing, or sweeping. The mowing terrain at LLNL is variable with typical reported conditions such as gopher holes, exposed tree roots, quick turns, normal and high growth, weeds, and mushy/wet ground that can result in jolts to the mower operator. The Larson Davis HVM200 seat pad with tri-axial accelerometer was used for monitoring workers. For the seated mowers, the HVM200 seat pad with tri-axial accelerometer was placed on the seat between the user and seat pan of the mower. For standing mowers, the seat pad with tri-axial accelerometer was wrapped with Dycem (for seat pad protection) and placed directly on top of the mower deck and underneath the end-user’s right foot. Monitoring for WBV was conducted while operators used different types, makes, and models of mowers. A total of sixteen samples were collected from the operators of seated mower and 6 samples were collected from operators of standing mower. This initial data sampling was part of a discovery phase to identify if WBV exposure met the 2016 ACGIH TLVs (Threshold Limit Values). This initial data identified that workers operating the mowers do not exceed either the TLV or Action Level (AL) for Weighted Acceleration Levels (A(8)) but typically exceed the Vibration Dose Value (VDV) lower boundary for the standing mowers.
Parameterized Neural Ordinary Differential Equations: Applications to Computational Physics Problems
This work proposes an extension of neural ordinary differential equations (NODEs) by introducing an additional set of ODE input parameters to NODEs. This extension allows NODEs to learn multiple dynamics specified by the input parameter instances. Our extension is inspired by the concept of parameterized ordinary differential equations, which are widely investigated in computational science and engineering contexts, where characteristics of the governing equations vary over the input parameters. We apply the proposed parameterized NODEs (PNODEs) for learning latent dynamics of complex dynamical processes that arise in computational physics, which is an essential component for enabling rapid numerical simulations for time-critical physics applications. For this, we propose an encoder-decoder-type framework, which models latent dynamics as PNODEs. We demonstrate the effectiveness of PNODEs with important benchmark problems from computational physics.
Assessment of ORIGEN Reactor Library Development for Pebble-Bed Reactors Based on the PBMR-400 Benchmark
This report provides an evaluation of present SCALE capabilities for modeling depletion of pebble-bed reactor systems, using the PBMR-400 benchmark as a test case. A specific aim of this work is to understand the system characteristics required to generate production-quality reactor data libraries for rapid depletion calculations with ORIGEN. This report includes a discussion of present SCALE capabilities for modeling doubly heterogeneous fuels, prior SCALE work modeling pebble bed–type reactors, and a detailed neutronic analysis of the PBMR-400 core for both fresh and equilibrium-composition core conditions.