Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lessons Learned Using Iridium to Communicate with a CubeSat in Low Earth Orbit

This paper presents the design and approval process for operating an Iridium transceiver on orbit and provide on-orbit performance data obtained from a CubeSat platform in Low Earth Orbit (LEO) (500 km orbit). On-orbit data demonstrates that use of a commercial, low-cost Iridium transceiver can serve as a valuable communication approach for low volume telemetry with less than a 30-minute lag for approximately 90% of the time. In this work, we also demonstrate that a radial differential velocity of 7 km/sec corresponding to about a 37.5kHz doppler shift and a distance of less than 2,000 km can be used for mission planning.

42 ENGINEERING↗
DOE OSTI2021

Investigation of Main Bearing Operating Conditions in a Three-Point Mount Wind Turbine Drivetrain

Wear related failures of spherical roller bearings in the main bearing position of three-point mount wind turbines have been higher than expected and can contribute to higher than anticipated operation-and-maintenance costs. In this paper, the operational conditions of such a main bearing, including measured axial displacement and velocity subject to the estimated axial loads are described for an instrumented, commercial wind turbine. The field measurements suggest a maximum axial speed between the bearing rings that is less than 2 millimeters per second. It is estimated that the axial speed between the ring and rollers is approximately 25% of this value. When compared to a speed of rolling from approximately 200-352 millimeters per second, the measured axial sliding is therefore significantly less than 1% of the speed of rolling. Previous numerical studies of lubricant film formation in rolling contacts have shown that the effect of axial sliding starts to be noticeable only when this ratio exceeds 10%; therefore, the axial velocity represents only a small disturbance to the nominal pure rolling case, and the influence on oil film building can be neglected. A simple analytic model of the main bearing motion was also developed and demonstrated similar displacement and velocity characteristics.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2022

X-ray Computed Tomography of Irradiated and Unirradiated AGR-3/4 Compacts

X-ray Computed Tomography (XCT) has been utilized to image and characterize compacts from the combined third and fourth irradiation of the Advanced Gas Reactor (AGR) Program, AGR-3/4, fuel. The experiment contained tristructural isotropic (TRISO)-coated fuel particles as well as designed-to-fail (DTF) fuel particles. Two irradiated compacts, representing the lower and higher range of AGR-3/4 burnup (4.85% and 14.92% fissions per initial heavy metal atom FIMA) were examined. These represent the first known highly irradiated TRISO fuel compacts to be examined via X-ray CT. Additionally, two unirradiated compacts from the same production batch as the examined irradiation compacts were also imaged for a baseline comparison. As XCT of irradiated TRISO compacts is not a commonly implemented characterization technique, a significant portion of the report focuses on developed methodology and imaging conditions. A specialized sample shielding device was developed and fabricated specifically to limit received dose to staff during sample preparation for XCT and to minimize excess gamma radiation dose to sensitive electronic components with the utilized X-ray system. Significant penetration through the uranium oxycarbide fuel kernels by significantly hardening the X-ray beam with specialized proprietary filters acquired from Carl Zeiss NTS Ltd. The filter utilized resulted in an average X ray photon energy of ~110 keV which approaches uranium’s K-edge (~115 keV), maximizing penetration for a microfocus X-ray source. The gamma-radiation emitted from the irradiated AGR-3/4 TRISO compacts, has the same properties and mechanisms for interaction with matter as X-rays, thus the detection of gamma-radiation by the utilized X-ray detectors was initially a concern. However, although ?-rays did produce an observable signal on the X-ray detector, its contribution to the overall imaging results appeared negligible upon 3D reconstruction. The neglibile impact on the resulting 3D reconstructed volumes were likely the result of: (1) a significantly lower detection efficiency for ?-rays relative to X-rays; (2) An X-ray flux at the detector several orders of magnitude higher than that of the impinging ?-rays from the irradiated compacts. These results suggest that irradiated compacts with significantly higher radiation fields can be examined in the future if an acceptable route for sample handling and preparation can be determined. Additionally, the 3D imaging results of XCT can provide a valuable means of assessing compacts. While in many ways complimentary to traditional post irradiation examination techniques such as optical ceramography, XCT can provide additional insight into compact features traditionally difficult to discern directly from cross-sectional imaging alone. Preliminary analyses on kernel size, morphology (aspect ratio and sphericity), and kernel orientation were presented. Sphericity, a simple morphological shape descriptor, was utilized to screen for kernel extrusions within the high burnup compact. The number of kernel extrusions identified via XCT represented an approximate two-fold increase from the quantity of extruded particles observed (via optical ceramography) in adjacent compacts from the same irradiation capsule. While numerical analysis of the compact datasets was highly preliminary, initial results show promise for providing complimentary metrics to current AGR-3/4 PIE and potentially additional insight into the processes driving TRISO fuel degradation during reactor operation. Additional analyses to be performed at a later date include a more detailed examination of kernel size, kernel sphericity (and observed kernel extrusions), and sphericity. Given all particles can be observed in a single data volume possible correlation of spatial position with observed kernel features will also be made at a later date.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Non-isometry, state dependence and holography

