Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

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

Machine learning-enabled discovery of ionic liquid–solvent electrolytes exhibiting high ionic conductivity

Ionic liquids (ILs), which are a class of materials with versatile nature and growing popularity, are facing impediments toward widespread usage as electrolytes due to various factors such as low ionic conductivity, high viscosity, high market price etc. One of the ways these limitations can be addressed is by mixing ILs with a molecular solvent. In a combinatorial sense, there exists an immense number of specific IL–solvent combinations. An exhaustive experimental or even simulation-based investigation of the chemical space spanned by such combinations can be extremely time-consuming, expensive, and nearly impossible. An alternative approach is to employ machine learning-based models developed from available databases. Although there exists prior literature that integrates machine learning to investigate mixtures of specific solvents with ILs, these models lack generalization necessitating development of a large number of ML models to handle various solvents. To remedy this shortcoming, as a part of designing green electrolytes with high ionic conductivity that can have potential applications in next-generation batteries and solar cells, this work aims to develop a unified machine learning model to predict ionic conductivity of any IL–solvent mixture system. In this regard, three models, namely, Random Forest, extreme gradient boosting (XGBoost), and artificial neural network (ANN) were formulated using the NIST ILThermo database. The dataset contained 549 unique ionic liquids from 16 cation families and 81 unique solvents, representing a total of 23 712 datapoints. SHAPLEY additive explanation (SHAP) method was used to assess the impact of various features on model prediction and their significance was compared with literature to gain physical insight about the model behavior. Finally, using the developed models, approximately 2.5 million IL–solvent mixtures at five different compositions were screened at room temperature. The high-throughput screening yielded nearly 19 000 IL–solvent mixtures for which ionic conductivity was found to exceed the ionic conductivity of conventional Li-ion battery electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗
DOE OSTI2024

Experimental Combustion and Flame Characterization of a Chemical Looping-Based Oxidative Dehydrogenation Byproduct Fuel Mixture Containing High CO2 Dilution

Abstract This study investigates the combustion performance of a CO2-rich fuel mixture containing ethane and methane as active species using a constant volume combustion chamber. This fuel is obtained as a byproduct of a chemical looping-based oxidative dehydrogenation (Cl-ODH) process ethylene production. The byproduct gas mixture has 40.79% CO2, 39.49% ethane, and 4.88% methane by weight with other minor compounds. Using this fuel for energy extraction would improve the process efficiency of the ethane to ethylene conversion. After initial combustion modeling, the gas fuel mixture was reduced to just the major species: CO2, ethane, and methane. The mixture was then tested for flammability limits and combustion performance under spark-ignition conditions. Effects of ambient conditions like temperatures between 300 and 400 K with initial pressures from 1 to 10 bar were tested. The effects of stoichiometry were tested to understand flame velocities and heat release. The fuel mixture showed an overall reduced flame velocity compared to gasoline. Instability in combustion was believed to be caused by the dissociation of ethane under elevated conditions. At higher pressures, the flame produces lower cumulative heat release. Simulations were also performed using a model tuned to replicate the operations of the combustion chamber used in the experiments. Heat release and unburnt fuel mass data were calculated to identify the discrepancies in the combustion completeness at elevated pressures. The effects of CO2 quenching the flame coupled with the increased dissociation of the fuel species can lead to up to more than 75% of the fuel mixture being unburnt. Data from this study were used to modify a small-scale spark-ignition engine to use this fuel and produce usable energy.

Energy & Fuels↗
DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

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

EXPERIMENTAL AND MODELING STUDY OF MOLECULAR TRANSPORT PHENOMENA IN COMPLEX LIQUID MIXTURES USING COVALENT ORGANIC FRAMEWORK MEMBRANES

Currently, 10-15% of the world’s energy is consumed by chemical separations, and more than 80% of that is used to purify organic liquids and recover rare earth elements and critical minerals. Membrane nanofiltration presents an energy-efficient, cost-effective, and eco-friendly alternative to current separation technologies. These complex liquid environments require high- performance, chemically-resistant membrane materials. Recently discovered Covalent Organic Frameworks (COFs) possess desirable properties as membrane materials for complex liquid separations. COFs are highly crystalline, chemically and thermally stable, with tunable size and charge. COFs, especially two-dimensional highly crystalline COFs, possess narrow pore-size distribution and controlled porous structures. Two-dimensional COFs are also resistant to swelling, another feature beneficial in complex liquid separations. Molecular interactions in complex liquid systems often dictate final membrane performance, resulting in a discrepancy between designed and apparent COF properties. Understanding and resolving this discrepancy is critical in utilizing COF membranes as a viable alternative to energy-intensive separations. The primary objective of this thesis was to use a commercially available COF TpPa-1 as a platform to investigate how mixed solvent environment affects COF membrane performance via experimental and molecular modeling studies. Specifically, the effects of solution pH, solvent- solvent-solute interactions in mixed solvents, and COF chemistry were carefully examined on apparent TpPa-1 pore size and permeability and targeted solute rejection of COF membranes for neat and mixed solvents. Experimental COF filtration performance was compared and corroborated with modeling predication. Findings from this study will provide guidance for future COF design, synthesis, and desired functional COF membrane performance.

