Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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 505 records · Page 28

Direct local parametrization of nuclear state densities using the back-shifted Bethe formula

Level densities are often parametrized using the back-shifted Bethe formula (BBF) for nuclei that possess experimental data for s-wave neutron resonance average spacings and a complete discrete level sequence at low excitation energies. However, these parametrizations require the additional modeling of the dependence of the spin-cutoff parameter on excitation energy. Here, in this work, we avoid the need to model the spin distribution of level densities by using the experimental data to parametrize directly the state densities, for which the BBF does not depend on the spin-cutoff parameter. This approach allows for a local parameterization of state densities that is independent of the spin-cutoff parameter. We provide these parameters in a tabulated form for applications in nuclear reaction calculations and for testing microscopic approaches to state densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Axion pulsarscope

Electromagnetic fields surrounding pulsars may source coherent ultralight axion signals at the known rotational frequencies of the neutron stars, which can be detected by laboratory experiments (e.g., pulsarscopes). As a promising case study, we model axion emission from the well-studied Crab pulsar, which would yield a prominent signal at 𝑓 ≈ 29.6 Hz regardless of whether the axion contributes to the dark matter abundance. We estimate the relevant sensitivity of future axion dark matter detection experiments such as DMRadio-GUT, Dark SRF, and CASPEr, assuming different magnetosphere models to bracket the uncertainty in astrophysical modeling. For example, depending on final experimental parameters, the Dark SRF experiment could probe axions with any mass 𝑚 𝑎 ≪ 10 −13 eV down to 𝑔 𝑎⁢𝛾⁢𝛾 ∼3 × 10 −13 GeV −1 with one year of data and assuming the vacuum magnetosphere model. These projected sensitivities may be degraded depending on the extent to which the magnetosphere is screened by charge-filled plasma. The promise of pulsar-sourced axions as a clean target for direct detection experiments motivates dedicated simulations of axion production in pulsar magnetospheres.

axions↗

Rates of Sea‐Level Rise Are Highly Sensitive to Ice Viscosity Parameters in Model Benchmarks

Glacier flow plays a major role in current and future rates of globally averaged sea-level rise. The viscosity of glacial ice, controlling the rate of flow, decreases as stress increases and is highly sensitive to the value of the stress exponent, $n$, in the constitutive equation for viscous flow. Glaciologists and climate modelers almost exclusively assume $n=3$ when modeling ice flow and projecting sea-level rise through forward modeling. However, recent work suggests that $n\approx 4$ better fits observations, prompting the question: How sensitive are projections of sea-level rise to the value of $n$? We use an established community ice flow model and standard benchmark experiments designed as an idealized representation of Pine Island Glacier, West Antarctica. While initializing an $n=3$ model to match observations of an $n=4$ ice sheet is possible, we find that incorrectly assuming $n=3$ when in fact $n=4$ dramatically underestimates rates of sea-level rise. The scale of this error grows nonlinearly with the magnitude of the climate forcing, acting to increase projection uncertainties. Additionally, we find that models often account for this stress-dependent rheology mismatch during model initialization in a way that masks this rheological effect in the short term while leaving model outputs vulnerable to larger biases in longer-term projections. Initializations to observations of Pine Island Glacier display similar rheology-mismatch fingerprints to our idealized example.

climate sensitivity↗

Association‐based activity coefficient model for electrolyte solutions

Abstract This work presents an association‐based activity coefficient model that explicitly considers the solution nonideality due to associations among ions and solvent species. Built upon the electrolyte nonrandom two‐liquid (eNRTL) model, the model greatly improves the accuracy of eNRTL model for strongly associating electrolyte solutions due to presence of ionic species with high surface charge density. The model successfully correlates mean ionic activity coefficients of 46 aqueous single‐salt systems from 10 cations and 5 anions at 298.15 K up to their solubility limits. With the ion‐specific association parameters identified, the model accurately predicts activity and osmotic coefficients for aqueous mixed‐salt systems at 298.15 K. The temperature dependence of the model results has also been examined at 273–373 K. With superior accuracy over a wide range of concentration and temperature, the model represents a major advancement over eNRTL model and has a great potential to be a next‐generation model for electrolyte solutions.

