Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic program analysis”

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 595 records · Page 33

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

CO Oxidation on Au/TiO 2 : Condition-Dependent Active Sites and Mechanistic Pathways

We present results of ab initio electronic structure and molecular dynamics simulations (AIMD), as well as microkinetic modeling of CO oxidation catalyzed by TiO2 supported Au nanocatalysts. A coverage-dependent microkinetic analysis, based on energetics obtained with density functional methods, shows that the dominant kinetic pathway, activated oxygen species, and catalytic active sites are all strongly depended on both temperature and oxygen partial pressure. A dynamic single atom catalytic (DSAC) mechanism for CO oxidation is proposed to be the dominant (but not the exclusive) pathway for this reaction under oxidative conditions below T < 400 K. This mechanism proceeds via an Au+-CO species that transfers from the Au-cluster to a surface oxygen ad-atom, where it reacts with the TiO2 support with a Mars van Krevelen mechanism to form CO2 and subsequently reintegrates back into the gold cluster to complete the catalytic cycle. At 300 = T = 600 K oxygen-bound single Au+CO sites and the perimeter sites of the Au-nanoparticle work in tandem to optimally catalyze the reaction. Above 600 K a variety of alternate pathways associated with the perimeter sites of the Au nanoparticle and/or the single atom site are found to be the most active pathways. Overall our results highlight the importance of the dynamic behavior of interfacial area between the catalysts and the support, which allows for a high degree of adaptation. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences and performed at the Pacific Northwest National Laboratory (PNNL). PNNL is a multi-program national laboratory operated for Department of Energy by Battelle. Computational resources were provided at W. R. Wiley Environmental Molecular Science Laboratory (EMSL), a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at PNNL and the National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory.

Wang, Yang-Gang↗

Field-Coupled Water Splitting with Metal-Free Donor–Acceptor Covalent Organic-Framework Junctions

Advancing metal-free electrocatalysts for hydrogen and oxygen evolution reactions (HER/OER) across acidic and alkaline media requires coordinated control of intermediate binding thermodynamics, interfacial charge delivery, and near-electrode transport dynamics. Here, we design amide-linked benzene–triazine covalent organic frameworks (BTA/TzTA-Hz COFs) and integrate them with carbon nanotubes (CNTs) to form COF–CNT junctions that establish a built-in interfacial electric field. Density functional theory (DFT) and electrostatic potential maps indicate complementary active motifs, with benzene-proximal fragments associated with HER and triazine-proximal motifs associated with OER. CNT integration shifts the contact-potential difference by ≈0.20 V, while operando electrochemical impedance spectroscopy suggests partially separable high-frequency junction-charging and lower-frequency Faradaic/transport responses. A 300 mT static magnetic field lowers the HER and OER overpotentials by tens of millivolts. Under anodic bias, the effective interfacial charging capacitance increases, and Mott–Schottky analysis shows an apparent ∼0.15 V flat-band shift with an essentially unchanged slope. Together, these observations are consistent with field-perturbed interfacial charging and altered bias partitioning. Field-dependent impedance and bubble imaging are consistent with magnetohydrodynamic convection that promotes bubble detachment and near-electrode mass transport for both half-reactions, and they reveal an OER-specific high-frequency perturbation under anodic bias. Under field, the heterostructure reaches an OER onset overpotential of ∼261 mV and requires an overpotential of 366 mV at 10 mA cm –2 in alkaline electrolyte. These results illustrate how reticular-framework chemistry, junction engineering, and both built-in and applied fields can program reactivity through interfacial electrostatics and near-electrode transport in organic-framework electrocatalysts.