Barnes, Anastasia↗
DOE OSTI2023

Monte Carlo goodness-of-fit tests for degree corrected and related stochastic blockmodels

Abstract We construct Bayesian and frequentist finite-sample goodness-of-fit tests for three different variants of the stochastic blockmodel for network data. Since all of the stochastic blockmodel variants are log-linear in form when block assignments are known, the tests for the latent block model versions combine a block membership estimator with the algebraic statistics machinery for testing goodness-of-fit in log-linear models. We describe Markov bases and marginal polytopes of the variants of the stochastic blockmodel and discuss how both facilitate the development of goodness-of-fit tests and understanding of model behaviour. The general testing methodology developed here extends to any finite mixture of log-linear models on discrete data, and as such is the first application of the algebraic statistics machinery for latent-variable models.

Karwa, Vishesh↗
DOE OSTI2021

A generalized distillation network synthesis model

In this report we propose a generalized superstructure-based distillation network synthesis model with improved modeling capabilities resulting in an extended solution space. The model can assign multiple mixtures to be separated (e.g., reactor network effluents) to different columns of the network while considering the interactions among different separation steps. In terms of outlets, products with general specifications, including pure components and multi-component mixtures, as well as streams without strict specifications (e.g., reactor recycle streams) can be readily handled. Furthermore, stream bypass is considered to avoid unnecessary separations, and thermal coupling is considered to yield novel solutions. The proposed model enables seamless integration with reactor network synthesis models thereby leading to solutions that are superior to the ones obtained by conventional sequential approaches.

42 ENGINEERING↗
DOE OSTI2024

Glass Design Using Machine Learning Property Models with Prediction Uncertainties: Nuclear Waste Glass Formulation

The United States Department of Energy is responsible for managing the legacy nuclear waste stored in underground tanks at the Hanford Site. The waste will be separately vitrified as low-activity waste and high-level waste fractions. Waste glass formulation algorithms have been traditionally developed using partial quadratic mixture property-composition models. Recently, machine learning (ML) techniques have been used to predict glass properties and discover new glass materials for nuclear waste vitrification, and these advancements can be utilized to improve waste glass composition design. In this proof-of-principle study, ML algorithms such as Gaussian process regression (GPR) were used to interpolate glass properties (e.g., viscosity, electrical conductivity, chemical durability). After selecting appropriate sets of GPR hyper-parameters for each property, an optimization program was developed to formulate glass compositions to maximize waste loading while simultaneously satisfying property within constraints. The results of the ML-based waste loadings and glass compositions were compared to those obtained using the traditional methods. Comparing to the previous glass design framework, the ML-based optimization methods offer improved glass designs and a streamlined approach to generation of optimally designed data and near real-time updates.

glass formulation, machine learning, constraints, ↗
DOE OSTI2021

Phase Separation in Ultramassive White Dwarfs

We report ultramassive white dwarfs are extreme endpoints of stellar evolution. Recent findings, such as a missing multi-Gyr cooling delay for a number of ultramassive white dwarfs and a white dwarf with a quasi-Chandrasekhar mass, motivate a better understanding of their evolution. A key process still subject to important uncertainties is the crystallization of their dense cores, which are generally assumed to be constituted of 16 O, 20 Ne, and a mixture of several trace elements (most notably 23 Na and 24 Mg). In this work, we use our recently developed Clapeyron integration technique to compute accurate phase diagrams of three-component mixtures relevant to the modeling of O/Ne ultramassive white dwarfs. We show that, unlike the phase separation of 22 Ne impurities in C/O cores, the phase separation of 23 Na impurities in O/Ne white dwarfs cannot lead to the enrichment of their cores in 23 Na via a distillation process. This severely limits the prospect of transporting large quantities of 23 Na toward the center of the star, as needed in the white dwarf core-collapse mechanism recently proposed by Caiazzo et al. We also show that despite representing ≈10% of the ionic mixture, 23 Na and 24 Mg impurities only have a negligible impact on the O/Ne phase diagram, and the two-component O/Ne phase diagram can be safely used in white dwarf evolution codes. We provide analytic fits to our high-accuracy O/Ne phase diagram for implementation in white dwarf models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Scrambling and operator entanglement in local non-Hermitian quantum systems

