Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fractional operators”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Search for the light hadron decay χ c 1 ( 3872 ) → π + π − η

With a data sample corresponding to an integrated luminosity of 11.5 fb − 1 collected with the BESIII detector operating at the BEPCII storage ring, for the first time the light hadron decay χ c 1 ( 3872 ) → π + π − η is searched for. While no significant signal is observed, the upper limits at the 90% confidence level for σ [ e + e − → γ χ c 1 ( 3872 ) ] B [ χ c 1 ( 3872 ) → π + π − η ] at center of mass energies from 4.13 GeV to 4.34 GeV are determined. By normalizing to the χ c 1 ( 3872 ) → π + π − J / ψ decay channel, a 90% confidence level upper limit for the branching fraction ratio R = B [ χ c 1 ( 3872 ) → π + π − η ] / B [ χ c 1 ( 3872 ) → π + π − J / ψ ] < 0.12 is given. These measurements provide important inputs for understanding the internal structure of the χ c 1 ( 3872 ) resonance. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

CO 2 to Bioplastics: Beneficial Re-use of Carbon Emissions from Coal-fired Power Plants using Microalgae

This project sought to address the technical and economic barriers to carbon dioxide (CO 2 ) capture and utilization using microalgae. Specifically, a dual photobioreactor (PBR)/open raceway pond (ORP) cultivation system was evaluated with respect to capital and operational costs, productivity, and culture health, and compared to an ORP-only cultivation system. In the dual (or “hybrid”) system, a cyclic flow PBR was used to provide inoculum for two open raceway ponds after each harvest, with the hypothesis that this PBR + pond system should result in improved performance compared to traditional raceway ponds. Two other raceway ponds were operated conventionally as a control, the experiments being performed at Duke Energy’s East Bend Station in northern KY using coal-derived flue gas as the CO 2 source. Although algae productivity achieved was mediocre due to the fact that the experiments had to be performed late in the growing season when climatic conditions were not optimal for high growth, the hybrid cultivation system showed a higher level of algae productivity (statistically significant) compared to the traditional ponds. In order to realize the maximum value of the algal biomass produced, fractionation of the biomass was examined. Bioplastic compounding and material characterization was subsequently performed using three feed stocks: whole biomass, lipid-extracted biomass and the proteinaceous residue from full fractionation. In general, the whole and lipid-extracted biomass gave similar results in terms of the properties of the resulting bioplastics, while the proteinaceous residue showed promise for the production of a polybutylene adipate terephthalate (PBAT) blend. Finally, sustainability assessment of a putative biorefinery system was conducted, including techno-economic and life cycle impact assessment. Nine different production scenarios were considered, comprising combinations of the three different biomass growth architectures (ORP, PBR and dual PBR-ORP systems) coupled with three different algae biomass processing pathways (drying only, lipid extraction and fractionation). Results show that the minimum selling price of the bioplastic feed stock (BPFS) is within the realm of economic competition with prices as low as $970 USD tonne -1 . Additionally, life cycle impact assessment results indicate drastic improvements in performance of the produced BPFS, with reductions in greenhouse gas emissions ranging between 67 and 116% compared to a petroleum based plastic feedstock.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Origin and characterization of single-electron events in Skipper-CCDs for light dark matter searches

This manuscript summarizes the work carried out within the SENSEI (Sub-Electron-Noise Skipper-CCD Experimental Instrument) collaboration, an effort focused on the search for light dark matter using Skipper-CCDs. Firstly, it presents the data acquisition and processing protocol developed to establish a selection criteria for events compatible with the dark matter signal. This protocol also enabled the structure and tools used for the study and characterization of phenomena that impact the quality of measurements made in the search for dark matter (measurement of the Fano factor, determination of the Compton background, and characterization of the high-energy surface background), which were carried out during this work. Secondly, it discusses the comprehensive study conducted on the origin of single-electron events in Skipper-CCDs, which allowed for the selection of detector operating parameters and improved the detection sensitivity in the search for dark matter. It is highlighted the characterization of the spatial correlation between single-electron events and high-energy events, among other developed event selection criteria, and the optimization of the output device to reduce sources of luminescence during its operation. As a result of this study, the lowest levels of dark current and spurious charge ever measured in a CCD are reported. Finally, the latest published results from SENSEI are presented. The obtained limits for the scattering of a light mediator are the best reported to date across the entire investigated mass range, while for the heavy mediator, they are the best for masses below 10 MeV. In terms of dark matter absorption, the most restringent limits below 10 eV are reported. These results, achieved with only a very small fraction of the total projected exposure for SENSEI, position Skipper-CCDs as leaders in the search for light dark matter worldwide.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Anyonic Membranes and Pontryagin Statistics

