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 1,225 records · Page 68

Spin-dependent Pomeron and Odderon in elastic proton–proton scattering

We introduce a new model of near-forward elastic proton-(anti)proton scattering at high energy based on the modern formulation of Pomeron and Odderon in terms of Wilson lines and generalized TMDs (GTMDs). We compute the helicity-dependent elastic amplitudes $\phi$1,2,3,4,5 in this model and study their energy dependence from the nonlinear small-x evolution equations. While both Pomeron and Odderon contribute to helicity-flip processes in general, in the forward limit t = 0 only the double helicity-flip amplitude $\phi$2, dominated by the spin-dependent Odderon, survives. This may effect the extraction of the ρ parameter as well as the total cross section in the LHC energy domain and beyond

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scale-dependent tipping points of bacterial colonization resistance

The commonly observed stability of natural microbiomes is important for their function, yet the ubiquity of microbiome stability remains enigmatic. The strongest form of stability, colonization resistance, protects residents against invaders and is often associated with specific porous structures, such as skin follicles or intestinal crypts. By systematically probing the colonization of fly gut–derived bacteria in microfluidic pores of varying sizes, we revealed that colonization patterns and invasion rates strongly depend on the pore size. Mathematical modeling shows that bacteria spontaneously tend to organize into a dense colonization-resistant state in pores exceeding a critical size. The scale dependence of stability and resilience could bias ecological filtering in microbiomes and should be considered in the design of microbial ecology experiments.

59 BASIC BIOLOGICAL SCIENCES↗

Imprints of high-momentum nucleons in nuclei on hard photons from heavy-ion collisions near the Fermi energy

The short-range correlation (SRC) induced by the tensor force in the isosinglet neutron-proton interaction channel leads to a high-momentum tail (HMT) in the single-nucleon momentum distributions n(k) in nuclei. Owing to the remaining uncertainties about the tensor force, the shape of the nucleon HMT may be significantly different from the dilute interacting Fermi gas model prediction n(k) ~ 1/k 4 similar to the HMT in cold atoms near the unitary limit. Within an isospin- and momentum-dependent Boltzmann-Uehling-Uhlenbeck transport model incorporating approximately the nucleon HMT, we investigate hard photon emissions in 14 N + 12 C and 48 Ca + 124 Sn reactions at beam energies around the Fermi energy. Imprints of different shapes of the HMT on the energy spectrum, angular distribution and transverse momentum spectrum of hard photons are studied. Furthermore, while the angular distribution does not carry any information about the shape of the nucleon HMT, the energy spectra and especially the mid-rapidity transverse momentum spectra of hard photons are found to bare strong imprints of the shapes of nucleon HMTs in the two colliding nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evidence of Mass Ordering of Charm and Bottom Quark Energy Loss in Au+Au Collisions at RHIC

Abstract Partons traversing the strongly interacting medium produced in heavy-ion collisions are expected to lose energy depending on their color charge and mass. We measure the nuclear modification factors for charm- and bottom-decay electrons, defined as the ratio of yields, divided by the number of binary nucleon–nucleon collisions, in $$\sqrt{s_{\textrm{NN}}}=200$$ s NN = 200 GeV Au+Au collisions to p + p collisions ( $$R_{\textrm{AA}}$$ R AA ), or in central to peripheral Au+Au collisions ( $$R_{\textrm{CP}}$$ R CP ). We find the bottom-decay electron $$R_{\textrm{AA}}$$ R AA and $$R_{\textrm{CP}}$$ R CP to be significantly higher than those of charm-decay electrons. Model calculations including mass-dependent parton energy loss in a strongly coupled medium are consistent with the measured data. These observations provide evidence of mass ordering of charm and bottom quark energy loss when traversing through the strongly coupled medium created in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sneutrino tribrid inflation, metastable cosmic strings and gravitational waves

