Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model reference control”

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 433 records · Page 24

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗

Verification of RESRAD-OFFSITE Code (V.4)

This report documents the verification of RESRAD-OFFSITE Version 4.0 and describes, where necessary, the verification of the following: • The data comprising the standard dose and risk coefficient libraries in the RESRAD database files Master_dcf_ICRP07.mdb and Master_dcf_2k.mdb. • The extraction and transfer of the data from the selected database file to the computational code by the RESRAD-OFFSITE 4.0 interface, ResOWin.exe. • The different processes that are modeled by the main computational code in RESRAD OFFSITE 4.0, ResOMain.exe. • The data displayed in the graphical and text reports. Many verifications were performed as part of the quality assurance quality control program associated with the development and release of RESRAD-OFFSITE 4.0, namely: • developer testing, • internal independent testing, and • release testing. Some were also performed in response to questions from users regarding the performance of the code. The main text of the report focuses on summarizing a subset of those tests, both independent and developer tests that verified the computations performed by the code. The verifications included in this report served as the basis for the development of the release tests of the computational executables and provided the quantitative results to be compared with the code output. The input and output interfaces and the data transfers between the various executables of the code were tested while performing the verification testing. They were tested intentionally during release testing. This report also provides some basic information to help in understanding the activities that were verified. The report: • outlines the components of RESRAD-OFFSITE 4.0 and the interconnections between these components, • outlines the processes modeled by the computational code, • provides summary figures and tables to offer confirmation of the verification of the computational components of the code, • reproduces the verifiers’ reports, if available, in individual appendices, • refers to the previous verification report (Yu et al. 2011) for more details about some of the verifications, and • reproduces the test cases and the testers’ reports from the release testing in individual appendices, when possible.

54 ENVIRONMENTAL SCIENCES↗

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗

Cybersecurity for Distance Relay Protection

This project is a DOE follow-up effort on the CREDC workshop held on September 13, 2018 in Cambridge, MA to discuss cybersecurity of distance relays, which considered the benefits, vulnerabilities and risk mitigations for the use of communication systems in power system protection. The objectives of this project are to define the taxonomy of relay protection and associated communications; define use cases describing approaches to reduce the cyber-attack surface on those protective relays; and evaluate the loss of operational functional capability from changes to communication coverage. Mitigating controls will also be evaluated to understand if there are other approaches to reduce attack surfaces while maintaining communications or partial communications. Distance relays are used to protect transmission lines of approximately 10 to 300 miles in length, by detecting short circuits (i.e., faults) on the lines and then tripping circuit breakers in the substation. Such protection systems are a subset of the power system and they incorporate sensing, logic and communication functions. Protection system exposure to cyberattack could be drastically limited by disconnecting relays from all vulnerable communication systems, but this may adversely impact overall power system performance in the absence of cyberattack. This project began with a use case analysis of protection systems with communications, as summarized in this report. It continued with modeling, testing and evaluation in a miniature power system (MPS), located in the Western Area Power Administration (WAPA) Electric Power Training Center (EPTC). The project also incorporated feedback from two industry meetings held in February and September 2019. The suggested next steps account for and complement the work already underway with DOE/CESER funding: 1. Study the performance of LCD and PC vs. PUTT, which is less reliant on communication system performance and GPS timing references. The PUTT scheme could prove to be more resilient to cyberattack or communications-related disruption. It could also be more tolerant of message re-routing with SDN/SDR communication systems. On the other hand, it will be more vulnerable to false tripping during dynamic events or to loss of the voltage signal. The optimum choice of scheme may depend on the specific power system and risk assessment. This study could provide a new template for evaluation based on business functions. 2. Research and develop new methods to detect and monitor distributed physical attacks, possibly using drones, video sensors, thermal sensors, machine learning and other advanced techniques. This will help mitigate the impact of cyberattack on the protection system, and will also help mitigate the impact of wild fires. 3. Implement a scalable PKI for use in electric utility protection systems. This will encourage widespread adoption of secure authentication methods that are already available, but not widely used at present. This will help secure engineering access to the relays. 4. Investigate the use of SDN in combination with SDR to achieve better cybersecurity and electromagnetic security of the network, incorporating path variability. This would help secure both engineering access and peer-to-peer GOOSE messaging. 5. Perform additional testing, with operator evaluation of “red button” scenarios, PUTT vs. LCD, relay mis-operations, and other cyberattacks in the EPTC. This is an important advantage of testing in the EPTC rather than by computer simulation or even hardware-in-the-loop simulation; the EPTC is already dedicated to managing the situational awareness, operator response times and other human impacts. One of the project objectives was to settle on a common nomenclature for this problem space. We have concluded that the OSI layer model, supplemented by ANSI device numbers and other IEEE standards, is already well-accepted by the industry. The IEEE PSRC knowledge base provides a great deal of public information

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimal Heliostat Assignment Strategy for Multiple-Receiver Systems