Anyons, unique to two spatial dimensions, underlie extraordinary phenomena such as the fractional quantum Hall effect, but their generalization to higher dimensions has remained elusive. The topology of Eilenberg-MacLane spaces constrains the loop statistics to be only bosonic or fermionic in any dimension. In this work, we introduce the novel anyonic statistics for membrane excitations in four dimensions. Analogous to the $\mathbb{Z}_N$-particle exhibiting $\mathbb{Z}_{N\times \gcd(2,N)}$ anyonic statistics in two dimensions, we show that the $\mathbb{Z}_N$-membrane possesses $\mathbb{Z}_{N\times \gcd(3,N)}$ anyonic statistics in four dimensions. Given unitary volume operators that create membrane excitations on the boundary, we propose an explicit 56-step unitary sequence that detects the membrane statistics. We further analyze the boundary theory of $(5{+}1)$D 1-form $\mathbb{Z}_N$ symmetry-protected topological phases and demonstrate that their domain walls realize all possible anyonic membrane statistics. We then show that the $\mathbb{Z}_3$ subgroup persists in all higher dimensions. In addition to the standard fermionic $\mathbb{Z}_2$ membrane statistics arising from Stiefel-Whitney classes, membranes also exhibit $\mathbb{Z}_3$ statistics associated with Pontryagin classes. We explicitly verify that the 56-step process detects the nontrivial $\mathbb{Z}_3$ statistics in 5, 6, and 7 spatial dimensions. Furthermore, in 7 and higher dimensions, the statistics of membrane excitations stabilize to $\mathbb{Z}_{2} \times \mathbb{Z}_{3}$, with the $\mathbb{Z}_3$ sector consistently captured by this process.

Abstract algebra↗
DOE OSTI2021

AGR-5/6/7 Irradiation Test Final As-run Report

