Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normalized”

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 361 records · Page 20

The influence of temperature history and flow mixing on the vapor-phase speciation of uranium oxide nanoparticles

The transport of radioactive material following a nuclear event depends heavily on fallout characteristics (e.g., composition, shape, and size), which in turn are controlled by the conditions under which the fallout forms. To this end, we use a novel plasma flow reactor (PFR) to explore how conditions including temperature history and gas mixing influence the chemical speciation and physical characteristics of uranium oxide nanoparticles. The PFR consists of a glass tube attached to an inductively coupled plasma (ICP) torch via an adaptor piece (ring flow injector) located downstream of the plasma-generating RF coil. We increase the gas flow rates from normal ICP operating conditions to higher values both upstream at the RF coil and downstream at the ring flow injector to create four distinct temperature profiles and two distinct gas flow mixing conditions (normal and higher Ar atomic ratios). Ex situ transmission electron microscopy (TEM) analysis is used to determine the crystal structure, morphology, and size distribution of the synthesized particles as a function of temperature and mixing. The particles display a striking dependence on these conditions in terms of the resulting speciation of uranium oxide: lower temperature (higher upstream) and slower cooling (normal downstream with normal Ar atomic ratio mixing) favor the formation of α-UO 3 particles, while faster cooling (higher downstream with higher Ar atomic ratio mixing) and higher temperature (normal upstream) favor the formation of UO 2 particles. The extent of gas mixing dictates particle aggregation: better mixing encourages the particles to aggregate more extensively. Here these results demonstrate the sensitivity of uranium to initial local conditions and can be used to improve our understanding of the kinetics that drive uranium speciation and ultimately inform fallout transport models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Permeability Evolution of Fractures in Shale in the Presence of Supercritical CO 2

We experimentally examined the impact of supercritical carbon dioxide (scCO 2 ) on fracture permeability and fracture surface characteristics in shales of various compositions. We measured permeability and fracture normal displacement at different effective stresses using both argon and scCO 2 as pore fluids. Both natural fractures and saw cuts in intact samples were used in our study. The permeability and fracture normal displacement decrease after multiple loading/unloading cycles, apparently as the result of inelastic compaction. In samples with high carbonate content, we observe an increase in permeability and fracture surface degradation after more than 3.5 days of exposure to scCO 2 , apparently due to carbonate dissolution. We demonstrate that the sensitivity of permeability to effective normal stress correlates well with the fracture normal displacement. The dependence of permeability on effective normal stress increases after exposure to scCO 2 as the fracture surface became more compliant.

58 GEOSCIENCES↗

Supernovae and their host galaxies – VII. The diversity of Type Ia supernova progenitors

ABSTRACT We present an analysis of the light curve (LC) decline rates (Δm15) of 407 normal and peculiar supernovae (SNe) Ia and global parameters of their host galaxies. As previously known, there is a significant correlation between the Δm15 of normal SNe Ia and global ages (morphologies, colours, and masses) of their hosts. On average, those normal SNe Ia that are in galaxies from the Red Sequence (early-type, massive, old hosts) have faster declining LCs in comparison with those from the Blue Cloud (late-type, less massive, younger hosts) of the colour–mass diagram. The observed correlations between the Δm15 of normal SNe Ia and hosts’ parameters appear to be due to the superposition of at least two distinct populations of faster and slower declining normal SNe Ia from older and younger stellar components. We show, for the first time, that the Δm15 of 91bg- and 91T-like SNe is independent of host morphology and colour. The distribution of hosts on the colour–mass diagram confirms the known tendency for 91bg-like SNe to occur in globally red/old galaxies, while 91T-like events prefer blue/younger hosts. On average, the youngest global ages of 02cx-like SNe hosts and their positions in the colour–mass diagram hint that these events likely originate from young population, but they differ from 91T-like events in the LC decline rate. Finally, we discuss the possible explosion channels and present our favoured SN Ia models that have the potential to explain the observed SN–host relations.

Hakobyan, A. A.↗

Lagrangian and Eulerian drag models that are consistent between Euler-Lagrange and Euler-Euler (two-fluid) approaches for homogeneous systems