Next-generation concentrating solar power concepts currently being developed under the U.S. Department of Energy's Gen3 program include several novel heliostat field and receiver arrangements. Among these concepts are solar fields utilizing multiple receivers or several distinct receiver flux control zones, and each of which has particular power or peak flux constraints. In this paper, we present a linear program that identifies optimal heliostat-to-receiver assignment strategies for multi-receiver or multi-zone solar fields based on detailed calculations made by an optics modeling package. Two distinct problems are addressed: namely, (1) a design problem that chooses the optimal set of heliostats to include in a final layout given receiver constraints and performance, and (2) an operations problem that assigns existing heliostats to a receiver, maximizing field power output while maintaining receiver power constraints. The optimization methodology is integrated into the concentrating solar power tower modeling package SolarPILOTTM, and results are presented showing that heliostat assignment strategies can impact both initial heliostat field layout characteristics and the cost of thermal energy produced by the field. The methodology succeeds in maintaining power requirements but at the expense of producing different flux profiles on each receiver surface. Consequently, the paper also discusses additional measures that are applied to ensure that desirable flux distribution properties are maintained, reducing local flux mean-absolute-deviation values from a baseline of 17% or greater to less than 5% compared to a desired reference flux profile.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Genomic-based biosurveillance for avian influenza: whole genome sequencing from wild mallards sampled during autumn migration in 2022–2023 reveals a high co-infection rate on migration stopover site in Georgia

The Caucasus region, including Georgia, is an important intersection for migratory waterbirds, offering potential for avian influenza virus (AIV) transmission between populations from different geographic areas. In 2022 and 2023, wild ducks were sampled during autumn migration events in Georgia to study the genetic relationships and molecular characteristics of influenza strains. Sequencing and phylogenetic analysis were used to compare the sampled strains to reference sequences from Africa, Asia, and Europe, allowing assessment of genetic relationships and virus transmission between migratory birds. Protein language modeling identified potential co-infections. Of 225 duck samples, 128 tested positive for the influenza M gene. 55 influenza-positive samples underwent whole-genome sequencing, revealing significant diversity. Analysis of the hemagglutinin (HA) segment showed notable differences among subtypes. Most samples were H6N1 and H6N6, but co-infections with combinations like H6H3, N8N1, N6H9, N2N6, and H9H6/N1N2 were also identified. These findings demonstrate the high variability of influenza viruses in migratory waterbirds in Georgia, including a notable rate of co-infections. Some samples exhibited uncommon genetic characteristics compared to other strains from the same year, suggesting Georgia’s role as a mixing vessel for influenza viruses. This facilitates reassortment during co-infections and contributes to the genetic diversity observed across flyways.

59 BASIC BIOLOGICAL SCIENCES↗

Investigation of Active Life Balancing to Recondition Li-ion Battery Packs for 2nd Life

An investigation is presented into the application of active life balancing battery management systems to recondition inhomogeneous Li-ion battery packs and achieve homogeneity in capacity for second life. This paper proposes a life balancing objective map, which controls the state of charge references for each cell in a battery pack based on its capacity relative to other cells. The impacts of life balancing objective maps on the reconditioning time and overall pack capacity fade are considered. Moreover, the proposed approach accounts for various load and environmental effects, including depth-of-discharge, battery pack temperature, number of cycles per day, and cells’ resting time at various states of charge. The study results in a battery pack with convergent capacities achieved in minimal time, while curtailing overall capacity fade. Simulation results for a battery pack with sixteen series connected 75 Ah Kokam NMC cells with 3.6% initial capacity imbalance are presented, indicating a reconditioning time of 1.34 months with a final capacity imbalance of 0.1% and an overall capacity fade of 0.005%.

