SEARCH · Engineering Papers

Results for “Statistical simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A superconducting quantum simulator based on a photonic-bandgap metamaterial

Synthesizing many-body quantum systems with various ranges of interactions facilitates the study of quantum chaotic dynamics. Such extended interaction range can be enabled by using nonlocal degrees of freedom such as photonic modes in an otherwise locally connected structure. Here, we present a superconducting quantum simulator in which qubits are connected through an extensible photonic-bandgap metamaterial, thus realizing a one-dimensional Bose-Hubbard model with tunable hopping range and on-site interaction. Using individual site control and readout, we characterize the statistics of measurement outcomes from many-body quench dynamics, which enables in situ Hamiltonian learning. Further, the outcome statistics reveal the effect of increased hopping range, showing the predicted crossover from integrability to ergodicity. Our work enables the study of emergent randomness from chaotic many-body evolution and, more broadly, expands the accessible Hamiltonians for quantum simulation using superconducting circuits.

Science & Technology - Other Topics↗

DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

DOE OSTI2021

Impact-induced twinning in a magnesium alloy under different stress conditions

We investigate deformation twinning in a highly textured magnesium alloy plate of Mg-3.1%Al-0.9%Zn-0.4%Mn under edge-on impact with in situ, synchrotron-based, ultrafast X-ray diffraction measurements, and corresponding stress states are simulated with the finite element method. Deformation twinning and its anisotropy under the triaxial stress condition are explained by a statistical analysis of resolved shear stress. The critical resolved shear stress criterion is applicable under complicated stress conditions induced by high strain rate impact loading. Three typical simple stress conditions are further explored as verification and application cases: uniaxial-stress, uniaxial-strain, and plane-stress. Extension twinning in the magnesium alloy is prone to occur for impact loading applied perpendicular to the crystallographic c-axis, regardless of the exact stress conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A statistical model for the product energy distribution in reactions leading to prompt dissociation

Direct dynamics calculations have been performed for three reactions: C 3 H 8 + H → i-C 3 H 7 + H 2 , C 3 H 8 + H → n-C 3 H 7 + H 2 , and C 2 H 3 + O 2 → HCO + CH 2 O. The fraction of the population for the radical products that promptly dissociates is computed. Here the results for C 3 H 8 + H are qualitatively similar to previous results for C 3 H 8 + OH, but the new results exhibit a slightly higher branching fraction for prompt dissociation products, owing to the fact that a greater fraction of the internal energy in the transition state ends up in the radical. For C 2 H 3 + O 2 → HCO + CH 2 O, the fraction of HCO that promptly dissociates is in excess of 99%. Consequently, the main product for C 2 H 3 + O 2 at lower temperatures should be written as H + CO + CH 2 O and not HCO + CH 2 O. These results are then compared with four previous systems: CH 2 O + H → HCO + H 2 , CH 2 O + OH → HCO + H 2 O, C 3 H 8 + OH → i-C 3 H 7 + H 2 O, and C 3 H 8 + OH → n-C 3 H 7 + H 2 O. Based upon these seven system, several statistical models are presented. The goal of these statistical models is to predict the fraction of the transition state energy that ends up in the rovibrationally excited radical. On average, these statistical models provide an excellent prediction of product energy distribution. Consequently, these models can be used instead of costly trajectory simulations for predicting prompt radical dissociation for larger species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗

DOE OSTI2023

Adaptation Strategies Strongly Reduce the Future Impacts of Climate Change on Simulated Crop Yields