The undisturbed flow of a particle is of fundamental importance since it controls both the undisturbed flow force and the perturbation force (which includes quasisteady, added-mass, and history forces). Here we use the pairwise interaction extended point particle framework to evaluate the undisturbed flow of each particle through superposition of the perturbation flow induced by all its neighbors. This approach allows calculation of various statistics related to undisturbed fluid velocity under conditions of both stationary and nonstationary particles. In a random distribution of stationary particles, while the macroscale undisturbed flow is slowly varying, the microscale undisturbed flow that arises due to the perturbation flow of neighbors varies substantially from one particle to another and this in turn leads to large variation in the hydrodynamic force exerted on the particles. The effect of particle motion is generally to increase the particle-to-particle variation in the undisturbed fluid velocity of the particles. We observe that this increase is greater for the transverse component than for the streamwise component. As a result, with increasing random particle motion, the distribution of undisturbed fluid-velocity fluctuation becomes isotropic. Three different normalized forces are defined: $Φ^L$ is the Lagrangian normalized force on an individual particle suitable for application in a microscale-informed Euler-Lagrange simulation, $Φ^E$ is the Eulerian normalized average force suitable for application in an Euler-Euler simulation, and $Φ^{LE}$ is the Lagrangian normalized force on an individual particle suitable for application in the standard Euler-Lagrange simulation. In this work, we establish precise relations between these different definitions. The drag laws developed based on particle-resolved direct numerical simulation results and experiments are appropriate for application only as the Eulerian normalized average force. We introduce the force consistency relation and use it to obtain an expression for $Φ^L$, which when applied to each particle and averaged over all the particles equals $Φ^E$. The results are first obtained in the limit of stationary particles and then extended to the general case of nonstationary particles.

42 ENGINEERING↗

Detecting Anomalies in Time Series Using Kernel Density Approaches

This paper introduces a novel anomaly detection approach tailored for time series data with exclusive reliance on normal events during training. Our key innovation lies in the application of kernel-density estimation (KDE) to scrutinize reconstruction errors, providing an empirically derived probability distribution for normal events post-reconstruction. This non-parametric density estimation technique offers a nuanced understanding of anomaly detection, differentiating it from prevalent threshold-based mechanisms in existing methodologies. In post-training, events are encoded, decoded, and evaluated against the estimated density, providing a comprehensive notion of normality. In addition, we propose a data augmentation strategy involving variational autoencoder-generated events and a smoothing step for enhanced model robustness. The significance of our autoencoder-based approach is evident in its capacity to learn normal representation without prior anomaly knowledge. Through the KDE step on reconstruction errors, our method addresses the versatility of anomalies, departing from assumptions tied to larger reconstruction errors for anomalous events. Our proposed likelihood measure then distinguishes normal from anomalous events, providing a concise yet comprehensive anomaly detection solution. The extensive experimental results support the feasibility of our proposed method, yielding significantly improved classification performance by nearly 10% on the UCR benchmark data.

Frehner, Robin↗

Impedance Response Influenced by Variability in the Random Distribution of Physical Properties of Coated Materials in Two-Dimensional Space

Heterogeneous physical characteristics of a system featuring a single-layer film on a metallic surface have been explored via its impedance response. The Nyquist plot showed a distorted semicircle, indicative of the system's unique distribution characteristics. Utilizing a copula-based probability method, a two-dimensional deterministic impedance model was successfully integrated, accounting for spatial physical properties such as permittivity and electrical conductivity. This strategy enabled in-depth exploration and mechanistic quantification of a broad spectrum of properties. A quantitative understanding of impedance signal alterations, characterized by normally or log-normally correlated variables, was achieved through the variation in aspect ratio and characteristic frequency of the impedance spectra. Log-normally distributed electrical properties provided a superior representation of the distorted impedance spectra. As coefficient of variation (CV) values fluctuated, the aspect ratio and characteristic frequency showed heightened sensitivity to log-normal permittivity compared to log-normal electrical conductivity. Notably, a marked positive linear correlation between electrical properties resulted in an impedance response that approximated perfect semicircular spectra. The variability in the electrical properties' distribution was demonstrated by considering the correlation coefficient between electrical conductivity and the z-direction position. Furthermore, the highest aspect ratio of the impedance spectra was observed when the electrical conductivity was randomly distributed across the z-direction space.

36 MATERIALS SCIENCE↗

Response patterns of simulated corn yield and soil nitrous oxide emission to precipitation change