47 OTHER INSTRUMENTATION↗

Moisture Performance of Unvented Attics With Vapor Diffusion Ports and Buried Ducts in Hot, Humid Climates

Energy efficiency measures, such as cool roofs, radiant barriers, interior radiative control coatings, and buried ducts are increasing in popularity and are promoted by energy codes because of their energy-saving potential. However, these strategies can also pose moisture risks in attics by lowering surface temperatures and increasing condensation potential and moisture accumulation. Of particular concern in hot-humid climates is dripping condensation on cold air-conditioning ductwork in the summer - commonly referred to as duct "sweating" - which threatens the attic floor with conditions conducive to mold growth and rot. One strategy to mitigate these moisture issues is to wrap ductwork in thicker duct-wrap insulation with an integrated exterior vapor barrier, but thick duct wrap can be difficult to come by, expensive, and unwieldy to work with. This study explores an alternative strategy of reducing moisture issues while embracing energy efficiency by using unvented attics with vapor diffusion ports and buried ductwork in hot-humid climates. Vapor diffusion ports have been studied so far in a wide range of U.S. climates, mostly in the context of conditioned attics. In this study, the strategy is implemented in the novel context of hot-humid climates with ductwork sitting atop blown-in attic floor insulation in unconditioned attics. Using a combination of field experiments and hygrothermal modeling, the findings of this project indicate that an unvented attic with vapor diffusion ports and buried ducts may be a key part of a successful low-cost method for reducing the attic moisture load by venting excess moisture out of the attic, keeping duct-jacket surfaces above dew point temperature, and keeping the roof deck safe from winter moisture accumulation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Transition Core Modeling for Extended Enrichment & Accident-Tolerant Fuels Using Polaris/PARCS

Commercial light water reactor (LWR) operators and fuel vendors are currently interested in increasing the low-enriched uranium (LEU) fuel enrichments from the current limit of 5.0 $^w/_o$ $^{235}U$ up to 10 $^ w/_o$ $^{235}U$ (referred to as "LEU+") in their current fleets; they are also interested in using accident-tolerant fuel (ATF) with both LEU and LEU+ fuel. This report aims to identify modeling challenges and accuracy concerns in transition core analysis using the SCALE Polaris lattice physics code and U.S. Nuclear Regulatory Commission core simulator PARCS. At the time this study was started, no publicly available LEU+ core designs existed for boiling water reactor (BWR) or pressurized water reactor (PWR) systems. Therefore, fuel lattices were shuffled within a multi-assembly model to mimic neutronically challenging lattice combinations seen in transition cores, such as a fresh LEU+ lattice next to depleted LEU lattices. In addition to multi-assembly models, whole-core BWR transition core calculations were performed for ATF and LEU+ fuel using an existing Hatch-1 Cycle 3 core model. A whole-core BWR model was chosen due to the more heterogeneous core designs compared to those for a PWR core. Since the original core is an old checkerboard core design and no core or fuel design optimization was performed for the modeled fuel types, these core calculations were intended only to provide: (1) Comparisons of core characteristics of interest, such as the pin power distributions and peaking factors, Doppler temperature coefficients (DTCs), and control blade worths (CBWs) under challenging core designs, (2) Identification of reactor physics challenges in modeling LEU+ and ATF cores, and (3) A stress test for the Polaris/PARCS two-step modeling approach, including characterization of the relative accuracy for predicting characteristics of interest such as pin power distributions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Adaptive anomaly detection for identifying attacks in cyber-physical systems: A systematic literature review

