Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bayes methods”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Periodicity significance testing with null-signal templates: reassessment of PTF’s SMBH binary candidates

Periodograms are widely employed for identifying periodicity in time series data, yet they often struggle to accurately quantify the statistical significance of detected periodic signals when the data complexity precludes reliable simulations. We develop a data-driven approach to address this challenge by introducing a null-signal template (NST). The NST is created by carefully randomizing the period of each cycle in the periodogram template, rendering it non-periodic. It has the same frequentist properties as a periodic signal template, and we show with simulations that the distribution of false positives is the same as with the original periodic template, regardless of the underlying data. Thus, performing a periodicity search with the NST acts as an effective simulation of the null (no-signal) hypothesis, without having to simulate the noise properties of the data. We apply the NST method to the supermassive black hole binaries (SMBHB) search in the Palomar Transient Factory (PTF), where Charisi et al. had previously proposed 33 high signal-to-noise candidates utilizing simulations to quantify their significance. Our approach reveals that these simulations do not capture the complexity of the real data. There are no statistically significant periodic signal detections above the non-periodic background. To improve the search sensitivity, we introduce a Gaussian quadrature based algorithm for the Bayes Factor with correlated noise as a test statistic. We show with simulations that this improves sensitivity to true signals by more than an order of magnitude. However, the Bayes Factor approach also results in no statistically significant detections in the PTF data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Comparison of modular analytical wake models to the Lillgrund wind plant

Efficient and accurate wake models are required for wind plant performance modeling and the suite of engineering processes that support wind plant layout, control, and monitoring. Although many analytical and engineering wake models with low computational costs have been proposed, their ability to represent the power production of large wind plants in a wide range of atmospheric conditions is not completely understood. The following validation study reviews the underlying theory for analytical wake models, outlines quality control procedures for observational data, and compares model results with observational data from the Lillgrund Wind Plant. Lillgrund makes a valuable case study for wake modeling because of its regular arrangement and the relatively close spacing of constituent wind turbines, which lead to regular and significant wake interactions within the wind plant and the development of deep array flow conditions. Formulations for the velocity deficit, wake-added turbulence, and wake superposition methods are considered in a modular sense, yielding many possible configurations to represent wind turbine wakes, of which seven are examined in detail. Velocity deficit models that account for flow conditions in the near wake are better able to reproduce power production for wind turbines in the transitional region of the wind plant, where wind turbines experience as many as five wakes from upstream turbines. Finally, in the deep array, where power production reaches asymptotic values and wake statistics become quasi-periodic, wake superposition schemes become the largest driver in error reduction; using the linear or maximum wake superposition methods can reduce the relative root mean square error by as much as 40% in the deep array.

17 WIND ENERGY↗
DOE OSTI2021

Characterizing the Non-linear Interactions Between Tide, Storm Surge, and River Flow in the Delaware Bay Estuary, United States

Low-lying coastal areas in the mid-Atlantic region are prone to compound flooding resulting from the co-occurrence of river floods and coastal storm surges. To better understand the contribution of non-linear tide-surge-river interactions to compound flooding, the unstructured-grid Finite Volume Community Ocean Model was applied to simulate coastal storm surge and flooding in the Delaware Bay Estuary in the United States. The model was validated with tide gauge data in the estuary for selected hurricane events. Non-linear interactions between tide-surge-river were investigated using a non-stationary tidal analysis method, which decomposes the interactions’ components at the frequency domain. Model results indicated that tide-river interactions damped semidiurnal tides, while the tide-surge interactions mainly influenced diurnal tides. Tide-river interactions suppressed the water level upstream while tide-surge interaction increased the water level downstream, which resulted in a transition zone of damping and enhancing effects where the tide-surge-river interaction was prominent. Evident compound flooding was observed as a result of non-linear tide-surge-river interactions. Furthermore, sensitivity analysis was carried out to evaluate the effect of river flooding on the non-linear interactions. The transition zone of damping and enhancing effects shifted downstream as the river flow rate increased.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Rapid wavefield forecasting for earthquake early warning via deep sequence to sequence learning

We propose a deep learning model, WaveCastNet, to forecast high-dimensional wavefields. WaveCastNet integrates a convolutional long expressive memory architecture into a sequence-to-sequence forecasting framework, enabling it to model long-term dependencies and multiscale patterns in both space and time. By sharing weights across spatial and temporal dimensions, WaveCastNet requires significantly fewer parameters than more resource-intensive models such as transformers, resulting in faster inference times. Crucially, WaveCastNet also generalizes better than transformers to rare and critical seismic scenarios, such as high-magnitude earthquakes. Here, we show the ability of the model to predict the intensity and timing of destructive ground motions in real time, using simulated data from the San Francisco Bay Area. Furthermore, we demonstrate its zero-shot capabilities by evaluating WaveCastNet on real earthquake data. Our approach does not require estimating earthquake magnitudes and epicenters, steps that are prone to error in conventional methods, nor does it rely on empirical ground-motion models, which often fail to capture strongly heterogeneous wave propagation effects.

