Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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 73 records · Page 4

Assessments of epistemic uncertainty using Gaussian stochastic weight averaging for fluid-flow regression

Here, we use Gaussian stochastic weight averaging (SWAG) to assess the epistemic uncertainty associated with neural-network-based function approximation relevant to fluid flows. SWAG approximates a posterior Gaussian distribution of each weight, given training data, and a constant learning rate. Having access to this distribution, it is able to create multiple models with various combinations of sampled weights, which can be used to obtain ensemble predictions. The average of such an ensemble can be regarded as the 'mean estimation', whereas its standard deviation can be used to construct 'confidence intervals', which enable us to perform uncertainty quantification (UQ) with regard to the training process of neural networks. We utilize representative neural-network-based function approximation tasks for the following cases: (i) a two-dimensional circular-cylinder wake; (ii) the DayMET dataset (maximum daily temperature in North America); (iii) a three-dimensional square-cylinder wake; and (iv) urban flow, to assess the generalizability of the present idea for a wide range of complex datasets. SWAG-based UQ can be applied regardless of the network architecture, and therefore, we demonstrate the applicability of the method for two types of neural networks: (i) global field reconstruction from sparse sensors by combining convolutional neural network (CNN) and multi-layer perceptron (MLP); and (ii) far-field state estimation from sectional data with two-dimensional CNN. We find that SWAG can obtain physically-interpretable confidence-interval estimates from the perspective of epistemic uncertainty. This capability supports its use for a wide range of problems in science and engineering.

97 MATHEMATICS AND COMPUTING↗

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anderson critical metal phase in trivial states protected by average magnetic crystalline symmetry

Transitions between distinct obstructed atomic insulators (OAIs) protected by crystalline symmetries, where electrons form molecular orbitals centering away from the atom positions, must go through an intermediate metallic phase. In this work, we find that the intermediate metals will become a scale-invariant critical metal phase (CMP) under certain types of quenched disorder that respect the magnetic crystalline symmetries on average. We explicitly construct models respecting average C 2z T, m, and C 4z T and show their scale-invariance under chemical potential disorder by the finite-size scaling method. Conventional theories, such as weak anti-localization and topological phase transition, cannot explain the underlying mechanism. A quantitative mapping between lattice and network models shows that the CMP can be understood through a semi-classical percolation problem. Ultimately, we systematically classify all the OAI transitions protected by (magnetic) groups Pm, P2', P4', and P6' with and without spin-orbit coupling, most of which can support CMP.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-averaged single-molecule tertiary structures reveal RNA self-folding through individual-particle cryo-electron tomography

Large-scale and continuous conformational changes in the RNA self-folding process present significant challenges for structural studies, often requiring trade-offs between resolution and observational scope. Here, we utilize individual-particle cryo-electron tomography (IPET) to examine the post-transcriptional self-folding process of designed RNA origami 6-helix bundle with a clasp helix (6HBC). By avoiding selection, classification, averaging, or chemical fixation and optimizing cryo-ET data acquisition parameters, we reconstruct 120 three-dimensional (3D) density maps from 120 individual particles at an electron dose of no more than 168 e - Å -2 , achieving averaged resolutions ranging from 23 to 35 Å, as estimated by Fourier shell correlation (FSC) at 0.5. Each map allows us to identify distinct RNA helices and determine a unique tertiary structure. Statistical analysis of these 120 structures confirms two reported conformations and reveals a range of kinetically trapped, intermediate, and highly compacted states, demonstrating a maturation folding landscape likely driven by helix-helix compaction interactions.

36 MATERIALS SCIENCE↗

Solution-phase sample-averaged single-particle spectroscopy of quantum emitters with femtosecond resolution

Here, the development of many optical quantum technologies depends on the availability of solid-state single quantum emitters with near-perfect optical coherence. However, a standing issue that limits systematic improvement is the significant sample heterogeneity and lack of mechanistic understanding of microscopic energy flow at the single emitter level and ultrafast timescales. Here we develop solution-phase single-particle pump-probe spectroscopy with photon correlation detection that captures sample-averaged dynamics in single molecules and/or defect states with unprecedented clarity at femtosecond resolution. We apply this technique to single quantum emitters in two-dimensional hexagonal boron nitride, which suffers from significant heterogeneity and low quantum efficiency. From millisecond to nanosecond timescales, the translation diffusion, metastable-state-related bunching shoulders, rotational dynamics, and antibunching features are disentangled by their distinct photon-correlation timescales, which collectively quantify the normalized two-photon emission quantum yield. Leveraging its femtosecond resolution, spectral selectivity and ultralow noise (two orders of magnitude improvement over solid-state methods), we visualize electron-phonon coupling in the time domain at the single defect level, and discover the acceleration of polaronic formation driven by multi-electron excitation. Corroborated with results from a theoretical polaron model, we show how this translates to sample-averaged photon fidelity characterization of cascaded emission efficiency and optical decoherence time. Our work provides a framework for ultrafast spectroscopy in single emitters, molecules, or defects prone to photoluminescence intermittency and heterogeneity, opening new avenues of extreme-scale characterization and synthetic improvements for quantum information applications.