We establish an equivalence between non-isometry of quantum codes and state dependence of operator reconstruction, and discuss implications of this equivalence for holographic duality. Specifically, we define quantitative measures of non-isometry and state dependence and describe bounds relating these quantities. In the context of holography we show that, assuming known gravitational path integral results for overlaps between semiclassical states, non-isometric bulk-to-boundary maps with a trivial kernel are approximately isometric and bulk reconstruction approximately state-independent. In contrast, non-isometric maps with a non-empty kernel always lead to state-dependent reconstruction. We also show that if a global bulk-to-boundary map is non-isometric, then there exists a region in the bulk which is causally disconnected from the boundary. Finally, we conjecture that, under certain physical assumptions for the definition of the Hilbert space of effective field theory in AdS space, the presence of a global horizon implies a non-isometric global bulk-to-boundary map.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Effects of knock intensity measurement technique and fuel chemical composition on the research octane number (RON) of FACE gasolines: Part 1 – Lambda and knock characterization

The Research and Motor Octane Number (RON and MON) rate the knock propensity of gasoline in the Cooperative Fuel Research (CFR) engine by comparing the knock intensity of sample fuels relative to that of primary reference fuels (PRF), a binary blend of iso-octane and n-heptane. Important differences exist between standard octane testing and automotive spark ignition (SI) engine knock testing including speed, load, air-to-fuel equivalence ratio (lambda), and knock characterization, which lead to a discrepancy between a fuel’s RON rating and its knock resistance characterized on an automotive SI engine based on knock-limited spark advance. This publication is the first of a set of three publications which modify operating parameters of the RON test method (ASTM D2699) to investigate the effects of these differences with automotive SI engine knock-limited spark advance testing. A fuel’s standard RON is evaluated at the lambda of the highest knock intensity, whereas automotive SI engines typically operate at stoichiometry. Here, we analyze the effects of a stoichiometric RON rating methodology. Additionally, the knock intensity response from the standard CFR knockmeter system is compared to a cylinder pressure oscillation-based knock intensity at several lambda settings. All experiments were performed with a set of seven Coordinating Research Council (CRC) Fuels for Advanced Combustion Engines (FACE) gasolines with approximately 95 RON. The fuel chemical composition impacted the lambda of the highest knock intensity, which resulted in fuel-specific offsets between the standard and stoichiometric RON ratings. The knock system comparison showed significant offsets between cylinder pressure-based and knockmeter-based knock intensity levels.

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

Mechanical characterization of fine-grain dispersion-strengthened tungsten as a plasma facing material

