Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 OSTI2022

Status Update on the Development of Transducers and Bonding Techniques for Enabling Acoustic Measurements of Damage in Microreactor Components

This report provides an overview of potential sensors and sensor-bonding techniques to enable the online acoustic interrogation of microreactor components and enhance structural health monitoring capabilities. The report focuses primarily on optical fiber–based acoustic sensors and describes initial experimental progress toward the deployment of these sensors for microreactor applications. The general approach is to monitor the resonant frequencies of microreactor components and search for evidence of structural defects that could indicate imminent failure. If properly identified, then the components could be repaired during the next reactor outage to prevent costly unplanned shutdowns. The ability to monitor the structural health of components could also reduce the need for time-consuming visual inspections and reduce staffing to improve microreactor economic viability. Increased sensor density is also one of the first steps to moving toward eventual semiautonomous operation. The expected microreactor conditions in which acoustic sensors must survive are characterized, including temperatures, neutron fluences, thermomechanical strains, and vibrational frequencies. Optical fiber–based acoustic sensors are identified as an attractive candidate for acoustic monitoring because of their high accuracy, immunity to electromagnetic interference, and resiliency in high-temperature, high-radiation environments. Optical fiber–based intrinsic sensors, such as type-II fiber-Bragg gratings and Fabry-Pérot Cavities (FPCs), are particularly attractive for a microreactor environment because of their high temperature stability, and FPCs also enable higher frequency interrogation with a lower sensitivity to radiation-induced drift. This report describes multiple interrogation systems, but the best interrogation system for a given situation will depend on the specific microreactor application, including the desired acoustic vibrational amplitudes, modes, and resonant frequencies. Initial experiments included fabricating three FPCs, tack-welding these FPCs to stainless-steel pipes or rods and performing room-temperature acoustic sensing tests to capture the vibrational frequency content. Peaks were identified in the measured frequency spectra and compared with the theoretical fundamental frequencies obtained from Euler-Bernoulli beam theory. Two of the three FPCs measured vibrational frequencies that generally matched those obtained from theory. Future work will include similar testing on pipes or other microreactor components with intentional flaws to evaluate the ability to determine changes in resonant frequencies. Finally, these tests will be repeated at high temperatures, potentially with an applied thermomechanical stress, to include environmental conditions similar to those for a microreactor application.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part I - Analysis with a Small Sample Size

Computer simulations, especially of complex phenomena, can be expensive, requiring high-performance computing resources. Often, to understand a phenomenon, multiple simulations are run, each with a different set of simulation input parameters. These data are then used to create an interpolant, or surrogate, relating the simulation outputs to the corresponding inputs. When the inputs and outputs are scalars, a simple machine learning model can suffice. However, when the simulation outputs are vector valued, available at locations in two or three spatial dimensions, often with a temporal component, creating a surrogate is more challenging. In this report, we use a two-dimensional problem of a jet interacting with high explosives to understand how we can build high-quality surrogates. The characteristics of our data set are unique - the vector-valued outputs from each simulation are available at over two million spatial locations; each simulation is run for a relatively small number of time steps; the size of the computational domain varies with each simulation; and resource constraints limit the number of simulations we can run. We show how we analyze these extremely large data-sets, set the parameters for the algorithms used in the analysis, and use simple ways to improve the accuracy of the spatio-temporal surrogates without substantially increasing the number of simulations required.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Observation and origin of multi-THz quantum beats in superfluorescent emission from argon

We report observation of multi-THz quantum beats in the amplitude of superfluorescence in pure argon gas, using a time delayed 780-nm probe beam that interrogates the relevant population densities. The superfluorescence is induced by three-photon absorption of an ultrashort 260.5-nm laser pulse that creates coherent and coupled populations in multiple closely spaced 3⁢d and 5⁢s Ar levels. Argon pressure and pump laser focusing parameters are seen to modify the frequency of the quantum beats, which is attributed to the absorption of the self-generated, nonresonant third harmonic of the pump from Ar near the zero of the real part of the overlapping linear susceptibilities of the closely spaced 3⁢d and 5⁢s excited states. The self-absorption of the phase-matched THG photons creates a coherent population of atoms which acts like an atomic excited state (3⁢ω-induced state). Furthermore, we have also observed a new two-step superfluorescent cascade that results from the population redistribution enabled by the 780-nm probe pulse between the 3⁢d'[5/2] 3 and the 3⁢ω-induced state and gives rise to a previously unknown transition that emits on 1587 nm.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