Geophysics↗
DOE OSTI2021

Comparing Emulation Methods for Computer Models with High Dimensional Output

This Master’s Capstone project will present a comparison of statistical models for computer simulation studies. The four models included in this comparison study were chosen for both their proven and diverse methodologies. We will present a case study on hurricane flood data in the Delaware bay which highlights the strengths and weaknesses of each model when applied to a very large spatial field. As computers have gotten faster, we have become interested in modeling increasingly large spatial fields in both size and resolution. Statistical algorithms that are able to efficiently handle these fields have never been more important. We therefore find this comparison to be extremely topical.

42 ENGINEERING↗
DOE OSTI2023

Storm Surge Modeling as an Application of Local Time-Stepping in MPAS-Ocean

This paper presents the first practical application of local time-stepping (LTS) schemes in the Model for Prediction Across Scales-Ocean (MPAS-O). We use LTS schemes in a single-layer, global ocean model that predicts the storm surge around the eastern coast of the United States during Hurricane Sandy. The variable-resolution meshes used are of unprecedentedly high resolution in MPAS-O, containing cells as small as 125 m wide in Delaware Bay. It is shown that a particular, third-order LTS scheme (LTS3) produces sea-surface height solutions that are of comparable quality to solutions produced by the classical four-stage, fourth-order Runge-Kutta method (RK4) with a uniform time step on the same meshes. Furthermore, LTS3 is up to 35% faster in the best cases considered, where the number of cells using the coarse time-step relative to those using the fine time-step is as low as 1:1. This shows that LTS schemes are viable for use in MPAS-O with the added benefit of substantially less computational cost. The results of these performance experiments inform us of the requirements for efficient mesh design and configuration of LTS regions for LTS schemes. In particular, we see that for LTS to be efficient on a given mesh, it is important to have enough cells using the coarse time-step relative to those using the fine time-step, typically at least 1:5 to see an increase in performance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermal-striping analysis methodology for sodium-cooled reactor design

Computational fluid dynamics simulations have been performed to study the applicability of engineering CFD methods for thermals striping analysis. Thermal striping is the fluctuating temperature profile in a solid caused by fluctuating fluid temperature, and the resulting fluctuating thermal stresses cause high cycle fatigue and eventual material failure. Here this study presents the methodology for thermal striping analysis that including a transient conjugate heat transfer RANS model and benchmark the results against a sodium triple jet experiment. The multiphysics calculations are implemented to analyze the thermal stresses in the solid domain caused by the coupled heat transfer between the fluid and solid domains. Validation data from a liquid sodium triple jet experiment include time-averaged temperature measurements and power spectra of the temperature signal. The numerical results agree well with these experimental measures, demonstrating key features such as the dominant frequency of temperature fluctuations. The applicability of a low-cost wall treatment method is demonstrated, enabling key computational savings. Finally, the performance of two finite element stress analysis software packages is compared, and the validity of the lower-cost method is confirmed. These results demonstrate the applicability of engineering methods for computational thermal striping calculations, enabling thermal striping estimations in large fluid systems such as the core of a Sodium-cooled nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Probabilistic Nanomagnetic Memories for Uncertain and Robust Machine Learning

This project evaluated the use of emerging spintronic memory devices for robust and efficient variational inference schemes. Variational inference (VI) schemes, which constrain the distribution for each weight to be a Gaussian distribution with a mean and standard deviation, are a tractable method for calculating posterior distributions of weights in a Bayesian neural network such that this neural network can also be trained using the powerful backpropagation algorithm. Our project focuses on domain-wall magnetic tunnel junctions (DW-MTJs), a powerful multi-functional spintronic synapse design that can achieve low power switching while also opening the pathway towards repeatable, analog operation using fabricated notches. Our initial efforts to employ DW-MTJs as an all-in-one stochastic synapse with both a mean and standard deviation didn’t end up meeting the quality metrics for hardware-friendly VI. In the future, new device stacks and methods for expressive anisotropy modification may make this idea still possible. However, as a fall back that immediately satisfies our requirements, we invented and detailed how the combination of a DW-MTJ synapse encoding the mean and a probabilistic Bayes-MTJ device, programmed via a ferroelectric or ionically modifiable layer, can robustly and expressively implement VI. This design includes a physics-informed small circuit model, that was scaled up to perform and demonstrate rigorous uncertainty quantification applications, up to and including small convolutional networks on a grayscale image classification task, and larger (Residual) networks implementing multi-channel image classification. Lastly, as these results and ideas all depend upon the idea of an inference application where weights (spintronic memory states) remain non-volatile, the retention of these synapses for the notched case was further interrogated. These investigations revealed and emphasized the importance of both notch geometry and anisotropy modification in order to further enhance the endurance of written spintronic states. In the near future, these results will be mapped to effective predictions for room temperature and elevated operation DW-MTJ memory retention, and experimentally verified when devices become available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

