Molten Salt Species Solubility Calculations Using Thermochimica's Gibbs Energy Minimizer (GEM)
Explore the source record for details and available documents.
Lee, Kyoung [ORNL]↗
SEARCH · Search NASA
Results for “Gibbs energy minimization”
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.
Acceleration of Thermochemistry Solves in MOOSE and Pronghorn
This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.
36 - MATERIALS SCIENCE↗
Equilipy: a python package for calculating phase equilibria
The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.
97 MATHEMATICS AND COMPUTING↗
Specification of Radionuclide Classes for MELCOR Molten Salt Reactor Simulations
MELCOR has been used extensively to facilitate virtual investigations into severe nuclear accidents for light-water reactors (LWRs). Non-light water reactors (non-LWRs) render some LWR-centric approaches potentially unsuitable. MELCOR has been instrumental in analyzing source terms for LWRs and has recently expanded its applicability to non-LWRs. To simplify radionuclide (RN) tracking, MELCOR currently groups elements into 17 classes, each containing representative species. This grouping, optimized for LWRs, is not appropriate for non-LWRs due to the different chemistry. This necessitates a reevaluation of radionuclide transport modeling. This report introduces a new class scheme for MELCOR tailored to MSR modeling, expanding the current 17 classes to 32 and are explained in the context of a UF 4 fueled FLiBe carrier MSR. It provides a discussion and justification for the new groupings and outlines a methodology for discovering and defining additional classes in MELCOR using a sample calculated RN inventory and a Gibbs energy minimizer (GEM).
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
High-Fidelity Multiphysics Assessment of the Molten Chloride Reactor Experiment (MCRE) Integrating Neutronics, Thermal-hydraulics and Thermochemistry
This work presents a MOOSE-based Multiphysics model of the LOTUS MSR reactor system. The model includes the coarse-mesh turbulent thermal-hydraulics including passive advection of delayed neutron precursors of the primary system through Pronghorn, the power density and neutronics calculation employing Griffin, and the chemical interactions through Gibbs energy minimization and redox potential with Thermochimica. Incorporating a plate-out model, the full multi-physics model is used to evaluate the effects of deposition and removal of solid compounds formed in the molten salt in the reactor surfaces, evaluating the concentrations and regions in which nickel, chromium and iron would be deposited. Results indicate an almost isothermal state of the fuel during operation at 25 kW that do not increase the plate-out effects due to low temperature gradients.
21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
Flux trapping in NbTiN strips
We use scanning superconducting quantum interference device (SQUID) microscopy to image vortices in superconducting structures fabricated from NbTiN thin films. We repeatedly cool superconducting strips with different width in an applied magnetic field and image the individual vortices. From these images we determine the threshold field at which the first vortex enters a strip and the number and configuration of vortices beyond this threshold field. We model the behavior of the vortices with and without considering the effect of pinning using numerical methods to minimize the Gibbs free energy of vortices in the strips. Our measurements provide a first benchmark to investigate the flux trapping properties of NbTiN thin films directly relevant to NbTiN-based superconducting digital circuits.
Bai, Ruiheng [Cornell University] (ORCID:000000025↗
Thermal Weight Determination and Interstate Coupling in State-Averaged ADAPT-VQE
Characterizing electronic thermal states at low temperatures is an important but challenging task in quantum chemistry and condensed matter physics, making it a prime candidate for a useful application in quantum computing. One of the most successful methods for state preparation on quantum computers is the Adaptive, Problem-Tailored (ADAPT) Variational Quantum Eigensolver (VQE), which has recently been generalized to treat excited states within a state-averaged framework as well as Gibbs states. In this work, we introduce Helmholtz-Optimized Thermal (HOT) ADAPT-VQE, an ancilla-free strategy for preparing Gibbs states that directly minimizes the Helmholtz free energy by targeting the dominant eigenstates of the thermal ensemble. We demonstrate the usefulness of HOT-ADAPT-VQE by predicting the free energy of two model systems with strongly correlated ground states: (1) the Fe 2+ cation in a magnetic field and (2) a [Cu 2 O 7 ] 10– fragment of the Mott insulator La 2 CuO 4 . Our results demonstrate that HOT-ADAPT-VQE significantly improves upon Gibbs-state estimates from multistate variants of ADAPT-VQE, often with substantially shallower quantum circuits, making it a promising candidate for thermal-state calculations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Electrically Powered High-Salinity Brine Separation Using Dimethyl Ether
Dewatering highly saline aqueous streams, from mining and geothermal leachates to industrial wastewater, is essential for effective resource recovery and safe disposal. Membraneless water extraction (MWE) uses a low-polarity solvent to separate water from concentrated aqueous solutions. In this study, we design a new MWE that uses dimethyl ether (DME) to selectively extract water from high-salinity brines, leveraging the volatility of DME to achieve rapid solvent recovery. By separating water and dissolved salts at a liquid–liquid interface, MWE minimizes the deleterious effects of scaling on vulnerable membrane and heat exchanger surfaces, reducing the need for extensive pretreatment and expensive materials. We begin by developing a computational framework for a multistage counterflow liquid–liquid contactor, which extracts water into DME, coupled with a multistage solvent regenerator that uses vapor compression to efficiently separate the desalinated water from the DME extractant. Excess Gibbs free energy and equation of state frameworks are used to model fluid phase equilibria in water–DME–sodium chloride (NaCl) mixtures, with interaction parameters estimated from experimental data. Incorporating equilibrium calculations into a system-scale computational model, we examine the performance of MWE using DME for the first time. Our analysis demonstrates that MWE can concentrate seawater desalination brine (>1.0 mol NaCl kg –1 ) to zero-liquid discharge salinities, with an energy consumption of under 50 kW h per m 3 of water extracted with a solvent recovery ratio greater than 99.9%. We highlight the importance of staging the vapor compression process to simultaneously minimize energy consumption while enabling brine concentration and product water solvent contamination. Finally, the thermodynamic framework developed here allows for the robust evaluation of new MWE solvents and systems for critical brine concentration and fractional precipitation applications.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Competitive Sorption of H 2 O and CO 2 in Clay Mineral Interlayers
The first-reported competitive CO 2 –H 2 O sorption isotherms on a naturally occurring expandable clay mineral (saponite) exchanged with Na, K, and Cs at near ambient CO 2 concentrations (1% CO 2 ) are reported by combining ATR-FTIR spectroscopy with gravimetric sorption methods. Little CO 2 sorption was observed at high relative humidity (RH) but increased significantly at low RH. Unlike prior work where CO 2 sorption was reported on samples that were previously subjected to heating and/or evacuation treatment to remove sorbed H 2 O, this study manipulated CO 2 sorption by controlling RH. CO 2 sorption increased in the order Cs- > K- > Na-exchanged saponite and was anticorrelated with the Gibbs energy of hydration of the exchangeable cation. As expected, H 2 O sorption followed the opposite trend, with increased H 2 O sorption occurring in the order Na > K > Cs. The amount of CO 2 sorbed using 1 atm CO 2 (p/p o = 1.0) ranged from 50 to 30 mg CO 2 /g clay and from 1.25 to 0.8 mg CO2 /g clay for 1% CO 2 (p/p o = 0.01) in N 2 . The sensitivity of the ATR-FTIR measurements is demonstrated by detecting the CO 2 stretch of adsorbed CO 2 from air containing 500 ppmv CO 2 . The spectral features of sorbed CO 2 were minimally affected by the nature of the exchangeable cation or variations in H 2 O content. However, the position of the HOH bending band shifted significantly from 1638 to 1610 cm –1 as H 2 O content decreased, reflecting changes in intermolecular hydrogen bonding between H 2 O molecules in the saponite interlayer. The shift in position of the HOH bending band to lower energy was coincident with increased CO 2 sorption. CO 2 sorption in clay mineral interlayers is consistent with sorption on weakly hydrated partially hydrophobic sites found on low charge density smectites and is enhanced by lower RH and by the presence of weakly hydrated exchangeable cations like Cs + or K + . In conclusion, the greatest CO 2 sorption occurs on Cs-saponite because this is the least hydrated cation of the three used in the study and indicates that CO2 sorption occurs on the neutral portion of the siloxane surface.
Johnston, Cliff T. [Purdue Univ., West Lafayette, ↗
Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics
Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Copper-Doped Tin Oxides Supported on Mesoporous Carbon Xerogel for Boosting the Electrochemical Reduction of CO 2 to Formate in Bicarbonate Solution Coupled with CO 2
The electrochemical reduction of CO 2 into valuable products at mild reaction conditions and using cheap renewable electrical energy are goals to sustain a low-carbon economy. Among the various CO 2 reduction reaction products, formic acid (FA) has received significant attention because of its low Gibbs free energy input requirement and the simple reduction reaction involving the transfer of 2 electrons and 2 protons. In this work, a copper-doped tin oxide catalyst supported on a mesoporous carbon xerogel was shown to enhance the electrochemical reduction of CO 2 to formate in a bicarbonate solution coupled with CO 2 . We observed that the synergistic SnCu oxides enhance the selectivity toward formate from 58.6% for Sn oxide and 28.7% for Cu oxide to over 71.2% for the SnCu oxides. The observed rate of formate production with SnCu oxide was 2.8 times higher compared to the rate of Cu oxide and about 1.5 times higher than with the Sn oxide catalyst. Our results reveal that selectivity for formate comes partly from the electrolysis of the bicarbonate solution and partly from continuous CO 2 gas purged into the solution. The contribution from the electrolysis of bicarbonate solution ranges from 15% to 40% when the concentration of bicarbonate solution ranges from 100 mM to 1 M. Chronoamperometric measurements for stability revealed that Cu oxide and Sn oxide showed stable current density for less than 30 h while under the same conditions, the stable current density was observed for more than 50 h with SnCu oxide catalyst. Additionally, the selectivity toward formate increased by 6% when the reactor pressure was increased from near ambient pressure to 4 psig. Our lab-scale electrochemical cell with SnCu oxide supported on the mesoporous carbon xerogel enhances the CO 2 solubility, minimizes the precipitation of salts that can degrade the catalytic performance, and suppresses the competitive hydrogen evolution reaction, demonstrating the feasibility of using our catalyst and system for the electrochemical conversion of CO 2 into formate with high selectivity, productivity, and stability. Furthermore, this could have significant implications for the mitigation of CO 2 emissions and the development of a sustainable chemical industry.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Integrated Framework of Multisource Data Fusion for Outage Location in Looped Distribution Systems
Accurate outage location is essential for expediting post-outage power restoration, minimizing outage duration, and enhancing the resilience of distribution networks. With the advent of advanced metering infrastructure, data-driven outage location methods have significantly advanced beyond traditional approaches that rely on manual inspections. However, existing methods still face critical challenges, like reliance on single-source data, limited ability to handle partially observable systems or difficulties with loop networks. To the best of our knowledge, no single approach has comprehensively addressed all of these challenges at once. To this end, this paper proposes a comprehensive multisource data fusion framework for outage locations via probabilistic graph networks. The framework consists of three key phases. First, a novel method for reconstituting distribution networks with loops is developed, transforming looped networks into multiple radial subnetworks that retain all outage causalities of the original network. Second, Bayesian network (BN) models are established for each subnetwork, integrating multiple data sources and network structures. Finally, a joint Gibbs sampling mechanism, featuring forward and backward information flow, is designed to merge data from separate BN models and maximize the utilization of limited evidence, ensuring accurate outage location identification. In conclusion, the framework was validated on two modified public test systems, and comparative studies confirmed its effectiveness.
24 POWER TRANSMISSION AND DISTRIBUTION↗