Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed Estimation”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Precise ROCOF estimation algorithm for low inertia power grids

The precise estimation of Rate Of Change Of Frequencies (ROCOFs) in a generation trip event can be helpful on power system inertia estimation, fast system response, and accurate event capturing. However, the ROCOF estimations from the existing Synchronized Measurement Devices (SMDs) are usually simply calculated by the finite difference between two adjacent frequency measurement points. The ineluctable noises, disturbances, and spikes from real-world frequency measurements can bring large dynamics to the ROCOF estimation and thus can result in an inaccurate estimation of the initial ROCOF. This issue becomes more serious when the target power grid has low inertia, and a large amount of distributed energy sources are deployed. To address this issue, a precise ROCOF estimation algorithm is designed based on least square method with flexible window size. In addition, a median filter is also designed and applied on the frequency measurements before using the proposed algorithm. The window size and thresholds in the proposed algorithm are determined with historical event data analysis. The proposed algorithm is deployed in a low cost, flexible, and distribution level universal grid analyzer (UGA) platform. In conclusion, multiple experiments are conducted in both a laboratory and the Hawaiian Islands to verify the effectiveness of the proposed algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Nuclear Data Uncertainty Quantification for Adjoint-Weighted Quantities via XGPT in Monte Carlo Codes

The Total Monte Carlo method (TMC) is an innovative approach to nuclear data uncertainty propagation. In TMC, a large number of independent ENDF files are randomly generated starting from nuclear models parameters and their uncertainties. These files are then processed with NJOY or similar codes to produce a set of formatted continuous-energy cross sections (i.e., ACE files), and independent Monte Carlo neutron transport simulations are run with these different ACE files as input. Finally the distributions of the considered response functions are obtained directly from the results of the Monte Carlo runs. Despite several improvements of the TMC method, this approach is often considered an inefficient way of propagating nuclear data uncertainties, as it requires several computationally expensive Monte Carlo runs. Covariance-based methods present important advantages for several applications, especially when different sources of uncertainties are to be investigated. For this reason, sensitivity analysis and uncertainty quantification capabilities are being implemented in several Monte Carlo neutron transport codes. On the other hand, the fact that no artificial multi-group covariance matrices are required and that higher moments of the responses distributions can be easily obtained, are definite advantages of the TMC method. In an effort to overcome the limitations of both approaches, a new method, called eXtended Generalized Perturbation Theory (XGPT) is under development. When performing nuclear data uncertainty propagation, such XGPT approach does not require generating discretized (i.e., multigroup) covariance matrices, and allows to estimate higher moments of the response probability distributions, via linear perturbation theory. One of the promising applications of the XGPT method involves its adoption in conjunction with the Proper Orthogonal Decomposition (POD) of nuclear data, to build a Reduced Order Model (ROM) for the approximation of thousands of independent TMC results in a single criticality source simulation. In the present work, this approach is employed for the first time to evaluate the impact of uncertainty in nuclear data on adjoint-weighted quantities. The method has been implemented in an extended version of the Serpent-2.1.24 Monte Carlo code and is demonstrated on two fast spectrum benchmarks, focusing on {sup 208}Pb nuclear data. A new method for propagating nuclear data uncertainties in criticality source Monte Carlo simulations has been developed. This eXtended Generalized Perturbation Theory (XGPT) approach allows to accurately reproduce reference Total Monte Carlo (TMC) uncertainty estimates in a single criticality source simulation. The present work aims to demonstrate that the proposed method can quantify nuclear data uncertainty based on adjoint-weighted quantities. The XGPT approach was applied to evaluate the impact of {sup 208}Pb nuclear data uncertainties on the effective prompt neutron lifetime (l{sub eff}) in two fast spectrum systems. The obtained XGPT results show a good agreement with the reference TMC results, obtained from independent Monte Carlo runs with {sup 208}Pb nuclear data from random TENDL-2013 evaluations. The XGPT method, implemented in a extended version of the Serpent Monte Carlo code and adopted along with the Proper Orthogonal Decomposition (POD) of the nuclear data, is able to provide an accurate reduced order model to approximate thousands TMC results within a single run, ensuring a large saving in computational requirements. The effect of {sup 208}Pb nuclear data uncertainties appears to be rather large. This is consistent with previous studies that adopted the same input random evaluations data from TENDL. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Multi-Fidelity Gaussian Process Regression Method for Probabilistic Wind Farm Power Curve Estimation

