Search NASA⌕ Search

SEARCH · Search NASA

Results for “tracking problem”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A low-latency graph computer to identify metastable particles at the Large Hadron Collider for real-time analysis of potential dark matter signatures

Abstract Image recognition is a pervasive task in many information-processing environments. We present a solution to a difficult pattern recognition problem that lies at the heart of experimental particle physics. Future experiments with very high-intensity beams will produce a spray of thousands of particles in each beam-target or beam-beam collision. Recognizing the trajectories of these particles as they traverse layers of electronic sensors is a massive image recognition task that has never been accomplished in real time. We present a real-time processing solution that is implemented in a commercial field-programmable gate array using high-level synthesis. It is an unsupervised learning algorithm that uses techniques of graph computing. A prime application is the low-latency analysis of dark-matter signatures involving metastable charged particles that manifest as disappearing tracks.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Advancing the Frontiers of Deep Learning for Low-Dose 3D Cone-Beam CT Reconstruction

X-ray computed tomography (CT) is an important noninvasive medical imaging modality for studying the structural details of internal organs. Image reconstruction in CT is an inverse problem of recovering an object's internal structure from the absorption profile of X-ray beams (sinogram) measured using a detector. The classical variational approach for CT reconstruction minimizes an energy functional using an appropriate iterative algorithm. Motivated by the success of deep learning (DL), researchers have begun to leverage training data and enhanced computing capabilities in recent years to produce high-fidelity reconstructed images. Nonetheless, much of the academic research in DL algorithms for CT has focused primarily on the two-dimensional setting (with simplified forward operators and noise model) for proofs-of-concept, and a comprehensive benchmarking of various classical and data-driven CT reconstruction approaches has not beenundertaken. The key objective of our CT reconstruction grand challenge was to promote methodological advancements for both classical and DL-based approaches for clinical CT with a reasonably accurately simulated 3D CT forward operator and noise model. We have utilized the publicly available LIDC-IDRI dataset and simulated sinograms and FDK images corresponding to two dose levels (clinical- and low-dose, constituting two tracks of the challenge) starting from the normal-dose images as the ground truth. In this paper, we summarize the motivation, context, and results of our challenge, and highlight the future research directions in DL for clinical CT.

X-ray tomography↗
DOE OSTI2024

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

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

Tracing the Rain Formation Pathways in Numerical Simulations of Deep Convection

Quantifying the microphysical process contributions to surface precipitation in numerical simulations can be challenging. This is due to the fact that many microphysical processes contribute to the formation and depletion of rain drops and there is almost always a spatial/temporal mismatch between where/when rain is formed and where/when it strikes the surface. In this work, we develop a tracing method that tracks the sources and sinks of raindrop mass and number as they are advected by the Weather Research and Forecasting model. Applying the method to an idealized squall line confirms that convective precipitation is dominated by warm rain processes (autoconversion and accretion) while stratiform precipitation is dominated by the melting of rimed and unrimed ice crystals. Sensitivity experiments in which the prescribed cloud drop number concentration is increased confirm the conventional wisdom that weakened autoconversion increases the fraction of raindrops originating from cold rain processes. The method also reveals that when applied to deep convection the Khairoutdinov and Kogan autoconversion scheme produces an excessive number of raindrops which are subsequently clipped in P3 microphysics to keep the rain size distribution within prescribed limits. This problem can mostly be mitigated by increasing the assumed radius for raindrops created by autoconversion.

58 GEOSCIENCES↗
DOE OSTI2021

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Updated SAM Model for the Molten Salt Reactor Experiment (MSRE)