This document presents the as-run analysis of the Advanced Gas Reactor (AGR) 5/6/7 irradiation experiment. AGR-5/6/7 is the last of a series of experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) in support of the development and qualification of tri-structural isotropic (TRISO) low enriched fuel for use in high-temperature gas-cooled reactors (HTGR). The test train contained five separate capsules that were independently controlled and monitored. Each capsule contained multiple 24.91-mm-long and 12.25-mm-dimeter compacts filled with low-enriched uranium carbide/oxide tri-structural isotropic fuel particles. The objectives of the AGR 5/6/7 experiment were to: 1. Irradiate reference-design fuel particles to support fuel qualification. 2. Establish operating margins for the fuel, beyond normal operating conditions. 3. Provide irradiated-fuel performance data and irradiated-fuel samples for post irradiation examination and safety testing. The primary objective of the AGR-5/6 test (Capsules 1, 2, 4, and 5) was to verify the successful performance of the reference-design fuel under normal operating conditions. The AGR-7 test (Capsule 3) was designed to explore fuel performance at higher temperatures. Its primary objective is to demonstrate the capability of the fuel to withstand conditions beyond normal operating conditions, in support of plant design and licensing. AGR 5/6/7 will also provide irradiated-fuel performance data based on the fission gas release from particles during irradiation. To achieve the test objectives, the AGR-5/6/7 experiment was irradiated in the northeast (NE) flux trap of the ATR with a planned duration of 500 effective full-power days (EFPDs). The NE flux trap was selected because its larger diameter provided greater flexibility for test-train design compared to the Large B positions used for the AGR-1 and AGR-2 irradiations, significantly enhancing test capabilities for the combined irradiation campaigns. Due to delays in the ATR schedule, the AGR-5/6/7 irradiation was significantly shorter than the originally planned 13-cycle schedule. Irradiation began on February 16, 2018 and ended on July 22, 2020, spanning nine ATR cycles (162B?168A) over two and a half years. Thus, the AGR-5/6/7 fuel compacts were irradiated for a total of approximately 360.9 EFPDs. Final burnup values, on a per-compact basis, ranged from 5.66% to 15.26% fissions per initial heavy metal atom, while fast fluence values ranged from 1.62 to 5.55 × 1025 n/m2 (E >0.18 MeV). Time-averaged volume-averaged fuel temperatures on a capsule basis at the end of irradiation ranged from 756°C in Capsule 5 to 1313°C in Capsule 3 excluding days with significantly lower temperature during the two short powered axial locator mechanism cycles, 163A and 167A. By the end of irradiation, 48 out of 54 installed thermocouples (TC) had failed (the bottom three capsules lost all TCs). During the first five cycles (162B ? 165A), fission-gas isotope release-rate-to-birthrate (R/B) ratios were stable in the 10-8?10-6 range, and no in-pile particle failures were observed based on the gross gamma counts. During this time, the high exposed kernel fraction and high fuel particle temperatures in Capsule 1 led to the maximum R/B value of around 2 ? 10-6 for Kr-85m. The fission gas release in all capsules started to increase from the second half of Cycle 166A, when a large number of in-pile particle failures occurred in Capsule 1 and a gas line problem in this capsule caused fission gas leakage at various degrees into the other four capsules. This gas line problem also prevented fission gas release measurement for Capsule 1 during the last three cycles due to gas flow isolation. By the end of irradiation, it is estimated that approximately 15 particle failures in Capsule 3, which was considered possible at the start of the irradiation, since the experiment was designed to operate beyond the high-temperature gas-cooled reactor (HTGR) normal operating temperature range. A few hundred in-pile particle failures were estimated for Capsule 1 by the end of Cycle 166A, but the total number of failures is unknown due to lack of fission gas release data in the later cycles. Additionally, four potential in-pile failures were identified for Capsule 2 during the last cycle, 168A. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on absence of the typical spikes in gross gamma counts and low failure estimates using the AGR-3/4 R/B per exposed kernel model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Analysis of Space-Conditioning Loads in Commercial Buildings

Space conditioning end-uses, which include heating, cooling, and ventilation, represent a significant fraction of commercial building energy use, with a wide variety of heating and cooling technology options available in the market. In the interest of improving the overall efficiency of heating, ventilation and air-conditioning (HVAC) technologies, governments, utilities and private sector entities have implemented a variety of market transformation policies that aim to influence consumer purchase decisions. To evaluate the costs and benefits of such programs, analysts typically postulate a hypothetical default equipment choice, and compare it to one that provides comparable service with lower energy and/or power use. The corresponding reduced operating cost provides a benefit that offsets the potential higher cost of improved efficiency. Typically, life-cycle cost or cash-flow analyses are used to quantify the net economic benefit. These analyses require the capability to assess how a given equipment design would perform across a broad range of characteristics, both of the building and of the local weather. While these assessments can be performed using customized building simulations, it is generally not practical to develop and validate detailed building simulation code to cover all the potential variations of equipment design and installation. An alternative, and somewhat simpler, approach is to solely use detailed building simulations to generate time series of heating and cooling loads in commercial buildings. These loads can then be used as input to more detailed, stand-alone engineering models that simulate HVAC system performance under different equipment designs. This approach was used to evaluate a range of high-efficiency commercial packaged air conditioner design options for the Department of Energy’s Appliance and Equipment Standards Program (DOE-EERE 2015). While there may be some loss of precision relative to full simulation, the accuracy of this approach is sufficient for practical applications of cost-benefit analysis. This report describes the development of a database of commercial building heating and cooling loads, generated using the EnergyPlus software package, a whole building energy use model supported by the Department of Energy (DOE-EERE 2020a). EnergyPlus takes as input a set of configuration files that describe the building itself (size, zoning, envelope characteristics, etc.) and the various systems within it (HVAC, lighting, water heating, etc.). This analysis uses a publicly available collection of commercial reference buildings (CRB), comprised of sixteen building types and three vintages (DOE-EERE 2020b; Deru et al 2011). Each building is simulated in eighteen different locations, covering a wide range of climatic conditions. The prototype building description files assign the type of HVAC equipment used, and capacities across climate zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Analysis of hematite attrition in a grid jet apparatus