Lin, Yu‐Jeng↗

Generation, dynamics, and correlations of the fission fragments' angular momenta

The generation of angular momentum in fissioning nuclei is not well understood. The predictions of different models disagree, particularly concerning the correlation between the fragments' angular momenta. In this article, a time-dependent collective Hamiltonian model is proposed to treat the generation of the angular momentum in the fission fragments due to the quantum uncertainty principle as well as the dynamics of the collective wave function during and after scission. The model is constructed in the framework of the frozen Hartree-Fock approximation using a Skyrme energy functional to extract deformations of the fission fragments as well as the interactions in a derived collective Hamiltonian. The fission reactions studied are $^{240}$Pu $\rightarrow$ $^{132}$Sn+$^{108}$Ru and $^{240}$Pu $\rightarrow$ $^{144}$Ba+$^{96}$Sr. The model can account for a large part of the angular momentum found in experimental data. The orientation of the angular momentum of each fragment is found to be mainly in the plane perpendicular to the fission axis, in agreement with the experiment. The magnitudes of the angular momenta in the two fragments are nearly uncorrelated, in agreement with the recent experimental data of Wilson et al., Nature (London) 590, 566 (2021). Some of the conclusions of the traditional collective vibration model are supported by the present model but some are not. Surprisingly, it is found that the angular momenta of the fragments are slightly correlated positively as in a wriggling mode. As a result, it is also found that the presence of an octupole deformation in a fragment can significantly increase the generated angular momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for high-mass resonances in final states with a τ-lepton and missing transverse momentum with the ATLAS detector

A search for high-mass resonances decaying into a τ-lepton and a neutrino using proton-proton collisions at a center-of-mass energy of $\sqrt{s}$ =13 TeV is presented. The full run 2 data sample corresponding to an integrated luminosity of 139 fb -1 recorded by the ATLAS experiment in the years 2015–2018 is analyzed. The τ-lepton is reconstructed in its hadronic decay modes and the total transverse momentum carried out by neutrinos is inferred from the reconstructed missing transverse momentum. The search for new physics is performed on the transverse mass between the τ-lepton and the missing transverse momentum. No excess of events above the Standard Model expectation is observed and upper exclusion limits are set on the W' →τ⁢v production cross section. Heavy W' vector bosons with masses up to 5.0 TeV are excluded at 95% confidence level, assuming that they have the same couplings as the Standard Model W boson. For nonuniversal couplings, W' bosons are excluded for masses less than 3.5–5.0 TeV, depending on the model parameters. In addition, model-independent limits on the visible cross section times branching ratio are determined as a function of the lower threshold on the transverse mass of the τ-lepton and missing transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constitutive modeling of δ-phase zircaloy hydride based on strain rate dependent nanoindentation and nano-scale impact dataset

In this work, a constitutive model for zircaloy-4 embedded with δ-phase hydride rim/blister structures (δ-phase zircaloy hydride) is presented. Nanoindentation and nano-scale impact experiments were performed on the δ-hydride blisters as a function of temperature in order to obtain strain rate dependent location specific constitutive response. A power law viscoplastic constitutive model was chosen to represent the strain rate dependent behavior. The strength coefficient and the strain hardening exponent in the power law relation were found to decrease with increase in temperature indicating softening behavior at higher temperatures. The strain rate exponent remained constant with increase in temperature indicating that the rate dependent ductility does not change significantly at higher temperatures. A finite element method based model was created to analyze the failure behavior and the ductile to brittle transition in zircaloy hydride. The model predicted that fracture strength of δ-phase zircaloy hydride decreases linearly with increase in temperature.

36 MATERIALS SCIENCE↗

Spin-dependent scattering of sub-GeV dark matter: Models and constraints

