MFIX-Exa: CFD-DEM Simulations of Thermodynamics and Chemical Reactions in Multiphase Flows
For presentation at the AIChE Annual Meeting, Orlando, FL, November 5-10, 2023.
Porcu, Roberto↗
SEARCH · Engineering Papers
Results for “ANN”
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.
Time-resolved measurements of OH during auto-ignition of syngas with trimethylsilanol and hexamethyldisiloxane
The effects of trimethylsilanol (TMSO) and hexamethyldisiloxane (HMDSO) addition on OH time histories during syngas (H 2 and CO) ignition were investigated using the University of Michigan rapid compression facility. Experiments spanned temperatures of 1010–1080 K, at a pressure of approximately 5 atm. Syngas mixtures of 1.2 % H 2 /2.8 % CO/20 % O 2 by volume (balance N 2 and Ar) provided a baseline for comparison with mixtures that included 100, 200, and 1000 ppm of the TMSO and 100 ppm of HMDSO. Narrow-line ultraviolet laser-absorption was used to measure OH mole-fraction during ignition. The addition of TMSO and HMDSO significantly shifted the OH time-histories earlier in time, by up to 51 %, compared with the baseline syngas mixture. The value of the maximum OH mole fraction was consistent between the 100 and 200 ppm TMSO mixtures and the 100 ppm HMDSO mixtures, but the maximum OH increased significantly with the 1000 ppm TMSO mixtures. Here, the OH data indicate TMSO and HMDSO were not direct sources of OH radicals. Analysis further indicates the TMSO and HMDSO decompose rapidly followed by reactions that enhance the production of H atoms, and the increased reactivity observed is via the H + O 2 = OH + O reaction.
Hexamethyldisiloxane↗
Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models
In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.
Crystal plasticity↗
Benchmarking of massively parallel phase-field codes for directional solidification
We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.
36 MATERIALS SCIENCE↗
Chemical timescale effects on detonation convergence
Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.
Adaptive mesh refinement↗
Nuclear microreactor transient and load-following control with deep reinforcement learning
The economic feasibility of nuclear microreactors will depend on minimizing operating costs through advancements in autonomous control, especially when these microreactors are operating alongside other types of energy systems (e.g., renewable energy). This study explores the application of deep reinforcement learning (RL) for real-time drum control in microreactors, exploring performance in regard to load-following scenarios. By leveraging a point kinetics model with thermal and xenon feedback, we first establish a baseline using a single-output RL agent, then compare it against a traditional proportional–integral–derivative (PID) controller. This study demonstrates that RL controllers, including both single- and multi-agent RL (MARL) frameworks, can achieve similar or even superior load-following performance as traditional PID control across a range of load-following scenarios. In short transients, the RL agent was able to reduce the tracking error rate in comparison to PID by one half to one third. Over extended 300-minute load-following scenarios in which xenon feedback becomes a dominant factor, PID maintained better accuracy, but RL still remained within a 1% error margin despite being trained only on short-duration scenarios. This highlights RL’s strong ability to generalize and extrapolate to longer, more complex transients, affording substantial reductions in training costs and reduced overfitting. Furthermore, when control was extended to multiple drums, MARL enabled independent drum control as well as maintained reactor symmetry constraints without sacrificing performance---an objective that standard single-agent RL could not learn. We also found that, as increasing levels of Gaussian noise were added to the power measurements, the RL controllers were able to maintain lower error rates than PID, and to do so with at least 10% and upwards of 150% less control effort. These findings illustrate RL's potential for autonomous nuclear reactor control, laying the groundwork for future integration into high-fidelity simulations and experimental validation efforts.
22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
Effect of cell compression on the performance and the structure of proton exchange membrane water electrolyzer (PEMWE) assembly
Here, in the field of water electrolysis, the proton exchange membrane water electrolyzer (PEMWE) is currently the most advanced technology for producing hydrogen without emitting CO 2 . Although PEMWE plants are already in operation, further research is needed to improve cell efficiency and reduce the use of rare materials, such as iridium oxide catalysts for the oxygen evolution reaction (OER). One of the main causes of performance loss in PEMWE is the relatively low electric conductivity of the porous transport layer (PTL) and of the anode catalyst layer, which results in ohmic losses and low catalyst utilization during high current density operation. The objective of this study is to investigate how optimization of the PTL and electrode interface can increase the cell performance. To this end, we tested different cell assemblies using fibrous and sintered PTLs, decreasing membrane thickness, reducing iridium loading, and inserting a microporous layer to increase contact surface area. Electrochemical characterization of each cell configuration was systematically performed at various compression levels as the pressure is a crucial parameter influencing the electrode/PTL contact area. In parallel, X-ray microcomputed tomography (micro-CT) was employed to investigate the effects of cell hydration and compression on the structure of PEMWE components. This study combining electrochemistry and micro-CT imaging presents how optimizing the electrode/PTL contact surface area, minimizes ohmic losses, and enables PEMWE operation with low iridium loading at high current densities.
Catalyst - PTL interface↗
What to expect when you're expecting engagement: Delivering procedural justice in large-scale solar energy deployment
Community engagement in the planning process to build large-scale solar (LSS) projects can win local support and advance procedural justice. However, an understanding of community engagement in current LSS development is lacking. Using responses from a U.S. nationwide survey (n = 979) of residential neighbors living within 3 miles (4.8 km) of completed LSS projects (i.e. “solar neighbors”) and project details from the U.S. Large-Scale Solar Photovoltaic Database (USPVDB), this study seeks to answer the following questions: How are solar neighbors' perceptions of community engagement associated with their attitudes toward their LSS projects? How do solar neighbors' perceptions of community engagement compare to their expectations? And, how do neighbors explain what they perceived about the planning process? We answer these questions using mixed methods, including regression modeling, a new gap analysis technique, and qualitative coding. We find that higher perceived engagement is associated with more positive attitudes toward the project, even when controlling for respondents who acted in opposition. Supporters and opponents alike expect more engagement than they perceived and information about projects both before construction and after operation is lacking. Awareness and engagement expectations increase at certain project size and proximity thresholds. However, most neighbors expect the public to offer input during engagement, but not make decisions. We contextualize these findings with explanatory comments from respondents.
14 SOLAR ENERGY↗
Multi-objective Decisions on Integrated Energy Systems Planning and Operation for Industrial Combined Heat and Power Supply
Unlike the power sector—which can transmit electricity over long distances via established grids—the industrial sector poses a unique challenge due to its geographically concentrated large-scale heat processes. Enhancing energy security in such industrial parks provides a dual benefit: reduced exposure to volatile fossil fuel prices and improved economic viability, largely driven by economies of scale in energy supply and distribution. This study presents a comprehensive technoeconomic analysis of a nuclear energy hub, employing a mechanism-focused approach to evaluate uncertainties in operational strategies and capacity optimization. Load profiles from three major energy intesive industries—chemical, refinery, and steelmaking—are examined, each presenting unique challenges and opportunities for nuclear energy integration. We adopt a multi-objective optimization framework, converting multiple objectives into a single objective function through the e-constraint method. The findings highlight clear trade-offs between system conditions and varying levels of energy independence. Overall, this analysis is granular enough to address industry-specific concerns yet sufficiently generalizable to provide actionable insights into the feasibility of nuclear-based clean heat solutions for the industrial sector.
25 - ENERGY STORAGE↗
Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics
A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of Zn and Ca on the deformation texture evolution of Mg–Zn–Ca alloys at elevated temperatures
This work utilized crystal plasticity simulations and experimental thermo-mechanical processing to examine how Zn and Ca content affects texture evolution in Mg–Zn–Ca alloys during simulated hot rolling. Four different compositions were studied: unalloyed Mg and three Mg alloys of ZX0p50 (Mg- 0.5 wt% Zn- 0.1 wt% Ca), ZX30 (Mg-3 wt% Zn- 0.1 wt% Ca), and ZX31 (Mg-3 wt% Zn- 0.3 wt% Ca). Multi-pass Gleeble experiments simulated hot rolling using plane strain conditions. Unalloyed Mg developed a strong basal texture, while the addition of Zn and Ca significantly weakened the texture in ZX30 and ZX31; this reduction was not seen in ZX0p50. Notably, ZX31 exhibited a split basal texture aligned with the rolling direction. Crystal plasticity simulations revealed that the weaker basal texture in ZX30 and ZX31 resulted from reduced activation of basal and twinning modes compared to unalloyed Mg and Zx0p50. Increased Zn and Ca content raised the ratios of basal to pyramidal critical resolved sheer stress (CRSS) and twin to pyramidal CRSS, enhancing pyramidal slip activity and weakening the basal texture. Extension twinning played a critical role in texture development during the first deformation pass. In unalloyed Mg and ZX0p50, twinning led to a strong basal texture, while ZX30 and ZX31 had weak basal textures that increased only slightly with further deformation. The split basal poles simulated in ZX31 were consistent with the experimental findings, highlighting the interplay between alloy composition and texture evolution.
Crystal plasticity finite element↗
Probing the elastic modulus and hardness of superionic boron cluster solid electrolytes
Polyhedral borane salts are highly tunable compounds with the potential to be used as solid electrolytes due to their high Li and Na superionic conductivity and relatively wide electrochemical stability window. In considering their application to all solid-state batteries, their mechanical properties play a critical role in their operation due to stresses the solid electrolyte is subjected to during battery operation. Density functional theory (DFT) calculations have provided an initial assessment of bulk and elastic moduli for some selected boron cluster solid electrolytes, but direct measurements are still scarce. Here, in this paper, we report the elastic moduli and hardnesses of LiCB 11 H 12 , LiCB 9 H 10 , NaCB 11 H 12 , and NaCB 9 H 10 for the first time using nanoindentation continuous stiffness method (CSM) measurements under inert atmosphere. Experimental modulus values ranging 8.8–12.3 GPa and hardness values ranging 0.15–0.38 GPa are lower than that of other inorganic solid-state electrolytes such as oxide- and sulfide-based solid electrolytes. DFT calculations on the expected moduli of these compounds are also presented and discussed. Analysis of the indentation plasticity index reveals that these compounds have higher plasticity index compared to other common oxide and sulfide solid electrolytes. According to the calculated Pugh's ratio, all compounds in this study except LiCB 11 H 12 are considered ductile.
25 ENERGY STORAGE↗
Angular-spectral filtering of recoil protons for optimization of fast neutron imaging employing proton converters
Fast neutron imaging is an important capability for diverse applications such as inertial confinement fusion diagnostics, cargo security, nuclear nonproliferation and arms control, and industrial inspection. Traditional phosphor image plates can be enhanced for fast neutron imaging using hydrogenous plastic converters which allow fast neutrons to scatter off hydrogen nuclei to produce energetic protons that can be recorded by the image plate. However, protons emitted by image plates are not constrained in their emission angle, which contributes to the blur of the resulting image. Here, we investigate two methods that can alter the spatial extent of converted protons that deposit energy in the image plate: reducing the converter thickness, and introducing a proton filter between the plastic converter and image plate to reduce the contribution of lower-energy, off-axis protons to the image. Here we determine the optimal plastic converter thickness for maximizing the signal intensity to be 2–3 mm through Monte Carlo simulations, and we benchmark this result against experimental measurements with a deuterium-tritium (DT) neutron generator. Next, we evaluate the image smearing and signal loss for various converters to show that solely reducing the converter thickness has the expected effect of reducing the blur from proton image smearing of the sharpness of an edge recorded on the image plate at the cost of reducing the signal intensity. The use of a proton filter is shown to achieve a similar improvement of edge sharpness as reducing the converter thickness while also sacrificing the signal intensity. We conclude that the use of proton energy filtering can improve the sharpness of fast neutron images in situations where the converter thickness cannot be reduced below some practical minimum. For more intense neutron sources, the signal intensity is of less concern, and optimizing the resolution of the image plate and therefore of the imaging system could have greater value. In these applications, proton filters may allow for improved fast neutron imaging measurements.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
In situ TEM annealing of neutron-irradiated Ti reveals a two-stage mechanism for elevated temperature radiation damage recovery
Understanding how irradiation-induced defects evolve at elevated temperatures is of critical importance to predicting materials' behavior under steady-state and accident scenarios. However, such mechanistic insight into microstructural evolution is limited by the nature of ex situ annealing and subsequent imaging. Here we show direct observation and quantification of defect recovery in neutron-irradiated Ti using in situ transmission electron microscopy (TEM) annealing experiments. In agreement with our prior work, and at temperatures below the irradiation temperature (T irr = 300 °C), dislocation loops are observed to glide. At elevated temperatures (>500 °C), dislocation lines become mobile and promote significant recovery of the microstructure. These mechanisms challenge the established electron irradiation-based model for radiation damage recovery, which originally suggests dissolution of static defect clusters, and demonstrates the importance of in situ characterization in understanding defect evolution in irradiated materials.
Hirst, Charles A. [Univ. of Michigan, Ann Arbor, M↗
Expanding market opportunities: cogeneration strategies for integrated PWR and thermal energy storage systems
We assess the economic viability of nuclear cogeneration by investigating three different modes—fixed dispatch, fully flexible dispatch, and flexible dispatch with minimum heat supply requirements. The analysis focuses on an existing pressurized water reactor (PWR) integrated with thermal energy storage (TES). Heat production costs are estimated under these modes for two U.S. electricity markets: the Electric Reliability Council of Texas (ERCOT) and the Pennsylvania–New Jersey–Maryland Interconnection (PJM). A sensitivity analysis examines profitability at varying heat market prices. Results indicate that fixed heat dispatch inflates heat production costs, often rendering projects economically feasible only at higher heat price levels. Fully-flexible dispatch lowers heat production costs by an average of 43 % compared to fixed dispatch. However, the current 30 % thermal dispatch limit may be insufficient to serve high baseline industrial demands cost‐effectively; higher maximum dispatch rates could enhance project economics. Markets with higher and more volatile electricity prices (e.g., ERCOT) offer greater total energy sales potential (i.e., heat and electricity), but also increase opportunity costs when heat production scheduling restrictions are imposed. In contrast, lower-price, less volatile markets (e.g., PJM) experience smaller impacts from such constraints and provide greater flexibility in accommodating varying cogeneration modes. In conclusion, these findings provide a framework to guide nuclear plant operators in aligning cogeneration strategies with industrial process requirements and electricity market conditions.
22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
Selective deuteration of an RNA:RNA complex for structural analysis using small-angle scattering
The structures of RNA:RNA complexes regulate many biological processes. Despite their importance, protein-free RNA:RNA complexes represent a tiny fraction of experimentally determined structures. Here, we describe a joint small-angle X-ray and neutron scattering (SAXS/SANS) approach to structurally interrogate conformational changes in a model RNA:RNA complex. Using SAXS, we measured the solution structures of the individual RNAs and of the overall RNA:RNA complex. With SANS, we demonstrate, as a proof of principle, that isotope labeling and contrast matching (CM) can be combined to probe the bound state structure of an RNA within a selectively deuterated RNA:RNA complex. Furthermore, we show that experimental scattering data can validate and improve predicted AlphaFold 3 RNA:RNA complex structures to reflect its solution structure. In conclusion, our work demonstrates that in silico modeling, SAXS, and CM-SANS can be used in concert to directly analyze conformational changes within RNAs when in complex, enhancing our understanding of RNA structure in functional assemblies.
HIV-1 dimerization initiation site↗
Microelectrode Arrays for Electrochemical Cycling of Individual Battery Particles
The electrochemical characterization of battery materials is usually performed on porous electrodes containing conductive additives, binders, and ensembles of redox-active particles at densities as high as 10 15 cm –3 . These complex structures often obscure individual battery particles’ intrinsic properties, hindering insights into the fundamentals of battery electrochemistry. To address this challenge, we developed microelectrode arrays to enable the charge and discharge of individual battery particles in liquid electrolytes. In this Methods and Protocols paper, we describe the fabrication of the microelectrode array chips, the assembly of particles onto microelectrodes, and the different types of electrochemical techniques that can be used, including galvanostatic cycling, cyclic voltammetry, potentiostatic intermittent titration, and electrochemical impedance spectroscopy. We use micron-sized lithium cobalt oxide as a model system to showcase the reliability and versatility of single-particle electrochemical experiments enabled by our microelectrode array. In conclusion, we aim to provide a detailed description to promote future uses of microelectrode arrays in fundamental research of batteries and other electrochemical systems.
25 ENERGY STORAGE↗