Engineering Papers⌕ Search

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.

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

At least 1,333 records · Page 74

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗

How Climate and Data Quality Impact Photovoltaic Performance Loss Rate Estimations

Different data pipelines and statistical methods are applied to photovoltaic (PV) performance datasets to quantify the performance loss rate (PLR). Since the real values of PLR are unknown, a variety of unvalidated values are reported. As such, the PV industry commonly assumes PLR based on statistically extracted ranges from the literature. However, the accuracy and uncertainty of PLR depend on several parameters including seasonality, local climatic conditions, and the response of a particular PV technology. In addition, the specific data pipeline and statistical method used affect the accuracy and uncertainty. To provide insights, a framework of (≈200 million) synthetic simulations of PV performance datasets using data from different climates is developed. Time series with known PLR and data quality are synthesized, and large parametric studies are conducted to examine the accuracy and uncertainty of different statistical approaches over the contiguous US, with an emphasis on the publicly available and “standardized” library, RdTools . In the results, it is confirmed that PLRs from RdTools are unbiased on average, but the accuracy and uncertainty of individual PLR estimates vary with climate zone, data quality, PV technology, and choice of analysis workflow. Best practices and improvement recommendations based on the findings of this study are provided.

14 SOLAR ENERGY↗

Improved thermonuclear rate of 42 Ti( p , γ ) 43 V and its astrophysical implication in the rp process

Context. Accurate 42 Ti(p,γ) 43 V reaction rates are crucial for understanding the nucleosynthesis path of the rapid capture process (rp process) that occurs in X-ray bursts. Aims. We aim to improve the thermonuclear rates of 42 Ti(p,γ) 43 V based on more complete resonance information and a more accurate direct component, together with the recently released nuclear masses data. We also explore the impact of the newly obtained rates on the rp process. Methods. We reevaluated the reaction rate of 42 Ti(p,γ) 43 V by the sum of the isolated resonance contribution instead of the Hauser-Feshbach statistical model. We used a Monte Carlo method to derive the associated uncertainties of new rates. The nucleosynthesis simulations were performed via the NuGrid post-processing code ppn. Results. The new rates differ from previous estimations due to the use of a series of updated resonance parameters and a direct S factor. Compared with the previous results from the Hauser-Feshbach statistical model, which assumes compound nucleus 43 V with a sufficiently high-level density in the energy region of astrophysical interest, large differences exist over the entire temperature region of rp-process interest, up to two orders of magnitude. We consistently calculated the photodisintegration rate using our new nuclear masses via the detailed balance principle, and found the discrepancies among the different reverse rates are much larger than those for the forward rate, up to ten orders of magnitude at the temperature of 10 8 K. Using a trajectory with a peak temperature of 1.95×10 9 K, we performed the rp-process nucleosynthesis simulations to investigate the impact of the new rates. Our calculations show that the adoption of the new forward and reverse rates result in abundance variations for Sc and Ca of 128% and 49%, respectively, compared to the variations for the statistical model rates. On the other hand, the overall abundance pattern is not significantly affected. The results of using new rates also confirm that the rp-process path does not bypass the isotope 43 V. Conclusions. Our study found that the Hauser-Feshbach statistical model is inappropriate to the reaction rate evaluation for 42 Ti(p,γ) 43 V. The adoption of the new rates confirms that the reaction path of 42 Ti(p,γ) 43 V(p,γ) 44 Cr(β + ) 44 V is a key branch of the rp process in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Multipoint Turbulence Analysis with HelioSwarm

Exploration of plasma dynamics in space, including turbulence, is entering a new era of multisatellite constellation measurements that will determine fundamental properties with unprecedented precision. Familiar but imprecise approximations will need to be abandoned and replaced with more-advanced approaches. We present a preparatory study of the evaluation of second- and third-order statistics, using simultaneous measurements at many points. Here, for specificity, the orbital configuration of the NASA Swarm mission is employed in conjunction with 3D magnetohydrodynamics numerical simulations of turbulence. The HelioSwarm nine-spacecraft constellation flies virtually through the turbulence to compare results with the exact numerical statistics. We demonstrate novel increment-based techniques for the computation of (1) the multidimensional spectra and (2) the turbulent energy flux. This latter increment-space estimate of the cascade rate, based on the third-order Yaglom–Politano–Pouquet theory, uses numerous increment-space tetrahedra. Our investigation reveals that HelioSwarm will provide crucial information on the nature of astrophysical turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving North American Wildfire Prediction by Integrating a Machine-Learning Fire Model in a Land Surface Model

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

54 ENVIRONMENTAL SCIENCES↗

Simulated wildfire burned area over the CONUS during 2001-2020

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM). A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

Liu, Ye↗

An intrinsic ductility parameter derived from anisotropic linear elasticity theory