Here, we present a successful realization of sneutrino tribrid inflation model based on a gauged U(1)_B-L extension of Minimal Supersymmetric Standard Model (MSSM). A single interaction term involving the B-L Higgs field and the right-handed neutrinos serves multiple purposes. These include the generation of heavy Majorana masses for the right-handed neutrinos to provide an explanation for the tiny neutrino masses via the seesaw mechanism, a realistic scenario for reheating and non-thermal leptogenesis with a reheat temperature as low as 10 6 GeV, and a successful realization of inflation with right-handed sneutrino as the inflaton. The matter parity which helps avoid rapid proton decay survives as a Z 2 subgroup of a U(1) R-symmetry. Depending on the choice of model parameters yields the following predicted range of the tensor to scalar ratio, 3 x 10 -11 ≲ r≲ 7x 10 -4 (6 x 10 -7 ≲ r ≲ 0.01), and the running of the scalar spectral index, -0.00022 ≲ dn_s/dln k ≲ -0.0026 (-0.00014 ≲ dn_s/dln k ≲ 0.005), along with the B-L breaking scale, 3 x 10 14 ≲ M/ GeV≲ 5 x 10 15 (6 x 10 15 ≲ M/ GeV≲ 2 x 10 16 ), calculated at the central value of the scalar spectral index, n_s =0.966, reported by Planck 2018. The possibility of realizing metastable cosmic strings in a grand unified theory setup is briefly discussed. The metastable cosmic string network admits string tension values in the range 10 -8 ≲ Gμ s ≲ 10 -6 , and predicts a stochastic gravitational wave background lying within the 2-σ bounds of the recent NANOGrav 12.5-yr data.

supersymmetry and cosmology↗

Constraining models for the origin of ultra-high-energy cosmic rays with a novel combined analysis of arrival directions, spectrum, and composition data measured at the Pierre Auger Observatory

The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. The model takes into account a rigidity-dependent magnetic field blurring and an energy-dependent evolution of the catalog contribution shaped by interactions during propagation. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around 20° for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of 4.5σ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the γ-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗

Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology

we report that despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow, which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13 C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13 C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13 C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13 C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13 C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P- 13 C MOMA and P- 13 C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗

Improved Weld Residual Stress Modeling System in BlackBear

This report presents enhancements to the MOOSE-based BlackBear application aimed at improving its capability to simulate welding and other thermo-mechanical manufacturing processes. Two primary avenues of improvement are pursued. First, to enhance user accessibility, we introduce a centralized default block restriction mechanism that ensures coverage checks are performed within user-specified default blocks. This default setting is applied consistently to all block-describable objects, such as variables, kernels, and more. In addition, we develop a modular action for moving heat source simulations, which integrates path file parsing, subdomain modification, and heat source kernel enforcement into a single, streamlined configuration. Second, to improve solver robustness, we implement an alternative method for assigning initial conditions to the updated active domain during the simulation, thereby enhancing convergence behavior. To validate the framework, we design and conduct several benchmark simulations, including heat conduction with progressive material addition, linear elasticity with time-dependent material deposition, and viscoplasticity model with isotropic hardening under similar conditions. Finally, we demonstrate the effectiveness of the proposed framework through large-scale thermo-mechanical welding simulations in both two and three dimensions.

42 ENGINEERING↗

Trion Formation Resolves Observed Peak Shifts in the Optical Spectra of Transition-Metal Dichalcogenides

Monolayer transition-metal dichalcogenides (ML-TMDs) have the potential to unlock novel photonic and chemical technologies if their optoelectronic properties can be understood and controlled. Yet, recent work has offered contradictory explanations for how TMD absorption spectra change with carrier concentration, fluence, and time. Here, we test our hypothesis that the large broadening and shifting of the strong band-edge features observed in optical spectra arise from the formation of negative trions. We do this by fitting an ab initio based, many-body model to our experimental electrochemical data. Our approach provides an excellent, global description of the potential-dependent linear absorption data. We further leverage our model to demonstrate that trion formation explains the nonmonotonic potential dependence of the transient absorption spectra, including through photoinduced derivative line shapes for the trion peak. Our results motivate the continued development of theoretical methods to describe cutting-edge experiments in a physically transparent way.