An implicit barotropic mode solver for MPAS-ocean using a modern Fortran solver interface

Here, we demonstrate use of a modern Fortran solver interface to manage solver algorithms for an implicit barotropic mode solver in the Model for Predictions Across Scales-Ocean (MPAS-O). ForTrilinos, a Fortran interface to Trilinos that contains a large collection of solver capabilities written in C++, has been implemented in MPAS-O to provide access to a suite of linear solver options. By virtue of the simplified wrapper and interface generator (SWIG) automation tool that generates modern Fortran interfaces to C++ code, we were able to implement the Fortran solver interface in MPAS-O using a familiar Fortran coding style while minimizing performance degradation. The ForTrilinos solver interface is written within MPAS-O’s time stepping modules as a subroutine in conjunction with MPAS-O code. Applied to an idealized ocean and a high-resolution realistic ocean test case, parallel performance of ForTrilinos solvers is examined. It is found that parallel scalability of the ForTrilinos solvers is highly dependent on the number of global synchronization points per solver iteration in each iterative solver algorithm. ForTrilinos solvers perform best compared to the Fortran hand-crafted (FHC) solver when the amount of work per processor is large enough. However, parallel scalability is better with the FHC solver and so when the work per core is modest FHC outperforms ForTrilinos. The intercomparison between the ForTrilinos and FHC solvers reveals that this performance hit in the ForTrilinos solver mostly comes from the global synchronization process, while suggesting that the matrix-vector multiplication process in the FHC solver needs to be optimized for better performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Measurement and calibration of non-linear shear terms in galaxy cluster fields

Introduction Galaxy cluster lensing is a powerful tool for measuring the mass of galaxy clusters, but accurate shear measurement and calibration are critical to obtaining reliable results. This study focuses on the measurement and calibration of weak lensing shears to improve mass estimates in cluster lensing. To deal with the problem, we first developed an image simulation pipeline, jedisim, which utilizes galaxy images extracted from the Hubble Space Telescope (HST) Ultra Deep Field (UDF) and the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Methods The simulations represent realistic galaxy distributions and morphologies as input sources. The foreground halo with a Navarro–Frenk–White (NFW) profile is constructed such that the lensing signals of background galaxies can be measured by the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) Science Pipelines. By comparing the measured reduced shear g meas and the true reduced shear g true , we observe non-linearity up to g ≲ 0.6 . We fit polynomials to the data with quadratic correction adequate to g ≲ 0.4 . Meanwhile, we conduct mass estimates using the pzmassfitter code on four different clusters. Results The mass estimate results are significantly improved after applying the shear calibration derived from the present work—from 4.954 ± 0.504 × 1 0 14 M ⊙ to 10.507 ± 0.498 × 1 0 14 M ⊙ after calibration for a simulated cluster with the mass of 10 × 1 0 14 M ⊙ . In multiple cases of validation, the estimated results are all consistent with true cluster mass. Discussion This study yields the first relationship between reality and shape measurement of the LSST Science Pipelines and serves as the first step toward the overall goal of mass calibration in cluster lensing. By addressing the challenges in shear measurement and calibration, we aim to enhance the accuracy and reliability of mass estimates in galaxy cluster lensing studies.

Astronomy & Astrophysics↗
DOE OSTI2022

Using probabilistic solar power forecasts to inform flexible ramp product procurement for the California ISO