The look-elsewhere effect from a unified Bayesian and frequentist perspective

When searching over a large parameter space for anomalies such as events, peaks, objects, or particles, there is a large probability that spurious signals with seemingly high signi ficance will be found. This is known as the look-elsewhere effect and is prevalent throughout cosmology, (astro)particle physics, and beyond. To avoid making false claims of detection, one must account for this effect when assigning the statistical significance of an anomaly. This is typically accomplished by considering the trials factor, which is generally computed numerically via potentially expensive simulations. In this paper we develop a continuous generalization of the Bonferroni and Sidak corrections by applying the Laplace approximation to evaluate the Bayes factor, and in turn relating the trials factor to the prior-to-posterior volume ratio. Here, we use this to define a test statistic whose frequentist properties have a simple interpretation in terms of the global p-value, or statistical significance. We apply this method to various physics-based examples and show it to work well for the full range of p-values, i.e. in both the asymptotic and non-asymptotic regimes. We also show that this method naturally accounts for other model complexities such as additional degrees of freedom, generalizing Wilks' theorem. This provides a fast way to quantify statistical significance in light of the look-elsewhere effect, without resorting to expensive simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Near Term Race for the Neutrino Mass Ordering

The muon neutrino disappearance experiments T2K and NOvA, measure a different combination of Delta m^2_31 and Delta m^2_32 then the electron neutrino disappearance experiments Daya Bay and JUNO. This can be used to determine the neutrino mass ordering, one of the outstanding questions of the neutrino sector. In this talk I will explain why this happens and what is the current state and future prospects for the mass ordering determination using this method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Utilizing Advanced Statistics to Determine Anomalistic Conditions in Pebble-Bed Reactors

Pebble-bed reactors (PBRs) utilize hundreds of thousands of fuel pebbles, which continuously circulate through the core, in lieu of traditional fuel assemblies to generate fissions and produce power. The use of unmarked fuel pebbles presents a challenge for international safeguards verification that nuclear material is not being diverted. To ensure pebble diversion is not taking place, new methods for accounting for and monitoring the pebbles should be examined to determine an appropriate methodology for performing bulk accountancy with pebbles. Here, this work examines the use of statistical methods for determining if the reactor is within a declared range of operation by examining the statistical distribution of pebble burnup as they are discharged from the core. Using this methodology, we created a model that detects diversion over 95% of the time, over multiple diversion pathways, if the reactor core maintains a constant power density during the diversion process and only falsely labels a diversion case nominal 2% of the time. For a diversion scenario where the reactor is maintained at a constant power, the statistical analysis can correctly identify if diversion is occurring over 80% of the time; however, nearly 20% of specific diversion pathways are mislabeled nominal. These results provide a basis and framework for exploring the further use of statistical methods to determine where these methods could be most useful and where additional methods, such as machine learning, could be used to capture if diversion is occurring in pebble-bed reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Non-standard interactions in SMEFT confronted with terrestrial neutrino experiments

The Standard Model Effective Field Theory (SMEFT) provides a systematic and model-independent framework to study neutrino non-standard interactions (NSIs). We study the constraining power of the on-going neutrino oscillation experiments T2K, NO ν A, Daya Bay, Double Chooz and RENO in the SMEFT framework. A full consideration of matching is provided between different effective field theories and the renormalization group running at different scales, filling the gap between the low-energy neutrino oscillation experiments and SMEFT at the UV scale. We first illustrate our method with a top- down approach in a simplified scalar leptoquark model, showing more stringent constraints from the neutrino oscillation experiments compared to collider studies. We then provide a bottom-up study on individual dimension-6 SMEFT operators and find NSIs in neutrino experiments already sensitive to new physics at ~20 TeV when the Wilson coefficients are fixed at unity. We also investigate the correlation among multiple operators at the UV scale and find it could change the constraints on SMEFT operators by several orders of magnitude compared with when only one operator is considered. Furthermore, we find that accelerator and reactor neutrino experiments are sensitive to different SMEFT operators, which highlights the complementarity of the two experiment types.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Analysis of Second Target Station Target Segment Removal Dose Rates

