Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “filter methods”

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

Fully Homomorphic Encryption

This code implements a Fully Homomorphic Encryption (FHE) system, enabling secure computation on encrypted data without requiring decryption. It supports encryption, decryption, and homomorphic operations like matrix multiplication and addition. This code is adaptable for integrating FHE into linear-time invariant (LTI) systems, including digital control and filtering. With proper configuration from subject matter expertise, encrypted system parameters and signals can be manipulated to perform tasks like state updates, output calculations, and convolution in the encrypted domain. By preserving the structure of LTI systems while ensuring privacy, the framework facilitates secure applications in areas such as autonomous systems, signal processing, and industrial automation. The code initializes the encryption system using parameters provided in the env dictionary. These parameters include the ciphertext modulus, key dimension, plaintext fixed-point scaling factor, and noise bound. During initialization, a secret key is generated, which is essential for encrypting and decrypting data securely. The modular design allows users to tailor these parameters to specific use cases or security requirements. The code implements multiple cryptographic schemes. The learning with errors (LWE) encryption method encodes cleartext message to their plaintext fixed-point representation then encrypted into ciphertext space with additive noise. This noise ensures the security of the scheme, relying on the computational hardness of the LWE problem. The code also includes the Gentry-Sahai-Waters (GSW) scheme based off the LWE problem. Homomorphic matrix multiplication is performed between the LWE and GSW to encrypted data. This is achieved using a decomposition function on the LWE ciphertext during the multiplication operation. For higher-dimensional data, the code includes a method to encrypt entire matrices (GSWMat) using GSW encryption. These encrypted matrices can then be used for homomorphic matrix multiplications (MatMult). The decryption function uses the secret key to recover the original plaintext, removing the added noise and scaling that was originally applied during encryption.