Particulate attrition is of interest for novel carbon-capture processes such as chemical looping combustion because the makeup cost of oxygen carrier is a significant portion of operating cost. As such, models to study and predict attrition of various oxygen carriers in fluidized bed systems are being developed. One of the regions of concern in fluidized bed systems is the high-velocity jet region near gas distributors in a fluid bed. This work studies the attrition of hematite particles using a modified ASTM apparatus to measure the particle size distribution throughout the experiment. Bed weight and gas velocity were varied. Hematite particles above the corresponding threshold value had a severe variation of particle size distributions which decreased with time. Weight fractions of the sieves over time were fit to a linear, time-variant population balance model to offer insight into particle attrition. The first-order rate constant was modified as a decaying exponential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Removal of Metal Oxides from Reactor Metals/Alloys by Laser Ablation

Storage of low and intermediate level radioactive waste in a deep geological repository (L and ILW DGR): - Large volumes of metallic waste = increased production of hydrogen gas. - Large fraction of the metallic wastes are out-of-core reactor components. - Carbon steel feeder pipes. - Radioactive contamination is concentrated within the surface oxides. - Removal of surface contaminants would lead to volume reduction of metallic waste that is stored in L and ILW DGR. - Current decontamination techniques generate large volumes of secondary waste. Objectives: Oxide removal from carbon steel using a ytterbium fiber laser (1064 nm) with a particular focus on identifying the key laser parameters that affect the removal efficiency. Advantages of laser ablation: Generate low volumes of waste, Flexibility of deployment, Remote operation. Grow an oxide film on carbon steel with the incorporation of cobalt. Optimize laser cleaning technology for the removal of surface oxides. Incorporate an oxide collection system for the ablated particles. Conclusions: Objective 1 - Simulated inactive oxides on carbon steel, using Co{sup 2+} as a surrogate for Co radionuclide. Objective 2 - Optimized laser parameters (200 MW/cm{sup 2}, 3 scans, 15 μm overlap); Complete removal of cobalt from surface and bulk metal with optimized cleaning settings. Objective 3 - Glass fiber media H14 HEPA filters had highest collection yield; High oxide collection yields (∼100%). Future Work: Further testing of the liquid trap will verify the efficacy of the HEPA filter and determine the amount of ablated oxides that pass through. Testing of the laser ablation system needs to be done on other reactor materials along with the incorporation of different transition metal cations. Improvements to the oxide collection system are required to improve collection yields.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Comparison of unirradiated and irradiated AGR-2 TRISO fuel particle oxidation response

The silicon carbide (SiC) coating in a tristructural isotropic (TRISO) particle acts as a barrier to fission product release during reactor operation and accident scenarios. Oxidation and subsequent failure of the SiC layer during a rare air ingress event is a proposed mechanism for fission product release in a high-temperature gas-cooled reactor (HTGR). Although previous oxidation studies have analyzed unirradiated TRISO particle response, this study compared the oxidation behavior of irradiated and unirradiated TRISO particles from the second Advanced Gas Reactor Fuel Development and Qualification Program irradiation experiment (AGR-2). Particles with exposed SiC were subjected to six varying oxidizing tests in the Furnace for Irradiated TRISO Testing (FITT), examined for failure fraction with the Irradiated Microsphere Gamma Analyzer (IMGA) and characterized with focused ion beam and scanning/transmission electron microscopy techniques to analyze the oxide layer. Uncorrelated unirradiated particle failures throughout the series of exposures suggests that external factors inherent to the experiment increased particle failure sensitivity. However, irradiated particle observations indicated an increased failure response at 400 h 1400 °C in both 2% and 21% O 2 atmospheres above failure associated with external factors. Oxide thickness measurements after 400 h at 1400 °C revealed a greater oxidation rate than predicted by parabolic growth, which was attributed to the increased complexity of the oxide structure at longer exposure times. Altering the atmosphere from 21% to 2% O 2 reduced the average oxide thickness by approximately 12%–14% in both irradiated and unirradiated particles at 400 h 1400 °C. Altogether, the minor variations observed between irradiated and unirradiated particles in this study led to the conclusion that unirradiated TRISO particles can be used to approximate irradiated TRISO oxidation kinetics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Thermal loading analysis of the ring injection dump for the Spallation Neutron Source facility