We calculate the scattering rate of sub-GeV dark matter in solid-state targets for spin-dependent dark matter–nucleon interactions. For dark matter particles with mass below 100 MeV, the scattering occurs predominantly through incoherent phonon production. For dark matter heavier than 100 MeV, we match onto the nuclear recoil calculation. To compare the sensitivity of future direct detection experiments with existing constraints, we consider three models with interactions that are mediated by spin-0 or spin-1 particles. This allows us to derive bounds on the cross section from searches for the mediating particle, including bounds from stellar cooling, beam dump experiments, meson factories, and dark matter self-interactions. The existing bounds are very stringent, though for 𝑚 𝜒 ≳100 MeV there is parameter space, which may be accessible with direct detection, depending on the exposure and background rates.

Gori, Stefania [University of California, Santa Cr↗

Growth instabilities shape morphology and genetic diversity of microbial colonies

Abstract Cellular populations assume an incredible variety of shapes ranging from circular molds to irregular tumors. While we understand many of the mechanisms responsible for these spatial patterns, little is known about how the shape of a population influences its ecology and evolution. Here, we investigate this relationship in the context of microbial colonies grown on hard agar plates. This a well-studied system that exhibits a transition from smooth circular disks to more irregular and rugged shapes as either the nutrient concentration or cellular motility is decreased. Starting from a mechanistic model of colony growth, we identify two dimensionless quantities that determine how morphology and genetic diversity of the population depend on the model parameters. Our simulations further reveal that population dynamics cannot be accurately described by the commonly-used surface growth models. Instead, one has to explicitly account for the emergent growth instabilities and demographic fluctuations. Overall, our work links together environmental conditions, colony morphology, and evolution. This link is essential for a rational design of concrete, biophysical perturbations to steer evolution in the desired direction.

Golden, Alexander↗

Efficient emulation of relativistic heavy ion collisions with transfer learning

Measurements from the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) can be used to study the properties of quark-gluon plasma. Systematic constraints on these properties must combine measurements from different collision systems and methodically account for experimental and theoretical uncertainties. Such studies require a vast number of costly numerical simulations. While computationally inexpensive surrogate models (“emulators”) can be used to efficiently approximate the predictions of heavy ion simulations across a broad range of model parameters, training a reliable emulator remains a computationally expensive task. We use transfer learning to map the parameter dependencies of one model emulator onto another, leveraging similarities between different simulations of heavy ion collisions. By limiting the need for large numbers of simulations to only one of the emulators, this technique reduces the numerical cost of comprehensive uncertainty quantification when studying multiple collision systems and exploring different models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimal Client Sampling in Federated Learning with Client-level Heterogeneous Differential Privacy

Federated Learning with client-level differential privacy (DP) provides a promising framework for collaboratively training models while rigorously protecting clients’ privacy. However, classic approaches like DP-FedAvg struggle when clients have heterogeneous privacy requirements, as they must uniformly enforce the strictest privacy level across all clients, leading to excessive DP noise and significant degradation in model utility. Existing methods to improve the model utility in such heterogeneous privacy settings often assume a trusted server and are largely heuristic, resulting in suboptimal performance and lacking strong theoretical foundations. Here, in this work, we address these challenges under a practical attack model where both clients and the server are honest-but-curious. We propose GDPFed, which partitions clients into groups based on their privacy budgets and achieves client-level DP within each group to reduce the privacy budget waste and hence improve the model utility. Based on the privacy and convergence analysis of GDPFed, we find that the magnitude of DP noise depends on both model dimensionality and the per-group client sampling ratios. To further improve the performance of GDPFed, we introduce GDPFed+, which integrates model sparsification to eliminate unnecessary noise and optimizes per-group client sampling ratios to minimize convergence error. Extensive empirical evaluations on multiple benchmark datasets demonstrate the effectiveness of GDPFed+, showing substantial performance gains compared with state-of-the-art methods.