36 MATERIALS SCIENCE↗

Producing Evaluation-quality 239 Pu Average Prompt Fission Neutron Multiplicities using a Correlated Fission Model

An evaluation of the average prompt fission neutron multiplicity, $\bar{ν}_p$, of 239 Pu(n,f) is shown. This evaluation includes (a) the correlated fission model CGMF, and (b) a detailed analysis of past and recently published experimental data. Using CGMF-calculated $\bar{ν}_p$ as prior enables to link, through the use of evaluated model input parameters, $\bar{ν}_p$ to other fission observables such as the prompt fission neutron spectrum (PFNS), preneutron emission fission yields as a function of mass, and the average total kinetic energy of the fragments. These evaluated parameters produce realistic predictions of many fission observables, while the evaluated $\bar{ν}_p$ agrees well (χ 2 ≈ 1) with data. Moreover, with the new evaluated $\bar{ν}_p$, the effective neutron multiplication factor of fast Pu ICSBEP critical assemblies are predicted with a mean bias of 58 pcm compared to 18 pcm with ENDF/B-VIII.0, when paired with a new 239 Pu PFNS and fission cross section. Due to these encouraging validation results, the evaluated $\bar{ν}_p$ is currently part of a release candidate for the 239 Pu ENDF/B-VIII.1 file. Hence, a correlated fission model was used for the first time for evaluating $\bar{ν}_p$ that is of evaluation quality. This is an important step towards consistent evaluations of prompt fission observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗

Efficiency enhancement of thermonuclear DD reaction in femtosecond laser plasma with the use of structured low-average-density targets

An increase in the yield of fast neutrons is experimentally demonstrated by exciting a nuclear DD reaction in the interaction of a relativistically intense (over 10{sup 18} W cm{sup −2}) ultrashort laser pulse with a deuterated low-average-density target volume-structured at the wavelength scale. It is shown that decreasing the average target density from 0.78 to 0.35 g cm{sup −3} doubles the neutron flux, which reaches 7 × 10{sup 4} particles per 1 J of input energy. The effect may be associated with an increase in the number of accelerated deuterium ions due to the three dimensional expansion of individual elements of the target structure. (paper)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Averaged cluster approach to including chemical short-range order in KKR-CPA

The single-site Korringa-Kohn-Rostoker coherent potential approximation (KKR-CPA) ignores short-range ordering present in disordered metallic systems. Here in this paper, we establish a technique to fix this shortcoming by embedding an averaged cluster that displays chemical short-range order (SRO). The degree of SRO can be tuned by externally defined order parameters. This averaged cluster can be embedded in the single-site CPA medium or a self-consistently obtained effective medium that contains SRO information. The validity of this method is demonstrated by applying it to two alloy systems—the CuZn body-centered cubic (BCC) solid solution, and AlCrTiV, a four-element BCC high entropy alloy. A comparison between the non-self-consistent and self-consistent modes is also provided for the two above-mentioned systems. We make the code available on the internet. Planned extensions to this paper are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Uncertainty quantification of mass models using ensemble Bayesian model averaging

Developments in the description of the masses of atomic nuclei have led to various nuclear mass models that provide predictions for masses across the whole chart of nuclides. These mass models play an important role in understanding the synthesis of heavy elements in the rapid neutron capture ( r ) process. However, it is still a challenging task to estimate the size of uncertainty associated with the predictions of each mass model. In this work, a method called ensemble Bayesian model averaging (EBMA) is introduced to quantify the uncertainty of one-neutron separation energies (S 1 n ) which are directly relevant in the calculations of r -process observables. Here, this Bayesian method provides a natural way to perform model averaging, selection, and uncertainty quantification, by combining the mass models as a mixture of normal distributions whose parameters are optimized against the experimental data, employing the Markov chain Monte Carlo method using the no-u-turn sampler. The EBMA model optimized with all the experimental S 1 n from the AME2003 nuclides are shown to provide reliable uncertainty estimates when tested with the new data in the AME2020.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Average nonlinear dynamics of particles in gravitational pulses: Effective Hamiltonian, secular acceleration, and gravitational susceptibility

Particles interacting with a prescribed quasimonochromatic gravitational wave (GW) exhibit secular (average) nonlinear dynamics that can be described by Hamilton’s equations. We derive the Hamiltonian of this “ponderomotive” dynamics to the second order in the GW amplitude for a general background metric. For the special case of vacuum GWs, we show that our Hamiltonian is equivalent to that of a free particle in an effective metric, which we calculate explicitly. We also show that already a linear GW pulse displaces a particle from its unperturbed trajectory by a finite distance that is independent of the GW phase and proportional to the integral of the pulse intensity.This effect is independent from the nonlinear memory effects that has been known. We calculate the particle displacement analytically and show that our result is in agreement with numerical simulations. We also show how the Hamiltonian of the nonlinear averaged dynamics naturally leads to the concept of the linear gravitational susceptibility of a particle gas with an arbitrary phase-space distribution. We calculate this susceptibility explicitly to apply it, in a follow-up paper, toward studying self-consistent GWs in inhomogeneous media within the geometrical-optics approximation