Abstract Background Precipitation plays an important role in crop production and soil greenhouse gas emissions. However, how crop yield and soil nitrous oxide (N 2 O) emission respond to precipitation change, particularly with different background precipitations (dry, normal, and wet years), has not been well investigated. In this study, we examined the impacts of precipitation changes on corn yield and soil N 2 O emission using a long-term (1981–2020, 40 years) climate dataset as well as seven manipulated precipitation treatments with different background precipitations using the DeNitrification-DeComposition (DNDC) model. Results Results showed large variations of corn yield and precipitation but small variation of soil N 2 O emission among 40 years. Both corn yield and soil N 2 O emission showed near linear relationships with precipitation based on the long-term precipitation data, but with different response patters of corn yield and soil N 2 O emission to precipitation manipulations. Corn yield showed a positive linear response to precipitation manipulations in the dry year, but no response to increases in precipitation in the normal year, and a trend of decrease in the wet year. The extreme drought treatments reduced corn yield sharply in both normal and wet years. In contrast, soil N 2 O emission mostly responded linearly to precipitation manipulations. Decreases in precipitation in the dry year reduced more soil N 2 O emission than those in the normal and wet years, while increases in precipitation increased more soil N 2 O emission in the normal and wet years than in the dry year. Conclusions This study revealed different response patterns of corn yield and soil N 2 O emission to precipitation and highlights that mitigation strategy for soil N 2 O emission reduction should consider different background climate conditions.

54 ENVIRONMENTAL SCIENCES↗

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗

Performance Comparison of Circular and Spherical Error Probable Estimators

This report compares the performance of three Circular Error Probable (CEP) estimators: the Grubbs-Patnaik estimator, a new, non-iterative, radial-integration estimator, and a median estimator. It also compares the performance of two Spherical Error Probable (SEP) estimators. The performance of each estimator is assessed in terms of bias, uncertainty, robustness, and computational complexity. Robustness is evaluated with respect to outliers, variations in the underlying statistical distribution characterizing munition impact positions, and impact-position measurement errors. The performance assessments indicate the radial-integration and Grubbs-Patnaik estimators perform nearly identically providing the statistical distribution of impact-position coordinates is jointly normal with zero means. In that case, both estimators outperform the median estimator by about 2% relative to the true CEP in terms of estimator uncertainty. The bias performance of the radial-integration and median estimators is close to zero for jointly normal impacts, however, the Grubbs-Patnaik estimator can be significantly biased for jointly normal impacts with non-zero means. When the statistical distribution characterizing impact positions is known, but not jointly normal, the radial-integration estimator is superior. In this case, the median estimator also outperforms the Grubbs-Patnaik estimator but is not quite as good as the radial-integration estimator. If the statistical distribution characterizing impacts is unknown and not jointly normal, or if distribution parameters are difficult or impractical to estimate, or if test data is corrupted with outliers, then the median estimator dramatically outperforms the other estimators, especially in terms of estimation bias. Unexpectedly, measurement noise did not significantly degrade the performance of any of the estimators, except for cases with signal to noise ratios less than five. Although the Grubbs-Patnaik estimator has remained the gold standard for CEP estimation for over half a century, the performance assessments indicate the new, non-iterative, radial-integration estimator and the median estimator offer significant advantages and, in most practical real-world conditions, are superior estimators. These estimators are also useful for SEP estimation whereas the Grubbs-Patnaik estimator does not extend to three dimensions.

97 MATHEMATICS AND COMPUTING↗

Hybrid Analytics Solution to Improve Coal Power Plant Operations