Modern cyberattacks in cyber-physical systems (CPS) rapidly evolve and cannot be deterred effectively with most current methods, which focus on characterizing past threats. Adaptive anomaly detection (AAD) is among the most promising techniques to detect evolving cyberattacks, with an emphasis on fast data processing and model adaptation. AAD has been researched extensively; however, to the best of our knowledge, our work is the first systematic literature review (SLR) on current research in this field. We present a comprehensive SLR, gathering 397 relevant papers and systematically analyzing 65 of them (47 research and 18 survey papers) on AAD in CPS from 2013 to November 2023. We introduce a novel taxonomy considering attack types, CPS application, learning paradigm, data management, and algorithms. Our findings show that most studies addressed either model adaptation or data processing, but rarely both simultaneously. This indicates a research gap in fully adaptive solutions. We also categorize algorithms, datasets, and attack characteristics, and summarize strengths and weaknesses across the literature. Our review provides a structured and accessible reference for researchers and practitioners, offering insights into key trends and highlighting limitations in current approaches. Finally, we outline several future research directions, including the need for integrated real-time processing and adaptive learning, explainability, and uncertainty quantification in AAD for CPS.

Adaptation↗

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗

MFANS 2024 - Formally Proving Characteristics of Cyber-Physical Systems

Cyber-physical systems (CPS) are engineered systems that rely on the smooth integration of computational algorithms and physical elements. This integration presents new challenges for verifying that systems will behave as expected. The goal of this presentation is to present current challenges and potential solutions for the formal verification of cyber-physical systems. For cyber systems, formal methods refer to systematically rigorous mathematical techniques employed in the specification, development, analysis, and verification of both software and hardware systems. Recent advancements in computer science have yielded sophisticated tools specifically designed to address challenges associated with formal methods in complex systems. These tools leverage various foundational concepts such as logic, formal languages, program semantics, type systems, type theory, and automata theory. A notable achievement in the application of formal methods is the seL4 microkernel, claimed to be the first general-purpose operating-system kernel to be verified. Its proof implies the absence of bugs and guarantees that the kernel meets specifications. For physical systems, dynamic and control theory has a history of using rigorous analytic techniques to prove functional correctness. Lyapunov, optimal, classical, modern, and robust control theories all provide rigorous mathematical methods both to analyze system performance and to design controller that can be guaranteed to meet certain objectives. Recent computational techniques like level set theory and reachability analysis provide assertions that a system's state will avoid unsafe regions. Even though success has been independently achieved for cyber systems and physical systems, the integration of such systems creates new challenges. In particular, there is an obvious discrepancy between finite-state machines and infinite-state systems, resulting in different approaches for modeling and analyzing these system. While it is possible to simulate hybrid systems, this provides only a demonstration of a performance and not proof. For hybrid systems, current formal methods and system analysis approaches typically require a workarounds to work on hybrid systems like CPS. This paper will outline the state of the art and limits of current practice for formally verifying CPS and will identify possible research directions that require attention.

97 MATHEMATICS AND COMPUTING↗

Performance of a CO 2 sorbent for indoor air cleaning applications: Effects of environmental conditions, sorbent aging, and adsorption of co-occurring formaldehyde

Indoor air cleaning systems that incorporate CO 2 sorbent materials enable HVAC load shifting and efficiency improvements. This study developed a bench-scale experimental system to evaluate the performance of a sorbent under controlled operation conditions. A thermostatic holder containing 3.15 g sorbent was connected to a manifold that delivered CO 2 -enriched air at a known temperature and relative humidity (RH). The air stream was also enriched with 0.8-2.1 ppm formaldehyde. The CO 2 concentration was monitored in real-time upstream and downstream of the sorbent, and integrated formaldehyde samples were collected at different times using DNPH-coated silica cartridges. Sorbent regeneration was carried out by circulating clean air in countercurrent. Almost 200 loading/regeneration cycles were performed in the span of 17 months, from which 104 were carried out at reference test conditions defined by loading with air at 25°C, 38% RH, and 1000 ppm CO 2 , and regenerating with air at 80°C, 3% RH and 400 ppm CO 2 . The working capacity decreased slightly from 43-44 mg CO 2 per g sorbent to 39-40 mg per g over the 17 months. The capacity increased with lower loading temperature (in the range 15-35°C) and higher regeneration temperature, between 40 and 80°C. The CO 2 capacity was not sensitive to the moisture content in the range 6-9 g/m 3 , and decreased slightly when dry air was used. Loading isothermal breakthrough curves were fitted to three simple adsorption models, verifying that pseudo-first-order kinetics appropriately describes the adsorption process. The model predicted that equilibrium capacities decreased with increasing temperature from 15 to 35°C, while adsorption rate constants slightly increased. The formaldehyde adsorption efficiency was 80%-99% in different cycles, corresponding to an average capacity of 86 ± 36 µg/g. Formaldehyde was not quantitatively released during regeneration, but its accumulation on the sorbent did not affect CO 2 adsorption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Assessment of Core Physics Characteristics of Extended Enrichment and Higher Burnup LWR Fuels using the Polaris/PARCS Two-Step Approach. Vol. I: PWR Fuel