Field-Assisted Sintering Technology (FAST) was used to produce fine-grained, dispersion-strengthened tungsten (W) materials. Investigated materials 4138, 4353, and 4355 composed of 3 wt% ZrC sintered at 1800 °C, 5 wt% ZrC sintered at 1800 °C, and 3 wt% ZrC sintered at 2000 °C, respectively. They were compared against ITER-grade W. A series of mechanical and thermal property testing and microstructure studies were conducted to study them as a potential plasma facing material (PFM) for fusion reactors. Hardness testing showed that manufacturing conditions substantially altered hardness. Material 4355 had an average HV10 value of 497.2 ± 16.8, slightly higher than ITER-grade at 378.5 ± 40.3. However, material 4353 was substantially higher with an HV10 value of 738.9 ± 31.7 over the investigated temperature range. Electron Backscatter Diffraction (EBSD) analysis showed that FAST produced substantially smaller grains than the hot-rolled ITER-grade W material, offering notable control over grain size. Materials 4353 and 4355 had grain sizes of 0.44 ± 0.20 µm and 3.67 ± 0.89 µm, respectively, whereas ITER-grade 27.14 ± 19.76 µm at room temperature. The fine grain structures showed no net coarsening after 1 hr. anneals up to 1800 °C, several hundred degrees above the 1100 – 1500 °C recrystallization range reported for conventional W. Inverse application of the Zener pinning relationship to the measured grain sizes indicates that these two FAST sintering conditions produce markedly different effective dispersoid populations, with effective particle diameters of approximately 90 nm at a peak sintering temperature of 1800 °C and approximately 460 nm at 2000 °C, respectively. This result demonstrates that the FAST thermal condition itself, and not the nominal ZrC content alone, governs the pinning effectiveness of the dispersion. Thermal diffusivity measurements support this finding independently. Materials of identical composition sintered at different temperatures differ by approximately 19% in measured thermal diffusivity with statistically indistinguishable density, while materials of different composition and sintering temperature converge to within approximately 2%. At a representative divertor heat flux of 10 MW/m², the lower thermal conductivity of the fine-grained materials corresponds to approximately 28 to 33 °C per millimeter of armor thickness relative to ITER-grade W, traded against a substantially larger margin to recrystallization-driven degradation. While high temperature tensile testing revealed likely contamination that motivates refinement of the manufacturing process, FAST-produced, fine-grained, dispersion-strengthened W offers process-controlled microstructural stability well above the operating temperatures of conventional W and supports its continued development as a PFM for economically viable commercial fusion power.

Parker, Gabe [ORNL] (ORCID:0000000190394100)↗
DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
DOE OSTI2021

The environmental footprint of data centers in the United States

Much of the world's data are stored, managed, and distributed by data centers. Data centers require a tremendous amount of energy to operate, accounting for around 1.8% of electricity use in the United States. Large amounts of water are also required to operate data centers, both directly for liquid cooling and indirectly to produce electricity. For the first time, we calculate spatially-detailed carbon and water footprints of data centers operating within the United States, which is home to around one-quarter of all data center servers globally. Our bottom-up approach reveals one-fifth of data center servers direct water footprint comes from moderately to highly water stressed watersheds, while nearly half of servers are fully or partially powered by power plants located within water stressed regions. Approximately 0.5% of total US greenhouse gas emissions are attributed to data centers. We investigate tradeoffs and synergies between data center's water and energy utilization by strategically locating data centers in areas of the country that will minimize one or more environmental footprints. Our study quantifies the environmental implications behind our data creation and storage and shows a path to decrease the environmental footprint of our increasing digital footprint.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Multi-scale impacts of climate change on hydropower for long-term water-energy planning in the contiguous United States