Accurate estimation of the power curve for wind turbines or wind farms is crucial to ensure their efficient operation and management. However, conventional methods for power curve estimation rely either on expensive and infrequent measurements or on low-quality numerical simulations. Moreover, the majority of previous studies on power curve estimation for wind turbines or wind farms focused on deterministic estimation, which provides a point estimate of the relationship between wind speed and power generation. Nevertheless, the deterministic approach fails to consider the inherent uncertainty associated with wind energy production resulting from varying turbine characteristics. This can lead to inaccurate power generation estimation and suboptimal decisions regarding energy management. In this paper, a kernel density estimation (KDE) based Multi-Fidelity Gaussian Process Regression (MFGPR) model is proposed to fuse theoretical power curve data and the ground true measurements to create a mapping of wind speed and wind power. By conducting a case study on an actual wind farm in China, the efficacy of the proposed MFGPR model was demonstrated in characterizing the variability of wind power. The probabilistic MFGPR model was also able to generate confidence intervals that encompassed the measured power, thereby improving the accuracy and confidence in wind power estimation or wind resource assessment. Overall, the proposed MFGPR model offers a reliable approach to integrate high-fidelity ground measurements and theoretical power curve data, resulting in precise wind resource assessment and power estimation.

Gaussian process regression↗
DOE OSTI2024

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗
DOE OSTI2022

Twist-three cross-sections in deeply virtual Compton scattering

We study the deeply virtual Compton scattering process with both twist-two and twist-three Compton form factors and present our cross-sections formulas with all polarization configurations. While the twist-three contributions are generally assumed to be negligible in the literature due to the kinematical suppression, we compare them with the twist-two ones at typical JLab 6 GeV and 12 GeV kinematics as well as EIC kinematics and show their kinematical suppression explicitly, justifying the leading-twist approximation made in the literature. In addition, we also estimate the twist-three Compton form factors using Wandzura-Wilczek relations and inputs of twist-two generalized parton distributions based on a reggeized spectator model. With those estimated Compton form factors, we analyze the kinematical behavior of twist-two and twist-three cross-sections in a wide range of kinematics, and discuss the optimal regions for separating the leading-twist effects from the higher-twist ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Application of the PHDS Germanium Gamma-ray Imager (GeGI) in Facility Decontamination Assessment - 19471

As part of a site wide investigation, U.S. EPA Region 5 had to assess the nature and extent of contamination in a former radioactive source manufacturing facility. This work involved verification of the radioisotope expected, ensuring no other radioisotopes were present, and quantification of the activity and areal extent of contamination throughout the facility. The device used to support this assessment work was the Germanium Gamma-ray Imager (GeGI) manufactured by PHDS Co. of Knoxville, Tennessee. GeGI is a germanium-semiconductor imaging detector system that provides identification, location, and quantification of radioactive materials within its field of view. GeGI provides the highest fidelity radiological assay information including the best gamma-ray energy resolution and spatial localization of any existing technology, including gamma-ray energy resolution of ∼0.3% full-width half-max (FWHM) at 662 keV and angular location resolution ∼ +/-1 degree. The GeGI is a mechanically cooled, portable system well-suited to field operations. PHDS Co. also provided on-site support during assessment operations and performed subsequent data analysis. Gamma-ray spectral data and associated images were collected with GeGI in six primary areas of the facility having significant residual contamination. A total of 61 gamma-ray spectrum/image sets were collected over the course of three working days (22 Compton images and 39 Pinhole/Compton images). During data collection, Compton imaging was used as the means by which to identify and prioritize hotspots. After contaminated regions were identified via Compton imaging, short-distance Pinhole spectrum/image sets were collected for additional analysis. All of the spectrum/image sets were analyzed to determine the radioisotope distribution within the GeGI field of view, as well as to provide an estimate of the activity associated with each hotspot. It was determined that the nature of residual contamination present was distributed across specific areas of elevated activity, where in some areas there was significant accumulation of contaminated material. It was also determined that no discrete sources of high-activity were present, as would be expected in the case of radioactive material product intended for sealed source manufacture. Results included the activity of residual contamination for each of the areas assessed, and an estimate of the physical area over which that contamination was distributed. It was anticipated by EPA that if decontamination were to be performed that it would involve some type of concrete scabbling or other means of physical removal of contaminated material. The GeGI data could be used to estimate the areal extent (square footage) of concrete removal that would be required for decontamination. Using the activity estimate data from the GeGI along with an estimated depth of concrete removal, planning could be performed for decontamination work (equipment and personnel), waste collection, packaging, transportation, and disposal, along with an associated cost estimate. With regard to health and safety, use of the GeGI reduced the radiation dose to the EPA R5 staff performing the assessment work. The GeGI is a non-destructive device that supports assessment of radioactively contaminated areas and sources from a distance, reducing the need for close contact to radiation, and the time necessary for personnel to operate near such hazards. Use of the GeGI in this facility assessment was our first attempt to apply this technology in what proved to be a quick and efficient decontamination assessment method. The GeGI can assess the nature and extent of residual contamination in a facility, first by quickly identifying the areas of elevated activity, followed by refined measurement and assessment of those elevated areas. The GeGI data can then be used to develop planning and a cost estimate of required decontamination and waste management activities. Since it's a non-destructive gamma identification and imaging device, it can reduce the distance of workers from radioactive sources and the dose incurred from assessment work. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Data from: Comparing traditional and Bayesian approaches to ecological meta-analysis