Xu, Jiahao [Univ. of Nevada, Reno, NV (United Stat↗

Natural Language Processing for Text Based Event Extraction: Identifying Events of Interest Related to Worldwide State-Sponsored Civil Nuclear Power

Beginning in FY20, SRNL was funded by the National Nuclear Security Administration’s Office of Defense Nuclear Non-Proliferation Research and Development to develop a prototype natural language processing/natural language understating machine learning-based modeling and analysis pipeline to extract and forecast events of interest from massive open data sources. The working hypothesis within the approach is that contextual shifts in key words and phrases act as indicators of events of interest over time. Therefore, by identifying points in time where contextual shifts occur, events of interest can be extracted along with explicit and implicit connections of entities and activities. The development of the preliminary prototype pipeline proved successful, meriting further testing of the pipeline on more broad topical domains and in a worldwide data environment. Therefore, SRNL, in collaboration with the Sanghani Center for Artificial Intelligence and Data Analytics at Virginia Tech, have continued development with a test case of identifying events of interest related to worldwide state-sponsored civil nuclear power in open data sources. In the first year of this follow-on effort, the team has curated domain-specific data corpuses using an automated scheme and applied the modeling and analysis pipeline. This robust, focused, and efficient approach consists of an ensemble of analyses applied to time dependent word embedding models that are trained on the data corpuses. In this report, the team has demonstrated the capability of the existing pipeline (as development has continued in parallel) by exploring several specific case-studies centered around Rosatom’s international activities regarding the planning, construction, operation, and/or shutdown of nuclear reactors. A basic timeline events has been generated by manually cataloging known “milestone” events that have occurred at reactors in Turkey, Finland, Hungary, and Egypt and compared with the output of the modeling pipeline. In this approach, the team has characterized the lead time using the prototype pipeline, as well as the ability to capture relevant information, which proved 100% successful. A deep dive example of the Akkuyu reactor (Turkey) is presented that shows the breadth of information that can be captured using the approach. In this case study, events were extracted pertaining to the planning/construction of Akkuyu including protests from the population, information campaigns in response to the protests, forged regulatory documents and lawsuits, budgetary/shareholder information, geopolitical tensions, and the various construction milestones. This has demonstrated the pipeline’s utility as a research aid or real-time event extraction tool, where summary-level information and detailed text extractions from millions of articles or Tweets across long time periods can be generated with significantly less effort than current techniques.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Molar-Mass-Dependent Partitioning of Polyethylene in Nanopores of Model Catalyst Supports from Small-Angle Neutron Scattering

Heterogeneous catalysis offers opportunities to enhance valorization of plastic waste via chemical recycling through control of the upcycled product distributions. Minimizing low-value light hydrocarbons is desired; however, fundamental insights into how to control selectivity are lacking. Here we use contrast variation with small-angle neutron scattering (SANS), model perdeuterated polyethylenes (dPEs), and a model liquid hydrocracking product (tetradecane) to quantify polymer partitioning within mesoporous silica (SBA-15). Polyethylene concentration within the mesopores is increased relative to the bulk solution, and this partitioning increases as the temperature increases. However, this polyethylene partitioning is maximized when the radius of gyration of the polymer chains is comparable to the SBA-15 pore size (10 nm). An increased partitioning at higher temperatures is attributed to entropically driven adsorption of PE within the mesopores. There is no observed preferential partitioning of hexatriacontane (a model oligomer) within the mesopores at the temperatures examined. Furthermore, these results suggest that pore size could promote the selective partitioning of polymer species into the mesopores by size. For plastic upcycling, pore-size-dependent partitioning should increase the probability for the reaction of long polymers over oligomeric and small-molecule polyolefin depolymerization products.

Adsorption↗

Examining the Impacts of Microphysical-Dynamical Feedbacks on Convective Clouds in Different Aerosol Environments Using Enhanced Observational and Modeling Strategies

This project investigated how microphysical–dynamical feedbacks influence deep convective clouds across a range of aerosol conditions, storm lifecycles, and meteorological regimes, and how these sensitivities depend on the modeling framework used to represent convection and microphysics. Using the Aerosol, Cloud, Precipitation, and Climate (ACPC) Working Group model intercomparison simulations of isolated deep convection in the Houston region, we found a robust warm-phase aerosol response in most cloud-resolving models. Increased aerosol loading tended to suppress warm-rain production, increase cloud water, and reduce surface rainfall and near-surface evaporation. These impacts were typically strongest early in the convective lifecycle, with the largest aerosol-driven differences occurring during the first half of storm evolution. The ice-phase response, on the other hand, varied widely across modeling frameworks, and differed in sign, timing, and vertical structure. Ice microphysical pathways and parameterizations are therefore leading sources of uncertainty in aerosol–deep convection interactions. Theoretical analyses further suggested that aerosol-driven invigoration through cold-phase processes was much weaker than previously hypothesized for cold-based storms, and in warm-based storms could even reduce updraft strength. These results imply that any invigoration signal is more likely linked to warm-phase processes and peaks early in storm development.

54 ENVIRONMENTAL SCIENCES↗

Mean-field model for the Curie-Weiss temperature dependence of coherence length in metallic liquids

The coherence length of the medium-range order (MRO) in metallic liquids is known to display a Curie-Weiss temperature dependence; its inverse is linearly related to temperature, and when extrapolated from temperatures above the glass transition, the coherence length diverges at a negative temperature with a critical exponent of unity. We propose a mean-field pseudospin model that explains this behavior. Specifically, we model the atoms and their local environment as Ising spins with antiferromagnetic exchange interactions. We further superimpose an exchange interaction between dynamical heterogeneities, or clusters of atoms undergoing cooperative motion. The coherence length in the metallic liquid is thus the correlation length between dynamical heterogeneities. Overall, our results reaffirm the idea that the MRO coherence length is a measure of point-to-set correlations, and that local frustrations in the interatomic interactions are prominent in metallic liquids.

36 MATERIALS SCIENCE↗

Pseudorapidity distributions of charged particles in pp($\bar p$), p(d)A and AA collisions using Tsallis thermodynamics

The pseudorapidity distributions of charged particles measured in p + p($\overline{\text{p}}$) collisions for energies ranging from $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=23.6$ GeV to 13 TeV and A + A collisions at RHIC and LHC are investigated in the fireball model with Tsallis thermodynamics. We assume that the rapidity axis is populated with fireballs following q-Gaussian distribution and the charged particles follow the Tsallis distribution in the fireball. We also extend the fireball model to asymmetric collision systems, i.e. d + Au collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=200$ GeV and p + Pb collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=5.02$ TeV, by taking into account the asymmetric geometry configuration. The model can fit the experimental data well for all the collision systems and centralities investigated. The collision energy and centrality dependence of the model parameters for the symmetric (asymmetric) collision system, i.e. the central position y 0 (y 0a , y 0A ) and its width σ (σ a , σ A ) of the fireball distribution, are also investigated and discussed. Furthermore, our results suggest that the fireball model with Tsallis thermodynamics can be used as a universal framework for the pseudorapidity distributions of charged particles in high energy collisions at RHIC and LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rapid QSTS Simulations for High-Resolution Comprehensive Assessment of Distributed PV

The rapid increase in penetration of distributed energy resources on the electric power distribution system has created a need for more comprehensive interconnection modeling and impact analysis. Unlike conventional scenario-based studies, quasi-static time-series (QSTS) simulations can realistically model time-dependent voltage controllers and the diversity of potential impacts that can occur at different times of year. However, to accurately model a distribution system with all its controllable devices, a yearlong simulation at 1-second resolution is often required, which could take conventional computers a computational time of 10 to 120 hours when an actual unbalanced distribution feeder is modeled. This computational burden is a clear limitation to the adoption of QSTS simulations in interconnection studies and for determining optimal control solutions for utility operations. The solutions we developed include accurate and computationally efficient QSTS methods that could be implemented in existing open-source and commercial software used by utilities and the development of methods to create high-resolution proxy data sets. This project demonstrated multiple pathways for speeding up the QSTS computation using new and innovative methods for advanced time-series analysis, faster power flow solvers, parallel processing of power flow solutions and circuit reduction. The target performance level for this project was achieved with year-long high-resolution time series solutions run in less than 5 minutes within an acceptable error.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