Climate change impacts on watersheds can potentially exacerbate water scarcity issues where water serves multiple purposes including hydropower. The long-term management of water and energy resources is still mostly approached in a siloed manner at different basins or watersheds, failing to consider the potential impacts that may concurrently affect many regions at once. There is a need for a large-scale hydropower modeling framework that can examine climate impacts across adjoining river basins and balancing authorities (BAs) and provide a periodic assessment at regional to national scales. Expanding from our prior assessment only for the United States (US) federal hydropower plants, we enhance and extend two regional hydropower models to cover over 85% of the total hydropower nameplate capacity and present the first contiguous US-wide assessment of future hydropower production under Coupled Model Intercomparison Project phase 6’s high-end Shared Socioeconomic Pathway 5-8.5 emission scenario using an uncertainty-aware multi-model ensemble approach. We present regional hydropower projections, using both BA regions and US Hydrologic Subregions (HUC4s), to consistently inform the energy and water communities for two future periods—the near-term (2020–2039) and the mid-term (2040–2059) relative to a historical baseline period (1980–2019). We find that the median projected changes in annual hydropower generation are typically positive—approximately 5% in the near-term, and 10% in the mid-term. However, since the risk of regional droughts is also projected to increase, future planning cannot overly rely on the ensemble median, as the potential of severe hydropower reductions could be overlooked. The assessment offers an ensemble of future hydropower generation projections, providing regional utilities and power system operators with consistent data to develop drought scenarios, design long duration storage and evaluate energy infrastructure reliability under intensified inter-annual and seasonal variability.

13 HYDRO ENERGY↗
DOE OSTI2021

Success of digital adiabatic simulation with large Trotter step

The simulation of adiabatic evolution has deep connections with adiabatic quantum computation, the quantum approximate optimization algorithm, and adiabatic state preparation. Here we address the error analysis problem in quantum simulation of adiabatic process using Trotter formulas. Here we show that with additional conditions, the circuit depth can be linear in simulation time T. The improvement comes from the observation that the fidelity error here can't be estimated by the norm distance between evolution operators. This phenomenon is termed the robustness of discretization in digital adiabatic simulation. It can be explained in three steps, from analytical and numerical evidence: (1) The fidelity error should be estimated by applying adiabatic theorem on the effective Hamiltonian instead. (2) Because of the specialty of Riemann-Lebesgue lemma, most adiabatic process is naturally robust against discretization. (3) As the Trotter step gets larger, the spectral gap of effective Hamiltonian tends to close, which results in the failure of digital adiabatic simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Distributed Software-Defined Network Architecture for Smart Grid Resilience to Denial-of-Service Attacks

An important challenge for smart grid security is designing a secure and robust smart grid communications architecture to protect against cyber-threats, such as Denial-of-Service (DoS) attacks, that can adversely impact the operation of the power grid. Researchers have proposed using Software Defined Network frameworks to enhance cybersecurity of the smart grid, but there is a lack of benchmarking and comparative analyses among the many techniques. In this work, a distributed three-controller software-defined networking (D3-SDN) architecture, benchmarking, and comparative analysis with other techniques is presented. The selected distributed flat SDN architecture divides the network horizontally into multiple areas or clusters, where each cluster is handled by a single Open Network Operating System (ONOS) controller. A case study using the IEEE 118-bus system is provided to compare the performance of the presented ONOS-managed D3-SDN, against the POX controller. In addition, the proposed architecture outperforms a single SDN controller framework by a tenfold increase in throughput; a reduction in latency of > 20%; and an increase in throughput of approximately 11% during the DoS attack scenarios.

Agnew Jr., Dennis↗
DOE OSTI2025

Hydraulic Testing of the Lost Lake Aquifer Near Recovery Well RWM001R

An aquifer pumping test was conducted on the Lost Lake Aquifer Zone (LLAZ) at the recently installed recovery well RWM001R in accordance with the approved test plan (Dixon, 2024). The objective of the testing was to determine baseline well performance parameters and aquifer hydraulic conductivity. This testing consisted of two step-drawdown tests to determine well performance properties and two constant pumping rate aquifer tests with recovery periods to determine aquifer hydraulic properties. Well performance parameters determined included specific capacity, well efficiency, and head loss coefficients. The results of the step-drawdown testing are presented in Tables ES1 and ES2 whereas the results of the constant pumping rate testing are presented in Table ES3. The average specific capacity of RWM001R operating at pump capacity (~54 gpm) was 2.2 gpm/ft using data from the final step of both tests. The average well efficiency determined from both step-drawdown tests was 56% using data from the final step of both tests. For both step-drawdown tests, well efficiency was inversely related to the pumping rate with the average efficiency decreasing from approximately 79% to 56%. Aquifer head loss coefficients from both tests were comparable (Table ES2).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Deriving cloud droplet number concentration from surface-based remote sensors with an emphasis on lidar measurements