How can independent system operators (ISOs) take advantage of probabilistic solar forecasts to lower generation costs and improve reliability of power systems? We discuss one three-step approach for doing so, focusing on how such forecasts might help the California Independent System Operator (CAISO) prepare unexpected net load ramps, where net load equals gross demand minus wind and solar production. First, we enhance an existing solar forecasting system to provide well-calibrated hours-ahead probabilistic forecasts. We then relate the degree of uncertainty reflected in the forecasted prediction intervals (independent variables) to error distributions for net load ramp forecasts for the CAISO real-time market (dependent variable) using machine learning and quantile regression. Projected ramp forecast errors conditioned on solar uncertainty are translated into flexible ramp requirements that therefore reflect real-time meteorological and solar conditions, improving on typical ISO procedures. Detailed descriptions are provided on the quantile regression and kth-nearest neighbor categorization methods for accomplishing that translation. Finally, a multiple time-scale look-ahead market simulation model is applied to a 118-bus IEEE Reliability Test System, modified to represent the CAISO generation mix and demand distributions. The model runs quantify how solar-conditioned ramp requirements can, first, decrease operating costs by reducing requirements compared to often conservative unconditional methods and, second, decrease generation scarcity events and consequently improve reliability by increasing flexibility requirements at times when unconditional forecast-based requirements understate actual ramp uncertainty. Solar-conditioned ramp requirements are found to reduce generation operating costs by about 2% for the test system (which would be equivalent to over $\$100$ million per year for a CAISO-size system).

14 SOLAR ENERGY↗
DOE OSTI2020

Quandary

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

Petersson, NilsA.↗
DOE OSTI2026

TDCOSMO XXIII. Measurement of the Hubble constant from the doubly lensed quasar HE 1104−1805

Time-delay cosmography leverages strongly lensed quasars to measure the Universe’s current expansion rate, H 0 , independently from other methods. The latest TDCOSMO milestone measurement primarily used quadruply lensed quasars for their mass profile constraints. However, doubly lensed quasars, being more abundant and offering precise time delays, could expand the sample by a factor of 5, significantly advancing towards a 1% precision measurement of H 0 . We present the first TDCOSMO analysis of a doubly imaged source, HE 1104−1805, including the measurement of the four necessary ingredients. First, by combining 17 years of data from the SMARTS, Euler, and WFI telescopes, we measured a time delay of 176.3 +11.4 −10.3 days. Second, using MUSE data, we extracted stellar velocity dispersion measurements in three radial bins with 5% to 13% precision. Third, employing F160W HST imaging for lens modelling and marginalising over various modelling choices, we measured the Fermat potential difference between the images. Fourth, using wide-field imaging, we measured the convergence added by objects not included in the lens modelling. By combining these four ingredients, we measured the time delay distance and the angular diameter distance to the deflector, favouring a power-law mass model over a baryonic and dark matter composite model. The measurement was performed blindly to prevent experimenter bias and resulted in a Hubble constant of H 0 = +5.8 −5.0 × λ int km s −1 Mpc −1 , where λ int is the internal mass sheet degeneracy parameter. This is in agreement with the TDCOSMO-2025 milestone and its precision for λ int = 1 is comparable to that obtained with the best-observed quadruply lensed quasars (4–6%). This work is a stepping stone towards a precise measurement of H 0 using a large sample of doubly lensed quasars, supplementing the current sample. The next TDCOSMO milestone paper will include this system in its hierarchical analysis, constraining λ int and H 0 jointly with multiple lenses.

cosmological parameters↗
DOE OSTI2021

Planets Across Space and Time (PAST). II. Catalog and Analyses of the LAMOST–Gaia–Kepler Stellar Kinematic Properties