A new indicator of intrinsic ductility (κ) is introduced based on linear elasticity. In the limit of elastic isotropy this parameter is equal to the Pugh ratio ( B / G ) plus a constant; but unlike the Pugh ratio, κ incorporates anisotropy and crystallography to improve predictive value. We identify a single ductile-to-brittle transition setpoint for κ that predicts crack-tip plasticity in atomistic simulations, experimental elongation to failure of polycrystalline elemental metals pulled in tension, and fracture energy measurements of glasses, suggesting commonality in the relationship between plasticity and fracture across all three of these cases. Statistical analysis supports the superiority of κ over B / G at predicting crack-tip plasticity in atomistic simulations and elongation to failure of polycrystalline elemental metals.

Winter, Ian S.↗

Accelerated statistical failure analysis of multifidelity TRISO fuel models

Statistical nuclear fuel failure analysis is critical for the design and development of advanced reactor technologies. Although Monte Carlo Sampling (MCS) is a standard method of statistical failure analysis for fuels, the low failure probabilities of some advanced fuel forms and the correspondingly large number of required model evaluations limit its application to low-fidelity (e.g., 1-D) fuel models. In this paper, we present four other statistical methods for fuel failure analysis in Bison, considering tri-structural isotropic (TRISO)-coated particle fuel as a case study. The statistical methods considered are Latin hypercube sampling (LHS), adaptive importance sampling (AIS), subset simulation (SS), and the Weibull theory. Using these methods, we analyzed both 1-D and 2-D representations of TRISO models to compute failure probabilities and the distributions of fuel properties that result in failures. The results of these methods compare well across all TRISO models considered. Overall, SS and the Weibull theory were deemed the most efficient, and can be applied to both 1-D and 2-D TRISO models to compute failure probabilities. Moreover, since SS also characterizes the distribution of parameters that cause TRISO failures, and can consider failure modes not described by the Weibull criterion, it may be preferred over the other methods. Finally, a discussion on the efficacy of different statistical methods of assessing nuclear fuel safety is provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory

Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Study of Pressure Fluctuations due to High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) are used to examine the pressure fluctuations generated by fully developed turbulence in supersonic turbulent boundary layers with an emphasis on both pressure fluctuations at the wall and the acoustic fluctuations radiated into the freestream. The wall and freestream pressure fields are first analyzed for a zero pressure gradient boundary layer with Mach 2.5 and Reynolds number based on momentum thickness of approximately 2835. The single and multi-point statistics reported include the wall pressure fluctuation intensities, frequency spectra, space-time correlations, and convection velocities. Single and multi-point statistics of surface pressure fluctuations show good agreement with measured data and previously published simulations of turbulent boundary layers under similar flow conditions. Spectral analysis shows that the acoustic fluctuations outside the boundary layer region have much lower energy content within the high-frequency region. The space-time correlations reflect the convective nature of the pressure field both at the wall and in the freestream, which is characterized by the downstream propagation of pressure-carrying eddies. Relative to those at the wall, the pressure-carrying eddies associated with the freestream signal are larger and convect at a significantly lower speed. The preliminary DNS results of a Mach 6 boundary layer show that the pressure rms in the freestream region is significantly higher than that of the lower Mach number case.

Duan, Lian↗

The Field Substellar Mass Function Based on the Full-sky 20 pc Census of 525 L, T, and Y Dwarfs

We present final Spitzer trigonometric parallaxes for 361 L, T, and Y dwarfs. We combine these with prior studies to build a list of 525 known L, T, and Y dwarfs within 20 pc of the Sun, 38 of which are presented here for the first time. Using published photometry and spectroscopy as well as our own follow-up, we present an array of color–magnitude and color–color diagrams to further characterize census members, and we provide polynomial fits to the bulk trends. Using these characterizations, we assign each object a T(eff) value and judge sample completeness over bins of T(eff) and spectral type. Except for types ≥T8 and T(eff) < 600 K, our census is statistically complete to the 20 pc limit. We compare our measured space densities to simulated density distributions and find that the best fit is a power law (dN/dM ∝ M^(-α) with α = 0.6 ± 0.1. We find that the evolutionary models of Saumon & Marley correctly predict the observed magnitude of the space density spike seen at 1200 K < T(eff) < 1350 K, believed to be caused by an increase in the cooling timescale across the L/T transition. Defining the low-mass terminus using this sample requires a more statistically robust and complete sample of dwarfs ≥Y0.5 and with T(eff) < 400 K. We conclude that such frigid objects must exist in substantial numbers, despite the fact that few have so far been identified, and we discuss possible reasons why they have largely eluded detection.

J. Davy Kirkpatrick↗

Impact of classical statistics on thermal conductivity predictions of BAs and diamond using machine learning molecular dynamics

Machine learning interatomic potentials (MLIPs) have greatly enhanced molecular dynamics (MD) simulations, achieving near-first-principles accuracy in thermal conductivity studies. In this work, we reveal that this accuracy, observed in BAs and diamond at sub-Debye temperatures, stems from an accidental error cancelation: classical statistics overestimates specific heat while underestimating phonon lifetimes, balancing out in thermal conductivity predictions. However, this balance is disrupted when isotopes are introduced, leading MLIP-based MD to significantly underpredict thermal conductivity compared to experiments and quantum statistics-based Boltzmann transport equation. This discrepancy arises not from classical statistics affecting phonon–isotope scattering rates but from its impact on the interplay between phonon–isotope and phonon–phonon scattering in the normal scattering-dominated BAs and diamond. In conclusion, this work underscores the limitations of MLIP-based MD for thermal conductivity studies at sub-Debye temperatures.

36 MATERIALS SCIENCE↗

ELM2.1-XGBfire1.0: improving wildfire prediction by integrating a machine learning fire model in a land surface model

Wildfires have shown increasing trends in both frequency and severity across the contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth system models (ESMs). Alternatively, fire models based on machine learning (ML), which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ELM2.1-XGBFire1.0) that integrates an eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran–C–Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001–2019, the ELM2.1-XGBFire1.0 outperforms process-based fire models in terms of spatial distribution and seasonal variations. The ELM2.1-XGBFire1.0 has proven to be a new tool for studying vegetation–fire interactions and, more importantly, enables seamless exploration of climate–fire feedback, working as an active component of E3SM.