Abstract Simulations of crop yield due to climate change vary widely between models, locations, species, management strategies, and Representative Concentration Pathways (RCPs). To understand how climate and adaptation affects yield change, we developed a meta‐model based on 8703 site‐level process‐model simulations of yield with different future adaptation strategies and climate scenarios for maize, rice, wheat and soybean. We tested 10 statistical models, including some machine learning models, to predict the percentage change in projected future yield relative to the baseline period (2000–2010) as a function of explanatory variables related to adaptation strategy and climate change. We used the best model to produce global maps of yield change for the RCP4.5 scenario and identify the most influential variables affecting yield change using Shapley additive explanations. For most locations, adaptation was the most influential factor determining the projected yield change for maize, rice and wheat. Without adaptation under RCP4.5, all crops are expected to experience average global yield losses of 6%–21%. Adaptation alleviates this average projected loss by 1–13 percentage points. Maize was most responsive to adaptive practices with a projected mean yield loss of −21% [range across locations: −63%, +3.7%] without adaptation and −7.5% [range: −46%, +13%] with adaptation. For maize and rice, irrigation method and cultivar choice were the adaptation types predicted to most prevent large yield losses, respectively. When adaptation practices are applied, some areas are predicted to experience yield gains, especially at northern high latitudes. These results reveal the critical importance of implementing adequate adaptation strategies to mitigate the impact of climate change on crop yields.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Simulating sticky particles: A Monte Carlo method to sample a stratification

Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗

DOE OSTI2023

Observing signals of spectral features in the cosmic-ray positrons and electrons from Milky Way pulsars

The Alpha Magnetic Spectrometer (AMS-02) has provided unprecedented precision measurements of the electron and positron cosmic-ray fluxes and the positron fraction spectrum. At the higher energies, sources as energetic local pulsars, may contribute to both cosmic-ray species. The discreteness of the source population, can result in features both on the positron fraction measurement and in the respective electron and positron spectra. For the latter, those would coincide in energy and would contrast predictions of smooth spectra as from particle dark matter. In this work, using a library of pulsar population models for the local part of the Milky Way, we perform a power-spectrum analysis on the cosmic-ray positron fraction. We also develop a technique to cross-correlate the electron and positron fluxes. We show that both such analyses, can be used to search statistically for the presence of spectral wiggles in the cosmic-ray data. For a significant fraction of our pulsar simulations, those techniques are already sensitive enough to give a signal for the presence of those features above the regular noise, with forthcoming observations making them even more sensitive. Lastly, by cross-correlating the AMS-02 electron and positron spectra, we find an intriguing first hint for a positive correlation between them, of the kind expected by a population of local pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

UrbanScaping: Community Spatial Data Visualization & Analytics

Evaluating the electrification potential of buildings through retrofitting is crucial for reducing carbon emissions and the carbon footprint of built environments. This study leverages the Automatic Building Energy Modeling (AutoBEM) software, integrating the Model America database to create an urban context-based spatial analysis platform for community engagement and development. We selected Camp Hill Borough, PA, as a case study to analyze building-specific energy performance and evaluate the electrification potential of each building by switching to different Heating, Ventilation, and Air Conditioning (HVAC) systems and measurement components. The simulation results generated by the workflow provide retrofitting suggestions to help mitigate the carbon footprint as well as energy saving statistics of buildings. Additionally, the developed web-based interface serves as a community engagement platform, allowing residents to provide feedback and further develop interactive communication protocols. The outcomes of this project offer a baseline for community electrification planning and contribute to the design of low-carbon communities.

Chowdhury, Shovan [ORNL]↗

DOE OSTI2020

PASSAT: particle accelerator helioScopes for Slim Axion-like-particle deTection