The ring injection dump (RID) is the largest beam dump in the Spallation Neutron Source facility at Oak Ridge National Laboratory and accepts a fraction of the beam from the ring that is not captured in the ring during injection. Thermo-fluid modeling of the RID components was performed during the initial design by using the software ANSYS CFX to characterize the power rating for the RID. From the calculations, the power rating was lowered from 200 kW to 150 kW, at 1.3 GeV beam energy due to concerns regarding the heating of the concrete structure. The Proton Power Upgrade project scope included reevaluating the power limit to the RID in anticipation that increasing it would provide more operational flexibility. The focus was on the shielding because it was the limiting factor, and data collected from thermocouples installed on the structure were used to benchmark the analysis. A three-stage process was adopted in the latest study. First, the already existing steady-state calculations were validated and extended with the latest Monte Carlo N-Particle source term calculations. Second, transient calculations were conducted to capture the dynamic state of the system in response to the energy dumped to the RID during the following operational period of the beam from 2005–2019 (14 years), and sensitivity analyses of crucial parameters were performed to benchmark the model with reasonable accuracy. Finally, an idealized transient cycle analysis was conducted with realistic duty factors to predict the temperature distribution at higher beam powers. Overall, the results confirmed the accuracy of the original steady state model but allowed the development of a new validated transient model with higher accuracy for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Constraining the Twomey effect from satellite observations: issues and perspectives

Abstract. The Twomey effect describes the radiative forcing associated with a change in cloud albedo due to an increase in anthropogenic aerosol emissions. It is driven by the perturbation in cloud droplet number concentration (ΔNd, ant) in liquid-water clouds and is currently understood to exert a cooling effect on climate. The Twomey effect is the key driver in the effective radiative forcing due to aerosol–cloud interactions, but rapid adjustments also contribute. These adjustments are essentially the responses of cloud fraction and liquid water path to ΔNd, ant and thus scale approximately with it. While the fundamental physics of the influence of added aerosol particles on the droplet concentration (Nd) is well described by established theory at the particle scale (micrometres), how this relationship is expressed at the large-scale (hundreds of kilometres) perturbation, ΔNd, ant, remains uncertain. The discrepancy between process understanding at particle scale and insufficient quantification at the climate-relevant large scale is caused by co-variability of aerosol particles and updraught velocity and by droplet sink processes. These operate at scales on the order of tens of metres at which only localised observations are available and at which no approach yet exists to quantify the anthropogenic perturbation. Different atmospheric models suggest diverse magnitudes of the Twomey effect even when applying the same anthropogenic aerosol emission perturbation. Thus, observational data are needed to quantify and constrain the Twomey effect. At the global scale, this means satellite data. There are four key uncertainties in determining ΔNd, ant, namely the quantification of (i) the cloud-active aerosol – the cloud condensation nuclei (CCN) concentrations at or above cloud base, (ii) Nd, (iii) the statistical approach for inferring the sensitivity of Nd to aerosol particles from the satellite data and (iv) uncertainty in the anthropogenic perturbation to CCN concentrations, which is not easily accessible from observational data. This review discusses deficiencies of current approaches for the different aspects of the problem and proposes several ways forward: in terms of CCN, retrievals of optical quantities such as aerosol optical depth suffer from a lack of vertical resolution, size and hygroscopicity information, non-direct relation to the concentration of aerosols, difficulty to quantify it within or below clouds, and the problem of insufficient sensitivity at low concentrations, in addition to retrieval errors. A future path forward can include utilising co-located polarimeter and lidar instruments, ideally including high-spectral-resolution lidar capability at two wavelengths to maximise vertically resolved size distribution information content. In terms of Nd, a key problem is the lack of operational retrievals of this quantity and the inaccuracy of the retrieval especially in broken-cloud regimes. As for the Nd-to-CCN sensitivity, key issues are the updraught distributions and the role of Nd sink processes, for which empirical assessments for specific cloud regimes are currently the best solutions. These considerations point to the conclusion that past studies using existing approaches have likely underestimated the true sensitivity and, thus, the radiative forcing due to the Twomey effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Renewable Diesel Production through Stand-Alone and Co-Hydrotreating of Catalytic Fast Pyrolysis Oil