Lois, Roberts [Idaho National Laboratory (INL), Id↗
DOE OSTI2021

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗
DOE OSTI2026

Fine-Grained Power and Energy Attribution on AMD GPU/APU-Based Exascale Nodes

Modern exascale GPU- and APU-based systems provide multiple power and energy sensors, but differences in scope, update rate, timing, and filtering complicate the attribution of short-lived accelerator activity. This paper presents a methodology to characterize and correct these effects on Cray EX systems with AMD Instinct MI250X GPUs (Frontier) and MI300A APUs (Portage). Using controlled square-wave workloads, we quantify update intervals, delay, aliasing, and variability across up to 512 GPUs and 480 APUs with on-chip (rocm-smi/amd-smi) and off-chip Cray Power Management sensors. We reconstruct power from cumulative energy counters to achieve faster response times, validate it against on-chip, off-chip, and node-level sensors, and integrate the resulting streams into a Score-P/PAPI-based tool for time-aligned, phase-level attribution. Applied to rocHPL, rocHPL-MxP, and HPG-MxP, the method separates energy savings due to reduced runtime from changes in power. Mixed precision reduces node energy on Frontier by 79% for rocHPL-MxP and 31% for HPG-MxP, with similar trends on Portage. These results provide portable guidance for sensor validation and power-aware optimization on current and future exascale systems.

Mcdaniel, Adam [ORNL] (ORCID:000000016926028X)↗
DOE OSTI2021

Electronic demodulation of optically phase demodulated signals

Device, methods and systems for the electronic demodulation of optically phase demodulated signals are described. An example optical local oscillator generator configured to generate a radio frequency (RF) tone at a desired RF frequency includes a first input configured to receive a broadband optical pulse train, a second input coupled to a delay line interferometer to receive a first control voltage for controlling a delay value of the interferometer and to produce an output optical pulse train, a dispersive element, coupled to the delay line interferometer, to map the output optical pulse train to a time-domain modulated optical pulse train, an optical-to-electrical converter, coupled to the dispersive element, to convert the time-domain modulated optical pulse train to an analog electrical signal, and an RF filter, coupled to the optical-to-electrical converter, to filter the analog electrical signal to generate the RF tone at the desired RF frequency.

Perlmutter, David Simon↗
DOE OSTI2023

Effects of turbulent diffusion and back-reaction on the dust distribution around two resonant planets

In evolved and dusty circumstellar discs, two planets with masses comparable to Jupiter and Saturn that migrate outwards while maintaining an orbital resonance can produce distinctive features in the dust distribution. Dust accumulates at the outer edge of the common gas gap, which behaves as a dust trap, where the local dust concentration is significantly enhanced by the planets’ outward motion. Concurrently, an expanding cavity forms in the dust distribution inside the planets’ orbits, because dust does not filter through the common gaseous gap and grain depletion in the region continues via inward drifting. There is no cavity in the gas distribution because gas can filter through the gap, although ongoing gas accretion on the planets can reduce the gas density in the inner disc. Such behaviour was demonstrated by means of simulations neglecting the effects of dust diffusion due to turbulence and of dust backreaction on the gas. Both effects may alter the formation of the dust peak at the gap outer edge and of the inner dust cavity, by letting grains filter through the dust trap. We performed high-resolution hydrodynamical simulations of the coupled evolution of gas and dust species, the latter treated as pressureless fluids, in the presence of two giant planets. We show that diffusion and backreaction can change some morphological aspects of the dust distribution but do not alter some main features, such as the outer peak and the expanding inner cavity. Furthermore, these findings are confirmed for different parametrizations of gas viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product Report

Objective analysis (OA) is a method of mapping unstructured data to a structured grid. In the context of scanning radar data, OA is used to interpolate data in antenna coordinates (range, azimuth, and elevation) onto a regularly spaced Cartesian grid (Trapp and Doswell 2000). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product (VAP) uses the Python ARM Radar Toolkit (Py-ART), a data model-driven interactive architecture for working with weather radar data, to map the data to a Cartesian grid (Helmus and Collis 2016). MMCG, with Py-ART built in, has the ability to take radar data in antenna coordinates and map the gates to a Cartesian grid using inverse distance weight functions such as Cressman (square) and Barnes (exponential), but also can filter the data during the interpolation. MMCG also allows arbitrary formulations for the radius of influence, which are matched to particular radar scanning strategies. This creates a complex parameter space for optimizing the retention of storm structure detail while minimizing artifacts. MMCG takes data processed with ARM’s Corrected Precipitation Radar Moments in Antenna Coordinates (CMAC) VAP and maps it to a Cartesian grid as the output product. A variety of fields that have been mapped to the Cartesian grid are then saved to plots to complement each grid file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Method for conducting in situ high-temperature digital image correlation with simultaneous synchrotron measurements under thermomechanical conditions

This work presents a novel method of obtaining in situ strain measurements at high temperature by simultaneous digital image correlation (DIC), which provides the total strain on the specimen surface, and synchrotron x-ray diffraction (XRD), which provides lattice strains of crystalline materials. DIC at high temperature requires specialized techniques to overcome the effects of increased blackbody radiation that would otherwise overexpose the images. The technique presented herein is unique in that it can be used with a sample enclosed in an infrared heater, remotely and simultaneously with synchrotron XRD measurements. The heater included a window for camera access, and the light of the heater lamps is used as illumination. High-temperature paint is used to apply a random speckle pattern to the sample to allow the tracking of displacements and the calculation of the DIC strains. An inexpensive blue theatrical gel filter is used to block interfering visible and infrared light at high temperatures. This technique successfully produces properly exposed images at 870 °C and is expected to perform similarly at higher temperatures. The average strains measured by DIC were validated by an analytical calculation of the theoretical strain. Simultaneous DIC and XRD strain measurements of Inconel 718 (IN718) tensile test specimens were performed under thermal and mechanical loads and evaluated. Finally, this approach uses the fact that with DIC, the total strain is measured, including plastic strain, while with XRD, only elastic strain is captured. The observed differences were discussed with respect to the effective deformation mechanisms.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Combining Remotely Sensed Evapotranspiration and an Agroecosystem Model to Estimate Center‐Pivot Irrigation Water Use at High Spatio‐Temporal Resolution

Abstract Estimating irrigation water use accurately is critical for sustainable irrigation and studying terrestrial water cycle in irrigated croplands. However, irrigation is not monitored in most places, and current estimations of irrigation water use has coarse spatial and/or temporal resolutions. This study aims to estimate irrigation water use at the daily and field scale through the proposed model‐data fusion framework, which is achieved by particle filtering with two configurations (concurrent, CON, and sequential, SEQ) by assimilating satellite‐based evapotranspiration (ET) observations into an advanced agroecosystem model, ecosys . Two types of experiments using synthetic and real ET observations were conducted to study the efficacy of the proposed framework for estimating irrigation water use at the irrigated fields in eastern and western Nebraska, United States. The experiments using synthetic ET observations indicated that, for two major sources of uncertainties of ET difference between observations and model simulations, which are bias and noise, noise had larger impacts on degrading the estimation performance of irrigation water use than bias. For the experiments using real ET observations, monthly and annual estimations of irrigation water use matched well with farmer irrigation records, with Pearson correlation coefficient ( r ) around 0.80 and 0.50, respectively. Although detecting daily irrigation records was very challenging, our method still gave a good performance with RMSE, BIAS, and r around 2.90, 0.03, and 0.4 mm/d, respectively. Our proposed model‐data fusion framework for estimating irrigation water use at high spatio‐temporal resolution could contribute to regional water management, sustainable irrigation, and better tracking terrestrial water cycle.

Environmental Sciences & Ecology↗
DOE OSTI2025

Characterization of Tank 9H Annulus Sample in Support of Residual Material Inventory Determinations

The Savannah River National Laboratory (SRNL) was requested by Savannah River Mission Completion (SRMC) to provide sample preparation and characterization of the Tank 9H annulus sample in support of Residual Material Inventory Determinations. One Tank 9H sample in three vials [HTF-9-25-13, HTF-9-25-14 and HTF-9-25-15], with each vial containing approximately 200 mL of the Tank 9H annulus salt solution, were delivered to the SRNL Shielded Cells for sample preparation and characterizations in February 2025. The density of the “as-received” solution contained in each of the three Tank 9H annulus sample vials were determined followed by a solid-liquid separation on each one using 0.45-micron Nalgene® nylon filter membranes. The resulting filtrates were combined to form the Tank 9H annulus sample with a total volume of about 600 mL. The combined wet solid fractions, about a total of 4.8 grams of salt material, remaining on the filter membranes were air-dried in the Shielded Cells for 72 hours. The total weight of the air-dried solids was 2.1 grams. These air-dried solids were washed with deionized water (DI water) at a phase ratio of 60 mL DI water/gram of solids to recover insoluble solids, if any. No visible or measurable quantity of insoluble solids were recovered after DI water washing of the air-dried solids because the air-dried solids completely dissolved in the DI water. The solid fraction-wash water was not combined with the 600 mL of the filtrate solution, and the resulting solution was not screened or analyzed for radionuclides. Aliquot sample volumes of the undiluted Tank 9H annulus sample were sent to the SRNL analytical services groups for radionuclides, elementals, anions and total mercury analysis by various methods including radiochemical separations/counting methods, inductively coupled plasma-atomic emission spectroscopy (ICP-AES), and Inductively Coupled Plasma Mass Spectroscopy (ICP-MS) and special preparations. All sample analyses were performed in triplicate. This report presents the analytical characterization results for the Tank 9H annulus sample. The results are also reported where analytical methods yielded additional analytes, other than those requested by SRMC. In the characterization of the Tank 9H annulus sample, the detection limits for all the analytes, as specified in the Technical Task Request (TTR) and Task Technical and Quality Assurance Plan (TTQAP), were met.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

D–MOPH–25: diverse MOF–molecule pairs for Henry’s constants prediction

Computational methods like grand-canonical Monte Carlo simulations and machine learning (ML) have accelerated metal–organic frameworks (MOF) exploration but are typically limited to a narrow range of adsorbates due to data availability and force field constraints. In this study, we introduce a dataset of diverse MOF–molecule pairs for Henry’s constant prediction, D–MOPH–25, which systematically explores a diverse chemical space by combining 113 molecular adsorbates with over 5000 MOF structures through an active learning process. D–MOPH–25 constitutes the most diverse adsorbate dataset used in any ML study of molecular adsorption in MOFs to date. Our workflow builds a benchmark for predicting Henry’s constants at 300 K, leveraging conformal prediction for uncertainty quantification. Assessment through Shannon entropy and uniform manifold approximation and projection confirms the comprehensiveness of D–MOPH–25 while highlighting the importance of robust classification to filter out unphysical data points in regression tasks. Although future enhancements in model architecture and sampling criteria could improve predictive performance, our dataset already spans the target space using only 2.31% of total possibilities. This comprehensive dataset facilitates assessment of model generalizability across adsorbate species and can establish a foundation for high-throughput MOF screening and ML-driven separation processes.

active learning↗
DOE OSTI2025

First Temperature Profile of a Stellar Flare Using Differential Chromatic Refraction

We present the first derivation of a stellar flare temperature profile from single-band photometry. Stellar flare DWF 030225.574−545707.45129 was detected in 2015 by the Dark Energy Camera as part of the Deeper, Wider, Faster program. The brightness (Δm g = −6.12) of this flare, combined with the high air mass (1.45 ≲ X ≲ 1.75) and blue filter (DES g, 398–548 nm) in which it was observed, provided ideal conditions to measure the zenithward apparent motion of the source due to differential chromatic refraction (DCR) and, from that, infer the effective temperature of the event. We model the flare’s spectral energy distribution as a blackbody to produce the constraints on flare temperature and geometric properties derived from single-band photometry. We additionally demonstrate how simplistic assumptions on the flaring spectrum, as well as on the evolution of flare geometry, can result in solutions that overestimate the effective temperature. Exploiting DCR enables studying chromatic phenomena with ground-based astrophysical surveys, and stellar flares on M dwarfs are a particularly enticing target for such studies due to their ubiquity across the sky and the heightened color contrast between their red quiescent photospheres and the blue flare emission. Our novel method will enable similar temperature constraints for a large sample of objects in upcoming photometric surveys like the Vera C. Rubin Legacy Survey of Space and Time.

Clarke, Riley W. [University of Delaware, Newark, ↗
DOE OSTI2025

Incident beam optics optimization for the single crystal neutron diffractometer Pioneer with a polarized beam option

Pioneer, a next-generation single-crystal neutron diffractometer, is under development for Oak Ridge National Laboratory’s Second Target Station. Designed to address a wide range of scientific questions, Pioneer will deliver homogeneous neutron beams with customizable size and divergence and provide a polarized beam option. This article introduces its incident beam optics, highlighting the optimization methodology and the simulated performance. Pioneer will utilize a modified elliptical-straight guide for neutron transport and deploy slit packages and insertable apertures to control beam size and divergence. The optimized guide geometry matches the optimal-and-full-sample-illumination condition, and the beam control system effectively filters out unwanted neutrons while preserving the desired ones. In addition, we have found that polygon-approximated guides provide satisfactory transport efficiency and beam homogeneity, eliminating the need for truly curved guides. Further, to enhance neutronics performance and reduce cost, the coatings of supermirror elements are individually optimized to the lowest half-integer m-values that are sufficient to deliver the desired neutrons. After evaluating polarizing V-cavities and 3 He spin filters over the default polarized wavelength band of 1.2–5.5 Å, we selected a translatable multichannel polarizing V-cavity as the incident beam polarizer. Strategically placed at a location where the beam divergence is low and a large in-guide gap has negligible impact on transport efficiency, the optimized V-cavity achieves an average P 2 T of ~35%.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A Fast Reactor Irradiation Experiment Design in the ATR

Modern modeling techniques were used to investigate a proposed method for fast neutron irradiations in an existing thermal-spectrum reactor, the Advanced Test Reactor (ATR). This method builds upon pre-existing ideas, where fast flux is increased by surrounding the specimens with fissionable “booster fuel” but diverges from historical approaches by using an already developed fuel element design used in the Belgian Reactor 2 (BR2) as the booster fuel while leveraging modern 3-D modeling and simulation techniques. Design evaluations and neutronics simulations were performed to evaluate the performance of a BR2 fuel element irradiated in an ATR flux trap with test pins in the central channel of the BR2 fuel element. These efforts have yielded promising results. Adding a BR2 fuel element in the northeast (NE) flux trap of ATR was predicted to result in a 150% increase to the incident fast neutron flux with a fast (>0.1 MeV) to thermal (<0.625 eV) neutron flux ratio ranging from approximately 50 to 150, dependent on the material used for thermal neutron filtering and volume of moderator within the central channel of the BR2 element. The predicted annual fast neutron fluence (>0.1 MeV) ranges from 7.9 × 1021 to 9.1 × 1021 n/cm2. Given the relatively large fast to thermal neutron flux ratio, the calculated radial power profiles within 4.3 mm outer diameter U10Zr fueled specimens irradiated within the BR2 booster fuel element are adequate representations of those within fast neutron reactors. The predicted radial power profiles are not flat, but they are more prototypic than those seen in advanced fuels tests which began in ATR in 2003. Another distinct advantage of this experiment design is that full-scale test pins can be irradiated to augment the ongoing series of reduced scale advanced fuels tests. The proposed experiment design irradiated within a BR2 fuel element in a flux trap of ATR offers an improved alternative to the current testing of advanced reactor fuels in ATR. Selected results from this design evaluation are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Next-Generation, High-temperature, High-frequency, High-efficiency, High-power-density Traction System

To meet performance and reliability requirements necessary for broader adoption of electric drive vehicles, the Electrical and Electronics Technical Team of the U.S. Drive partnership has established aggressive design goals for next-generation electric vehicle drivetrains. Specifically, the 2025 roadmap stipulates a 100 kW/L power density target and a $\$$2.7/kW cost target for power electronics, in addition to high-voltage operation (i.e., greater than 800 VDC). The additional targets for traction motor and the overall system performance impose further challenges on the power electronics design. For example, many high specific power machines have reduced iron content, and therefore reduced intrinsic filtering, thus requiring the inverter to supply a low-distortion drive current. These machines also typically have a high pole count, thus requiring drive current at a higher electrical frequency. Other motors, such as brush-less dc and switch reluctance machines, require a carefully-shaped, non-sinusoidal drive current (Yang, Shang, Brown, & Krishnamurthy, 2015), (Zhang, Bowman, O'Connel, & Haran, 2018), (Anderson, et al., 2018). Two- and three-level inverter topologies are the conventional framework for the power electronics design of the drivetrain, and some demonstrations have shown recent progress towards addressing cost, power density and efficiency goals (Gurpinar & Ozpineci, 2018), (Zhu, Kim, Chen, Erickson, & Maksimović, 2018), (Deshpande, Chen, Narayanasamy, Sathyanarayanan, & Luo, 2018), (Alizadeh, et al., 2019). However, an unconventional approach may be necessary to take the dramatic leap in power density necessitated by the roadmap—while simultaneously addressing the other system needs. Therefore, this project leverages the flying capacitor multilevel (FCML) topology, together with a scalable, modular approach, to address these needs. This type of hybrid converter has several advantages: lower voltage (i.e., less than 300 V) transistors can be used, energy-dense capacitors process most of the power, and the output current waveform is multilevel and exhibits a frequency multiplying effect—in other words, the output has reduced dv/dt and filtering requirements for the same high voltage dc bus. For example, in an electric vehicle with an 800 V bus, a 10-level FCML could leverage 100 V, commercially available GaN devices switching at 115 kHz to produce a ~1 MHz switching waveform (modulated according to the motor drive requirements) with one ninth of the dv/dt of a two-level converter. Prior work has already demonstrated promising performance and gravimetric power density figures for more electric aircraft applications (Pallo, Foulkes, Modeer, Coday, & Pilawa-Podgurski, 2018). This project leverages lessons learned to achieve the volumetric power density of 100 kW/L by employing advanced liquid cooling, address the 300,000 mile reliability challenge with redundant design, topology failure studies and online health monitoring, and reduce costs to $\$$2.7/kW through the use of low-cost GaN devices, modular converter assemblies, and modest modifications to traditional manufacturing methods. The project involved several hardware designs, each achieving increasing performance. At the conclusion of the project, a volumetric power density of 380 kW/L was achieved, in a 800V dc-ac converter, greatly surpassing even the aggressive target goal.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Thermal Model Heat Rate Predictions of the AGR-5/6/7 Experiment

Several fuel and material irradiation experiments have been planned for the U.S. Department of Energy Advanced Gas Reactor Fuel Development and Qualification Program, which supports the development and qualification of tristructural isotropic (TRISO) coated particle fuel for use in high-temperature gas-cooled reactors. The goals of these experiments are to provide irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, support development of fuel performance models and codes, and provide irradiated fuel and materials for post-irradiation examination and safety testing. Originally planned and named as separate fuel experiments, but subsequently combined into a single test train. AGR-5/6/7 will test low-enriched uranium oxycarbide (UCO – a heterogeneous mixture of uranium oxide and uranium carbide) TRISO fuel. The AGR-5/6 portion of the experiment will provide data to support qualification of the selected reference fuel design, while the AGR-7 portion will serve as a margin test, irradiating the fuel beyond normal operating conditions. During irradiation, the temperature in each capsule is controlled by varying the helium / neon gas mixture flowing through the capsule until the desired thermocouple temperature is reached. Further adjustments can be made by changing the lobe power and the neutron filter. Three neutron filters are available to use as the experiment burns up the fissionable fuel. A thermal finite element model has been created for the five capsules comprising the AGR-5/6/7 experiment. Heat rates from a physics analysis are imported into the model along with fast neutron fluence. Graphite shrinkage due to the fast neutron fluence is incorporated into the model. Gas gaps change as a function of fast neutron fluence. Radiation heat transfer is a major contributor in this model. This is a large model with more than 1 million finite element brick elements. More than 150 different parts are modeled in the finite element model and communicate with each other from a heat transfer sense. More than 50 thermocouples are used in the experiment and are compared to actual measurements. The experiment is composed of five separate stainless steel capsules all welded together. This paper presents a method used to predict the gas mixture in the event that all of the thermocouples fail in a particular capsule. Curve fitting of heat rates from previous cycles and projecting them into the future is the basis of the method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗
DOE OSTI2022

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