Despite the wide application of meta-analysis in ecology, some of the traditional methods used for meta-analysis may not perform well given the type of data characteristic of ecological meta-analyses. We reviewed published meta-analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies (ni) and the number of studies or records (k) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta-analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that ni and k were highly variable, distributions were right-skewed, and were generally small (median ni =5, median k=44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval based on the t-distribution that uses an adjusted standard error (the Hartung-Knapp-Sidik-Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z-distribution confidence intervals had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta-analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random-effects meta-analysis, the performance of the best frequentist and Bayesian methods were similar for the same combinations of factors (k and mean replication), though the Bayesian approaches had higher than nominal (>95%) coverage for the mean effect when k was very low (k<15). Our literature review suggests that many meta-analyses that used z-distribution or bootstrapping confidence intervals may have over-estimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta-analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Machine Learning-Based Method to Estimate Transformer Primary-Side Voltages with Limited Customer-Side AMI Measurements

Distribution control applications such as volt/var optimization, network reconfiguration, and distribution automation require accurate knowledge of the distribution system state. The lack of sufficient sensors on the primary side of distribution networks often limits the accuracy of the control decisions by these applications. The deployment of advanced metering infrastructure (AMI) provides utilities an opportunity to translate the AMI data on the secondary onto the primary so that it can be used as pseudo-measurements to augment the limited existing measurements on the primary. This paper develops an approach for estimating service transformer primary-side voltages by using limited secondary-side AMI measurements. The estimated primary-side voltages can be used by utilities as pseudo-measurements for distribution control applications. The detailed secondary model topology, which is an essential input data for many existing algorithms, is not required for the proposed method. The performance of the proposed method is validated by using AMI measurements from the field and an actual distribution feeder model of San Diego Gas & Electric Company.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Quantifying concordance of correlated cosmological data sets

We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the 1% probability threshold. Furthermore, this shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A Bayesian Approach to Real-Time Dynamic Parameter Estimation Using Phasor Measurement Unit Measurement

Here, we develop a polynomial-chaos-expansion (PCE)-based approach for decentralized dynamic parameter estimation through Bayesian inference. Using this approach, the non-Gaussian distribution of the inverted parameters is obtained. More specifically, we first represent the decentralized generator model with the PCE-based surrogate. This surrogate allows us to efficiently evaluate the time-consuming dynamic solver at parameter values through Metropolis-Hastings (M-H)-based Markov chain Monte Carlo (MCMC). Then, we propose a two-stage hybrid Markov chain Monte Carlo (MCMC) to recover a posteriori distribution of the decentralized generator model parameters. In the first stage, we use the gradient-enhanced Langevin MCMC algorithm to characterize an intermediate posterior parameter distribution. This algorithm is computationally scalable to the high-dimensional parameter space. Based on the intermediate posterior distribution, during the second stage, we use the adaptive MCMC algorithm to fine-tune the strong correlations between the parameters. Finally, the fully recovered a posterior distribution is obtained in the end. The simulation results show that the proposed PCE-based hybrid MCMC algorithm can accurately and efficiently estimate the high-dimensional generator dynamic model parameters with full probabilistic distribution provided.

42 ENGINEERING↗
DOE OSTI2016

ASSESSING ASTROPHYSICAL UNCERTAINTIES IN DIRECT DETECTION WITH GALAXY SIMULATIONS

We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions that are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Furthermore, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulations. We show that the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way’s estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models that are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deep Learning At Depth: Estimating subsurface parameters from geophysical monitoring data

Geophysical imaging techniques are a non-invasive way to image the subsurface and understand both subsurface solid (rock/soil) and fluid property distributions and their evolution in time. Inversions of the geophysical data, such as Electrical Resistance Tomography (ERT) data, are solved to estimate the subsurface property distributions, such as conductivity, and many inversion techniques smooth out sharp gradients in rock or fluid property distributions. Sharp gradients in subsurface properties tend to be present in situations with complex subsurface structures, which are common in many subsurface applications. We have successfully demonstrated that it is possible to inform, or constrain, inversions with neural networks trained on synthetic data with complex subsurface structures. Initial results suggest this process may be optimizable to yield property distributions that better represent the true property distributions than the same inversion process without the neural network constraint. Future work would optimize the neural network performance for this application and then apply the synthetic-data trained neural network to real data to understand the utility and performance of this technique for real data sets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Demonstration of machine-learning-enhanced Bayesian quantum state estimation

Machine learning (ML) has found broad applicability in quantum information science in topics as diverse as experimental design, state classification, and even studies on quantum foundations. Here, we experimentally realize an approach for defining custom prior distributions that are automatically tuned using ML for Bayesian quantum state estimation methods that generally better conform to the physical properties of the underlying system than standard fixed prior distributions. Previously, researchers have looked to Bayesian quantum state tomography for advantages like uncertainty quantification, the return of reliable estimates under any measurement condition, and minimal mean-squared error. However, practical challenges related to long computation times and conceptual issues concerning how to incorporate prior knowledge most suitably can overshadow these benefits. Using both simulated and experimental measurement results, we demonstrate that ML-defined prior distributions reduce net convergence times and provide a natural way to incorporate both implicit and explicit information directly into the prior distribution. These results constitute a promising path toward practical implementations of Bayesian quantum state tomography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Shielding of the azimuthal magnetic field by the anode plasma in a relativistic self-magnetic-pinch diode

In relativistic electron beam diodes, the self-generated magnetic field causes electron-beam focusing at the center of the anode. Generally, plasma is formed all over the anode surface during and after the process of the beam focusing. In this work, we use visible-light Zeeman-effect spectroscopy for the determination of the magnetic field in the anode plasma in the Sandia 10 MV, 200 kA (RITS-6) electron beam diode. The magnetic field is determined from the Zeeman-dominated shapes of the Al III 4 s –4 p and C IV 3 s –3 p doublet emissions from various radial positions. Near the anode surface, due to the high plasma density, the spectral line-shapes are Stark-dominated, and only an upper limit of the magnetic field can be determined. The line-shape analysis also yields the plasma density. The data yield quantitatively the magnetic-field shielding in the plasma. In conclusion, the magnetic-field distribution in the plasma is compared to the field-diffusion prediction and found to be consistent with the Spitzer resistivity, estimated using the electron temperature and charge-state distribution determined from line intensity ratios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Logarithmic Resilience Risk Metrics That Address the Huge Variations in Blackout Cost

Resilience risk metrics must address the customer cost of the largest blackouts of greatest impact. However, there are huge variations in blackout cost in observed distribution utility data that make it impractical to properly estimate the mean large blackout cost and the corresponding risk. These problems are caused by the heavy tail observed in the distribution of customer costs. To solve these problems, we propose resilience metrics that describe large blackout risk using the mean of the logarithm of the cost of large-cost blackouts, the slope index of the heavy tail, and the frequency of large-cost blackouts.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Research Trends and Applications of PMUs

This work is a survey of current trends in applications of PMUs. PMUs have the potential to solve major problems in the areas of power system estimation, protection, and stability. A variety of methods are being used for these purposes, including statistical techniques, mathematical transformations, probability, and AI. The results produced by the techniques reviewed in this work are promising, but there is work to be performed in the context of implementation and standardization. As the smart grid initiative continues to advance, the number of intelligent devices monitoring the power grid continues to increase. PMUs are at the center of this initiative, and as a result, each year more PMUs are deployed across the grid. Since their introduction, myriad solutions based on PMU-technology have been suggested. The high sampling rates and synchronized measurements provided by PMUs are expected to drive significant advancements across multiple fields, such as the protection, estimation, and control of the power grid. This work offers a review of contemporary research trends and applications of PMU technology. Most solutions presented in this work were published in the last five years, and techniques showing potential for significant impact are highlighted in greater detail. Being a relatively new technology, there are several issues that must be addressed before PMU-based solutions can be successfully implemented. This survey found that key areas where improvements are needed include the establishment of PMU-observability, data processing algorithms, the handling of heterogeneous sampling rates, and the minimization of the investment in infrastructure for PMU communication. Solutions based on Bayesian estimation, as well as those having a distributed architectures, show great promise. The material presented in this document is tailored to both new researchers entering this field and experienced researchers wishing to become acquainted with emerging trends.

42 ENGINEERING↗
DOE OSTI2016

Turbulence Model Study for Pressure and Velocity Distributions for Molten Salts in a Crossflow Tube Bundle

Empirical pressure drop correlations for various compositions of molten salts in crossflow across tube bundles/banks are limited due to the scarcity of experimental data for these fluids. This lack of experimental data presents a limitation when designing heat exchangers for Molten Salt Reactors (MSRs) that utilize cross-flow configurations with turbulent molten salts flowing in the shell side. In an effort to provide a starting point for future design and experimental work, a computational fluid dynamics (CFD) campaign is being conducted to estimate pressure drops over a range of Reynolds numbers, pitch to diameter ratios, and for both staggered and in-line tube banks. In an effort to continue this work, the study of various turbulence models were conducted to determine the applicability for steady RANS simulations. An in-line tube bank with FLiBe (LiF-BeF{sub 2}), a frequently selected molten salt in MSRs, was selected for these simulations. The estimated pressure drops and pressure and velocity distributions for each turbulence model was compared to determine the applicability. The k-ε turbulence model family was found to be the most appropriate model for steady Reynolds averaged Navier Stokes within this type of domain. Difference turbulence models were investigated for simulation of molten salt experiencing crossflow within an in-line tube bundle. The flow was simulated using a steady RANS formulation with different families of turbulence models. The different turbulence models provide a view of this specific flow behavior potential requirements and applicable for steady RANS. Of particular note, the k-ε models were able to converge easily within this domain. The predicted behavior for both models were observed to be mostly similar for both velocity and pressure distributions. Whereas, k-ω SST and SA models did not converge well and had observable differences for the predictions of pressure and velocity distributions. The global behaviors were generally similar, but the variations occurring suggest k-ω SST and SA turbulence models are not useful in the steady RANS formulation of this domain. When comparing to the empirical pressure drop, each turbulence model resulted in overpredictions ranging from 8.7% and 19.75% difference. The standard low-Re k-ε model having the lowest difference is likely the better model for simplified predictions for pressure drops. Further investigations will be conducted to determine the convergence issues for the k-ω SST and SA turbulence models. These might be due to the flow geometry or potentially a revision of the meshing is needed. The flow behavior may not be resolved using these models in a steady formulation due to inherent unsteadiness of the flow. (authors)

42 ENGINEERING↗
DOE OSTI2016

Regional oxygen reduction and denitrification rates in groundwater from multi-model residence time distributions, San Joaquin Valley, USA

Rates of oxygen and nitrate reduction are key factors in determining the chemical evolution of groundwater. Little is known about how these rates vary and covary in regional groundwater settings, as few studies have focused on regional datasets with multiple tracers and methods of analysis that account for effects of mixed residence times on apparent reaction rates. This study provides insight into the characteristics of residence times and rates of O 2 reduction and denitrification (NO 3 – reduction) by comparing reaction rates using multi-model analytical residence time distributions (RTDs) applied to a data set of atmospheric tracers of groundwater age and geochemical data from 141 well samples in the Central Eastern San Joaquin Valley, CA. The RTD approach accounts for mixtures of residence times in a single sample to provide estimates of in-situ rates. Tracers included SF 6 , CFCs, 3 H, He from 3 H (tritiogenic He), 14 C, and terrigenic He. Parameter estimation and multi-model averaging were used to establish RTDs with lower error variances than those produced by individual RTD models. The set of multi-model RTDs was used in combination with NO 3 – and dissolved gas data to estimate zero order and first order rates of O 2 reduction and denitrification. Results indicated that O 2 reduction and denitrification rates followed approximately log-normal distributions. Rates of O 2 and NO 3 – reduction were correlated and, on an electron milliequivalent basis, denitrification rates tended to exceed O 2 reduction rates. Estimated historical NO 3 – trends were similar to historical measurements. Here, results show that the multi-model approach can improve estimation of age distributions, and that relatively easily measured O 2 rates can provide information about trends in denitrification rates, which are more difficult to estimate.

54 ENVIRONMENTAL SCIENCES↗