Nuclear fuel with extended enrichment (235U enrichment within 5-8 wt%) is one of the evolutionary changes that have been pursued in recent years by commercial light water reactor operators and fuel vendors to improve the fuel cycle economy and operation performance of a nuclear plant. This work assesses the performance of the Polaris/PARCS two-step approach in core physics modeling of the pressurized water reactor cores with extended enrichment fuel, referred to as “LEU+” in this report. A representative LEU+ core with a 24-month fuel cycle developed by Southern Nuclear Company (SNC) was modeled using this two-step approach. A representative LEU core with an 18-month fuel cycle was also modeled to provide a reference for the LEU+ core. As expected, significantly more burnable poison absorbers were used in the LEU+ core to accomodate its higher fuel enrichment. Nine different fuel assembly types were modeled using Polaris for each core to generate the assembly cross sections, which were then processed by GenPMAX to prepare the cross-section data for PARCS. The average specific powers of each fuel batch in each core were derived from VERA results and higher specific powers in fresh assemblies were found in the LEU core due to its less total uranium loading included in the VERA LEU model, given that the total core power was assumed to be the same for both cores. PARCS models were developed to simulate the steady-state operations of both cores. PARCS results on the LEU+ core were first compared with the VERA results for verification purpose; good agreements were seen in soluble boron and burnup distribution results, indicating that the Polaris/PARCS modeling and simulation were correctly implemented. Core physics parameters calculated by PARCS, at zero power physics tests, beginning of cycle (BOC), and end of cycle conditions (EOC), were compared between the LEU+ core and the LEU core, including soluble boron concentration, burnup distributions, assembly and pin power peaking factors, fuel temperature reactivity coefficients, moderator temperature and density reactivity coefficients, control rod worth, and shut down margin. The main differences in PARCS results between the LEU+ and the LEU cores are summarized below: 1)The critical boron concentrations were found to be much higher in the LEU+ core than in the LEU core (1582 vs. 1335 ppm for peak values). 2)Higher assembly radial power peaking factors (1.4 vs. 1.3 for peak values), 2D pin peaking factors (1.53 vs. 1.42 for peak values), and 3D pin peaking factors (1.89 vs. 1.81 for peak values) were found in the LEU+ core than in the LEU core. 3)Significantly higher reactivity coefficients of moderator temperature (and density) were found in LEU+ than in LEU.4)Significantly lower control rod worth at EOC were found in LEU+ than LEU for all but one control banks.5)Significantly lower shut down margins were found in the LEU+ core than in the LEU core. The Polaris/GenPMAX/PARCS code suite was found to be capable of modeling the LEU+ PWR core for steady-state operations and no unexpected results in core physics parameters were observed, in spite of that a) several bugs in PARCS were identified and workarounds were used; b) several features were found lacking in the current version of PARCS that would be useful for core modeling. A list of requests for bug fixes and feature upgrades for PARCS originated from this work were transmitted to the code developers. The assessments on the performance of the Polaris/PARCS two-step approach in core modeling for boiling water reactor with LEU+ fuel is ongoing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A large-scale comparison of Artificial Intelligence and Data Mining (AI&DM) techniques in simulating reservoir releases over the Upper Colorado Region