54 ENVIRONMENTAL SCIENCES↗

Large-eddy simulations of the Northeastern US coastal marine boundary layer

In this study, large eddy simulations (LES) of offshore boundary layers near the Nantucket coast are performed using Nalu-Wind. The marine boundary layer conditions are chosen to match the predominantly unstable and neutral conditions measured by the Cape Wind platform. The appropriate domain, resolution, and boundary condition settings required for the LES are established through this work. Differences between stable and unstable cases are found in the wind speed profiles, averaged statistics, and wind spectra, and explained in terms of stratification effects.

54 ENVIRONMENTAL SCIENCES↗

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗

The Roles of Large- and Small-Scale Structures in Turbulent Free Shear Layers

The occurrence of large-scale coherent structures in turbulent free shear flows (especially the planar mixing layer) has been recognized for some time. Indeed, the observation of such structures in mixing layers did much to promote interest in the study of coherent structures in turbulence. It has been widely assumed that the large-scale structures in these flows are responsible for the entrainment of free-stream fluid and the overall growth of the layer, while the small-scale structures provide mixing and dissipation. A model of scalar mixing based on these ideas was proposed for these flows. However, recent experimental and computational evidence suggests that the dominance of the large-scale structures in turbulent mixing layers is not universal. In addition, there is a substantial variation among experiments in several statistical measures of self-similar mixing layers, for example growth rate and velocity variances. To investigate the importance of large-scale structures, several free shear flows (mixing layers and wakes) have been simulated via direct numerical simulation. The simulations are designed to mimic experimental mixing layers in which the splitter plate boundary layers are turbulent. Different levels of two-dimensional forcing are included resulting in large-scale structures of differing strength and importance. These simulations are used to investigate the role of large-scale coherent structures in free shear layers and the effect of these structures on relevant turbulence statistics and scalar mixing. It is found that the statistics and structures in several experiments involving turbulent mixing layers are in better agreement with simulations that do not exhibit dominant large-scale structures than those in which the common mixing layer structures do dominate. It is also found that the level of forcing can have a profound effect on the qualitative and quantitative features of these shear layer, even when they are nominally self-similar.

Moser, Robert D.↗

Monte Carlo turbulence simulation

The paper describes turbulence simulation experiments based on the principles of control system theory, that is, the construction of a system characterized by a system function such that upon exciting the system with prescribed noise processes the output of the system is a realization of a random processing the desired statistical attributes of turbulence. An experimental autocorrelation of Jimsphere measurements of wind velocity was approximated to simulate turbulent wind. From the approximate autocorrelation function, the required system function is obtained, and a discrete time system is designed. Another method of simulation is to solve the convolution integral by filter techniques. Other methods include discrete Fourier simulation and self-similar simulation.

Fichtl, G. H.↗

An Investigation of Time Distributions for Task Primitives to Support the HUNTER Dynamic Human Reliability Analysis

As an effort to support data collection for dynamic human reliability analysis (HRA), this study investigates time distributions for task primitives defined in the Goals, Operators, Methods, and Selection rules (GOMS)–Human Reliability Analysis (HRA) method. GOMS-HRA was developed to provide cognition-based time and human error probability information for dynamic HRA calculation in the Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) framework. HUNTER is a framework to support the dynamic modelling of human error in conjunction with other modelling tools. In this paper, we investigate time distributions using experimental data collected from the Simplified Human Error Experimental Program (SHEEP) study, which suggests an HRA data collection framework to complement full-scope simulator research as well as collect input data for dynamic HRA using simplified simulators such as the Rancor Microworld Simulator. In this study, time required for GOMS-HRA task primitives to satisfy thirteen statistical distributions is investigated. Then, the time distributions from student operators and professional operators are compared and discussed. As a result, this study identified several time distributions on five GOMS-HRA task primitives at a statistically significant level. According to analyses to date, a greater number of significant time distributions was found in abnormal or emergency operating procedures rather than standard operating procedures. In the future, it is expected that the result of this study can provide objective reference on elapsed time data for task primitives as well as help to realistically simulate scenarios within dynamic HRA.

99 GENERAL AND MISCELLANEOUS↗