Catalytic Fast Pyrolysis (CFP) of biomass followed by hydrotreating can be used for the production of renewable transportation fuels, such as sustainable aviation fuel and renewable diesel, which can reduce the dependence on fossil fuels and help achieve greenhouse gas reduction goals. Stand-alone hydrotreating of CFP oil and co-hydrotreating of the CFP oil with straight-run diesel (SRD) showed potential to produce high-quality diesel products in this study. Stand-alone hydrotreating of CFP oil produced a diesel fraction with oxygen content < 0.1 wt.% and an indicated cetane number (ICN) of 45. Co-hydrotreating of CFP oil and SRD using a 1:4 volumetric ratio produced diesel products with 0.1 wt.% of oxygen and ICNs of 42-45. Compared to CoMo/Al2O3, NiMo/Al2O3 was a more attractive catalyst for co-hydrotreating of two streams, leading to a higher carbon efficiency, a higher selectivity of cycloalkanes, a higher ICN, and enhanced desulfurization. Compared to stand-alone hydrotreating of CFP oil, co-hydrotreating with SRD could be achieved at a less severe operating condition and it produced hydrotreated products with a higher carbon efficiency. Synergy between CFP oil and SRD during co-hydrotreating process was suggested by the better experimental results in comparison with theoretical interpolation results. Co-hydrotreating gave high biogenic carbon incorporations of 91%-97% as determined by C-14 analysis.

biomass↗
DOE OSTI2020

Detection of DoS Attacks Using ARFIMA Modeling of GOOSE Communication in IEC 61850 Substations

Integration of Information and Communication Technology (ICT) in modern smart grids (SGs) offers many advantages including the use of renewables and an effective way to protect, control and monitor the energy transmission and distribution. To reach an optimal operation of future energy systems, availability, integrity and confidentiality of data should be guaranteed. Research on the cyber-physical security of electrical substations based on IEC 61850 is still at an early stage. In the present work, we first model the network traffic data in electrical substations, then, we present a statistical Anomaly Detection (AD) method to detect Denial of Service (DoS) attacks against the Generic Object Oriented Substation Event (GOOSE) network communication. According to interpretations on the self-similarity and the Long-Range Dependency (LRD) of the data, an Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) model was shown to describe well the GOOSE communication in the substation process network. Based on this ARFIMA-model and in view of cyber-physical security, an effective model-based AD method is developed and analyzed. Two variants of the statistical AD considering statistical hypothesis testing based on the Generalized Likelihood Ratio Test (GLRT) and the cumulative sum (CUSUM) are presented to detect flooding attacks that might affect the availability of the data. Our work presents a novel AD method, with two different variants, tailored to the specific features of the GOOSE traffic in IEC 61850 substations. The statistical AD is capable of detecting anomalies at unknown change times under the realistic assumption of unknown model parameters. The performance of both variants of the AD method is validated and assessed using data collected from a simulation case study. We perform several Monte-Carlo simulations under different noise variances. The detection delay is provided for each detector and it represents the number of discrete time samples after which an anomaly is detected. In fact, our statistical AD method with both variants (CUSUM and GLRT) has around half the false positive rate and a smaller detection delay when compared with two of the closest works found in the literature. Our AD approach based on the GLRT detector has the smallest false positive rate among all considered approaches. Whereas, our AD approach based on the CUSUM test has the lowest false negative rate thus the best detection rate. Depending on the requirements as well as the costs of false alarms or missed anomalies, both variants of our statistical detection method can be used and are further analyzed using composite detection metrics.