The breakdown of Lieb-Robinson bounds in local, non-Hermitian quantum systems opens up the possibility for a rich landscape of quantum many-body phenomenology. We elucidate this by studying information scrambling and quantum chaos in non-Hermitian variants of paradigmatic local quantum spin-chain models. We utilize a mixture of exact diagonalization and tensor network techniques for our numerical results and focus on three dynamical quantities: (i) out-of-time-ordered correlators (OTOCs), (ii) operator entanglement of the dynamics, and (iii) entanglement growth following a quench from product initial states. We show that while OTOCs fail to capture information scrambling in a simple, local, non-Hermitian transverse-field Ising model, the closely related operator entanglement is a robust measure of dynamical properties of interest. Moreover, we show that the short-time growth of operator entanglement can generically detect entanglement phase transitions in these systems while its long-time average is shown to be a reliable indicator of quantum chaos and entanglement phases. This allows us to extend operator entanglement-based diagnostics from previous works on closed and open quantum systems, to the new arena of monitored quantum dynamics. Finally, we remark on the efficacy of these dynamical quantities in detecting integrability/chaos in the presence of continuous monitoring.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Measurements and kinetic modeling of O 2 vibrational kinetics in O 2 –Ar mixtures partially dissociated by a Ns pulse discharge

Vibrational kinetics of O 2 is studied during the O atom recombination in an O 2 –Ar mixture, partially dissociated by a burst of ns discharge pulses in a heated plasma flow reactor. The time-resolved temperature in the discharge afterglow is determined by Rayleigh scattering. Time-resolved O atom number density is measured by ps Two-Photon absorption Laser Induced Fluorescence, calibrated in xenon. Time-resolved vibrational level populations of molecular oxygen, O 2 (v= 8–20), are measured by ps Laser Induced Fluorescence (LIF), with the absolute calibration by NO LIF. Time-resolved ozone number density is monitored by broadband UV absorption. The results are compared with the predictions of a state-specific kinetic model. The experimental data indicate a rapid initial decay of O 2 (v) populations generated by electron impact in the discharge, due to the vibration-translation (V–T) relaxation by O atoms. This is followed by a slower population reduction, on the time scale much longer compared to that for V–T relaxation or vibration-vibration (V–V) exchange. Both O atoms and the O 2 (v) populations decay on the same time scale, indicating that chemical reactions initiated by the O atom recombination result in the generation of vibrationally excited O 2 molecules. These trends are reproduced by the kinetic model, which shows that the reaction of O atoms with ozone is the dominant pathway of O 2 (v) generation at the present conditions. The predicted relative O 2 (v) populations are close to the experimental results, but absolute number densities differ from the experimental data. This is likely due to uncertainties in the absolute calibration of LIF measurements and in the spectroscopic model used in the data reduction. The present work demonstrates the capability for the absolute, time-resolved measurements of vibrationally excited O 2 in recombining gas flows, to quantify the energy partition in the recombination reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Polychlorinated biphenyls, polychlorinated dibenzo- p -dioxins, polychlorinated dibenzofurans, pesticides, and diabetes in the Anniston Community Health Survey follow-up (ACHS II): single exposure and mixture analysis approaches

Dioxins and dioxin-like compounds measurements were added to polychlorinated biphenyls (PCBs) and organochlorine pesticides to expand the exposure profile in a follow-up to the Anniston Community Health Survey (ACHS II, 2014) and to study diabetes associations. Participants of ACHS I (2005–2007) still living within the study area were eligible to participate in ACHS II. Diabetes status (type-2) was determined by a doctor's diagnosis, fasting glucose ≥125 mg/dL, or being on any glycemic control medication. Incident diabetes cases were identified in ACHS II among those who did not have diabetes in ACHS I, using the same criteria. Thirty-five ortho-substituted PCBs, 6 pesticides, 7 polychlorinated dibenzo-p-dioxins (PCDD), 10 furans (PCDF), and 3 non-ortho PCBs were measured in 338 ACHS II participants. Dioxin toxic equivalents (TEQs) were calculated for all dioxin-like compounds. Main analyses used logistic regression models to calculate odds ratios (OR) and 95 % confidence intervals (CI). In models adjusted for age, race, sex, BMI, total lipids, family history of diabetes, and taking lipid lowering medication, the highest ORs for diabetes were observed for PCDD TEQ: 3.61 (95 % CI: 1.04, 12.46), dichloro-diphenyl dichloroethylene (p,p’-DDE): 2.07 (95 % CI 1.08, 3.97), and trans-Nonachlor: 2.55 (95 % CI 0.93, 7.02). The OR for sum 35 PCBs was 1.22 (95 % CI: 0.58–2.57). To complement the main analyses, we used BKMR and g-computation models to evaluate 12 mixture components including 4 TEQs, 2 PCB subsets and 6 pesticides; suggestive positive associations for the joint effect of the mixture analyses resulted in ORs of 1.40 (95% CI: -1.13, 3.93) for BKMR and 1.32 (95% CI: -1.12, 3.76) for g-computation. Furthermore, the mixture analyses provide further support to previously observed associations of trans-Nonachlor, p,p’- DDE, PCDD TEQ and some PCB groups with diabetes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Behavior of astatine and bismuth in non-conventional solvents: Extraction into imidazolium-based ionic liquid and methyl anthranilate with active pharmaceuticals binary mixtures from nitric acid media