The Kepler telescope has discovered over 4000 planets (candidates) by searching ~200,000 stars over a wide range of distance (order of kpc) in our Galaxy. Characterizing the kinematic properties (e.g., Galactic component membership and kinematic age) of these Kepler targets (including the planet candidate hosts) is the first step toward studying Kepler planets in the Galactic context, which will reveal fresh insights into planet formation and evolution. In this paper, the second part of the Planets Across the Space and Time (PAST) series, by combining the data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Gaia and then applying the revised kinematic methods from PAST I, we present a catalog of kinematic properties (i.e., Galactic positions, velocities, and the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 35,835 Kepler stars. Further analyses of the LAMOST–Gaia–Kepler catalog demonstrate that our derived kinematic age reveals the expected stellar activity-age trend. Furthermore, we find that the fraction of thin (thick) disk stars increases (decreases) with the transiting planet multiplicity (N {sub p} = 0, 1, 2 and 3+) and the kinematic age decreases with N {sub p}, which could be a consequence of the dynamical evolution of planetary architecture with age. The LAMOST–Gaia–Kepler catalog will be useful for future studies on the correlations between the exoplanet distributions and the stellar Galactic environments as well as ages.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

HarDWR - Harmonized Water Rights Records

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - Recalculated based on data sourced from WestDAAT - Changed using a Site ID column to identify unique records to using aa combination of Site ID and Allocation ID - Removed the Water Management Area (WMA) column from the harmonized records. The replacement is a separate file which stores the relationship between allocations and WMAs. This allows for allocations to contribute to water right amounts to multiple WMAs during the subsequent cumulative process. - Added a column describing a water rights legal status - Added "Unspecified" was a water source category - Added an acre-foot (AF) column - Added a column for the classification of the right's owner v1.02 - Added a .RData file to the dataset as a convenience for anyone exploring our code. This is an internal file, and the one referenced in analysis scripts as the data objects are already in R data objects. v1.01 - Updated the names of each file with an ID number less than 3 digits to include leading 0s v1.0 - Initial public release Description Here we present an updated database of Western U.S. water right records. This database provides consistent unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of seven broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of inter-sectoral water allocation changes though water markets (Grogan et al., *in review*). Specifically, the data were formatted for use as input to a process-based hydrologic model, Water Balance Model (WBM), with a water rights module (Grogan et al., *in review*). While this specific study motivated the development of the database presented here, water management in the U.S. West is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al., 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. West. We produced the water rights database presented here in four main steps: (1) data collection, (2) data quality control, (3) data harmonization, and (4) generation of cumulative water rights curves. Each of steps (1)-(3) had to be completed in order to produce (4), the final product that was used in the modeling exercise in Grogan et al. (*in review*). All data in each step is associated with a spatial unit called a Water Management Area (WMA), which is the unit of water right administration utilized by the state in which the right came from. Steps (2) and (3) required use to make assumptions and interpretation, and to remove records from the raw data collection. We describe each of these assumptions and interpretations below so that other researchers can choose to implement alternative assumptions an interpretation as fits their research aims. Motivation for Changing Data Sources The most significant change has been a switch from collecting the raw water rights directly from each state to using the water rights records presented in WestDAAT, a product of the Water Data Exchange (WaDE) Program under the Western States Water Council (WSWC). One of the main reasons for this is that each state of interest is a member of the WSWC, meaning that WaDE is partially funded by these states, as well as many universities. As WestDAAT is also a database with consistent categorization, it has allowed us to spend less time on data collection and quality control and more time on answering research questions. This has included records from water right sources we had previously not known about when creating v1.0 of this database. The only major downside to utilizing the WestDAAT records as our raw data is that further updates are tied to when WestDAAT is updated, as some states update their public water right records daily. However, as our focus is on cumulative water amounts at the regional scale, it is unlikely most records updates would have a significant effect on our results. The structure of WestDAAT led to several important changes to how HarWR is formatted. The most significant change is that WaDE has calculated a field known as `SiteUUID`, which is a unique identifier for the Point of Diversion (POD), or where the water is drawn from. This separate from `AllocationNativeID`, which is the identifier for the allocation of water, or the amount of water associated with the water right. It should be noted that it is possible for a single site to have multiple allocations associated with it and for an allocation to be able to be extracted from multiple sites. The site-allocation structure has allowed us to adapt a more consistent, and hopefully more realistic, approach in organizing the water right records than we had with HarDWR v1.0. This was incredibly helpful as the raw data from many states had multiple water uses within a single field within a single row of their raw data, and it was not always clear if the first water use was the most important, or simply first alphabetically. WestDAAT has already addressed this data quality issue. Furthermore, with v1.0, when there were multiple records with the same water right ID, we selected the largest volume or flow amount and disregarded the rest. As WestDAAT was already a common structure for disparate data formats, we were better able to identify sites with multiple allocations and, perhaps more importantly, allocations with multiple sites. This is particularly helpful when an allocation has sites which cross WMA boundaries, instead of just assigning the full water amount to a single WMA we are now able to divide the amount of water between the number of relevant WMAs. As it is now possible to identify allocations with water used in multiple WMAs, it is no longer practical to store this information within a single column. Instead the stAllocationToWMATab.csv file was created, which is an allocation by WMA matrix containing the percent Place of Use area overlap with each WMA. We then use this percentage to divide the allocation's flow amount between the given WMAs during the cumulation process to hopefully provide more realistic totals of water use in each area. However, not every state provides areas of water use, so like HarDWR v1.0, a hierarchical decision tree was used to assign each allocation to a WMA. First, if a WMA could be identified based on the allocation ID, then that WMA was used; typically, when available, this applied to the entire state and no further steps were needed. Second was the spatial analysis of Place of Use to WMAs. Third was a spatial analysis of the POD locations to WMAs, with the assumption that allocation's POD is within the WMA it should belong to; if an allocation still had multiple WMAs based on its POD locations, then the allocation's flow amount would be divided equally between all WMAs. The fourth, and final, process was to include water allocations which spatially fell outside of the state WMA boundaries. This could be due to several reasons, such as coordinate errors / imprecision in the POD location, imprecision in the WMA boundaries, or rights attached with features, such as a reservoir, which crosses state boundaries. To include these records, we decided for any POD which was within one kilometer of the state's edge would be assigned to the nearest WMA. Other Changes WestDAAT has Allowed In addition to a more nuanced and consistent method of assigning water right's data to WMAs, there are other benefits gained from using the WestDAAT dataset. Among those is a consistent categorization of a water right's legal status. In HarDWR v1.0, legal status was effectively ignored, which led to many valid concerns about the quality of the database related to the amounts of water the rights allowed to be claimed. The main issue was that rights with legal status' such as "application withdrawn", "non-active", or "cancelled" were included within HarDWR v1.0. These, and other water rights status' which were deemed to not be in use have been removed from this version of the database. Another major change has been the addition of the "unspecified water source category. This is water that can come from either surface water or groundwater, or the source of which is unknown. The addition of this source category brings the total number of categories to three. Due to reviewer feedback, we decided to add the acre-foot (AF) column so that the data may be more applicable to a wider audience. We added the ownerClassification column so that the data may be more applicable to a wider audience. File Descriptions The dataset is a series of various files organized by state sub-directories. In addition, each file begins with the state's name, in case the file is separate from its sub-directory for some reason. After the state name is the text which describes the contents of the file. Here is each file described in detail. Note that st is a placeholder for the state's name. stFullRecords_HarmonizedRights.csv: A file of the complete water records for each state. The column headers for each of this type of file are: state - The name of the state to which the allocations belong to. FIPS - The two digit numeric state ID code. siteID - The site location ID for POD locations. A site may have multiple allocations, which are the actual amount of water which can be drawn. In a simplified hypothetical, a farm stead may have an allocation for "irrigation" and an allocation for "domestic" water use, but the water is drawn from the same pumping equipment. It should be noted that many of the site ID appear to have been added by WaDE, and therefore may not be recognized by a given state's water rights database. allocationID - The allocation ID for the water right. For most states this is the water right ID, and what is recommended to use should a right be looked up on a given state's water rights database. The water amounts associated with these IDs tend to be finer scaled than those associated with siteID. It should be noted that some allocations may be extracted from multiple sites, particularly for larger Places of Use. ownerClassification - A classification of the types of owners for water rights. The most common is `Private` which incorporates a wide range of entities. Several classifications would be grouped into a government category, most of which are for the U.S. Federal Government. These allocations could be listed as "Federal", "United States of America", or as the names of any number of federal agencies. The last major grouping of entities is for "Native American"s. priorityDate - The date we use as the water right priority date for our modeling analysis. This is the legal priority date when it is available. However, for some rights, specifically from California and New Mexico, we used a pseudo priority date (e.g. well completion date or start of well drilling date) when a legal priority date was not available. The most questionable dates come from New Mexico, where the only date associated with certain water right records was the date the allocation was recorded in the database. As the allocation record creation tended to be within a few months of the filing of the application of the water right, from manually double checking the water rights, and our analysis focuses on aggregating water rights on the timescale of years, we determined it was acceptable to use such dates to include as many records as possible. primaryBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories WestDAAT. This column is the original WaDE category for the primary water use at the PoD site. allocationBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories for WestDAAT. This column is the original WaDE category

Economics↗
DOE OSTI2023

Implementation of a feature selection algorithm in FARM to identify important state variables and time-invariant matrices

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM aids the HERON software module in the evaluation of the optimal dispatch for the different IES components. Set-point trajectories are required to meet limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To evaluate the feasibility of HERON generated set-points and to do so in an acceptable time, FARM employs reduced order models to represent the dynamic behavior of the systems to be dispatched. These surrogate models take the form of a linear dynamic system with sets of Linear Parameter Varying (LPV) matrices that are mapped to the system operating space. These matrices are derived from the trajectories of system state variables and system output variables during transients. The accuracy of LPV matrices depends on the selection of state variables. In previous reports, state variables were selected by adopting a complicated workflow requiring multiple software licenses and an advanced level of user expertise. In this report, a new workflow that automates the state variable selection process is presented. It significantly reduces the frequency of user interventions and does not require multiple software licenses. Each module in the new workflow is described in detail, and the input / output examples in each step of the workflow are provided. It was demonstrated that this workflow can greatly reduce the complexity of the state variable selection process, and that the updated FARM-Gamma and FARM-Delta validators can benefit from this workflow when solving the power dispatch problem of a representative IES test case. Finally, some code improvements that can further enhance the efficiency are suggested.

42 ENGINEERING↗
DOE OSTI2025

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

The CXSFIT spectral fitting code: Past, present and future

Magnetically confined plasma experiments generate a wealth of spectroscopic data. The first step toward extracting physical parameters is to fit a spectral model to the often complex spectra. The CXSFIT (Charge eXchange Spectroscopy FITting) spectral fitting code was originally developed for fitting charge exchange spectra on JET from the late 1980s onward and has been further developed over decades to keep up with the needs of the users. The primary use is to efficiently fit a large number of spectra with many constrained Gaussian spectral lines of which the physical parameters can be coupled in a user-friendly manner. More recent additions to the code include time-dependent couplings between parameters, flexible background subtraction, and a non-linear coupling scheme between fit parameters. The latter was a pre-requisite for implementing Zeeman and motional Stark effect multiplets in the library of spectral features. The ability to save and replay “fit recipes,” even when multiple iterations are required, has ensured the traceability of the results and is one of the keys to the longevity and success of the code. The code is also in use on other tokamaks (AUG, ST-40) and to fit data from other spectroscopic diagnostics on JET. In this paper, we document the current capabilities and philosophy behind the structure of the code, including some of the algorithms used to calculate spectral features numerically efficiently. We also provide an outline of how CXSFIT could be transferred into a framework that would be able to meet the spectral fitting requirements of future devices, such as ITER.

Delabie, Ephrem G.↗
DOE OSTI2024

Model-Based Detection of Coordinated Attacks (DCA) in Distribution Systems

The fast-paced growth in digitization of smart grid components enhances system observability and remote-control capabilities through efficient communication. However, enhanced connectivity results in heightened system vulnerability towards cybersecurity risks in the cyber-physical power system. Coordinated cyber-attacks (CCA), when undetected, lead to system-wide impact in terms of large disturbances or widespread outages. Detecting CCA in the cyber layer is critical to thwart cyber-attacks in real-time before the attack impacts the physical system. The challenge of locating CCA stems from the complex grid dynamics, making it difficult to distinguish between normal operational variations and cyber-attack impact. CCA often employs multiple attack vectors targeting geographically distributed components, further complicating CCA identification. Existing research in intrusion detection is primarily focused on the transmission network and limited to detecting individual attacks. In this paper, a novel proactive DCA strategy is proposed for early detection of CCA by establishing correlations among distinct attack events through model-based reinforcement learning that utilizes abductive reasoning to conclude the attacker goal. The solution includes understanding the system model, learning the system dynamics, and correlating individual cyber-attacks to extract the attacker’s objective. The developed learning algorithm identifies the most probable attack path to reach the attacker’s objective by predicting the next attack steps. A DNP3-based cyber-physical co-simulation testbed is developed to test the proposed algorithm using the IEEE 13-node test feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Time-resolved momentum imaging of UV photodynamics in structural isomers of iodopropane probed by site-selective XUV ionization

The photodynamics of 1- and 2-iodopropane (1 and 2-IP) were studied in a time-resolved scheme incorporating ultraviolet (UV) excitation and extreme ultraviolet (XUV) probing, which initiates photoionization selectively from the I 4d core orbital. UV absorption in the A-band of both isomers leads to prompt C–I bond fission, with significant disposal of internal energy into the propyl radical product. Site-selective ionization enables a range of charge transfer (CT) processes between the nascent highly charged iodine ions and neutral propyl radicals, dependent on the interfragment distance at the instant of ionization. Subtle differences in the dynamics of these CT processes between the two isomers are observed. In 1-IP, the kinetic energies of iodine ions produced by UV photodissociation and subsequent XUV multiple ionization increased notably over the first few hundred femtoseconds, which could be understood in terms of differing gradients along the photodissociation coordinates of the neutral and polycationic states involved in the pump and probe steps, respectively. Led by a recent report of HI elimination in UV photoexcited 2-IP [Todt et al., Phys. Chem. Chem. Phys., 22(46), 27338 (2020)], we also model the most likely signatures of this process in the present experiment, and can identify signal in the 2-IP data (that is absent or significantly weaker in the data from the unbranched 1-IP isomer) that is consistent with such a process occurring on ultrafast timescales.

Allum, Felix [Deutsches Elektronen-Synchrotron (DE↗
DOE OSTI2024

Fiber loop quantum buffer for photonic qubits

Abstract We report a fiber loop quantum buffer based on a low-loss 2 × 2 switch and a unit delay made of a fiber delay line. We characterize the device by using a two-photon polarization entangled state in which one photon of the entangled photon pair is stored and retrieved at a repetition rate up to 78 kHz. The device, which enables integer multiples of a unit delay, can store the qubit state in a unit of fiber delay line up to 5.4 km and the number of loop round-trips up to 3. Furthermore, we configure the device with other active elements to realize integer multiplier and divider of a unit delay of a qubit. The quantum state tomography is performed on the retrieved photon and its entangled photon. We obtain a state fidelity > 94 % with a maximum storage time of 52 s with an insertion loss of 5.56 dB. To further characterize the storing and retrieving processes of the device, we perform entanglement-assisted quantum process tomography on the buffered qubit state. The process fidelity of the device is > 0.98. Our result implies that the device preserves the superposition and entanglement of a qubit state from a two-photon polarization-entangled state. This is a significant step towards facilitating applications in optical asynchronous transfer mode based quantum networks.

Fook Lee, Kim (ORCID:0000000290901640)↗
DOE OSTI2021

Critical and Subcritical Californium Source Driven Noise Analysis Experiments with Fresh PWR Fuels Pins

These experiments were performed in 1983 between August 1 and September 16 and used 35 operational days of critical facility time at the Babcock & Wilcox critical facility in Lynchburg, Va. These measurements were to assess the capability of the Cf source driven noise analysis method to measure the subcriticality of light water reactor fresh fuel configurations. This report documents the experiments of 1983 that were not reported at that time. These measurement with 2.459 wt. % 235 U enriched uranium oxide fresh fuel pins started at the delayed critical configuration of the fully assembled pressurized water reactor (PWR) core configuration with 1511 parts per million (ppm) natural boron (19.8 wt. % 10 B) in the water. This was a cylindrical configuration of 4961 fuel pins. After the initial measurements at delayed criticality, subcritical measurements by the Cf source driven noise analysis method were performed in which the number of detectors and Cf source locations were varied. Additional boron was added to the water moderator and reflector to reduce the neutron multiplication factor, with resulting boron concentrations of 1561, 1613, 1765, 1880, 2104, 2384, 2975, 3606, and 4303 ppm. After the completion of these measurements the increased, the boron concentration was returned to that at delayed critical (1511 ppm) and the reactivity was decreased by reducing the number of fuel pins in steps to 3713, 2553, 1281, 749, and 333 fuel pins and then to that (289) of a s ingle square 17 x 17 pin array, corresponding to a single PW fuel element. For all measurement that started at delayed criticality and then reduced the neutron multiplication factor to various subcritical states, inverse kinetics rod drop measurements at the higher reactivities were also performed but not reported here. At all subcritical states the break frequency noise analysis data was fitted to obtain the prompt neutron decay constant. Some data presented in this report are from notes not in the logbook. The purpose of this report is to document the experimental information for the measurements performed so that a later date someone could perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Program (INCSBEP) or Nuclear Energy Agency (NEA) benchmark. The data from these measurements is available from the Laboratory Records Department of Oak Ridge National Laboratory and the logbook is available from INCSBEP at Idaho National Laboratory.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