We propose a novel method to search for axion-like particles (ALPs) at particle accelerator experiments. ALPs produced at the target via the Primakoff effect subsequently enter a region with a magnetic field, where they are converted to photons that are then detected. Dubbed Particle Accelerator helioScopes for Slim Axion-like-particle deTection (PASSAT), our proposal uses the principle of the axion helioscope but replaces ALPs produced in the Sun with those produced in a target material. Since we rely on ALP-photon conversions , our proposal probes light (slim) ALPs that are otherwise inaccessible to laboratory-based experiments which rely on ALP decay, and complements astrophysical probes that are more model-dependent. As a first application, we reinterpret existing data from the NOMAD experiment in light of PASSAT, and constrain the parameter space for ALPs lighter than $$\sim 100~\mathrm{eV}$$ ~ 100 eV and ALP-photon coupling larger than $$\sim 10^{-4}~\mathrm{GeV}^{-1}$$ ~ 10 - 4 GeV - 1 . As benchmarks of feasible low-cost experiments improving over the NOMAD limits, we study the possibility of re-using the magnets of the CAST and the proposed BabyIAXO experiments and placing them at the proposed BDF facility at CERN, together with some new detectors. We find that these realizations of PASSAT allow for a direct probe of the parameter space for ALPs lighter than $$\sim 100~\mathrm{eV}$$ ~ 100 eV and ALP-photon coupling larger than $$\sim 4\times 10^{-6}~\mathrm{GeV}^{-1}$$ ~ 4 × 10 - 6 GeV - 1 , which are regions that have not been probed yet by experiments with laboratory-produced ALPs. In contrast to other proposals aiming at detecting single or two-photon only events in hadronic beam dump environments, that rely heavily on Monte Carlo simulations, the background in our proposal can be directly measured in-situ , its suppression optimized, and the irreducible background statistically subtracted. Sensitivity evaluations with other beams will be the subject of a future paper. The measurements suggested in this paper represent an additional physics case for the BDF at CERN beyond those already proposed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

The Role of Subcloud Mesoscale Convergence in Sculpting Convective Updraft Width and Depth

The initiation of deep moist convection is governed in part by the horizontal width of updrafts near cloud base, which limits the deleterious effects of entrainment-driven dilution on buoyant thermals ascending through the free troposphere. However, the factors controlling cloud-base updraft width, which in turn dictates cloud depth, are not well understood. We track the evolving three-dimensional structure of the mesoscale subcloud forcing for vertical motion and near-cloud thermodynamic ingredients within a high-resolution ensemble of simulations of seven realistic daytime orographic convection initiation events to determine their relative roles in controlling cloud width and depth. Statistical analysis of approximately 5000 cloudy updraft samples indicates that the most important contributors to the width of cloudy updrafts across the ensemble are the depth and magnitude of the subcloud mesoscale ascent. However, the depth achieved by clouds is more consistently predicted by the near-cloud ambient relative humidity within the lower to middle free troposphere and convective available potential energy. Therefore, although the width of cloudy updrafts may be partly set at low levels by the mesoscale vertical mass and moisture flux, the likelihood of deep moist convection is governed by the generation of positive buoyancy within cumulus thermals and entrainment-driven dilution that reduces it. The persistence of the low-level mesoscale vertical forcing locally consolidates and vertically transports boundary layer moisture, helping to reduce updraft dilution. However, these factors vary in relative impacts on cloudy updrafts across individual cases, indicating multiple pathways for deep convection initiation.

Convective storms↗

DOE OSTI2026

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Detectability of Varied Hybridization Scenarios Using Genome-Scale Hybrid Detection Methods

Hybridization events complicate the accurate reconstruction of phylogenies, as they lead to patterns of genetic heritability that are unexpected under traditional, bifurcating models of species trees. This phenomenon has led to the development of methods to infer these varied hybridization events, both methods that reconstruct networks directly, as well as summary methods that predict individual hybridization events from a subset of taxa. However, a lack of empirical comparisons between methods – especially those pertaining to large networks with varied hybridization scenarios – hinders their practical use. Here, we provide a comprehensive review of popular summary methods: TICR, MSCquartets, HyDe, Patterson’s D-Statistic (ABBA-BABA), D3, and Dp. TICR and MSCquartets are based on quartet concordance factors gathered from gene tree topologies and HyDe, Patterson’s D-Statistic, D3, and Dp use site pattern frequencies to identify hybridization events between sets of three taxa. We then use simulated data to address questions of method accuracy and ideal use scenarios by testing methods against complex networks which depict gene flow events that differ in depth (timing), quantity (single vs. multiple, overlapping hybridizations), and rate of gene flow (γ). We find that deeper or multiple hybridization events may introduce noise and weaken the signal of hybridization, leading to higher relative false negative rates across all methods. Despite some forms of hybridization eluding quartet-based detection methods, MSCquartets displays high precision in most scenarios. While HyDe results in high false negative rates when tested on hybridizations involving extinct or unsampled ghost lineages, HyDe is the only method able to identify the direction of hybridization, distinguishing the source parental lineages from recipient hybrid lineages. Lastly, we test the methods on a dataset of ultraconserved elements from the bee subfamily Nomiinae, finding possible hybridization events between clades which correspond to regions of poor support in the species tree estimated in a previous study.