Abstract. Given the importance of constraining cloud droplet number concentrations (Nd) in low-level clouds, we explore two methods for retrieving Nd from surface-based remote sensing that emphasize the information content in lidar measurements. Because Nd is the zeroth moment of the droplet size distribution (DSD), and all remote sensing approaches respond to DSD moments that are at least 2 orders of magnitude greater than the zeroth moment, deriving Nd from remote sensing measurements has significant uncertainty. At minimum, such algorithms require the extrapolation of information from two other measurements that respond to different moments of the DSD. Lidar, for instance, is sensitive to the second moment (cross-sectional area) of the DSD, while other measures from microwave sensors respond to higher-order moments. We develop methods using a simple lidar forward model that demonstrates that the depth to the maximum in lidar-attenuated backscatter (Rmax⁡) is strongly sensitive to Nd when some measure of the liquid water content vertical profile is given or assumed. Knowledge of Rmax⁡ to within 5 m can constrain Nd to within several tens of percent. However, operational lidar networks provide vertical resolutions of > 15 m, making a direct calculation of Nd from Rmax⁡ very uncertain. Therefore, we develop a Bayesian optimal estimation algorithm that brings additional information to the inversion such as lidar-derived extinction and radar reflectivity near the cloud top. This statistical approach provides reasonable characterizations of Nd and effective radius (re) to within approximately a factor of 2 and 30 %, respectively. By comparing surface-derived cloud properties with MODIS satellite and aircraft data collected during the MARCUS and CAPRICORN II campaigns, we demonstrate the utility of the methodology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Measurement of Fission Product Concentration Profiles in AGR-3/4 TRISO Fuel Graphitic Matrix and Nuclear Graphites

