Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “qualitative evaluation”

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 325 records · Page 18

Fabrication of Black Body Grids by Thick Film Printing for Quantitative Neutron Imaging

Neutron imaging offers deep penetration through many high-Z materials while also having high sensitivity to certain low-Z isotopes such as 1H, 6Li, and 10B. This unique combination of properties has made neutron imaging an attractive tool for a wide range of material science and engineering applications. However, measurements made by neutron imaging or tomography are generally qualitative in nature due to the inability of detectors to discriminate between neutrons which have been transmitted through the sample and neutrons which are scattered by the sample or within the detector. Recent works have demonstrated that deploying a grid of small black bodies (BBs) in front of the sample can allow for the scattered neutrons to be measured at the BB locations and subsequently subtracted from the total measured intensity to yield a quantitative transmission measurement. While this method can be very effective, factors such as the scale and composition of the sample, the beam divergence, and the resolution and construction of the detector may require optimization of the grid design to remove all measurement biases within a given experimental setup. Therefore, it is desirable to have a method by which BB grids may be rapidly and inexpensively produced such that they can easily be tailored to specific applications. In this work, we present a method for fabricating BB patterns by thick film printing of Gd2O3 and evaluate the performance with variation in feature size and number of print layers with cold and thermal neutrons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computed Tomography Scanning and Geophysical Measurements of the CarbonSAFE Seal Integrity Wells in the Illinois Basin