IEC 61850 electrical substations↗
DOE OSTI2024

GCAM–GLORY v1.0: representing global reservoir water storage in a multi-sector human–Earth system model

Abstract. Reservoirs play a significant role in modifying the spatiotemporal availability of surface water to meet multi-sector human demands, despite representing a relatively small fraction of the global water budget. Yet the integrated modeling frameworks that explore the interactions among climate, land, energy, water, and socioeconomic systems at a global scale often contain limited representations of water storage dynamics that incorporate feedbacks from other systems. In this study, we implement a representation of water storage in the Global Change Analysis Model (GCAM) to enable the exploration of the future role (e.g., expansion) of reservoir water storage globally in meeting demands for, and evolving in response to interactions with, the climate, land, and energy systems. GCAM represents 235 global water basins, operates at 5-year time steps, and uses supply curves to capture economic competition among renewable water (now including reservoirs), non-renewable groundwater, and desalination. Our approach consists of developing the GLObal Reservoir Yield (GLORY) model, which uses a linear programming (LP)-based optimization algorithm and dynamically linking GLORY with GCAM. The new coupled GCAM–GLORY approach improves the representation of reservoir water storage in GCAM in several ways. First, the GLORY model identifies the cost of supplying increasing levels of water supply from reservoir storage by considering regional physical and economic factors, such as evolving monthly reservoir inflows and demands, and the leveled cost of constructing additional reservoir storage capacity. Second, by passing those costs to GCAM, GLORY enables the exploration of future regional reservoir expansion pathways and their response to climate and socioeconomic drivers. To guide the model toward reasonable reservoir expansion pathways, GLORY applies a diverse array of feasibility constraints related to protected land, population, water sources, and cropland. Finally, the GLORY–GCAM feedback loop allows evolving water demands from GCAM to inform GLORY, resulting in an updated supply curve at each time step, thus enabling GCAM to establish a more meaningful economic value of water. This study improves our understanding of the sensitivity of reservoir water supply to multiple physical and economic dimensions, such as sub-annual variations in climate conditions and human water demands, especially for basins experiencing socioeconomic droughts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Light new physics in B → K ( * ) ν ν ¯ ?

The study of the rare decays B → K ( * ) ν ν ¯ offers a window into the dynamics operating at the electroweak scale, allowing studies of the Standard Model and searches for heavy new physics. However, the analysis of these decays is also potentially sensitive to the on shell production of new light bosons X through the process B → K ( * ) X . In particular, Belle II has recently measured B + → K + ν ν ¯ , finding a 2.8 σ excess under the assumption of heavy new physics. Since this excess is rather localized in the kaon energy, a fit that includes the decay mode B + → K + X to the kinematic distributions prefers m X ≈ 2 GeV with branching fraction Br [ B → K X ] = ( 8.8 ± 2.5 ) × 10 − 6 and a significance of ≈ 3.6 σ . However, no excess was found in the measurements of B → K ( * ) ν ν ¯ , and a global analysis of the Belle II and data leads to Br [ B → K X ] = ( 5.1 ± 2.1 ) × 10 − 6 with a reduced significance of ≈ 2.4 σ . We then study various simplified dark-flavored models and present a possible UV completion based on a gauged B 3 − L 3 symmetry, highlighting the discovery potential of dedicated searches for B → K ( * ) X at Belle II. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2022

Liquefaction of pelleted corn cobs

Liquid slurries at high loadings (pumpable slurries) are critical to handle the streams between the unit operations in a biorefinery. Their liquefaction potential was initially tested with a small-scale method using 3FPU Ctec-2 /g dry biomass, and results were confirmed in 1L bioreactors experiments. Liquid slurries of cobs were observed in small scale tests (based on visual appearance, and reducing sugars released). Results were confirmed in 1L bioreactors under similar experimental conditions. Yield stress value in bioreactor was 34 ± 6Pa at 96 h, and glucan and xylan conversion to glucose and xylose, respectively, were 41 and 21%. In addition, work has been initiated with liquefaction of other pelleted corn stover fractions: Husk/Leaf (A),Husk (B), Stems (D) and controls (whole corn stover; F and G). Compositional analysis of the pellets showed glucan content (%) in the samples ranged from 31.9 to 36.1, xylan from 19.1 to 20.9, total lignin 13.1 to 18.8,and total ash 1.73 to 10.19. The moisture content ranged from 6.21 to 11.66%. Given these similar characteristics, we also examined liquefaction of these other samples, and they also showed that liquid slurries are possible. Work is continuing to define favourable conditions for liquid slurries formation.