The third Advanced Gas Reactor (AGR) irradiation experiment, AGR-3/4, was designed to investigate the migration of fission products in fuel compact graphitic matrix and reactor graphite components. Using destructive methods, radial fission product concentration profiles were measured for gamma-emitting fission products (e.g., Ag-110m, Cs-134, and Eu-154) and beta-emitting Sr-90 in irradiated graphitic and graphite components from six different AGR-3/4 irradiation capsules. These new measured concentration profiles can now be compared to non-destructive measurements and fission product transport simulations and will be used to derive new diffusivities and sorptivities to support refinement of fission product transport models and high-temperature gas-cooled reactor (HTGR) source-term analyses. Each capsule in the AGR-3/4 experiment had four fuel compacts in the middle of two concentric rings of graphitic matrix material, PCEA graphite, or IG-110 graphite. In addition to the approximately 1898 tristructural isotropic (TRISO) coated particles in each compact, there were 20 designed to fail (DTF) particles coated only in pyrocarbon so that they released fission products into the surrounding cylindrical rings of carbonaceous materials. Destructive sampling of the rings involved machining/milling material from around the circumference of the rings, collecting that material, and performing radiochemical analyses on it. Milling operations were performed in multiple steps or segments, and each segment was generally 0.508 mm (0.020 in) thick. Knowing the radial position at which each segment was milled, the volume of the milled material at each segment, and the fission product content in each segment, the radial fission product concentrations were constructed for select isotopes in each ring. Ag-110m profiles had the most variation. Some profiles were peaked at an inner or outer surface. Some were peaked at the middle of the ring wall thickness. Some increased radially outward, and some decreased radially outward. These types of variations and the fact that the measured profiles do not generally compare favorably with the transport model employed for AGR-3/4 may adversely impact the ability to extract reasonable transport parameters for this isotope. In many cases, the Cs-134 profiles decreased somewhat linearly in the outward radial direction, and in cursory comparisons, the shapes of these profiles appeared similar to those from model predictions. The step changes in concentration across the inner-outer ring gap were generally consistent with the model predictions as well. In some cases, there were local maxima in concentration at the outer surface of the rings. This suggests that fission products could have transported in the small gaps between the inner ring and the outer ring and between the outer ring and the sink ring such that some portion of a given fission product can bypass diffusion through the ring itself. The analysis of the small nubs on the outer surfaces of some of the outer rings revealed fission product concentrations in the nubs that were often higher than in the outermost segments of the rings. This further supports the hypothesis that short-circuit, gap transport occurred, causing relatively high surface concentrations on the outer surfaces of the rings. Eu-154 and Sr-90 profiles tend to have very similar shapes, suggesting that they transport via the same mechanisms. The observed profiles were indicative of a transport process where the isotopes are sorbed on the inner surface of the ring, but diffusion into the ring from that surface is quite slow. Some elevated concentrations of Sr-90 (relative to Eu-154) on the outer surface of a ring suggested that rapid, gas-gap transport of gaseous precursor Kr-90 and volatile Rb-90 could have occurred prior to their decaying to Sr-90. Overall, the Eu 154 and Sr-90 profiles were still very similar, which indicates that the transport of short-lived Sr-90 precursors is not a major effect. In some capsules, the qualitative Sr-90 behavior across the ring gaps was consistent with the model (using the available legacy Sr-90 transport parameters), but in other capsules the model was inconsistent with the measurements and seems to underestimate the amount of Sr-90 in the outer rings. The total ring Sr-90 inventories were estimated for all the rings that were subject to physical sampling. These results will be used to adjust the predicted particle and/or compact releases used in the AGR-3/4 fission product transport model. Given the different irradiation temperatures among the capsules and the rings, it was not possible to discern fundamental differences in the transport of isotopes within the different carbon materials, i.e., graphitic matrix, IG-110, or PCEA. It may be possible to do this in the course of determining transport from the concentration profiles in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Spectral deferred correction methods for high-order accuracy in poroelastic problems

In this work, we investigate high-order accuracy in time integration by examining two operator splitting methods for poroelastic problems: the two-pass and the spectral deferred correction (SDC) methods. To enhance the order of accuracy, the two-pass method partitions a coupled operator symmetrically, whereas the SDC method corrects truncation errors by establishing an error equation. These high-order methods are applied to underlying solution strategies, i.e., monolithic, fixed-stress sequential, and undrained sequential methods. We observe that semi-discretized systems from spatial discretization have forms similar to those of index-1 differential algebraic equations (DAEs), causing order reduction against the two-pass method when it is used in conjunction with either the monolithic or sequential method. On the other hand, the SDC in conjunction with the monolithic method exhibits the desired second-order accuracy in poroelastic problems while increasing the order of accuracy for index-1 DAEs. However, the SDC in conjunction with either of the two sequential methods does not achieve the desired order of accuracy, and maintains first order because the flow equation for poroelasticity has an additional approximation associated with the volumetric strain rate term, which does not yield exactly the same forms as those of conventional DAEs. Thus, the monolithic SDC method can achieve higher-order accuracy, but may require higher computational costs because it involves solving matrix systems larger than those for the sequential methods.

02 PETROLEUM↗
DOE OSTI2023

Demonstration of Scaled-Production of Rare Earth Oxides and Critical Materials from U. S. Coal-Based Sources (Final Report)