This project focused on developing advanced methods for thermal performance monitoring of a coal-fueled power plant. The specific goal was to develop and demonstrate a new thermal performance monitoring approach using a hybrid model that integrates a physics-based heat balance model with a machine learning-based pattern recognition model. The hybrid model enables increased accuracy and scope of the thermal analysis and an improved ability to monitor and detect changes in plant operation. This new approach takes full advantage of the individual model capabilities and creates an important new set of capabilities not previously possible using the two types of models separately. Using the heat balance model, a rich set of derived parameters (virtual sensors) are calculated from the measured plant operating data at each time point. The combined measured and derived data values are used by machine learning algorithms to create pattern recognition models over the range of normal unit operation. To create the monitoring models, historical data from normal operation of the plant is first processed by the heat balance model to compute the derived parameter data. The result is a greatly expanded set of normal operating data that can be used as input to create the pattern recognition model. Once the models are calibrated for normal operation, the hybrid model is suitable for use in continuous online monitoring. During online monitoring, new plant operating data is processed first by the heat balance model and then by the pattern recognition model. Results from the pattern recognition model quantify the deviation of each measured or derived parameter from its expected value in normal operation. The hybrid models can detect abnormal changes in plant operating data with very high accuracy and sensitivity. When abnormal behavior is detected, alerts are generated automatically for evaluation by the plant monitoring staff. The new hybrid solution product was developed and verified in the performance of the project. The hybrid solution was tested first in a simulation environment that mimicked the plant data systems and infrastructure used by U.S. power generating plants and utilities. The hybrid solution was then deployed for real-time, online monitoring of an operating coal-fueled power plant at a field test site. Field testing demonstrated that all hybrid solution development objectives were accomplished. The project work was based on combining the capabilities of two existing software products to create the new hybrid solution product. One of these was the existing MapEx® heat balance product and the other was the existing SureSense® advanced pattern recognition product. Each of these separate products was assessed to be at a Technology Readiness Level (TRL) of 9 at the start of the effort. The hybrid solution product was assessed to be at a TRL of 2 at the start of the project based on early feasibility work by the project team. At completion of the field testing performed in the project, the hybrid solution product was assessed to be at a TRL of 7. The project team expects that the hybrid solution product will be deployed commercially and will achieve a TRL of 9 within one year after completion of the project.

01 COAL, LIGNITE, AND PEAT↗

SPARC 3D Field Physics and Support of the Non-Axisymmetric Coil Assessment

Commonwealth Fusion Systems (CFS) is deploying SPARC, a net-energy tokamak, by 2025. A notable challenge facing all tokamak approaches to fusion energy production is maintaining a stable plasma and thereby steady energy production. “Disruptions” to plasma operation are often encountered in devices when they operate at high normalized pressure, high normalized density, or high normalized current. SPARC is unique in that it is expected to demonstrate Q>=2 in a plasma with low normalized pressure and density, but with a more modest current limit buffer sufficient to avoid disruptions. Due to both the high normalized current and high magnetic field of the SPARC design, it is expected to decrease its relative resilience to instabilities driven by non-axisymmetry in the externally applied magnetic field, termed ‘error fields’. To prevent error field driven instabilities, a primary line of defense is strict engineering tolerances in magnet fabrication and installation. A secondary line of defense is a purpose-built set of magnetic coils to correct these errors, termed error field correction coils (EFCCs). The overarching goals of this INFUSE project were to evaluate the effectiveness of the EFCCs in the SPARC design and provide guidance as to how to improve this effectiveness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced Component Performance Study: Motor-Driven Pumps 1998-2020

This report presents an enhanced performance evaluation of motor-driven pumps (MDPs) at U.S. commercial nuclear power plants. The data used in this study are based on the operating experience failure reports from calendar year 1998 through 2020 as reported in the Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS). The MDP failure modes considered for standby systems are fail to start (FTS), fail to run (FTR) for one hour of operation (FTR≤1H), FTR after one hour of operation (FTR>1H), and for normally running systems FTS and FTR. An eight-hour unreliability estimate is also calculated and trended. The component reliability estimates and the reliability data are trended for the most recent 10-year period while yearly estimates for reliability are provided for the entire study period. The following increasing trends were identified for MDPs for the most recent 10-year period: • Standby MDP frequency of start demands (demands per reactor year) • Standby MDP frequency of FTR≤1H hours (hours per reactor year) • Standby MDP frequency of FTR>1H hours • Normally running MDP frequency of run hours. The following decreasing trends were identified for MDPs for the most recent 10-year period: • Standby MDP FTR≤1H failure probability • Normally running MDP FTR failure rate • Standby MDP unavailability • Normally running MDP total unreliability (8-hour mission) • Standby MDP frequency of FTR≤1H events (failures per reactor year) • Normally running MDP frequency of FTR events.

99 GENERAL AND MISCELLANEOUS↗

High Solids Performance Testing in a Scaled TSCR System