The U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) researchers in Morgantown, West Virginia, utilized computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) to evaluate the integrity of three cores (Geophysical Monitor #2 (GM2), Verification Well #1 (VW #1), and Wabash #1 for the Illinois Basin Carbon Storage Assurance Facility Enterprise (CarbonSAFE) efforts. This work was carried out in collaboration with the Illinois State Geological Survey (ISGS) as part of their efforts to analyze cores from two field locations, including the Illinois Basin-Decatur Project (IBDP) and second stage Illinois Industrial Carbon Capture and Storage (IL-ICCS) sites at Decatur, Illinois, and the Wabash CarbonSAFE project at the Wabash Valley Resources Integrated Gasification Combined Cycle plant in Terre Haut, Indiana. The results of this study are presented in several formats and are available online on the Energy Data eXchange (EDX). The rock characterization was conducted using non-destructive techniques, allowing for future analysis of the cores. While the equipment used did not provide direct visualization of shale pore space, it allowed for detection of fractures and discontinuities. Low resolution CT imagery with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with x-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were utilized to identify areas of interest for further study as well as fractured zones. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provided a multi-scale analysis of this core and provided both a macro and micro description of the core that is relevant for many subsurface energy-related examinations traditionally performed at NETL.

58 GEOSCIENCES↗

Evaluation of BISON's Transient Fission Gas Release Model on High-Burnup UO 2 MiniFuel

BISON's fission gas release (FGR) model has been evaluated and tested over several hypothetical thermal transient cases at high burnup and validated against FGR experimental data from an annealing test at high burnup. Diffusion-controlled FGR is accounted for in the model by including several physical mechanisms such as the diffusion of fission gas to grain faces, grain boundary sweeping, grain growth, and intergranular bubble growth. Under temperature variations, FGR is increased due to the development of micro-cracks at grain faces. This existing transient capability, however, was developed based on transient fuel behavior during power ramps and the associated micro-cracking process. Nevertheless, a fraction of the FGR during temperature transients such as LOCA can ensue from different mechanisms (in particular, fuel fragmentation) and present an important dependency on local burnup and high burnup structure formation. It follows that the performance and necessary developments of the existing BISON fission gas model for the analysis of LOCA-type transients need to be investigated through a systematic assessment study. Such a study is performed in this work. The model shows an overall good qualitative agreement with experimental observations on FGR during annealing tests with different temperature conditions. It also accounts well for microstructural effects on FGR. When quantitatively compared to FGR data from previously irradiated 103 MWd/kgU UO 2 discs under thermal annealing, the model shows a less satisfactory agreement with the experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Cyber risk assessment and investment optimization using game theory and ML-based anomaly detection and mitigation for wide-area control in smart grids

The electric power grid is increasingly becoming susceptible to cyber attacks that exploit vulnerabilities in the smart grid control, information, and physical layers. Successful cyber attacks can have catastrophic impacts on the social and economic well-being of any nation all over the globe. It has, thus, become imperative to secure the smart grid against such adversarial actions to ensure stable, secure, and reliable operation of the grid. The existing research and industry practices prove to be inadequate in terms of providing pragmatic and effective defense methodologies and measures for long-term cybersecurity planning and real-time cybersecurity for grid operation. For example, existing works lack models that incorporate uncertain behavior of cyber-attackers and pragmatic defense measures for cyber risk assessment and cybersecurity investment optimization which often provide unreliable and strictly qualitative solutions to these problems. At the same time, with the growing number of cyber incidents in the grid, there still exists a need to develop attack-resilient algorithms for wide-area monitoring, protection, and control (WAMPAC) applications like the wide-area voltage control systems (WAVCS) for Flexible AC Transmissions Systems (FACTS) that lack in scalable and feasible solutions from the cybersecurity perspective. This dissertation proposes novel models and methodologies for: (1) Cybersecurity planning, and (2) Cybersecurity for system operation. The cybersecurity planning is achieved through cyber risk assessment and cybersecurity resource investment optimization for long-term cybersecurity of the grid using game theory and attack-defense trees. Cybersecurity for system operation consists of development of cyber anomaly detection and mitigation algorithms for flexible AC transmission system (FACTS) controller-based wide-area voltage control systems (WAVCS) using machine learning (ML), and software defined networking-based moving target defense network routing for achieving real-time cyber-physical security for grid operations. This is followed by hardware-in-the-loop (HIL) implementation and evaluation of these attack prevention, detection, and mitigation algorithms and methodologies showcasing their feasibility in a close to real-world environment. For cybersecurity planning, a novel approach involving a combination of game theory and attack defense trees (ADT) for optimal cybersecurity resource allocation in the smart grid is proposed. This methodology involves modeling of the cyber-physical smart grid substations as ADTs, defining attacker costs, defense costs, and attack probabilities for attack access points. Using game theoretical formulation, optimal defense strategies for the defender of the system to invest cybersecurity resources in the grid are obtained. Additionally, a game-theoretic framework is developed for quantitative cyber-physical risk assessment of the grid under a dynamically changing cyber threat space and uncertain behavior of cyber attackers which is further used to optimize investments in the smart grid's cybersecurity resources. The attacker, defender, and the smart grid system are modeled while incorporating attacker-stochasticity and federal guidelines for smart grid cybersecurity. This allows quantification of threat, vulnerabilities, and attack impact of the grid for quantitative risk assessment. The defender's budget to invest in the security resources in the grid is optimized based on the strategies leading to minimum system risk. The evaluation of the proposed solutions highlight the feasibility for practical implementation of these methodologies and algorithms in the smart grid, while taking the federal requirements and guidelines for smart grid security into consideration. For achieving cybersecurity for system operation, attack prevention, detection, and mitigation algorithms and methodologies are developed specifically for FACTS-based WAVCS. Anomaly detection and mitigation in the WAVCS are achieved using algorithms based on machine learning which involves offline training and testing of ML models with CPS datasets incorporating physics-based features that allow accurate distinction between system faults and cyber attacks. For attack prevention, a methodology based on software defined network (SDN)-based moving target defense (MTD) network routing is proposed that enables prevention of Denial of Service (DoS) type attacks on the smart grid communication system. Subsequently, these methodologies and algorithms are implemented and evaluated on an HIL testbed that allows for real-time attack prevention, detection, and mitigation of emulated cyber attacks on the WAVCS in a close to real-world environment. The results show highly accurate and efficient performance of the implemented algorithms and methodologies with the smart grid system operating within the NERC's system operation limits even in the presence of DoS and data integrity cyber attacks. This work opens up future research opportunities in other directions such as (1) Expanding cybersecurity planning methodologies to real-time cyber contingency analysis with different game formulations; and (2) Applying the cybersecurity for system operation algorithms to broader categories of wide-area control applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An Evaluation of Tropical Cyclone Rainfall Structures in the HighResMIP Simulations against Satellite Observations

This study evaluates tropical cyclone (TC) rainfall structures in the CMIP6 HighResMIP global climate model (GCM) simulations against satellite rainfall retrievals. Here, we specifically focus on TCs within the deep tropics (25°S–25°N). Analysis of TC rain rate composites indicates that in comparison to the satellite observations at the same intensity, many HighResMIP simulations tend to overproduce rain rates around TCs, in terms of both maximum rain rate magnitude and area-averaged rain rates. In addition, as model horizontal resolution increases, the magnitude of the peak rain rate appears to increase. However, the area-averaged rain rates decrease with increasing horizontal resolution, partly due to the TC eyewall being located closer to the TC center, thus occupying a smaller area and contributing less to the area-averaged rain rates. The effect of ocean coupling is to lower the TC rain rates, bringing them closer to the satellite observations, due to reduced horizontal moisture flux convergence and surface latent heat flux beneath TCs. Examination of horizontal rain rate distributions indicates that vertical wind shear–induced rainfall asymmetries in HighResMIP-simulated TCs are qualitatively consistent with the observations. In addition, a positive relationship is observed between the area-averaged inner-core rainfall and TC intensification likelihoods across the HighResMIP simulations, as GCM simulations producing stronger TCs more frequently have the greater rainfall close to the center, in agreement with previous theoretical and GCM simulation results.

54 ENVIRONMENTAL SCIENCES↗

Relationships between Immersion Freezing and Crystal Habit for Arctic Mixed-Phase Clouds—A Numerical Study

Abstract The number concentration of ice particles in Arctic mixed-phase clouds is a major controlling factor of cloud lifetime. The relationships between ice nucleation mode and ice crystal habit development are not yet constrained by observations. This study uses a habit-predicting microphysical scheme within a 3D large-eddy simulation model to evaluate the relationship between immersion freezing and ice habit in a simulated Arctic mixed-phase cloud case. Three immersion freezing parameterizations are considered: a volume-dependent freezing scheme (VF), a parameterization limited to activated droplets (C-AC), and a parameterization limited to coarse aerosol particles (C-CM). Both C-AC and C-CM are based on classical nucleation theory. The freezing rate with VF is found to be greater in downdraft regions than in updraft regions due to the downdraft having a higher number concentration of large droplets. The C-AC cases show active freezing of small droplets near cloud top, whereas in the C-CM cases, mainly the 8–32- μ m-sized droplets freeze in updraft regions near the cloud base. Because the initial crystal size is assumed to affect the axis ratio of hexagonal plates, the VF cases produce crystals with larger axis ratios, resulting in smaller mode radii than the C-AC cases. In all cases, irregular polycrystals dominate near cloud top and a band-like structure develops within the cloud, which qualitatively agrees with previous observations. In the VF and C-CM cases, unactivated large droplets arising from coarse-mode aerosol particles contributed significantly to the freezing rate, producing an important influence on crystal habit.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Impurities in Neptunium Metal Using Double-Focusing SF-ICP-MS

In the transmutation process, minor actinides such as neptunium, americium and curium obtained from spent nuclear fuel are incorporated into new fuel that can be placed into a fast reactor, where these elements fission into products with shorter half-lives which reduces the health and safety hazards. In order to better understand the all the involved processes, fundamental physical properties of actinides should be investigated, and for that reason, high-purity materials need to be produced. One of the main challenges in characterizing these materials is the lack of pure metal available for study, and lack of analytical methods to determine the purity. Trans-uranic elements are challenging to study because of the need to be handled in specialized facilities, safeguards concerns and they are not easily purchased and/or produced. High-purity Np metal is being produced at the Materials and Fuels Complex (MFC) of Idaho National Laboratory (INL), in a process of direct oxide reduction developed by Squires and Lessing, with a target purity of 99.999%. The goal of this research is to develop methods to determine the impurities in the Np metal, with detection limits low enough to achieve this purity level. In this work, we have used several ICP-MS and ICP OES techniques, but the focus was the method development using a recently acquired Double-Focusing Sector Field Inductively Coupled Plasma Mass Spectrometer (SF-ICP-MS). To produce the Np metal, approximately 20 g of the neptunium oxide powder was calcined at 800 oC for 6 h in a magnesium oxide crucible, which is placed inside a stainless steel crucible, that goes into the Hot Uniaxial Press (HUP) furnace insert. After the calcining step, the Np oxide was mixed with approximately 9 g of calcium metal and 100 g of ultra-dry calcium chloride. The HUP furnace was heated in a ramp rate of 50 oC min -1 up to 1000 oC. When the reaction mixture reached 850 oC, a tantalum stirrer was activated at 100 rpm. After the tantalum stirrer was removed and the mixture cooled overnight, the MgO crucible was opened with a hammer and chisel to retrieve the Np metal button at the bottom of the reaction mixture. Prior to the dissolution process, the surface contaminants of the sample were removed. The Np button was washed in 10 mL of ultra-pure water for 1 min, with a faint reaction observed on the surface, indicating the presence of Ca metal. The water-washed sample was placed in a tube with 2 M HCl for only 3 s, as a vigorous reaction in the button surface was observed. The leachates were saved for further analysis. The sample was then rinsed with ultra-pure water and air-dried prior to dissolution. The cleaned sample was placed in a PTFE beaker and 25 mL of 6 M HCl and 150 µL of 24 M HF were added. After 10 minutes of reaction, the beaker was placed on a hot plate for 1 h and 30 min, followed by the addition of 1 mL of 16 M HNO3 and allowed to react for a further 1 h and 20 min. The sample was removed from heat and 10 mL of HNO3 16 M was added to react with the HCl. The beaker was covered and left for overnight reaction, producing a solution with approximately 20 mL of final volume. A Quadrupole Inductively Coupled Plasma Mass Spectrometer (Q-ICP-MS, Perkin Elmer ELAN DRC II Q-ICP-MS) was used to obtain qualitative mass spectra for the dissolution blank, each of the leachates and the dissolved sample. In the water leachate spectrum mainly Ca and Mg were observed, while the acid wash had a similar peak pattern to the one observed for the Np metal sample, indicating that the process not only removed surface contamination, but also reacted with the sample itself. In the Np material spectrum, Ta and actinides were identified as candidates for major contaminants, as the most intense peaks corresponded to the mass-to-charge ratios of these elements. Minor contaminants included Al, Au, Ca, Cr, Fe, I, Na, Nb, Pb, Sn, Ti, Tl, V, Y. The presence of impurities in the material could have originated from several sources, such as the materials used for the casting process, the use of impure NpO2 or contamination of reagents and materials used for the analysis. Some of the isotopes listed as potential contaminants are susceptible to spectral interferences, so Double-Focusing SF-ICP-MS (Nu Instruments Attom ES HR-ICP-MS) was used for the quantitative analysis. For each isotope determined, several items were evaluated: raw signal of a blank solution, sensitivity and determination coefficient of calibration curves, possibility of interferences and comparison of results obtained for at least two resolving powers (RP): 300 and 4000. The evaluation of the signal in a blank solution can indicate if a spectral interference from the plasma gas and/or the solvent was successfully resolved (commonly observed for 52Cr+/40Ar12C+, for example). It is important to highlight that a signal reduction is expected when a higher RP is used, as the ion beam transmission is decreased by the use of narrower slits

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanocomposites Based on Coil-Comb Diblock Copolymers

The polymer chain architecture is an important factor determining the phase behavior of nanoparticle (NP) assembly in polymer matrices. Block copolymers (BCPs) containing a random copolymer (RCP) block present a convenient variation on traditional BCPs to tune the interaction parameters between the polymer blocks and the nanofillers as well as to evaluate the effect of the chain architecture on the NP arrangements within BCP microdomains. Here, we synthesized BCPs with a coil polystyrene (PS) block and a comb RCP block through reversible addition-fragmentation chain transfer polymerization. The RCP block consists of methyl- and lauryl acrylates, the latter of which confers a long-chain alkyl moiety to favorably interact with alkyl-passivated NPs. BCPs showing lamellar, cylindrical, and mixed morphologies were obtained by varying the volume fractions of the RCP block (f RCP ). In comparison to coil-coil BCP, the coil-comb BCPs show highly asymmetric phase behavior with respect to fRCP, where lamellar morphologies were observed at fRCP from 0.31 to 0.51. NPs in the size of 4-5 nm were successfully incorporated in the RCP block of the BCPs with periodicities of 30-60 nm. An order-to-order phase transition from lamellae to PS cylinders was observed after the addition of only 1-2 vol % 5 nm NPs into the BCP with a periodicity of 25 nm and f RCP of 0.51. Self-consistent field theory-density functional theory simulations qualitatively described the observed morphologies and phase transitions in the nanocomposites. Here, the current study presents a platform to fabricate nanocomposites with NP assemblies in coil-comb BCPs that contain a random copolymer block and provides insight into how polymer chain architectures can affect the phase behavior of BCPs and nanocomposites.

36 MATERIALS SCIENCE↗

Deactivation trends of Pd/SSZ-13 under the simultaneous presence of NO, CO, hydrocarbons and water for passive NO x adsorption

Pd-functionalized chabazite (Pd/SSZ-13) was evaluated for passive NO x adsorption (PNA) using low temperature combustion with diesel (LTC-D) reaction feed of the “United States Driving Research and Innovation for Vehicle efficiency and Energy sustainability” (U.S.DRIVE) protocol. Notably as per the protocol conditions, 12% O 2 + 6% CO 2 + 6% H 2 O was flown in all cases. NO-uptake studies in the presence of the LTC-D feed showed a systematic decline with (NO:Pd)molar changing from 0.5 to 0.4 after 10 trials. Control experiments showed more pronounced decline for NO + CO case, with appreciable intial (NO:Pd)molar value of 0.4 yet comparable-to-the-LTC-D-feed decline after 8 trials. CO-induced particle formation and larger extent of particle sintering was also evident from TEM analysis. Other controls did not exhibit trial-dependent deactivation as(NO:Pd)molar values were constant, at 0.3, 0.2, 0.5 and 0.2 for the NO, NO+H 2 , NO + unsaturated hydrocarbon (C 2 H 4 , C 3 H 6 ) and NO + saturated HC (C 3 H 8 , C 10 H 22 ) feeds respectively. In addition to these quantitative differences, desorption behaviors are qualitatively different. While only one major desorption event is observed for the full LTC-D, NO+ CO and NO+ unsaturated HC-controls, desorption occurs in two distinct stages for NO, NO + saturated HCs and initial-NO+H 2 trials. This arises due to inherent differences in Pd sites while exposed to these chemically distinct feeds. Furthermore, presence of reducing agents such as CO and unsaturated HCs in the feed result in almost complete elimination of lower, sub-200 °C desorption peak. The higher temperature desorption peaks, at > 300 °C is associated with NO strongly bound to ionic Pd sites and are prevalent under reducing conditions. Using DRIFTS also, three complexes leading to PNA are assigned, [O = N–Pd 2 +(OH)–Z], [O = N–Pd 2 +(Z2)] and [O = N–Pd 2 +(H 2 O)y–Z] with the latter being clearly observed upon water exposure. Pd/SSZ-13 showed higher hydrocarbon trapping than the SSZ-13 counterpart.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microjet printing of metal salt or oxide targets for nuclear reaction studies on radionuclides

Radioactive targets for direct measurements of neutron-induced reactions are required to improve evaluated cross-section data and, ultimately, the fidelity of neutron reaction network simulations. Electrodeposition and molecular plating are the current state-of-the-art radioactive target production techniques, but not all metals can be electrodeposited or molecular plated with high yields. Alternative techniques can be expensive or may produce targets that are unsatisfactory in terms of thickness, yield, or purity. Microjet printing is a new, rather inexpensive technique that utilizes equipment with a small footprint and has the potential to produce thin, highly radioactive targets with good uniformity and minimal impurities from the target fabrication process. This study involved optimizing the process of producing microjet printed targets to allow for the fabrication of a stable vanadium(V) oxide (V 2 O 5 ) target that was compared with an analogous electrodeposition V 2 O 5 target manufactured via the application of a vanadium chemical conversion coating on aluminum (Al) foil. Finally, the results from these studies suggest microjet printing could be used to produce relatively uniform target layers with adequate film thicknesses (< 15 μm). V 2 O 5 microjet printed targets, when compared with chemical conversion coating V 2 O 5 targets, appeared to be qualitatively less uniform and quantitatively larger in thickness (11.7(27) μm vs. 5.3(18) μm). However, the conversion coating target contained more impurities and the Al backing had a higher background contribution to measurements of neutron-induced reactions as opposed to targets produced via microjet printing. Overall, the results from this study suggest microjet printing has the capability to be an excellent alternative target production technique to the electrodeposition and molecular plating methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dynamical instability in multi-orbiter systems with gas friction

ABSTRACT Closely packed multiplanet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in active galactic nuclei discs around supermassive black holes. We use N-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of K (the difference in the planetary semimajor axes in units of the mutual Hill radius) and (unequal) planet masses. We find that, in general, the stable region (large K) and unstable region (small K) are separated by a “grey zone”, where the (in)stability is not guaranteed. We report the numerical values of the critical spacing for stability Kcrit and the “grey zone” range in different systems, and provide fitting formulae for arbitrary frictional forcing strength. We show that the stability of a system depends on the damping time-scale τ relative to the zero-friction instability growth time-scale tinst: two-planet systems are stable if tinst ≳ τ; three-planet systems require tinst ≳ 10τ−100τ. When K is sufficiently small, tinst can be less than the synodic period between the planets, which makes frictional stabilization unlikely to occur. As K increases, tinst tends to grow exponentially, but can also fluctuate by a few orders of magnitude. We also devise a linear map to analyse the dynamical instability of the “planet + test mass” system, and find qualitative agreement with N-body simulations.

Li, Jiaru (ORCID:0000000155507421)↗

Business Models for Coal Plant Decommissioning

Aging coal-fired power plants are retiring across the United States. Researchers at Pacific Northwest National Laboratory studied coal-fired power plant decommissioning business models to support program development for the Department of Energy’s Office of Energy Efficiency & Renewable Energy (EERE) to engage with communities affected by power plant retirement processes. Evaluating and reducing the impact on coal-dependent communities from phasing out coal production is a crucial piece of nationwide economic development and community-directed engagement in the clean energy future. This report describes the steps in the typical coal-fired power plant decommissioning process, including analysis of the following key factors: Drivers including policy and regulations, competition with other fuels, and corporate responsibility goals; Types of coal plant decommissioning including plant retirement without full decommissioning, repurposing with fuel switch, redeveloping to utilize existing transmission, and decommissioning and repurposing with other commercial activities; and, Costs of coal plant decommissioning, including variations across regulated and deregulated markets and the funding resources to support decommissioning activities. The details described in this report relied heavily on the availability of local news reports and communications regarding plant retirements, as the researchers identified a lack of formal literature surrounding business models for coal plant decommissioning and a need for a more robust source material collection. In addition, a lack of oversight entities responsible for decommissioning processes may have led to this shortfall of information. Future research regarding coal plant decommissioning can include quantitative studies on power plant retirement drivers, procedures, and costs, qualitative assessments on affects decommissioning has on surrounding communities, investors, and other stakeholders, and identifying best practices to engage community stakeholders in the decommissioning process.

01 COAL, LIGNITE, AND PEAT↗

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

Discrete modeling of a transformer with ALEGRA

We report progress on a task to model transformers in ALEGRA using the “Transient Magnetics” option. We specifically evaluate limits of the approach resolving individual coil wires. There are practical limits to the number of turns in a coil that can be numerically modeled, but calculated inductance can be scaled to the correct number of turns in a simple way. Our testing essentially confirmed this “turns scaling” hypothesis. We developed a conceptual transformer design, representative of practical designs of interest, and that focused our analysis. That design includes three coils wrapped around a rectangular ferromagnetic core. The secondary and tertiary coils have multiple layers. The tertiary has three layers of 13 turns each; the secondary has five layers of 44 turns; the primary has one layer of 20 turns. We validated the turns scaling of inductance for simple (one-layer) coils in air (no core) by comparison to available independent calculations for simple rectangular coils. These comparisons quantified the errors versus reduced number of turns modeled. For more than 3 turns, the errors are <5%. The magnetic field solver failed to converge (within 5000 iterations) for >10 turns. Including the core introduced some complications. It was necessary to capture the core surfaces in thin grid sheaths to minimize errors in computed magnetic energy. We do not yet have quantitative benchmarks with which to compare, but calculated results are qualitatively reasonable.

42 ENGINEERING↗

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Field Test Report Neutron Scintillator Array Dry Storage Cask Scanner FY2024

During two weeks of Field Testing at the Idaho National Laboratory INTEC Cask Farm in July and August 2024, the LLNL Dry Storage Cask Scanner Array was lifted on top of an MC-10 dry storage fuel cask and operated to acquire neutron and gamma-ray data from the 24 fuel bundle positions. Neutron and gamma-ray data acquisition scans across the top of the cask of varying dwell times were performed July 15-18, 2024 and August 19-22, 2024 to evaluate the ability of the scanner data to reveal asymmetries in the fuel positions that reflect asymmetries in the MC-10 cask fuel bundle loading. The MC-10 cask 24 position fuel bundle loading at the INTEC Cask Farm is well documented, including the locations of six empty fuel bundle positions. This loading presents an opportunity to test the ability of the scanner system to detect diversion of spent fuel bundles as well as to validate the MC-10 cask MCNP modeling. The cask scanner array consists of six Stilbene crystal scintillator detectors and a linear actuator frame that moves the six detectors across the MC-10 dry storage cask to obtain data above each of the 24 fuel bundle positions. The detectors are connected to a pulse-shape discrimination data acquisition system capable of generating separate neutron and gamma-ray spectra for each detector and for each scan position. From the prior single detector Field Test in 2021 and iteration with MCNP modeling, the neutron and gamma-ray data were analyzed in multiple energy regions to identify an analysis method that would provide the strongest and most consistent signature of the asymmetric MC-10 cask fuel loading1 . From both the 2021 Field Test and the current Field Test results, the neutron capture gamma-ray count rate around 2.2 MeV provides the strongest signature of the asymmetric MC-10 cask fuel loading and has qualitative agreement with MCNP calculations. Counting all gamma-rays produces a similar signature. Neutrons emerging from the cask top are moderated and captured by the hydrogen in the polyethylene moderator and scintillator detector, producing a 2.2 MeV gamma ray which is seen in the scintillator gamma-ray spectrum. The count rate in the 2.2 MeV gamma-ray region is ~50 c/s, which is ~1000x higher than the ~0.05 n/s rate in the > 4MeV neutron region, and ~50x greater than the ~1 n/s rate in the neutrons > 500 keV region. Analysis of the 2.2 MeV neutron-capture Compton-scattered gamma-rays produces a statistically significant signature of the INTEC Cask Farm MC-10 asymmetric fuel loading. MCNP simulations indicate that the average neutron energy spectrum offers the potential to detect a large asymmetry from several missing bundles as well as individual missing fuel bundles. Testing this feature will require measurements on a cask with single missing elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An uncertainty-focused database approach to extract spatiotemporal trends from qualitative and discontinuous lake-status histories

Changes in lake status are often interpreted as palaeoclimate indicators due to their dependence on precipitation and evaporation. The Global Lake Status Database (GLSDB) has since long provided a standardised synopsis of qualitative lake status over the last 30,000 14C years. Potential sources of uncertainty however are not recorded in the GLSDB. Here we present an updated and improved relational-database framework that incorporates uncertainty in both chronology and the interpretation of palaeoenvironmental data. The database uses peer-reviewed palaeolimnological studies to produce a consensus on qualitative lake-status histories, whose chronologies are revised and standardized through the recalibration of radiocarbon dates and the application of Bayesian age-depth modelling for stratigraphic archives. Quantitative information on absolute water-level elevation is preserved if available from geomorphological sources. We also propose a new probabilistic analytical framework that accounts for these uncertainties to reconstruct synoptic, integrated environmental signals. The process is based on a Monte Carlo algorithm that iteratively samples individual lake-status histories within the limits of their uncertainties to produce many possible scenarios. We then use Recursively-Subtracted Empirical Orthogonal Function analysis to extract dominant patterns of lake-status variability from these scenarios. As a proof of concept, we apply this framework to 67 sites in eastern and southern Africa whose lake-status histories cover part of the late Pleistocene and/or Holocene. We show that, despite the sometimes large temporal and interpretation uncertainties, and the inclusion of highly discontinuous lake-status time series, identifying the major known millennial-scale climatic phases during the last 20,000 years is possible. Our framework was also able to identify an antiphased response between the lake basins in eastern and interior southern Africa to these changes. Here, we propose that our new database and methodology framework serves as a template for efficient lake-status data synthesis, encourages the incorporation of lake-status data in palaeoclimate syntheses, and expands the possibilities for the use of such data in the evaluation of climate models.

58 GEOSCIENCES↗