In recent years, the Artificial Intelligence and Data Mining (AI&DM) models have become popular tools in assisting various aspects of reservoir operation. However, the practical uses are still rarely reported. Comparison experiment of many AI&DM models over a large number of reservoir cases is particularly valuable to help reservoir operators first examine the usefulness and transferability of different AI&DM models, and then identify the most stable and reliable AI&DM model in assist of various decision-making processes. In this study, a total of 12 AI&DM models with different parameterizations and simulation scenarios are comprehensively tested out and compared in simulating the controlled reservoir outflows of 33 reservoir cases over the Upper Colorado Region, United States. Results show that the Random Forecast and the Long-Short-Term-Memory model could consistently derive the best statistical performance than other models under the baseline simulation scenario. The employed AI&DM models could obtain satisfactory statistical interquartile ranges (25–75%) between [0.6–0.9], [0.3–0.8], and [0.2–0.8], for CORR, NSE, and KGE measurements, respectively, and [1.5–6.5], [–15 to 20], and [0.5–8.5] for the normalized RMSE, PBIAS and RSR measurements, respectively. Results also show Multi-Layer Perceptron model and Extreme Gradient Boosting Tree Algorithm produced more stable and superior performance than other models under more complex input scenarios. We also found that the performance of different AI&DM models are closely relevant to the reservoir elevations, sizes, and functionalities. Discussions were made about the sensitivity of AI&DM models’ parameterizations and the key advantages of AI&DM models over the rule-based reservoir models. We further identify that the main advantage of AI&DM models is the flexibility in designing input structures, whereas the rule-based simulation model is rather limited. Future studies were suggested regarding the best way reservoir operators and researchers could use, select, and apply different AI&DM models in simulating reservoir releases under different natural and modeling environments. Finally, this comparison study also serves as a reference and a piece of groundwork for further promoting the practical uses of AI&DM models in assisting reservoir operation.

54 ENVIRONMENTAL SCIENCES↗

Design and Optimization of Structured Multi-Functional Trapping Catalysts for Conversion of Hydrocarbons and NOx from Diesel and Advanced Combustion Engines

Oxides of nitrogen in the form of nitric oxide (NO) and nitrogen dioxide (NO 2 ) commonly referred to as NOx, is one of the two chemical precursors that lead to ground-level ozone, a ubiquitous air pollutant in urban areas. A major source of NOx is generated by equipment and vehicles powered by diesel engines, which have a combustion exhaust that contains NOx in the presence of excess O 2 . Vehicular emission control catalysts are ineffective in eliminating CO, hydrocarbons, and NOx during engine cold-start when exhaust temperatures are below 200°C. The objective of the project was to develop and demonstrate a multi-functional, catalyzed trap that enables vehicles with advanced combustion strategies to meet Tier 3 emissions standards while achieving the 150 °C challenge for sustained co-oxidation of HCs and CO and ≥90% NO trapping and release during warmup. Specifically, the multi-functional Lean HC+NOx (LHCNT) was developed for application in the exhaust aftertreatment of conventional diesel engines and engines having low temperature combustion (LTC) regimes. Activities included the design and synthesis of adsorbents and catalysts, screening and evaluation. Passive NOx absorbers (PNA), hydrocarbon (HC) traps, and oxidation catalysts (OC) were evaluated for use in series or as integrated devices. Predictive tools were developed utilizing the characterization and analysis of these materials, and an emission system was designed and optimized utilizing the catalyst systems. Microkinetic models were developed for the PNA for the simple NO-only feed and complex feed containing CO, H 2 , and model hydrocarbons (ethylene and dodecane). A first-principles, mechanistic-based model of the PNA was developed which utilized molecular-scale estimates (density functional theory) of energy barriers, mechanistic-based kinetics and realistic treatments of the flow and transport processes. Two new oxidation catalysts were developed (PdCu alloy, mixed copper-ceria-cobalt oxide), both of which significantly lessened the detrimental inhibition by CO on hydrocarbon and NO oxidation. A method for lessening the detrimental impact of CO on PNA activity was developed that involves use of an oxidation catalyst upstream of the PNA. The SwRI Ectolab TM burner system was applied to evaluate the baseline PNA material and confirmed performance comparable to the benchflow PNA studies using simulated exhaust. Spatially-resolved mass spectrometry (SpaciMS) was used to measure the transient spatial profiles of reacting species spanning the length of a three-function LHCNT containing PNA, HCT, and OC. The findings from this study provide diesel vehicle and catalyst companies valuable information to develop more cost effective emission control catalysts which helps to expand the use of more fuel efficient diesel power. The fundamental modeling and experimental tools and findings from this project can be applied to catalyst technologies used in the energy and chemical industries. The project led to 14 publications in the peer-reviewed literature with 2 additional currently under review. Finally, the project also led to training of several doctoral students who were placed in research jobs in industry and academia. Specifically, Mugdha Ambast (UH) has joined Cummins, Kevin Gu (UVa) has joined GM, and Abhay Gupta (UH) is to join Caterpillar.