The development of reference standard problems based on prototypical reactor designs is of particular importance to verify the adequacy of computer codes and evaluation models for specific reactor types. To support the multiphysics coupled simulation of molten-salt-fueled reactor (MSR) using SAM and Griffin computer codes for safety and licensing analysis, much efforts have been put into enhancing code capabilities and developing reference models for the MSR primary loop in SAM. In this work, a previously developed Molten Salt Reactor Experiment (MSRE) primary loop model was updated to include a two-dimensional (2-D) core region and external core components in one-dimension (1-D) or zero-dimension (0-D). To ensure accurate feedback calculation in multi-physics simulations, the delayed neutron precursor tracking model and solid graphite model were added in the SAM model. In addition, the 2-D and 1-D domains are tightly coupled using the recently developed single-solve approach in SAM. The updated model has been tested under both steady-state and transient scenarios to demonstrate its potential for the multi-physics simulation of MSRE with coupled SAM and Griffin.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Quantifying microbial roles in environmental iron oxidation via an integrated kinetics, `omics and metabolic modeling study (Final Report)

Iron oxyhydroxides are extremely reactive components of environmental systems, and therefore exert a strong influence on biogeochemical cycles. These oxyhydroxides strongly adsorb many biologically-relevant elements, including organic carbon and phosphate, as well as a wide range of metals including uranium and actinide species. Thus, the formation mechanism of iron oxyhydroxides are key to understanding both nutrient and contaminant cycling. Microorganisms can catalyze iron oxidation and promote the formation of Fe biominerals and thus are increasingly recognized as important players in biogeochemical cycling. However, it is completely unknown how much of environmental iron oxidation is biologically mediated versus abiotic, and various challenges in studying microbial iron oxidation have hindered accurate incorporation into hydrobiogeochemical models. The overarching goal of our work was to quantify and constrain microbial iron oxidation rates and use ‘omics to gain insight into the controls on this process, while developing tools to enable integration of biotic iron oxidation into hydrobiogeochemical models. Our work focused on the Savannah River Site (SRS) in South Carolina, where extensive microbial iron oxidation has been observed. At Tims Branch, part of the Argonne National Laboratory Wetland Hydrobiogeochemistry Science Focus Area (Argonne SFA), where groundwater discharges into a stream, iron-oxidizing microbial mats form and appear to be a major sink of uranium. In the wetlands that surround Tims Branch, there are wide swaths of iron microbial mats and flocs (mobilized mat). We measured biotic and abiotic iron oxidation rates using mats and water sampled from these sites and found that iron oxidation is primarily carried out by chemolithotrophic microorganisms. The resulting rate constants can be incorporated into models. These mats were characterized by metagenomics and metatranscriptomics, which showed that aerobic chemolithotrophs were the dominant iron-oxidizing bacteria (FeOB), and these included Gallionellaceae and Leptothrix, and possibly Rhodoferax, which is known as an Fe-reducer but may also oxidize Fe(II). This demonstrated that diverse FeOB can coexist and suggests that there are a range of niches and therefore drivers of chemolithotrophic iron oxidation. Analysis of reconstructed genomes strongly suggests that a major factor in diversity is carbon source, as genomes contained varied pathways for autotrophy and heterotrophy. We performed an in-depth analysis of Leptothrix ochracea genomes, since this sheath-former is one of the primary mat builders, yet its physiology remained unresolved. A combination of genomics, transcriptomics, and metabolic modeling suggest that L. ochracea grows mixotrophically using a combination of Fe(II) and organics for energy and both inorganic and organic carbon to create biomass. This contrasts with the largely autotrophic Gallionellaceae (Gallionella, Sideroxydans, and Ferriphaselus) also present in the mats and flocs. Remarkably, multiple FeOB, both Leptothrix and Gallionellaceae, showed activity in response to Fe(II) in live mat incubations. We tracked the gene expression of individual MAGs to Fe(II) and found that various autotrophic and heterotrophic FeOB responded to Fe(II), increasing expression of both carbon fixation and organic utilization genes. The results of the integrated field, kinetics, and omics studies give detailed insight into 1) the taxa that oxidize Fe, and 2) how they connect Fe, C, and N cycles. Towards the goal of connecting omics data to hydrobiogeochemical models, we worked with the KBase team to create a template metabolic model for chemolithotrophic iron oxidation. We initially modeled the well-characterized isolate Gallionellaceae Sideroxydans lithotrophicus, and also applied the model to the mixotroph L. ochracea. In all, we have characterized diverse FeOB in a representative wetland system and solved key problems that enable better incorporation of iron-oxidizing microbes into hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Limiter-aware adaptive excitation control for hydropower generators with legacy voltage regulation systems

Hydropower units often rely on legacy excitation systems with fixed automatic voltage regulator (AVR) settings, but parameter drift, changing operating conditions, and evolving grid-service demands can degrade terminal-voltage regulation over time. This problem should be addressed because poor voltage regulation affects reactive-power support, system stability, and reliable plant operation, and replacing installed AVR and limiter logic is costly and disruptive in practice. To address this issue, this paper proposes a retrofit-friendly enhancement for hydropower generator excitation systems that combines adaptive tuning of the AVR middle lead–lag branch with a small proportional-integral correction for removing steady bias and slow drift. The adaptive laws are normalized, projection-bounded, and equipped with leakage, while adaptation is enabled only in the middle branch. Thus, the intent of the underexcitation and overexcitation limiters, as well as legacy current and rate limits, remains unchanged. Stability is analyzed under incremental performance assumptions, yielding boundedness and asymptotic tracking in the absence of saturation. The method is evaluated using real operational data from a utility-scale hydropower unit at Rocky Reach Dam, Unit C-8. Compared with the fixed AVR baseline, the proposed method reduces terminal-voltage mean squared error by about 70%, removes long periods of steady-state bias, and improves transient tracking while coexisting with limiter actions. These results show that substantial voltage-control improvement can be achieved on existing hydropower units with minimal integration effort and without replacing the installed AVR structure.

ESST5B model↗
DOE OSTI2024

Preserving Tracer Correlations in Moment-Based Atmospheric Transport Models

A linear non-diffusive algorithm for advective transport is developed that greatly improves the detail at which aerosols and clouds can be represented in atmospheric models. Linear advection schemes preserve tracer correlations but the most basic linear scheme is rarely used by atmospheric modelers on account of its excessive numerical diffusion. Higher-order schemes are in widespread use, but these present new problems as nonlinear adjustments are required to avoid occurrences of negative concentrations, spurious oscillations, and other non-physical effects. Generally successful at reducing numerical diffusion during the advection of individual tracers, for example, particle number or mass, the higher-order schemes fail to preserve even the simplest of correlations between interrelated tracers. As a result, important attributes of aerosol and cloud populations including radial moments of particle size distributions, molecular precursors related through chemical equilibria, aerosol mixing state, and distribution of cloud phase are poorly represented. We introduce a new transport scheme, minVAR, that is both non-diffusive and preservative of tracer correlations, thereby combining the best features of the basic and higher-order schemes while enabling new features such as the tracking of sub-grid information at arbitrarily fine scales with high computational efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Optimizing Vehicle Fleet and Assignment for Concentrating Solar Power Plant Heliostat Washing

Concentrating solar power central-receiver plants use thousands of sun-tracking mirrors, i.e., heliostats, to reflect sunlight to a central receiver, which collects and uses the heat to generate electricity. Over time, soiling reduces the reflectivity of the heliostats and, therefore, the efficiency of the system. Current industry practice sends vehicles to wash heliostats in an ad hoc fashion. We present a mixed-integer nonlinear program that determines wash vehicle fleet size, mix, and assignment of wash crews to heliostats to minimize the sum of (i) the revenues lost due to heliostat soiling, (ii) the costs of hiring wash crews and operating the vehicles, and (iii) the costs of purchasing wash vehicles. We establish conditions for convexity of the objective function, and then propose a decomposition method that enables near-optimal solutions to the wash vehicle fleet sizing and assignment problem on the order of a couple of minutes. Furthermore, these solutions yield hundreds of thousands of dollars in savings per year over current industry practices.

14 SOLAR ENERGY↗
DOE OSTI2021

Realization of a high luminosity muon beam at the Jefferson Lab

The proton is the major component of the visible matter of our universe, so being able to determine its fundamental properties, such as its charge radius, is of vital importance. In recent years, several experiments have been performed for the accurate measurement of the proton radius with different experimental techniques. The main and most common are based on two different techniques: ep scattering and Lamb shift spectroscopy of regular hydrogen atoms. The latter can also be performed on muonic hydrogen atoms with a significant improvement in accuracy. Recent measurements that use this method have highlighted a discrepancy between the results obtained with the diffusion of electrons and those obtained with the spectroscopic technique leading to the so-called ?proton radius puzzle?. The explanation of why this difference between the measured values exists is still a matter of discussion. The last attempt to address the puzzle was the PRad Experiment. It was performed with a unique experimental setup that gave the scientists an excellent control on the systematic errors that are usually present in electron-proton scattering experiments. This has significantly increased the precision of the measurement, bringing it in the same area as the spectroscopic ones done on muonic hydrogen atoms, which are considered the most precise. However, to solve the puzzle once and for all, even more precise electron scattering measurements are needed. Alternatively, new methods can be used, such as eg. muon-proton scattering at low transferred momenta in order to lower the systematic errors that comes with the measurement. This thesis tries to address this problem by studying the set up of an intense muon beam at the Jefferson Laboratory (JLab) that can be used to carry out a high-precision m ? p scattering experiment. Using Montecarlo simulations I studied the muon production resulting from the interaction of the primary electron beam with the beam-dump at JLab Hall-A. The simulations were performed using the FLUKA package via the FLAIR graphical interface. In this study, I modified an existing geometry (built for the ?Beam Dump eXperiment?) and performed several studies in order to optimize the muon beam parameters. In order to estimate the beam intensity, I studied the muon generation and transportation through the concrete bunker surrounding the Hall-A beam-dump. Detailed studies on muon attenuation and dispersion at different depth allowed me to track the beam profile in space and momentum as a function of the distance from the beam-dump. Last, but not least, I studied a possible tagging system for the muon beam. The tagging system will allow us to know with good precision (O(1=1000)) the energy of each muon of the beam. Based on the results obtained, I created an event generator that incorporates the detailed description of the muon beam (energy, position, emittance) useful to define and optimize the experimental setup for a future measurement of the proton radius at the JLab.

Fulci, Antonino↗
DOE OSTI2024

Avoiding fusion plasma tearing instability with deep reinforcement learning

For stable and efficient fusion energy production using a tokamak reactor, it is essential to maintain a high-pressure hydrogenic plasma without plasma disruption. Therefore, it is necessary to actively control the tokamak based on the observed plasma state, to manoeuvre high-pressure plasma while avoiding tearing instability, the leading cause of disruptions. This presents an obstacle-avoidance problem for which artificial intelligence based on reinforcement learning has recently shown remarkable performance. However, the obstacle here, the tearing instability, is difficult to forecast and is highly prone to terminating plasma operations, especially in the ITER baseline scenario. Previously, we developed a multimodal dynamic model that estimates the likelihood of future tearing instability based on signals from multiple diagnostics and actuators. Here we harness this dynamic model as a training environment for reinforcement-learning artificial intelligence, facilitating automated instability prevention. We demonstrate artificial intelligence control to lower the possibility of disruptive tearing instabilities in DIII-D, the largest magnetic fusion facility in the United States. The controller maintained the tearing likelihood under a given threshold, even under relatively unfavourable conditions of low safety factor and low torque. In particular, it allowed the plasma to actively track the stable path within the time-varying operational space while maintaining H-mode performance, which was challenging with traditional preprogrammed control. This controller paves the path to developing stable high-performance operational scenarios for future use in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

One-shot omnidirectional pressure integration through matrix inversion

In this work, we present a method to perform 2D and 3D omnidirectional pressure integration from velocity measurements with a single-iteration matrix inversion approach. This work builds upon our previous work, where the rotating parallel ray approach was extended to the limit of infinite rays by taking continuous projection integrals of the ray paths and recasting the problem as an iterative matrix inversion problem. This iterative matrix equation is now 'fast-forwarded' to the 'infinity' iteration, leading to a different matrix equation that can be solved in a single step, thereby presenting the same computational complexity as the Poisson equation. We observe computational speedups of ~10 6 when compared to brute-force omnidirectional integration methods, enabling the treatment of grids of ~10 9 points and potentially even larger in a desktop setup at the time of publication. Further examination of the boundary conditions of our one-shot method shows that omnidirectional pressure integration implements a boundary condition where the boundary points are treated as interior points to the extent that information is available. Finally, we show how the method can be extended from the regular grids typical of particle image velocimetry to the unstructured meshes characteristic of particle tracking velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Fault diagnosis and fault tolerant control for T-S fuzzy stochastic distribution systems subject to sensor and actuator faults

The problem of fault diagnosis (FD) and fault tolerant control (FTC) for a class of Takagi-Sugeno (T-S) fuzzy stochastic distribution control (SDC) systems subject to sensor and actuator faults is discussed in this paper. First, fuzzy logic models are used to approximate the output probability density function (PDF). Next, an adaptive augmented state/fault diagnosis observer is proposed to estimate the system state, sensor and the actuator faults simultaneously. New expected weights based on the sensor fault estimation information and a PI-type fuzzy feedback fault tolerant (FT) controller are designed to compensate the effect of sensor fault and actuator fault simultaneously. When the sensor fault occurs, the expected objective is redesigned to compensate the sensor fault. Meanwhile, the PI controller can compensate the effect of actuator fault, and the output PDF of the system can still track the desired PDF after the fault occurs. Finally, an example of quality distribution control in chemical reaction process is given to confirm the effectiveness of the algorithm.

42 ENGINEERING↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Cybersecurity for Distributed Energy Resources: All Hazards Approach for Grid Modernization

Distributed energy resources (DER) sit at the intersection of IoT and critical infrastructure. As we work towards clean energy and decarbonization targets, high rates of growth of DER are expected, making them an important component of generation and grid services. They have been a part of the explosion of internet connected devices developed to solve real-world problems and make life easier by providing clean, local energy and enable the smart management of grid services. However, the quick-to-market, low-cost drivers for DER, combined with a historically relative low-impact has meant that cybersecurity has not been a top priority. So as DER penetration increases, will DER be a part of the growing challenge of securing the grid, or part of the solution for necessary grid modernization? In this talk, we will discuss the factors that have led to our current position and what makes cybersecurity for DER unique. A wide range of research is conducted at national labs, in partnership with industry and academia, to inform everything from standards and regulation, to the development of cyber-physical anomaly detection through use of digital twins, and even advanced threat analytics to SBOM and HBOM tracking. Even with this multi-faceted approach to the challenges, the question remains: is this too little too late or are we doing enough to secure the grid for the next 30 years?

14 SOLAR ENERGY↗
DOE OSTI2022

Bayesian State-Space Modeling Framework for Understanding and Predicting Golden Eagle Movements Using Telemetry Data

Predicting raptor movements through a wind power plant under given atmospheric and topographical conditions is a crucial first step in the overall goal of quantifying the risk of turbine-related collisions and mortalities. Extracting behavioral traits of golden eagles (Aquila chrysaetos) from telemetry data requires the fusion of noisy and sparse movement data (location, heading, velocity) with a stochastic mathematical representation of the eagles' decision-making processes. In this study, we framed this problem in a Bayesian state-space framework where both observations and decision-making are assumed to be stochastic processes connected through hidden states (mode of flight, intent), and the unknown model parameters are assumed to be random variables that are calibrated using the available telemetry data. This framework allowed for rigorous consideration of underlying uncertainties while allowing for both data and prior biological knowledge to contribute to a probabilistic and predictive agent-based movement model. We implemented and applied the Bayesian framework to understand movement behavior of 23 GPS-tagged golden eagles travelling in the western US for years 2019 and 2020. Our preliminary findings show that the Bayesian state-space framework provides a robust inverse modeling apparatus to decode eagle behavioral characteristics from telemetry data. This study was primarily aimed at verifying and validating the framework with selected golden eagle tracks (both long- and short-ranged), with future research aimed at extending the framework to include multi-mode flight, consideration of atmospheric data and uplift mechanisms, eagle-to-eagle interaction, and eagle-to-turbine interaction.

Bayesian modeling↗
DOE OSTI2022

The “Fresnel Equations” for Diffuse radiation on Inclined photovoltaic Surfaces (FEDIS)

Here the well-known Fresnel equations solve for the reflection and transmission of light for precise incident angles. The transmission of diffuse radiation incident on a planar or domed surface is often needed for real-world applications. Due to the complexity of the Fresnel equations, the analytical solution of the integration has hitherto been unobtainable over the last centuries. Therefore, this problem was numerically solved by integrating the angular transmittances in space often leading to substantial computing burden and bias in the results. To efficiently estimate the solar energy resource for a glass-covered photovoltaic (PV) module, we derive an analytical solution of diffuse transmission based on the rigorous integration of an alternate form of the Fresnel equations. The approach leads to a simple yet accurate relative transmittance model that reconciles the solar energy sensed by pyranometers and PV panels. With limited and clearly stated approximations, the complex mathematical derivation resulted in an elegant solution. An experiment using 1-year of data at the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) shows that the new model dramatically decreases the disparity between the solar radiation measurements by a Kipp and Zonen CM Pyranometer 22 (CMP22) and an IMT reference cell on a 1-axis tracking system. The solution in this paper can be widely used in scientific and engineering research, development, and applications wherever the Fresnel equations are used.

14 SOLAR ENERGY↗