Garcia-Enriquez, Lissette [Univ. of Texas at El Pa↗

TEAMER: Biofouling Analysis for Wave Energy Piston Design - Load Cell Data

Biofouling and corrosion are a major concern for all ocean-deployed components, especially when mechanical motion is involved. Triton has developed the concept of a biofouling mitigation seal as part of the piston sealing assembly for the Triton Wave Energy Converter (TSI-WEC). This mitigation seal has the purpose of preventing the formation of a biofilm on the inside of the piston cylinder. It is hypothesized that the prevention of a biofilm will reduce the amount of macro-biofouling that can occur in the piston assembly. The mitigation seal can also reduce the wear on the main dynamic seal, helping to maintain smooth operation and water-tightness. The cylinder is made from a thermoset composite epoxy, which is resistant to corrosion. However, no studies have researched the material's performance with biofouling. Triton placed two prototype Power Take-Off (PTO) assemblies in a PNNL biofouling tank, one with a biofouling mitigation seal and one without, allowing for an evaluation of seal effectiveness at the prevention of biofouling. In actual WEC operation, wave action would react against the piston, which would drive the linear actuator and electric generator, providing electrical power. In the test setup, this was reversed; a linear actuator was powered to drive the piston in a consistent motion within the cylinder. There are two assemblies: one has a biofouling mitigation seal, the other (control) does not. The following data encompasses a 4 month test period, with load cells being used to monitor piston friction force. Results from this testing will be used to improve seal design and material selection, mitigating risk of premature failure during open water testing and evaluation. This project is part of the TEAMER RFTS 3 (request for technical support) program.

16 TIDAL AND WAVE POWER↗

Behavior of U 3 Si 2 Fuel and FeCrAl Cladding under Normal Operating and Accident Reactor Conditions

As part of the Department of Energy's Nuclear Energy Advanced Modeling and Simulation program, an Accident Tolerant Fuel High Impact Problem was initiated at the beginning of fiscal year 2015 to investigate the behavior of \usi~fuel and iron-chromium-aluminum (FeCrAl) claddings under normal operating and accident reactor conditions. The High Impact Problem was created in response to the United States Department of Energy's renewed interest in accident tolerant materials after the events that occurred at the Fukushima Daiichi Nuclear Power Plant in 2011. The High Impact Problem is a multinational laboratory and university collaborative research effort between Idaho National Laboratory, Los Alamos National Laboratory, Argonne National Laboratory, and the University of Tennessee, Knoxville. This report primarily focuses on the engineering scale research in fiscal year 2016 with brief summaries of the lower length scale developments in the areas of density functional theory, cluster dynamics, rate theory, and phase field being presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Projected precipitation changes within the Great Lakes and Western Lake Erie Basin: a multi‐model analysis of intensity and seasonality

ABSTRACT The Great Lakes region encompasses the largest freshwater lake network in the world and supports a diverse network of agriculture, transportation, and tourism. Recently, Lake Erie has experienced increased hypoxia events, which have been attributed to agricultural practices and changes in run‐off. Here we examine the projected changes in extreme precipitation events to address concerns regarding regional agriculture, surface run‐off, and subsequent water quality. Precipitation projections within the overall Great Lakes Basin and the Western Lake Erie Basin subregion are examined using climate model simulations of varying spatial resolutions to understand historical precipitation and projected future precipitation. We develop three model ensembles for the historical period (1980–1999) and the mid‐century (2041–2060) that cover a range of spatial resolutions and future emissions scenarios, including: (1) 12 global model members from the fifth Climate Model Intercomparison Project ( CMIP5 ) using Representative Concentration Pathway ( RCP ) 8.5, (2) ten regional climate model ( RCM ) members from the North American Regional Climate Change Assessment Program driven by CMIP3 global models using the A2 emissions scenario, and (3) two high resolution RCM simulations ( RCM4 ) driven by CMIP5 global models using the RCP 8.5 scenario. For the historical period, all model ensembles overestimate winter and spring precipitation, and many of the models simulate a summer drying that is not observed. At mid‐century, most of the models predict a 10–20% increase in precipitation depending on the time of year. Daily probability distribution functions from three model ensembles reveal spring seasonal increases in high precipitation event probabilities when compared to the historical period, suggesting an increase in the frequency of high intensity precipitation at mid‐century. Overall, the presence of lakes or higher spatial resolution does not ensure improved representation of historical processes, and more complex interactions between large‐scale dynamics, local feedbacks, and physical parameterizations drive the model spread.

Basile, Samantha J.↗

Exploring the Nucleon Structure via Deep Electroproduction Processes

Understanding the internal structure of the nucleon is a fundamental goal of modern physics, which aims at a comprehensive framework describing the internal dynamics of quarks and gluons. Among other structure functions, Generalized Parton Distributions (GPDs) offer a powerful framework for describing the nucleon dynamics by correlating the longitudinal momentum and the transverse position of its internal partons. Such a correlation provides a three-dimensional picture of the nucleon and enables access to fundamental properties, including the internal pressure distributions and the parton’s angular momentum contribution to the nucleon's total spin, thereby playing a central role in resolving the nucleon spin puzzle. At the Thomas Jefferson National Facility (JLab), polarized electron beam experiments allow for probing GPDs through the measurement of hard exclusive processes. Among the cleanest experimental channels, we find the electro-production of a real photon through the Deeply Virtual Compton Scattering (DVCS) mechanism. The first data-taking period of the CLAS12 program, taking place in 2018, allowed for unique DVCS Beam Spin Asymmetry (BSA) measurements in the phase space covered by a 10.6 GeV polarized electron beam impinging on an unpolarized liquid hydrogen target. Although detecting all final-state particles ensures exclusivity of the process, conservation laws indicate that it is not mandatory. I adopt an approach omitting the direct detection of the recoil proton, providing a simplified yet effective event selection strategy that boosts statistics and gives access to a larger phase space sensitive to the underlying GPD dynamics through BSA and cross section measurements. The Double DVCS (DDVCS) process promises a dedicated mapping of GPDs. Contributing to the electro-production of a lepton pair cross-section, the DDVCS reaction extends DVCS by allowing the final-state photon to be virtual, enriching the kinematic phase space and providing unique access to the internal correlations encoded by GPDs. A feasibility study is conducted to assess the potential of future DDVCS measurements at Jefferson Lab and the future Electron-Ion Collider (EIC). While Jefferson Lab will provide DDVCS measurements in the valence region through the SoLID$\mu$ and $\mu$CLAS12 experimental projects, in the long term, the EIC will provide complementary measurements in the sea region, both accessing unprecedented information about GPDs in a phase space region otherwise inaccessible. Taken together, these investigations demonstrate both the current capabilities and future opportunities for probing GPDs through exclusive processes. The experimental analysis of DVCS at CLAS12 provides precise measurements within an established framework, while the phenomenological study of DDVCS opens the door to richer and more comprehensive explorations with future detectors and facilities.

Alvarado, Juan [Université Paris-Saclay: Gif-sur-Y↗

Modular Controllable Transformers (MCT)

Large Power Transformers (LPTs) are at the heart of the electricity infrastructure and are critical pieces in today’s grid. Aging assets, long turn-around times, transportation delays and dependence on foreign manufacturing for these components have created conditions for points of failure across the system. Moreover, there is no available dynamic control on these passive sections of the grid. For all the automation embedded on the low voltage side of the low voltage side of the system there is absolutely no flexibility on the high voltage side. Changes in power flow patterns are brought about by altering dispatch patterns for generators which requires a coordinated effort and is often time consuming. Traditional FACTS based approaches such as UPFC or HVDC light involve massive modifications to the existing infrastructure. Recent instances of physical attacks on LPTs and the slow recovery process associated with fixing or replacing these damaged units has created an increased awareness for this issue. LPT designs are highly customized owing to specific field requirements. Long lead times coupled with transportation delays make procuring and installing LPTs a lengthy process. In an attempt to minimize the recovery time following an LPT failure, transformer manufacturers have developed mobile transformer units. These units typically consist of either a single-phase or three-phase transformer mounted on a truck and characterized by compact form factors. By reducing the assembly time and transport time these units offer a fast and temporary solution to LPT failures. Several efforts led by DOE, Edison Electric Company (EEI) and North American Electric Reliability Corporation (NERC) have led to the RecX program. RecX transformers consist of single phase transformers transported on specialized assemblies to enable fast restoration and minimal installation time. NERC and EEI have instituted transformer sharing programs like Spare Equipment Database (SED) as well as Spare Transformer Equipment Database (STEP) which encourage utilities to maintain and share a fleet of LPTs. Mobile transformers are typically limited to less than 100 MVA to enable transportation on a trailer. Also, the efficiency of mobile transformer units and RecX transformers is lower since they are temporary solutions. The ideal replacement solution would be one that could utilize a set of standardized transformers while incorporating flexibility to influence metrics like power flow, apparent impedance and voltage profiles. In view of these issues the team at Georgia Tech proposed a concept called Modular Controllable Transformers (MCTs). This project was commissioned as a part of DOE’s Transformer Resiliency and Advance Component (TRAC) program. The approach consists of splitting one large unit rated upwards of 100 MVA with smaller modular, standardized units. Each of these smaller units are augmented with a converter to add power flow, voltage and impedance control. Splitting the power level into multiple units has two important implications. By making smaller units, the transportation times can be reduced drastically. Moreover, even though these units have very specific impedance and other parameters, the augmented converters can make it emulate the unit being replaced. Thus, faster replacement times and added controls are achieved making this an extremely resilient approach. In the event of one modular units failure, the other two controllable units can operate in the most resilient fashion keeping the net system overloads to a minimum. Moreover, for the same probability of failure only a fraction of the system capacity is lost. The converter is augmented on to the tertiary winding of the transformer making it retrofittable in the field without massive modifications to the system. The material proposed in this report highlights design considerations to minimize recovery time as well as results displaying overall improvement in resiliency. The cost analysis of this approach shows significant savings compared to UPFC and HVDC lite based approaches. Extensive analysis was conducted on numerous aspects of this project over the course of this project. Simulation studies were conducted to prove the 5MVA converter design. Challenging BIL constraints and transient issues were addressed over the course of this project. Simulations on EMTP platform presented in consequent sections were used to finely tune the control for the same. A final design for a back-to-back 5MVA neutral-point-clamped (NPC) converter was completed in this project. Furthermore, switching components like high power IGBTs from IXYS Westcode Press-Pack series and ABB Stackpack series were identified for this application. These devices exhibit a blocking voltage of 6.5 kV and a nominal collector current of 900 A. Potential housing options as well as existing commercial option were identified for this converter design. A 56 MVA modular transformer unit was designed with Delta Star Inc. with consultation and inputs from all project members. The transformer cooling systems were chosen to be Oil Natural Air Forced (ONAF) to optimize ease of transport. A full Product Requirement Document (PRD) for the converter and the transformer was generated. The team also investigated key aspects of resiliency and how to quantify the same. A metric was developed to analyze resiliency and numerous simulation cases were presented. Increased (N-1) resiliency and outage resistance was showcased through simulation on the Texas Large scale power system. Dynamic control capabilities were also showcased on numerous system level simulations. An IEEE 30 bus system was used to show the ability of the MCT approach to reduce line congestion. Power flow control; a part of the MCTs abilities, showcased a 29% reduction in transformer loading on the same system. Full control over tie line flows was demonstrated on a modified IEEE 13 bus system. Dynamic LTC capabilities were shown through simulation and were compared to traditional LTCs. Furthermore, the ability of the MCTs to alter apparent transformer impedance was simulated to promote standardization of modular transformer units while retaining the ability to emulate a varied range of impedance characteristics. Thus granular control over power flow, voltages, impedances and was presented. Regulatory analysis conducted by Oak Ridge National Laboratory (ORNL) further bolstered the efficacy of the MCT approach and verified the added resiliency achieved due to this approach. As a conclusion, it was proved that the resilience improvement obtained from using the MCT approach is justified by just one failure in a 30-year period. This means that the benefits obtained from using the MCT approach far outweigh the cost of setting this approach up. As a part of the future efforts, the team is proposing building a (12.47-24) kV/400 kVA scaled down unit to control 5MVA of power flows. Emphasis would be laid on using commercial components to speed up the build of these units. Further the actual transport of the modular transformer unit would be a part of this operation and would provide insight into transportation issues. The approach shows much promise and adds key dynamic control on passive sub transmission and transmission sectors.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Development and Application of a Data-Driven Methodology for Validation of Risk-Informed Safety Margin Characterization Models

The document is the Final Technical Progress Report for the Nuclear Energy University Program’s Integrated Research Project (IRP) on “Development and Application of a Data-Driven Methodology for Validation of Risk-Informed Safety Margin Characterization Models” that supports the LWR Sustainability Program’s RISMC R&D Pathway. The project goal is to develop and demonstrate a data-driven methodology for validation of advanced computer models used in nuclear power plant safety analysis. Specifically, the advanced computer models are those in the toolkit developed to support risk-informed safety margin characterization (RISMC), an integrated deterministic/probabilistic safety analysis methodology developed in the Department of Energy’s Light Water Reactor Sustainability (LWR-S) program. The report reflects the progress made towards the project’s stated goal by contributions by researchers and graduate students from universities, with support from researchers from national laboratories and industry companies. The project organization, effort coordination and technical implementation are summarized, followed by discussion of main findings, issues, and path forward. Selected chapters provide a more detailed description of tasks, approaches and respective findings and recommendations. Noteworthy are contributions that serve as guidelines for methodology development. It is also noted that this report is complemented by other milestone reports (as stand-alone deliverables) that provide detailed discussion of the technical developments. The project results have been documented in a number (12) dissertations and these, 50+ peer-reviewed publications in technical journals and conference proceedings.

42 ENGINEERING↗

ESnet Requirements Review Program Through the IRI Lens: A Meta-Analysis of Workflow Patterns Across DOE Office of Science Programs (Final Report)

The Department of Energy (DOE) ensures America’s security and prosperity by addressing its energy, environmental, and nuclear challenges through transformative science and technology solutions. The DOE’s Office of Science (SC) delivers groundbreaking scientific discoveries and major scientific tools that transform our understanding of nature and advance the energy, economic, and national security of the United States. The SC’s programs advance DOE mission science across a wide range of disciplines and have developed the research infrastructure needed to remain at the forefront of scientific discovery. The DOE SC’s world-class research infrastructure — exemplified by the 28 SC scientific user facilities — provides the research community with premier observational, experimental, computational, and network capabilities. Each user facility is designed to provide unique capabilities to advance core DOE mission science for its sponsor SC program and to stimulate a rich discovery and innovation ecosystem. Research communities gather and flourish around each user facility, bringing together diverse perspectives. A hallmark of many facilities is the large population of students, postdoctoral researchers, and early-career scientists who contribute as full-fledged users. These facility staff and users collaborate over years to devise new approaches to utilizing the user facility’s core capabilities. The history of the SC user facilities has many examples of wildly inventive researchers challenging operational orthodoxy to pioneer new vistas of discovery; for example, the use of the synchrotron X-ray light sources for study of proteins and other large biological molecules. This continual reinvention of the practice of science — as users and staff forge novel approaches expressed in research workflows — unlocks new discoveries and propels scientific progress. Within this research ecosystem, the high-performance computing (HPC) and networking user facilities stewarded by SC’s Advanced Scientific Computing Research (ASCR) program play a dynamic cross-cutting role, enabling complex workflows demanding high performance data, networking, and computing solutions. The DOE SC’s three HPC user facilities and the Energy Sciences Network (ESnet) high-performance research network serve all of the SC’s programs as well as the global research community. Argonne Leadership Computing Facility (ALCF), the National Energy Research Scientific Computing Center (NERSC), and Oak Ridge Leadership Computing Facility (OLCF) conceive, build, and provide access to a range of supercomputing, advanced computing, and large-scale data-infrastructure platforms, while ESnet interconnects DOE SC research infrastructure and enables seamless exchange of scientific data. All four facilities operate testbeds to expand the frontiers of computing and networking research. Together, the ASCR facilities enterprise seeks to understand and meet the needs and requirements across SC and DOE domain science programs and priority efforts, highlighted by the formal requirements reviews (RRs) methodology. In recent years, the research communities around the SC user facilities have begun experimenting with and demanding solutions integrated with HPC and data infrastructure. This rise of integrated-science approaches is documented in many community and high-level government reports. At the dawn of the era of exascale science and the acceleration of artificial intelligence (AI) innovation, there is a broad need for integrated computational, data, and networking solutions. In response to these drivers, DOE has developed a vision for an Integrated Research Infrastructure (IRI): To empower researchers to meld DOE’s world-class research tools, infrastructure, and user facilities seamlessly and securely in novel ways to radically accelerate discovery and innovation.

42 ENGINEERING↗

Comparing event generator predictions and ab-initio calculations of $\nu$-$^{12}$C neutral current quasi-elastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event Generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs use to predict kinematic observables for a range of neutrino energies. These simulations may lack physics captured by more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the physics contained in EG simulations. Neutral current quasi-elastic (NCQE) scattering events for neutrinos and anti-neutrinos on a $^{12}$C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the Relativistic Fermi Gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Phasor-Measurement-Unit-Based Data Analytics Using Digital Twin and PhasorAnalytics Software

A major objective of this project was to apply GE’s commercial machine learning and data analytics toolsets to large-scale, real-world, anonymized Phasor Measurement Unit (PMU) datasets in order to extract signatures, correlated and/or causal factors, and precursor patterns associated with significant power system phenomena. The project had a particular emphasis on extraction of insights relevant to asset health monitoring, real-time load modeling and cybersecurity monitoring. Additionally, the team was directed to undertake a comprehensive data quality analysis for the provided datasets and encouraged to estimate the ‘machine-learning readiness’ of the datasets by documenting any major obstacles to the application of commercial machine learning algorithms. To accomplish the aforementioned objectives, the project team’s work centered around the identification of key event signatures and application of the identified event signatures for event detection and event classification. The industry-validated, semi-supervised machine learning strategy employed for event signature identification involved several major tasks, including data-preprocessing, generation of an overabundance of features, normal data identification, normality modeling, and event signature identification through a methodical, quantitative ranking of features in order of relevance to each studied event type. Throughout the project, data quality issues and mitigation techniques were investigated. In this report, insights are provided regarding the readiness of the provided synchrophasor datasets for application of machine learning and data analytics. The methodologies employed for this technical strategy are summarized in this report. With regards to data preprocessing and feature generation, the provided Training and Test Datasets were ingested into GE’s big data environment. Subsequently, the team applied bad data cleansing and data imputation scripts, event detection scripts, and application programming interfaces (APIs) to the datasets for convenient data access. The project team completed development and validation of dozens of physics-based, statistics-based and transformation-based feature functions used for the extraction of over 60 synchrophasor features. Using a new parallel feature generation technology developed on this project, over 60 features have been rapidly generated for the full two years’ worth of Training and Test Dataset data associated with both the Eastern and Western interconnects. Even accommodating for temporal down-sampling inherent to the feature extraction procedure, this parallel feature generation activity resulted in a massive feature set with a storage requirement approximately equal to that of the raw training dataset itself. With regards to normal data identification and normality modeling, a normality model was built using the feature data extracted from the Training Dataset and iteratively refined subsequent to incremental adjustments and expansions of the Training Dataset feature data. With respect to event characterization and signature identification, an event signature identification pipeline was developed and used in conjunction with the normality model to identify over 15 event signatures for key event categories within the Training Dataset. The identified event signatures were used to characterize hundreds of key events in terms of relative severity, duration, and location of the event. An investigation was undertaken to identify correlated and causal factors involved in transformer events. A separate investigation into temporal trends in ring-down analysis results was undertaken to determine possible associations between system dynamics and various other factors such as loading, season or year. To validate the identified event signatures, additional work was undertaken to develop signature-based anomaly detection and classification tools suitable for convenient application to the synchrophasor datasets. The anomaly detection and classification tools, suitable for online application, were then applied to the entirety of the Eastern Interconnect Training and Test Datasets. Performance of the event detection and classification tools was evaluated upon receipt of the Test Dataset event logs (i.e., the labels for events contained in the Test Dataset), and promising results were obtained despite several challenges (documented herein) associated with application of supervised or semi-supervised machine learning methods to large-scale, anonymized datasets. Finally, the detection and classification tools were used to detect, classify, and characterize thousands of new events not included in the original event logs provided by the DOE within both the Training and Test Datasets.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An EvoDevo Study of Salmonid Visual Opsin Dynamics and Photopigment Spectral Sensitivity

Salmonids are ideal models as many species follow a distinct developmental program from demersal eggs and a large yolk sac to hatching at an advanced developmental stage. Further, these economically important teleosts inhabit both marine- and freshwaters and experience diverse light environments during their life histories. At a genome level, salmonids have undergone a salmonid-specific fourth whole genome duplication event (Ss4R) compared to other teleosts that are already more genetically diverse compared to many non-teleost vertebrates. Thus, salmonids display phenotypically plastic visual systems that appear to be closely related to their anadromous migration patterns. This is most likely due to a complex interplay between their larger, more gene-rich genomes and broad spectrally enriched habitats; however, the molecular basis and functional consequences for such diversity is not fully understood. This study used advances in genome sequencing to identify the repertoire and genome organization of visual opsin genes (those primarily expressed in retinal photoreceptors) from six different salmonids [Atlantic salmon ( Salmo salar ), brown trout ( Salmo trutta ), Chinook salmon ( Oncorhynchus tshawytcha ), coho salmon ( Oncorhynchus kisutch ), rainbow trout ( Oncorhynchus mykiss ), and sockeye salmon ( Oncorhynchus nerka )] compared to the northern pike ( Esox lucius ), a closely related non-salmonid species. Results identified multiple orthologues for all five visual opsin classes, except for presence of a single short-wavelength-sensitive-2 opsin gene. Several visual opsin genes were not retained after the Ss4R duplication event, which is consistent with the concept of salmonid rediploidization. Developmentally, transcriptomic analyzes of Atlantic salmon revealed differential expression within each opsin class, with two of the long-wavelength-sensitive opsins not being expressed before first feeding. Also, early opsin expression in the retina was located centrally, expanding dorsally and ventrally as eye development progressed, with rod opsin being the dominant visual opsin post-hatching. Modeling by spectral tuning analysis and atomistic molecular simulation, predicted the greatest variation in the spectral peak of absorbance to be within the Rh2 class, with a ∼40 nm difference in λ max values between the four medium-wavelength-sensitive photopigments. Overall, it appears that opsin duplication and expression, and their respective spectral tuning profiles, evolved to maximize specialist color vision throughout an anadromous lifecycle, with some visual opsin genes being lost to tailor marine-based vision.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling and analysis of heat emissions from buildings to ambient air

Heat emissions from buildings is a significant source of anthropogenic heat influencing the urban microclimate; however, they are usually oversimplified in urban climate and microclimate modeling. This study developed a bottom-up physics-based approach to calculate heat emissions from buildings to the ambient air and implemented the approach in EnergyPlus. In this work, a simple result verification was conducted by comparing the EnergyPlus simulated results against the spreadsheet calculations. Simulations covering 16 commercial building types, four climates, and two energy efficiency levels were conducted to understand and evaluate the building heat emissions and their temporal patterns as well as three major components: (1) building envelope (convective heat transfer to ambient air), (2) zones (air exfiltration and exhaust air), and (3) HVAC systems (relief air and heat rejection from condensers or cooling towers). The main findings are: (1) heat emissions are usually higher than the site energy use (about 2.5 times), and their dynamics should be considered; (2) building characteristics and their energy systems lead to differences in heat emission contributions from the three components, and their dynamics, for example, in the warehouse models, the envelope component accounts for 90.4%, while it is 12.7% for the large office models; (3) for most building typologies, the climate has a strong impact on heat emissions, for example, buildings with dominant heat emissions from the zone exhaust air and/or the HVAC reject heat, a general decrease in heat emissions in hotter climates is observed, while envelope-dominated buildings show the opposite; and (4) building technologies that reduce energy use in buildings may perform differently in reducing heat emissions. The developed heat emissions calculation method can be adopted in EnergyPlus and most other building energy modeling programs. It can provide dynamic building heat emissions as an input to urban climate computational fluid dynamics (CFD) models at a higher spatial and temporal resolution than is currently available, to improve the simulation accuracy of the urban microclimate and capture the urban heat island effect and urban overheating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Legacy mercury: A watershed scale cleanup challenge - 15559

Over the next three decades, the US Department of Energy (DOE) Oak Ridge Office of Environmental Management (OREM) will focus on reducing the flux of mercury contamination that exits at the Y-12 National Security Complex (Y-12 NSC) [1, 2]. OREM, Tennessee Department of Environment and Conservation (TDEC), and US Environmental Protection Agency (EPA) have identified mercury contamination within and released from the Y-12 NSC as the highest priority cleanup risk on the Oak Ridge Reservation (ORR). Furthermore, mercury contamination is a worldwide issue with more than 3,000 mercury-contaminated sites currently identified [3], and the United Nations Environment Programme has recently highlighted the risk of this contamination to human health [4]. Unlike other contaminants-metals, radionuclides, and organic solvents-that impact the DOE Environmental Management (DOE-EM) cleanup program at the ORR and other DOE sites, mercury has several unique characteristics that make environmental remediation of the Y-12 NSC one of the most complex efforts ever encountered. These distinctive physicochemical properties for mercury include the following: it is a dense non-aqueous liquid at ambient-temperature-pressure; it is the only metal that bio-magnifies; and it is the only contaminant transported as both a cation and/or a dissolved or gaseous elemental metal (similar to an organic solvent) under environmental conditions. Most importantly, mercury also undergoes biogeochemical transformation processes, which include aqueous and surface complexation, redox reactions, and atypical methylation reactions, producing the potent neurotoxin methylmercury [5, 6]. Each of the aforementioned distinctive characteristics for mercury is observed within this heterogeneous and dynamic watershed-scale system-which includes the Y-12 NSC boundary and 23 km of contaminated creek and flood plain downstream. Because of these complexities, implementing cost effective and sustainable solutions that reduce mercury flux from various primary and secondary contamination sources will require the linking of basic science understanding and applied research advancements into OREM's cleanup process. In 2009, the DOE-EM Office of Site Restoration's Soil and Groundwater Remediation Program and Oak Ridge National Laboratory (ORNL) initiated an applied field research initiative with the goal of developing a systems-based approach to control mercury flux in soil and water to protect surface water, groundwater, and ecological receptors. The Remediation of Mercury and Industrial Contaminants Applied Field Research Initiative (ROMIC-AFRI) represents the vehicle that enables the linking of science to remedial applications. Here we discuss the progress made to date in addressing key applied research and technology development needs to address this watershed scale conundrum. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Measurement of the Neutron Electromagnetic Form Factor Ratio from a Rosenbluth Technique with Simultaneous Detection of Neutrons and Protons

The internal structure of protons and neutrons provides insight into both the dynamical behavior of the constitute quarks and gluons, and emergent properties of the nucleons (such as mass, spin, and electromagnetic distributions). Elastic electron-nucleon scattering can probe the elastic electromagnetic form factors of the nucleon. The electric and magnetic form factors, respectively, encode information about the internal charge and magnetization distributions within the nucleon. Precision data for these form factors, over a broad range of the four-momentum transfer squared, Q^2, can benchmark theoretical models describing the strong interaction of nuclear physics. The Super BigBite Spectrometer (SBS) program in Hall A at Jefferson Lab, is a series of high-precision experiments which seek to significantly extend the Q^2 reach of previous data for the nucleon electromagnetic form factors. The first two experiments of this program are known as G_M^n and the neutron Two Photon Exchange (nTPE) and the data were collected from October 2021 to February 2022. Both experiments were conducted with the simultaneous measurement of D(e,e'n) and D(e,e'p) reactions for quasi-elastic electron-deuteron scattering. The scattered electrons were detected in the BigBite Spectrometer, which features multiple large-acceptance Gas Electron Multiplier (GEM) detectors. The Super BigBite Spectrometer provided simultaneous detection of scattered nucleons, and utilized a large acceptance dipole magnet and Hadron Calorimeter (HCal). The G_M^n experiment provides precision measurements of the neutron magnetic form factor, via the ratio method, over a Q^2 range of 3.0 to 13.5 (GeV/c)2. From this data analysis, preliminary values for G_M^n/µ_n G_D are extracted. For Q^2=4.48 (GeV/c)2 we find G_M^n/µ_n G_D=0.9546±0.0132 and for Q^2=4.476 (GeV/c)2 we find G_M^n/µ_n G_D=0.9563±0.0110. These preliminary G_M^n/µ_n G_D values are more precise than existing world data in this Q^2 regime and are consistent with the most recent parameterization of the G_M^n/µ_n G_D world data. The nTPE experiment provides a first measurement of the neutron Rosenbluth Slope and seeks to quantify the two-photon exchange(TPE) contribution to elastic electron-neutron scattering at a fixed Q^2=4.5 (GeV/c)2 with two different beam energies and scattering angle values. For data of the proton form factor ratio, µ_p G_E^p/G_M^p, significant discrepancies exist between values obtained from Rosenbluth Separation and polarization transfer measurement, particularly at large Q^2, and TPE contributions are thought to resolve this discrepancy. The impacts of TPE contributions have not yet been experimentally established for the neutron. From the data analysis presented in this dissertation, a preliminary result for the neutron Rosenbluth Slope is found as S^n=(G_E^n )^2/t_n (G_M^n )^2=0.0916±0.0476 for Q^2=4.48 (GeV/c)2. This value of the neutron Rosenbluth Slope is consistent with the world data extrapolation and the absence of large TPE corrections.

Wertz, Ezekiel [Thomas Jefferson National Accelera↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

Comparison of two software systems for quantification of myocardial blood flow in patients with hypertrophic cardiomyopathy

Backgorund: Quantification of myocardial blood flow (MBF) by positron emission tomography (PET) is important for investigation of angina in hypertrophic cardiomyopathy (HCM). Several software programs exist for MBF quantification, but they have been mostly evaluated in patients (with normal cardiac geometry), referred for evaluation of coronary artery disease (CAD). Software performance has not been evaluated in HCM patients who frequently have hyperdynamic LV function, LV outflow tract (LVOT) obstruction, small LV cavity size, and variation in the degree/location of LV hypertrophy. Aim: We compared results of MBF obtained using PMod, which permits manual segmentation, to those obtained by FDA-approved QPET software which has an automated segmentation algorithm. Methods: {sup 13}N-ammonia PET perfusion data were acquired in list mode at rest and during pharmacologic vasodilation, in 76 HCM patients and 10 non-HCM patients referred for evaluation of CAD (CAD group.) Data were resampled to create static, ECG-gated and 36-frame-dynamic images. Myocardial flow reserve (MFR) and MBF (in ml/min/g) were calculated using QPET and PMod softwares. Results: All HCM patients had asymmetric septal hypertrophy, and 50% had evidence of LVOT obstruction, whereas non-HCM patients (CAD group) had normal wall thickness and ejection fraction. PMod yielded significantly higher values for global and regional stress-MBF and MFR than for QPET in HCM. Reasonably fair correlation was observed for global rest-MBF, stress-MBF, and MFR using these two softwares (rest-MBF: r = 0.78; stress-MBF: r = 0.66.; MFR: r = 0.7) in HCM patients. Agreement between global MBF and MFR values improved when HCM patients with high spillover fractions (> 0.65) were excluded from the analysis (rest-MBF: r = 0.84; stress-MBF: r = 0.72; MFR: r = 0.8.) Regionally, the highest agreement between PMod and QPET was observed in the LAD territory (rest-MBF: r = 0.82, Stress-MBF: r = 0.68) where spillover fraction was the lowest. Unlike HCM patients, the non-HCM patients (CAD group) demonstrated excellent agreement in MBF/MFR values, obtained by the two softwares, when patients with high spillover fractions were excluded (rest-MBF: r = 0.95; stress-MBF: r = 0.92; MFR: r = 0.95). Conclusions: Anatomic characteristics specific to HCM hearts contribute to lower correlations between MBF/MFR values obtained by PMod and QPET, compared with non-HCM patients. These differences indicate that PMod and QPET cannot be used interchangeably for MBF/MFR analyses in HCM patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