The Tank Side Cesium Removal (TSCR) project is a technology demonstration that will pretreat Hanford tank waste supernatant in support of the Direct Feed Low-Activity Waste (DFLAW) mission. The TSCR system employs two key separation technologies: dead-end filtration (DEF) and ion exchange (IX) using crystalline silicotitanate (CST) media. DEF will be used to remove undissolved solids from tank waste to protect the functionality of the IX columns and the IX system will remove Cs-137 from tank waste. The separation technologies (DEF and IX) used in TSCR are technically mature and have also been successfully deployed at the Savannah River site in a similar facility known as the Tank Closure Cesium Removal (TCCR) system. While testing with simulants and real waste has been successfully performed under conditions expected during the initial operation of TSCR, test data is absent for assessing off normal high solids loading that may be in the TSCR waste feed. Normal TSCR treatment operations are expected to handle wastes with solids content on the order of 200 ppm, and off normal solids loading could be much larger than the nominal level. The testing program described in this report was conducted to understand the consequence of operating the TSCR system at elevated solids loadings up to the high-solids limit of 15,000 ppm [i.e., 1.5-wt%] identified in the TSCR design basis. Although the system is not required to make throughput above the nominal solids loading, the testing was intended to provide important information related to potential off normal operations. At off normal levels near the high-solids limit, there are potential implications for TSCR performance in the areas of throughput, DEF pressure drop, filter backflush frequency, and IX column pressure drop. In addition, intrusion of solids into the IX column was postulated to impact the Cs-137 loading behavior by promoting channeling or flow maldistribution in the column; since the magnitude of the postulated effect was unknown, assessing it was also of interest. The testing was performed using representative waste simulants and a prototypic, integrated TSCR system designed and assembled specifically to conduct the high solids performance assessment. Overall, the scaled TSCR testing demonstrated that full-scale unit operations can succeed in fulfilling their processing objectives in the presence of solids up to 3,000 ppm, but there are potential performance challenges to filtration operations at solids loadings as low as ~500 ppm. The severity of the challenge is likely to be dependent on the type and size distribution of solids, of which the current testing only examined a single type and size distribution. To provide some flexibility for future full-scale operations, the results of the testing suggest two possible risk reduction strategies that can be implemented without any changes in TSCR design or configuration. One option would be to enact an administrative limit on the solids loading to protect TSCR from feeds that are likely to require a high DEF swap frequency. Another option is to permit operation of the DEFs at differential pressures greater than 2 psid before swapping filters. The selection of a higher differential pressure target is not anticipated to adversely impact DEF backflushing efficacy and would reduce both swap frequency and the amount of waste sent to AP-108.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Enhanced Component Performance Study: Turbine Driven Pumps 1998–2022

This report presents an enhanced performance evaluation of turbine driven pumps (TDPs) at U.S. commercial nuclear power plants. The data used in this study are based on the operating experience failure reports from calendar year 1998 through 2022 as reported in the Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS). The TDP failure modes considered for standby systems are fail to start (FTS), fail to run (FTR) for one hour of operation (FTR=1H), FTR after one hour of operation (FTR>1H), and for normally running systems FTS and FTR. An eight hour unreliability estimate is also calculated and trended. The component reliability estimates and the reliability data are trended for the most recent 10 year period while yearly estimates for reliability are provided for the entire study period. No increasing trends were identified for TDPs for the most recent 10 year period: The following decreasing trends were identified for TDPs for the most recent 10 year period: • Standby TDP FTR>1H failure rate • Normally running TDP FTR failure rate • Standby TDP unavailability • Standby MDP total unreliability (8-hour mission) • Standby TDP frequency of start demands (demands per reactor year) • Standby TDP frequency of FTR=1H hours (hours per reactor year) • Standby TDP frequency of FTR>1H events (failures per reactor year) • Normally running TDP frequency of start demands • Normally running TDP frequency of run hours • Normally running TDP frequency of FTR events.

99 GENERAL AND MISCELLANEOUS↗

Process Anomaly Detection for Sparsely Labeled Events in Nuclear Power Plants