The project objective was to demonstrate scaled production of high purity rare earth oxides (REO), nominally exceeding 90% grade, from coal refuse sources using innovative technologies that reduce cost and improve environmental outcomes relative to traditional rare earth processing technologies. The project utilized a critical material pilot plant constructed and tested as part of a previous U.S. Department of Energy project. Target performance criteria was a 50% reduction in production cost based on previous optimum values, 150% increase in recovery and greater than 2% concentrates of rare earth oxides, cobalt and manganese. Concentrate production goals were to produce a rare earth mix oxide product at a rate of 200 grams per day having a minimum purity of 50% as well as products of cobalt and manganese having a minimum purity of 2%. A previous pilot plant investigation identified acid cost as the major contributor to an operating cost that made the recovery of rare earth and other critical metals from bituminous coal sources economically challenging. As such, acid cost reduction was a major target using bio-oxidation reactors to produce sulfuric acid from naturally occurring coal pyrite. Based on laboratory data, a bio-oxidation circuit was designed for the pilot plant to produce 7.5 l/min of acid using two 11-m3 (3000 gallon) reactors equipped with 40 hp aerators for air dispersion. The pilot-scale tests revealed that acid concentration equivalent to as high as 1.0 M sulfuric acid could be continuously produced. However, the bio-acid contained exceptionally high iron concentrations, which complicated downstream processing of the pregnant leach solution (PLS). A TEA of the bio-oxidation circuit showed that production cost was approximately $0.13 per kg acid equivalent if produced using a four-day retention time in the reactors. This value represents a 48% reduction from that of purchased bulk sulfuric acid. Calcination (or roasting) studies were conducted on coarse refuse from West Kentucky No. 13 and Fire Clay coal seam sources. The test results revealed the potential to increase recovery by nearly 100% using temperatures between 500°C to 700°C with light REE recovery value being the most improved. Acid baking of the calcined products using sulfuric acid at 250°C increased heavy rare earth recovery from around 40% to 80% while decreasing the acid requirements by over 50%. The existing pilot plant was upgraded for the pilot scale demonstration of REE and CM recovery. The primary feedstocks were West Kentucky No.13 and heap leach pregnant leach solution (PLS) while a secondary feedstock was a lignite waste material from a construction sand operation. The pilot scale operation started with PLS generation through leaching followed by iron and aluminum removal, nominally at 3.3 and 4.5 pH, respectively. Leaching lixiviants used for the test were industrial grade sulfuric acid or bio-acid generated at the pilot scale facility. For most of the tests, the solid feed rate was 200 lb/hr whereas lixiviant was added at 2 gpm to provide an optimal residence time of 45 minutes. Following the contaminant removal step, several different process schematics were tested with the goal of maximizing REE recovery and purity. In the first test, direct processing of aluminum precipitation raffinate for REE recovery using oxalic acid at pH 1.5 was investigated. Overall REE recovery was approximately 45%. Unfortunately, elevated calcium content in the PLS significantly impacted the RE-Oxide product grade. Similarly, high calcium content decreased both the product purity and grades of CM products. As such, a new flowsheet was tested with the same initial process schematic but different precipitation stages for REEs and CMs at pH 6.0 and 9.0, respectively. It was noted that the overlapping precipitation behavior of Co, Ni and Zn with REEs limited the applicability of this process schematic. While this change increased the RE-Oxide grade from 36% in the first test to 87%, the loss of critical metals to the REE cake and bypass of the REEs to the CM cake significantly impacted the recovery of both REEs and CMs. Therefore, the modified process flowsheet combined oxalic acid precipitation stage raffinate and redissolved CM cake filtrate to maximize both the recovery and purity of the products. Consequently, a RE-Oxide product with 85% purity and CM cakes with over 19% Co, 38% Ni, 14% Zn and 9% Mn content were generated with significantly higher recoveries. While the modified process flowsheet improved recoveries and grades, elemental losses observed in separate precipitation and redissolution losses inspired the adaptation of a single precipitation stage at pH 9.0 for both REEs and CMs. This change was anticipated to maximize the REE recovery while minimizing the costs associated with separated redissolution and processing stages. As expected, REE recovery in this new circuit arrangement was over 56% with a product grade of over 87% RE-Oxide content. Similarly, Co, Ni, Mn, and Zn recoveries of 54%, 40%, 67%, and 66%, respectively, were achieved. Unfortunately, the elevated calcium content present in the solution due to its precipitation at pH 9.0 caused a decrease in the CM cake quality. Therefore, the final process flowsheet involved the addition of calcium oxalate precipitation following the oxalic acid precipitation stage, which effectively eliminated calcium contamination of the CM products. Finally, the pilot scale experiments conducted using bio-acid achieved comparable REE leaching recoveries to the conventional sulfuric acid leaching. Elevated iron concentration in the solution caused the co-precipitation of REEs with the iron and aluminum cake, resulting in the REE losses. Furthermore, elevated iron content bypassing the iron and aluminum precipitation stages contaminated the metal sulfide and manganese cake, respectively. A techno-economic analysis was performed based on a commercial facility capable of treating 500 tph of coal-based material. The production cost for West Kentucky No. 13 coarse refuse material ranged from approximately $500-$700/kg of total rare earth oxide whereas the lignite source had significantly lower production cost of $100-$300/kg. The significant difference in cost was due to the easier leaching characteristics of the lignite material and the higher feed concentrations. All process scenarios resulted in a negative net present value (NPV). For the lignite feedstock, laboratory REE leach recovery values were about 30% higher than the pilot plant data. Using the lab leach results, a positive net present value was achieved and the production cost decreased from $100-$300 $/kg to less than $150/kg of total REO.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Precursor Analysis Report: SQL Slammer Worm Infection of Davis-Besse Nuclear Power Plant 2003