This report documents the analysis of the dose rate fields in two configurations of the target system of the Second Target Station (STS) at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) in support of the remote handling target segment removal operation. The two configurations detailed in this report are both assuming 8 hours since beam was incident on the target. The first configuration assumes all shielding is in place just as it was during operation and allows for confirmation of hands-on maintenance capability of the target drive components. The second configuration assumes that a hatch in the target drive room ceiling has been opened to the high bay, the target segment removable core vessel shield block has been removed, and there is a direct line-of-sight to the target segment some 5 m below the target drive room elevation. Dose rate fields in both configurations are documented in this report along with details of the methods of computation. The dose rate fields included in this report are intended to support the engineering of remote handling special equipment and the planning of remote handing operations. The models and analysis detailed in this report are in support of the preliminary design of the STS and future analysis is needed to support the final design of the STS.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Identification of contamination in the pulse-compression chamber of the OMEGA EP Laser System originating from clean room garments

The contamination of optical components with microscopic particles strongly impacts their ability to handle high-power or high-intensity laser pulses via an array of mechanisms. This work explores the nature and origin of a subset of contamination particles found inside the pulse-compression chamber of the OMEGA EP Laser System. Using Raman microscopy accompanied by other analytical methods, it was shown that sodium nitrate constitutes a significant fraction of the particles that possess diameters on the order of a few micrometers. Further investigation aiming to reveal the origin of these particles suggests that these sodium nitrate particles are deposited via shedding from the reusable clean room garments worn in the Laser Bay and compression chamber. Raman microscopy also revealed various plastic particles such as polystyrene, polyethylene, and poly(diallyl isophthalate).

Raman spectroscopy↗
DOE OSTI2024

Implementing Superresolution of Nonstationary Tides with Wavelets: An Introduction to CWT_Multi

Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.

Lobo, Matthew↗
DOE OSTI2022

Optimized Machine Learning Model for Predicting Groundwater Contamination

The use of physical models to predict groundwater contaminant movement remains technically challenging due to the complexity of the phenomena, the heterogeneity of key parameters in nature, and the presence of poorly defined interactive and feedback processes. New approaches to address these challenges are needed. In this study, we evaluate various Artificial Intelligence (AI)-based approaches to understand a hexavalent chromium (Cr(VI)) plumes located on the U.S. Department of Energy’s (DOE) Hanford Site in Richland, WA. The groundwater monitoring dataset used in this study included data from the 100 Area along the Columbia River and included data collected between 2010 to 2019. This study investigates the most prominent contaminant, Cr(VI), with the Extreme Gradient Boosting (XGBoost) machine learning model. The XGBoost models were compared with optimized versions using an Empirical Bayes Search Cross-Validation technique for better prediction. The optimized XGBoost model yielded an R^2 value of 0.99 on the training set and 0.85 on the testing set, whereas XGBoost without optimization yielded a value of 0.83 on the training set and 0.85 on the testing set. This paper provides an overview of a computational method for groundwater contamination modeling that shows promise for improving current remediation efforts.

Mazumdar, Hirak↗
DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A comparison of eight optimization methods applied to a wind farm layout optimization problem

Abstract. Selecting a wind farm layout optimization method is difficult. Comparisons between optimization methods in different papers can be uncertain due to the difficulty of exactly reproducing the objective function. Comparisons by just a few authors in one paper can be uncertain if the authors do not have experience using each algorithm. In this work we provide an algorithm comparison for a wind farm layout optimization case study between eight optimization methods applied, or directed, by researchers who developed those algorithms or who had other experience using them. We provided the objective function to each researcher to avoid ambiguity about relative performance due to a difference in objective function. While these comparisons are not perfect, we try to treat each algorithm more fairly by having researchers with experience using each algorithm apply each algorithm and by having a common objective function provided for analysis. The case study is from the International Energy Association (IEA) Wind Task 37, based on the Borssele III and IV wind farms with 81 turbines. Of particular interest in this case study is the presence of disconnected boundary regions and concave boundary features. The optimization methods studied represent a wide range of approaches, including gradient-free, gradient-based, and hybrid methods; discrete and continuous problem formulations; single-run and multi-start approaches; and mathematical and heuristic algorithms. We provide descriptions and references (where applicable) for each optimization method, as well as lists of pros and cons, to help readers determine an appropriate method for their use case. All the optimization methods perform similarly, with optimized wake loss values between 15.48 % and 15.70 % as compared to 17.28 % for the unoptimized provided layout. Each of the layouts found were different, but all layouts exhibited similar characteristics. Strong similarities across all the layouts include tightly packing wind turbines along the outer borders, loosely spacing turbines in the internal regions, and allocating similar numbers of turbines to each discrete boundary region. The best layout by annual energy production (AEP) was found using a new sequential allocation method, discrete exploration-based optimization (DEBO). Based on the results in this study, it appears that using an optimization algorithm can significantly improve wind farm performance, but there are many optimization methods that can perform well on the wind farm layout optimization problem, given that they are applied correctly.

17 WIND ENERGY↗