An essential aspect of online monitoring, subtle anomaly detection increases the detection lead time for equipment failure and enables a nuclear power plant (NPP) to mitigate unexpected partial or full outages, resulting in significant cost saving to the plant. Once an anomaly is detected by plant staff, its cause and severity are investigated. Because the vast majority of anomalies require some level of investigation, including some that require time-consuming examination, before they are passed over to the engineering organization for further analysis, plants are often equipped with tools to assist the staff in performing anomaly detection. Those tools operate as a black box and are often based on statistical methods that establish sensor correlations using preconfigured mathematical models and flag correlation deviations as anomalies. Due to the number of anomalies detected at a given NPP on a daily basis, a significant number of flagged anomalies usually await examination for days or weeks. A primary cause of this backlog is that the methods used by the tools generate many false positives. Though this is usually attributed to oversensitive model settings due to very narrow normal operation bands, it can also be associated with the model development being inadequate for the process being monitored, or with missing model inputs that could have explained misclassified positives. The performance of anomaly detection tools impacts their plant acceptance and utilization, especially when the effort to address false positives generated by the tool depletes the value or cost saved by using that tool. Thus, means to advance anomaly detection performance have been investigated by the Department of Energy’s Light Water Reactor Sustainability program. Previous and ongoing efforts have targeted unsupervised machine-learning (ML) methods, which do not require the labeling of any data fed into the ML model. By contrast, in supervised anomaly detection methods, every data point is labeled as either a normal or abnormal process condition, and the model is trained to replicate the classification process. Supervised methods usually outperform unsupervised methods, due to the added value in differentiating normal from anomalous states of the monitored process. An NPP’s corrective action program requires it to track and document, via a dedicated report, the resolution of any issues that occur within the plant. Once created, each report is reviewed by a plant screening committee, and several classifications and decisions are made. Recently, a collaborating NPP developed an artificial intelligence and ML-based classifier to categorize a condition report (CR) into classes that can serve to label the data as normal or anomalous. Applying CRs as labels represents a semi-supervised use case. Semi-supervised ML assumes that labels exist for some data points (i.e., labeled anomalies, in this case) but not for the rest. In this effort, semi-supervised ML methods were used to fuse data from CRs with anomaly detection methods in order to test the hypothesis that partially labeled anomalies would improve the accuracy of the anomaly detection methods. Specifically, two methods were used. The first is the deep Semi-supervised Anomaly Detection (deep SAD) method, which can handle labels ranging from fully unsupervised to fully supervised cases. The second is a newly designed ML method developed specifically for this effort and referred to as the high-order feature (HOF)-based method. To evaluate these two methods in controlled environments, synthetic data generators were developed and used. The first datasets used a spring-mass-damper (SMD) system simulator commonly found in mechanical engineering references. This was used to create two use cases: a one- and a three-mass system. Anomalies were introduced by changing the spring and damper coefficients while the system was actuated by random forces. The second datasets used the commercial Dymola-Modelica software to build a simplified nuclear reactor model. Anomalies were added in the form of corrupted sensor readings and/or control commands. The deep SAD method was tested using the SMD system, while the HOF method was tested using both datasets. Application of the deep SAD semi-supervised ML method demonstrated that labels can generate increased confidence in detecting true anomalies. This helped increase the number of true positives and decrease the number of false negatives—something that would aid in addressing the backlog of possible anomalies. Application of the HOF method demonstrated that labels can aid in down selecting from a candidate set of features to a more optimal subset in order to better differentiate between normal and anomalous conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗

Deep Learning Based Frequency Stability Assessment in Power Grid with High Renewables

Frequency stability assessment is one critical aspect of power system security assessment. Traditional N-1 screening method is based on the simulations of a few typical daily and seasonal operation scenarios. However, the increasing integration of inverter-based renewables and the retirement of conventional synchronous generators result in decreasing system inertia and growing complexity of system operating conditions. Selecting a few typical operation scenarios cannot cover all operating conditions, and the time-domain simulation of all operation conditions requires tremendous time. This paper proposes a more efficient frequency stability assessment method based on deep learning. The affinity propagation clustering algorithm is used to divide the dataset into different clusters, so the selected dataset for training can cover the diversified operating conditions as much as possible. Also, feature normalization is applied to both the training dataset and testing dataset in order to remove any unnecessary bias. Especially, trained model based on full dataset normalization has bounded error in the prediction. The case study on the reduced 240-bus WECC system demonstrates that the proposed method can predict accurate frequency nadir with limited training dataset. The deep learning model using the revised feature normalization can predict more accurate frequency nadir than that using the traditional feature normalization and has very small maximum prediction error.

affinity propagation↗

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility. Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson’s R and Spearman’s ρ coefficients as well as Kendall’s τ and Top–Down correlation. We also test the behavior of association measures, including the coefficient of determination R 2 , Kendall’s W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spear man’s ρ, Kendall’s τ, and Kendall’s W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R 2 coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships. Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

59 BASIC BIOLOGICAL SCIENCES↗