The SQL Slammer Worm Infection of Davis-Besse Nuclear Power Plant 2003 Precursor Analysis Report leverages publicly available information about Davis-Besse’s 2003 cyber attack and catalogs anomalous observables for each technique employed in the attack. This analysis is based upon the methodology of the Cybersecurity for the Operational Technology Environment (CyOTE) program. On 25 January 2003, the SQL Slammer worm infected more than 90% of vulnerable hosts and crashed the internet in 10 to 15 minutes, making it one of the fastest spreading worms in history. SQL Slammer is a fileless, memory-resident worm that remotely exploits a stack-based buffer overflow vulnerability on local hosts to intensively scan and rapidly self-propagate across the internet. The worm infected approximately 300,000 unpatched hosts running Microsoft Structured Query Language (SQL) Server 2000 or Microsoft Desktop Engine (MSDE) 2000 with SQL Server Resolution Service. The SQL Slammer worm indirectly infected FirstEnergy’s Davis-Besse nuclear power plant by first infecting a consultant’s company network server and then propagating through an external misconfigured connection into Davis-Besse’s site network. The infection caused major network congestion, slow performance, data overloads, and the inability of local hosts to communicate with each other, which eventually caused a loss of availability and a loss of view when the Safety Parameter Display System (SPDS) and Plant Process Computer (PPC) crashed. At the time of the infection, the plant was already offline, the digital monitoring systems had redundant analog backups, and the plant control and safety functions were not affected, so there were no concerns of a safety breach. However, this incident resulted in many lessons learned and spawned important discussions about cybersecurity’s role in nuclear safety and electric power reliability regulation, policy, and guidance. Researchers and analysts identified 10 unique techniques utilized during the attack with a total of 640 observables using MITRE ATT&CK® for Industrial Control Systems. The CyOTE program assesses observables accompanying techniques used prior to the triggering event to identify opportunities to detect malicious activity. If observables accompanying the attack techniques are perceived and investigated prior to the triggering event, earlier comprehension of malicious activity can take place. Eight of the identified techniques used during Davis-Besse cyber attack were precursors to the triggering event. Analysis identified 596 observables associated with these precursor techniques, 428 of which were assessed to have an increased likelihood of being perceived in the 331 days preceding the triggering event. The response and comprehension time could have been reduced if the observables had been identified earlier. The information gathered in this report contributes to a library of observables tied to a repository of artifacts, data sources, and technique detection references for practitioners and developers to support the comprehension of indicators of attack. Asset owners and operators can use these products if they experience similar observables or to prepare for comparable scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