36 MATERIALS SCIENCE↗

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Non-adiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (wherethe electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase- space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Furthermore, the F-PSSH algorithm sets the stage for modeling nonadiabatic dynamics under strong externally pumped circular polarization in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Light transport with weak angular dependence in fog

Random scattering and absorption of light by tiny particles in aerosols, like fog, reduce situational awareness and cause unacceptable down-time for critical systems or operations. Computationally efficient light transport models are desired for computational imaging to improve remote sensing capabilities in degraded optical environments. To this end, we have developed a model based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes, thereby covering the applicability domain of both the small angle and diffusion approximations. An analytic solution was derived and validated using experimental data acquired at the Sandia National Laboratory Fog Chamber facility. The evolution of the fog particle density and size distribution were measured and used to determine macroscopic absorption and scattering properties using Mie theory. A three-band (0.532, 1.55, and 9.68 µ m) transmissometer with lock-in amplifiers enabled changes in fog density of over an order of magnitude to be measured due to the increased transmission at higher wavelengths, covering both the moderate and highly scattering regimes. The meteorological optical range parameter is shown to be about 0.6 times the transport mean free path length, suggesting an improved physical interpretation of this parameter.

74 ATOMIC AND MOLECULAR PHYSICS↗

Determinants of cognitive and brain resilience to tau pathology: a longitudinal analysis

Mechanisms of resilience against tau pathology in individuals across the Alzheimer’s disease spectrum are insufficiently understood. Longitudinal data are necessary to reveal which factors relate to preserved cognition (i.e. cognitive resilience) and brain structure (i.e. brain resilience) despite abundant tau pathology, and to clarify whether these associations are cross-sectional or longitudinal. We used a longitudinal study design to investigate the role of several demographic, biological and brain structural factors in yielding cognitive and brain resilience to tau pathology as measured with PET. In this multicentre study, we included 366 amyloid-β-positive individuals with mild cognitive impairment or Alzheimer’s disease dementia with baseline 18 F-flortaucipir-PET and longitudinal cognitive assessments. A subset (n = 200) additionally underwent longitudinal structural MRI. We used linear mixed-effects models with global cognition and cortical thickness as dependent variables to investigate determinants of cognitive resilience and brain resilience, respectively. Models assessed whether age, sex, years of education, APOE-ε4 status, intracranial volume (and cortical thickness for cognitive resilience models) modified the association of tau pathology with cognitive decline or cortical thinning. We found that the association between higher baseline tau-PET levels (quantified in a temporal meta-region of interest) and rate of cognitive decline (measured with repeated Mini-Mental State Examination) was adversely modified by older age (Stβ interaction = -0.062, P = 0.032), higher education level (Stβ interaction = -0.072, P = 0.011) and higher intracranial volume (Stβ interaction = -0.07, P = 0.016). Younger age, higher education and greater cortical thickness were associated with better cognitive performance at baseline. Greater cortical thickness was furthermore associated with slower cognitive decline independent of tau burden. Higher education also modified the negative impact of tau-PET on cortical thinning, while older age was associated with higher baseline cortical thickness and slower rate of cortical thinning independent of tau. Our analyses revealed no (cross-sectional or longitudinal) associations for sex and APOE-ε4 status on cognition and cortical thickness. In this longitudinal study of clinically impaired individuals with underlying Alzheimer’s disease neuropathological changes, we identified education as the most robust determinant of both cognitive and brain resilience against tau pathology. The observed interaction with tau burden on cognitive decline suggests that education may be protective against cognitive decline and brain atrophy at lower levels of tau pathology, with a potential depletion of resilience resources with advancing pathology. Finally, we did not find major contributions of sex to brain nor cognitive resilience, suggesting that previous links between sex and resilience might be mainly driven by cross-sectional differences.

59 BASIC BIOLOGICAL SCIENCES↗

Time-dependent orbital-free density functional theory: Background and Pauli kernel approximations

Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. In this work, we formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{ω}{qk_F}$ and $\frac{q}{2k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TRINIDI: Time-of-Flight Resonance Imaging With Neutrons for Isotopic Density Inference

Accurate reconstruction of 2D and 3D isotope densities is a desired capability with great potential impact in applications such as evaluation and development of next-generation nuclear fuels. Neutron time-of-flight (TOF) resonance imaging offers a potential approach by exploiting the characteristic neutron absorption spectra of each isotope. However, it is a major challenge to compute quantitatively accurate images due to a variety of confounding effects such as severe Poisson noise, background scatter, beam non-uniformity, absorption non-linearity, and extended source pulse duration. We present the TRINIDI algorithm which is based on a two-step process in which we first estimate the neutron flux and background counts, and then reconstruct the areal densities of each isotope and pixel. Both components are based on the inversion of a forward model that accounts for the highly non-linear absorption, energy-dependent emission profile, and Poisson noise, while also modeling the substantial spatio-temporal variation of the background and flux. Further, to do this, we formulate the non-linear inverse problem as two optimization problems that are solved in sequence. We demonstrate on both synthetic and measured data that TRINIDI can reconstruct quantitatively accurate 2D views of isotopic areal density that can then be reconstructed into quantitatively accurate 3D volumes of isotopic volumetric density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