Astatine is one of the least chemically studied elements and its behavior in the presence of non- conventional solvents has not been investigated before. This work considers both hydrophobic ionic liquids and binary mixtures as alternatives to conventional solvents. The study is based on the extraction of astatine and bismuth (target material required to produce astatine) into imidazolium-based ionic liquid, and binary mixtures formed by active pharmaceuticals (ibuprofen, lidocaine) and a food grade ingredient (methyl anthranilate). It is shown that both the ionic liquid and binary mixtures can successfully extract At from nitric acid media, but extraction of Bi into the ionic liquid is very inefficient, resulting in a good separation factor for these two elements in the entire studied acidity range. Extraction of At into binary mixtures is very efficient, having distribution ratio values as high as 1000, while the behavior of Bi under these conditions depends on the composition of the mixture. Furthermore, a mathematical model has been developed to fit both At and Bi experimental data and applied to determine corresponding thermodynamic extraction constants.

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

Polymer principles behind solubilizing lignin with organic cosolvents for bioenergy

Lignin solubilization is key to a viable biorefinery because its removal leads to facile deconstruction of biomass and because the isolated lignin can serve to derive precursors of novel high-value materials. The mixing of organic solvents with water has been shown to improve biomass fractionation and lignin conversion reactions. However, generally-applicable solubilization strategies are lacking because of the remarkable variability of lignin across plant feedstocks. Here, to obtain a predictive understanding of lignin solvation, we perform molecular dynamics simulations of model lignin polymers in two mixtures of water with polar aprotic solvents: tetrahydrofuran (THF) : water and γ-valerolactone (GVL) : water. The model lignins include H-, G- and S-only homopolymers and lignins with an S : G 1 : 1 ratio. We find that a well-established theory of self-avoiding polymers in a “good” solvent describes accurately the physical conformations of all types of lignin in both solvents. As the degree of methoxy substitution increases in the homopolymers, the distributions of the lignin radius of gyration and the Flory exponent ν , which describes the lignin-solvent interactions, do not change in THF : water, while ν shifts to slightly higher values in GVL : water. We attribute this increase to the interaction between the methyl group of GVL with the lignin methoxy groups. We also find that the reduction in the lignin radius of gyration due to branching is accurately described by the Zimm–Stockmayer theory for both THF : water and GVL : water. The above findings validate the applicability of polymer physics concepts to lignin and suggest that GVL : water may have the most favorable interaction with S-lignin, whereas the interactions of THF : water with lignin are independent of lignin monomeric content.

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

Modeling and simulation of gas separations with spiral–wound membranes

Models for gas separations with spiral–wound membranes are developed and found to exhibit good agreement with experiments performed on N 2 /O 2 mixtures. The two–dimensional (2D) model can be accurately approximated by a one–dimensional (1D) surrogate model when the spacer widths are chosen to make the channel pressure drops small. Subsequently, the separation of propane/propylene mixtures from the recycle purge stream of a polypropylene reactor is investigated. Assuming ideal gas is found to lead to significant overestimations in membrane stage cuts (sometimes more than 10%), an extent comparable to that associated with extrapolating constant olefin permeance from a low–pressure condition. While olefin permeance can change significantly with pressure, using a constant–permeance formulation can result in a small (< 2.5%) underprediction in stage cut if the value for the permeance is taken from the feed condition. Finally, membrane properties and costs necessary for a viable separation process are discussed.

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

Screening of metal complexes and organic solvents using the COSMOSAC-LANL model to enhance the energy density in a non-aqueous redox flow cell: an insight into the solubility

In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.

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

SAM Enhancements for Air-Ingress Event Modeling in HTGRs

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes recent efforts and progress in addressing the code capability gaps in SAM for the modeling of the air-ingress phenomenon in High-Temperature Gas-Cooled Reactors (HTGRs). The capability enhancements implemented to the code include: a multi-component flow model to capture the air-helium mixture during air ingress, a 0-D graphite oxidation model to capture the reaction of graphite with oxygen in the air, and an isentropic critical flow model to more accurately predict the de-pressurization of the reactor system due to a small break in the primary loop. Verification cases and demonstrations are provided to showcase these capabilities in SAM. Existing gaps in the code’s capability are also identified.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