79 ASTRONOMY AND ASTROPHYSICS↗

Improved information criteria for Bayesian model averaging in lattice field theory

Bayesian model averaging is a practical method for dealing with uncertainty due to model specification. Use of this technique requires the estimation of model probability weights. Here, we revisit the derivation of estimators for these model weights. Use of the Kullback-Leibler divergence as a starting point leads naturally to a number of alternative information criteria suitable for Bayesian model weight estimation. We explore three such criteria, known to the statistics literature before, in detail: a Bayesian analog of the Akaike information criterion which we call the BAIC, the Bayesian predictive information criterion, and the posterior predictive information criterion (PPIC). We compare the use of these information criteria in numerical analysis problems common in lattice field theory calculations. We find that the PPIC has the most appealing theoretical properties and can give the best performance in terms of model-averaging uncertainty, particularly in the presence of noisy data, while the BAIC is a simple and reliable alternative.

97 MATHEMATICS AND COMPUTING↗

Average-atom model with Siegert states

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. Furthermore, this is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Average Rényi entanglement entropy in Gaussian boson sampling

Recently, many experiments have been conducted with the goal of demonstrating a quantum advantage over classical computation. One popular framework for these experiments is Gaussian boson sampling, where quadratic photonic input states are interfered via a linear optical unitary and subsequently measured in the Fock basis. In this paper, we study the modal entanglement of the output states in this framework just before the measurement stage. Specifically, we compute Page curves as measured by various Rényi- α entropies, where the Page curve describes the entanglement between two partitioned groups of output modes averaged over all linear optical unitaries. We derive these formulas for α = 1 (i.e., the von Neumann entropy) and, more generally, for all positive integer α , in the asymptotic limit of infinite number of modes and for input states that are composed of single-mode-squeezed-vacuum state with equal squeezing strength. We then analyze the limiting behaviors when the squeezing is small and large. Having determined the averages, we then explicitly calculate the Rényi- α variance for integers α > 1 and are able to show that these entropies are weakly typical. Published by the American Physical Society 2025

Youm, Jason (ORCID:0009000057597782)↗

Quantifying uncertainty in physics-based predictions of rare-isotope production cross sections via Bayesian-inspired model averaging across nuclear mass tables

Accurate prediction of fragmentation cross sections is essential for rare-isotope beam production, planning new-isotope searches, and designing experiments to study the most exotic regions of the nuclear chart. However, existing reaction models and phenomenological cross-section parametrizations often exhibit significant deviations over broad regions of mass and charge. In this work, a Bayesian-inspired model-averaging framework is developed to combine abrasion-ablation (AA) calculations based on multiple nuclear mass tables into a single statistically weighted estimate. For the calibrated systems, the model weights are assigned empirically according to the relative quality of fit to measured cross sections, thereby reducing systematic model bias while preserving the underlying physics content of the AA description. The weights are constrained using proton-rich fragmentation data for the 78 Kr and 124 Xe projectiles. The resulting parameter trends are then propagated to the 92 Mo and 144 Sm systems through a controlled scaling procedure. In the present implementation, the excitation-energy prescription is fixed, while the averaging is performed across nuclear-mass inputs; the framework provides both weighted cross sections and associated uncertainty estimates. Applied to proton-rich fragmentation, the present approach provides a practical basis for interpolation and limited extrapolation in regions relevant to rare-isotope production. The resulting predictions are used to assess the production of very proton-rich nuclei, and candidate new isotopes are discussed.

Bayesian methods↗

Shifting and tilting towards enhanced piezoelectricity in high-temperature ceramics: an average structure study of (1 − x )BiFe 2/8 Ti 3/8 Mg 3/8 O 3 – x PbTiO 3 through X-ray and neutron diffraction

A comprehensive report is provided of the average structure of a high Curie temperature ferroelectric ceramic, (1 − x)BiFe 2/8 Ti 3/8 Mg 3/8 O 3 –xPbTiO 3 (BFTM–xPT) when x = 0.25–0.375 at room temperature. Both neutron and synchrotron X-ray powder diffraction data were collected and combined Rietveld refinements were completed. These data revealed that, of the compositions studied, BFTM–xPT crystallizes as mixed phase material. At x = 0.25–0.275, the average structure of each phase is best described by space groups R3c and Cc. At 0.30 ≥ x ≥ 0.35, a morphotropic phase boundary is present and is best modelled by space groups P4mm and Cc, where space group Cc is the dominant phase. As the concentration of PbTiO 3 increases, the amount of octahedral tilting decreases and x = 0.375 exhibits a loss of octahedral tilting and crystallizes in space groups P4mm and Cm. In conclusion, an updated phase diagram at room temperature of this novel system is given and structural intricacies responsible for its piezoelectric properties are revealed.

high Curie temperature↗