Cruz, Antonio↗
DOE OSTI2023

Measurement of neutrino oscillations using neutrino and antineutrino beams in the NOvA experiment

NOvA is a long-baseline accelerator neutrino oscillation experiment using the NuMIneutrino beam from Fermilab. Its main physics goals are to probe the 3-flavour oscillationparameters: neutrino mass hierarchy, CP-violating phase dcp and octant of .23 mixingangle by observing electron neutrino appearance and muon neutrino disappearance. Twofunctionally identical detectors are placed off-axis from the centre of the NuMI beam.The near detector at Fermilab is 100 m underground, and the far detector is locatedon the surface at Ash River, 810 km away from the beam source. The initial neutrinobeam spectra are measured using the near detector data and the oscillation parametersare extracted by fitting the observed data to the predicted neutrino spectrum in the fardetector.This thesis is centered around how to improve the sensitivity of |.m232| and .23 measurementsin the muon neutrino disappearance analysis. NOvA will take data for about12 years. The operation of the NOvA experiment for each year costs tens of millions ofdollars, thus it is valuable to maximise the sensitivity of the analysis. Three samples ofmuon neutrino events are studied in this thesis to improve the analysis sensitivities. First,higher energy muon neutrinos are investigated by extending the energy range in NOvA’scurrent standard analysis. Second, for the sample of events used in NOvA’s existing analysis,a new energy estimator which has been developed to improve the neutrino energy resolution is considered. Furthermore, in addition to binning the events as function ofenergy and hadronic energy fraction, three particle identifiers are introduced to separateneutrino events by signal purity to reduce the effects from backgrounds. Third, an additionallower purity sample of muon neutrino charged current (CC) events that look similarto neutral current events and have not been included in NOvA’s existing analyses havebeen studied.This thesis reanalyses NOvA’s data used in the 2020 analysis, corresponding to anexposure of 13.60×1020 protons on target (POT) in the neutrino beam mode recordedfrom February 6, 2014 to March 20, 2020, and 12.50×1020 protons on target in theantineutrino beam mode recorded between June 29, 2016 to February 26, 2019. Thisthesis has implemented a fit to Asimov fake data, generated where sin2 .23 = 0.59 and.m232 = 2.5 × 10-3 eV2. These sensitivity studies show that the uncertainty range of|.m232| at 1 s in the new analysis is reduced by 5.5% and the significance of maximaldisappearance rejection improves by 7.7%, compared to the standard analysis. This isequivalent to adding 11-16% more data. The best fit values of the oscillation parametersfrom fitting to the far detector (FD) data with the new analysis are found to besin2 .23 = 0.568+0.025-0.043 (sin2 .23 = 0.454+0.046-0.026) and .m232 = 2.399+0.055-0.070 × 10-3 eV2 (.m232= -2.427+0.055-0.067 × 10-3 eV2) for the normal (inverted) hierarchy. The difference in thebest fit for sin2 .23 (.m232) between the new analysis and NOvA’s 2020 analysis is around2% (1.4%). The uncertainty range at 1 s for .m232 decrease by 8% (4%) for the normalhierarchy (inverted hierarchy) compared to the standard analysis. The uncertainty rangefor sin2 .23 is close to the standard analysis. This thesis also implements the fit from combiningelectron neutrino appearance and muon neutrino disappearance. The combinedanalysis shows that the best fit values are very close to the standard analysis. However,the uncertainty range of .m232 at 1 s is reduced by 3.7% using the new analysis. The maximaldisappearance significance is not improved in the new analysis, but the new analysisslightly improves the rejection of the disfavoured octant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