Bjorner, Marianne B.↗

DOE OSTI2020

The Dark Energy Spectrographic Instrument (DESI) Guide-Focus Assembly

The specific research of the DESI project is to study in detail the expansion history of the Universe over the past 10 billion years. In order to do this, the project designed, fabricated, tested and commissioned the DESI instrument which has been deployed at the Mayall Telescope at the Kitt Peak National Observatory near Tucson, Arizona. The DESI project was conducted by approximately 30 US and foreign national laboratories. Our CRADA with the The Institut de Física d’Altes Energies (IFAE) in Barcelona (Spain) involved both research planning for eventual observing on the DESI telescope; and in providing hardware jointly designed by LBNL and IFAE; ten guide and focus sensor systems (GFAs) that either provide feedback to the Mayall telescope tracking system or provide focus position information to the DESI corrector hexapod. The key major goals of this CRADA were achieved.For survey planning, Drs. Fernandez, Fosalba, Garcia-Bellido, and Gaztañaga were all important contributors to the working groups to which they were assigned. Their application of redshift techniques, development of simulations to support target selection, exploration of the use of alternative observables to extract cosmological information like AP and void statistics, and the use of PAU data to assist in the optimization of targeting selections have been vital in the early planning for targets during commissioning. The major role in the project involved Drs. Castander, Miquel, Padilla, and Sànchez. They were all directly involved with the development of the Guide-Focus Assembly (GFA) hardware, software and filter validation. In order to accomplish the goals of the DESI Project, the Mayall telescope must track DESI targeted galaxies for up to 15 minutes and the DESI corrector must be repositioned on its hexapod system for different telescope orientations. The guide and focus assemblies (GFA) provide in this information. There are ten GFA modules on the focal plane, six providing tracking information and four providing focus information. A common GFA design supports both functions by a simple switching of its operating mode. IFAE designed, fabricated, tested and delivered twelve GFAs. Once received at LBNL, the GFAs were integrated with the focal plane. This was an iterative process as the GFA’s were delivered as they were made. Our working relationship with the IFAE scientists and engineers was remarkable given the distance and time zone issues. The GFAs have performed extremely well during the commissioning phase of the instrument.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Multifidelity Approach to Sensitivity Estimation in Large-Eddy Simulation

A novel approach to compute affordable approximate sensitivities in a large-eddy simulation (LES) is proposed and assessed. The approach is based on solving a Reynolds-averaged Navier–Stokes (RANS) problem that has been linearized around the mean LES solution, with closure modeling required for the linearized changes in turbulent Reynolds stresses. In the present study, the closure modeling is based on a linearized algebraic turbulence model. The method is assessed for the flow over a NACA0012 airfoil at a fixed angle of attack, with the Reynolds number as the variable parameter. The results show that, provided an accurate linearized closure model, the method predicts the correct sensitivity of the skin friction coefficient and the mean velocity field at the cost of a linearized RANS, which provides an important proof-of-concept for this approach. Here, the linearized algebraic turbulence model with standard model coefficients produces a reasonably accurate sensitivity, but the results also suggest that accuracy could be gained from recalibrating the model coefficients for this new use of the model.

42 ENGINEERING↗