02 PETROLEUM↗

Design and Optimization of Structured Multi-Functional Trapping Catalysts for Conversion of Hydrocarbons and NOx from Diesel and Advanced Combustion Engines

Oxides of nitrogen in the form of nitric oxide (NO) and nitrogen dioxide (NO 2 ) commonly referred to as NOx, is one of the two chemical precursors that lead to ground-level ozone, a ubiquitous air pollutant in urban areas. A major source of NOx is generated by equipment and vehicles powered by diesel engines, which have a combustion exhaust that contains NOx in the presence of excess O 2 . Vehicular emission control catalysts are ineffective in eliminating CO, hydrocarbons, and NO x during engine cold-start when exhaust temperatures are below 200°C. The objective of the project was to develop and demonstrate a multi-functional, catalyzed trap that enables vehicles with advanced combustion strategies to meet Tier 3 emissions standards while achieving the 150 °C challenge for sustained co-oxidation of HCs and CO and ≥90% NO trapping and release during warmup. Specifically, the multi-functional Lean HC+NOx (LHCNT) was developed for application in the exhaust aftertreatment of conventional diesel engines and engines having low temperature combustion (LTC) regimes. Activities included the design and synthesis of adsorbents and catalysts, screening and evaluation. Passive NOx absorbers (PNA), hydrocarbon (HC) traps, and oxidation catalysts (OC) were evaluated for use in series or as integrated devices. Predictive tools were developed utilizing the characterization and analysis of these materials, and an emission system was designed and optimized utilizing the catalyst systems. Microkinetic models were developed for the PNA for the simple NO-only feed and complex feed containing CO, H 2 , and model hydrocarbons (ethylene and dodecane). A first-principles, mechanistic-based model of the PNA was developed which utilized molecular-scale estimates (density functional theory) of energy barriers, mechanistic-based kinetics and realistic treatments of the flow and transport processes. Two new oxidation catalysts were developed (PdCu alloy, mixed copper-ceria-cobalt oxide), both of which significantly lessened the detrimental inhibition by CO on hydrocarbon and NO oxidation. A method for lessening the detrimental impact of CO on PNA activity was developed that involves use of an oxidation catalyst upstream of the PNA. The SwRI Ectolab TM burner system was applied to evaluate the baseline PNA material and confirmed performance comparable to the benchflow PNA studies using simulated exhaust. Spatially-resolved mass spectrometry (SpaciMS) was used to measure the transient spatial profiles of reacting species spanning the length of a three-function LHCNT containing PNA, HCT, and OC. The findings from this study provide diesel vehicle and catalyst companies valuable information to develop more cost effective emission control catalysts which helps to expand the use of more fuel efficient diesel power. The fundamental modeling and experimental tools and findings from this project can be applied to catalyst technologies used in the energy and chemical industries. The project led to 14 publications in the peer-reviewed literature with 2 additional currently under review. Finally, the project also led to training of several doctoral students who were placed in research jobs in industry and academia. Specifically, Mugdha Ambast (UH) has joined Cummins, Kevin Gu (UVa) has joined GM, and Abhay Gupta (UH) is to join Caterpillar.

42 ENGINEERING↗