Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical solution”

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

Chain conformations and phase separation in polymer solutions with varying solvent quality

Molecular dynamics simulations are used to investigate the conformations of a single polymer chain, represented by the Kremer-Grest bead-spring model, in a solution with a Lennard-Jones liquid as the solvent when the interaction strength between the polymer and solvent is varied. Results show that when the polymer-solvent interaction is unfavorable, the chain collapses as one would expect in a poor solvent. For more attractive polymer-solvent interactions, the solvent quality improves and the chain is increasingly solvated and exhibits ideal and then swollen conformations. However, as the polymer-solvent interaction strength is increased further to be more than about twice the strength of the polymer-polymer and solvent-solvent interactions, the chain exhibits an unexpected collapsing behavior. Correspondingly, for strong polymer-solvent attractions, phase separation is observed in the solutions of multiple chains. These results indicate that the solvent becomes effectively poor again at very attractive polymer-solvent interactions. Nonetheless, the mechanism of chain collapsing and phase separation in this limit differs from the case with a poor solvent rendered by unfavorable polymer-solvent interactions. In the latter, the solvent is excluded from the domain of the collapsed chains while in the former, the solvent is still present in the pervaded volume of a collapsed chain or in the polymer-rich domain that phase separates from the pure solvent. Finally, in the limit of strong polymer-solvent attractions, the solvent behaves as a glue to stick monomers together, causing a single chain to collapse and multiple chains to aggregate and phase separate.

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

Subtle Molecular Changes Largely Modulate Chiral Helical Assemblies of Achiral Conjugated Polymers by Tuning Solution-State Aggregation

Understanding the solution-state aggregate structure and the consequent hierarchical assembly of conjugated polymers is crucial for controlling multiscale morphologies during solid thin-film deposition and the resultant electronic properties. However, it remains challenging to comprehend detailed solution aggregate structures of conjugated polymers, let alone their chiral assembly due to the complex aggregation behavior. Herein, we present solution-state aggregate structures and their impact on hierarchical chiral helical assembly using an achiral diketopyrrolopyrrole-quaterthiophene (DPP-T4) copolymer and its two close structural analogues wherein the bithiophene is functionalized with methyl groups (DPP-T2M2) or fluorine atoms (DPP-T2F2). Combining in-depth small-angle X-ray scattering analysis with various microscopic solution imaging techniques, we find distinct aggregate in each DPP solution: (i) semicrystalline 1D fiber aggregates of DPP-T2F2 with a strongly bound internal structure, (ii) semicrystalline 1D fiber aggregates of DPP-T2M2 with a weakly bound internal structure, and (iii) highly crystalline 2D sheet aggregates of DPP-T4. These nanoscopic aggregates develop into lyotropic chiral helical liquid crystal (LC) mesophases at high solution concentrations. Intriguingly, the dimensionality of solution aggregates largely modulates hierarchical chiral helical pitches across nanoscopic to micrometer scales, with the more rigid 2D sheet aggregate of DPP-T4 creating much larger pitch length than the more flexible 1D fiber aggregates. Combining relatively small helical pitch with long-range order, the striped twist-bent mesophase of DPP-T2F2 composed of highly ordered, more rigid 1D fiber aggregate exhibits an anisotropic dissymmetry factor (g-factor) as high as 0.09. This study can be a prominent addition to our knowledge on a solution-state hierarchical assembly of conjugated polymers and, in particular, chiral helical assembly of achiral organic semiconductors that can catalyze an emerging field of chiral (opto)electronics.

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

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A General Binary Isotherm Model for Amines Interacting with CO 2 and H 2 O

CO 2 capture by primary or secondary amines has been of great research interests for a century because of its industrial importance. Interest has grown even more, because of the need to eliminate the CO 2 emissions that drive global warming. Experimental evidence shows that CO 2 sorption in primary or secondary amines is accompanied by co-absorption of H 2 O. A quantitative analysis of such CO 2 -H 2 O co-absorption behavior is important for practical process design and theoretical understanding. Even though there is almost an experimental consensus that water enhances CO 2 uptake capacity, an analytic model to explain this phenomenon is not well established. Instead, some empirical models such as the Toth model are used to describe the isotherm without accounting for the presence of water. Recently, we have demonstrated that the isotherm equation of CO 2 sorption into strong-base anion exchange materials with quaternary ammonium can be derived from that of strong-base aqueous alkaline solutions by correcting for the drastic change in water activity and by including an appropriate parameterization of the water activity terms. In this paper, we generalize this model from quaternary ammonium to primary, secondary and tertiary amines either in solutions or as functional groups in polymer resins. For primary, secondary and tertiary amines, the isotherm equation can be derived by extending that of a weak-base aqueous alkaline solution such as aqueous ammonia. The model has been validated using experimental data of aqueous ammonia in literature. Furthermore, this general model even includes quaternary ammonium as a special limit. Hence, this general model offers a platform that can treat the isotherms of solid amines, aqueous amines and aqueous alkaline solutions in a unified way.

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

Nonlinear Anti-(Parity-Time) Symmetric Dimer

In the present work we propose a nonlinear anti- P T -symmetric dimer, that at the linear level has been experimentally created in the realm of electric circuit resonators. We find four families of solutions, the so-called upper and lower branches, both in a symmetric and in an asymmetric (symmetry-broken) form. We unveil analytically and confirm numerically the critical thresholds for the existence of such branches and explore the bifurcations (such as saddle-node ones) that delimit their existence, as well as transcritical ones that lead to their potential exchange of stability. We find that out of the four relevant branches, only one, the upper symmetric branch, corresponds to a spectrally and dynamically robust solution. We subsequently leverage detailed direct numerical computations in order to explore the dynamics of the different states, corroborating our spectral analysis results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

FY25 Electric Grid Security Annual Report

Sandia’s Electric Grid Security program advances a national vision of energy dominance and accessibility, while applying our national security -emphasis on ensuring of a secure, resilient, and affordable electric system for all users. Our achievements reflect a strategic approach combining technology development; modeling, simulation, and data analytics; and partnered demonstrations and outreach to further the adoption of advanced grid and storage technologies. Our FY25 efforts leverage the strengths of our partnerships—spanning Sandia’s core science and technology competencies as well as external technology leaders—to develop the solutions today which enable the grid of tomorrow. Key accomplishments in this report that support our strategy span our technical program areas and include: • New open-source analytical tools for systems -level planning and optimization, including significant advances to the QuESt analytical environment; • Further advancement of artificial intelligence and machine learning to enhanced grid operations and planning as we rise to the challenge of new large loads; • Development of solid-state power conversion technologies and a new medium-voltage research lab; • New technologies to assess wildfire vulnerabilities and mitigate potential impacts; • Advanced applications of new cybersecurity technologies with industry partners; • Contributions to understanding the impacts of electromagnetic pulses and geomagnetic disturbances on grid components; and • Digital twin development for hybrid microgrids with multiple generators, storage, and loads. This report indicates key areas of research and engagement and summarizes the impact of Sandia’s contributions through notable accomplishments, journal publications, patents, and technical conferences and presentations. It is provided with the hope that readers discover ways we can further team to create our modern grid and apply the outcomes of our efforts. The bulk of work described herein is funded by several offices within the U.S. Department of Energy (USDOE), including the Office of Electricity (OE); Cybersecurity, Energy Security, and Emergency Response (CESER); former offices such as the Office of Energy Efficiency and Renewable Energy (EERE), the Grid Deployment Office (GDO), the Office of Clean Energy Demonstrations (OCED), and other key programs at USDOE. As we continue to state in these annual reports, the contributors to our successes are too numerous to name here, though our team wishes to express our deep gratitude to the numerous program and project sponsors at the US Department of Energy, who often function equally as technical collaborators; our many partners in industry, academia, utilities, and other national labs; and fellow researchers and business partners at Sandia whose leadership and creativity have enabled the accomplishments described herein.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Out-of-time-ordered-correlators for the pure inverted quartic oscillator: classical chaos meets quantum stability

Out-of-time-ordered-correlators (OTOCs) have been suggested as a means to diagnose chaotic behavior in quantum mechanical systems. Recently, it was found that OTOCs display exponential growth for the inverted quantum harmonic oscillator, mirroring the fact that this system is classically and quantum mechanically unstable. In this work, I study OTOCs for the inverted anharmonic (pure quartic) oscillator in quantum mechanics, finding only oscillatory behavior despite the classically unstable nature of the system. For higher temperature, OTOCs seem to exhibit saturation consistent with a value of –2< x 2 > T < p 2 > T at late times. I provide analytic evidence from the spectral zeta-function and the WKB method as well as direct numerical solutions of the Schrödinger equation that the inverted quartic oscillator possesses a real and positive energy eigenspectrum, and normalizable wave-functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Detection of Br and I as atomic anions using liquid sampling – atmospheric pressure glow discharge/Orbitrap mass spectrometry

The quantification and determination of halogen isotope composition is essential in many fields including geochemistry and nuclear forensics. Detection of the halogen elements is commonly attempted with inductively coupled plasma mass spectrometry (ICP-MS) however, there are multiple challenges faced. Most significantly this includes very poor ionization efficiency and the potential for isobaric interferences from either matrix or plasma species. Thus, sensitivity and accuracy are often limiting issues. Here a liquid sampling-atmospheric pressure glow discharge (LS-APGD) ionization source is coupled with an ultrahigh resolution Orbitrap MS to perform halogen detection of bromine and iodine as initial analytes. This facilitates simple and sensitive detection of these elements as atomic anions from simple aqueous salt solutions. Collision-induced dissociation (CID) and higher-energy collisional dissociation (HCD) energies were optimized to produce the maximum response of Br − and I − . The LS-APGD conditions were optimized using a design of experiments (DOE) approach for Br − and I − concurrently. Response curves for Br − and I − solutions determined limit of detection (LOD) values of 50 pg and 5 pg, respectively, in 20 μL aliquots. The curves indicated response factors of 0.67 for 79 Br − and 0.90 for 127 I − . Br isotope ratios were determined with precision between 0.7 and 7.3% RSD, with the isotope ratios determined with precision 0.8% RSD at concentrations above the limit of quantification. Preliminary tests were conducted to evaluate the effect of different cations (Na + , K + , and Mg 2+ ) on Br − responses, with no discernible impacts observed on the halogen signal responses. This study demonstrates the potential use of the LS-APGD-Orbitrap-MS for the detection of multiple halogens while avoiding interferences, minimizing sample preparation, and overcoming ionization barriers.

Lang, Dehlia A. [Clemson University, SC (United St↗
DOE OSTI2021

ITER cold VDEs in the limit of perfectly conducting walls

Recently, it has been shown that a vertical displacement event (VDE) can occur in ITER even when the walls are perfect conductors, as a consequence of the current quench (CQ) [A. H. Boozer, Phys. Plasmas 26, 114501 (2019)]. We used the extended-MHD code M3D-C1 with an ITER-like equilibrium and induced a CQ to explore cold VDEs in the limit of perfectly conducting walls, using different wall geometries. In the case of a rectangular first wall with the side walls far away from the plasma, we obtained very good agreement with the analytical model developed by Boozer that considers a top/bottom flat-plates wall. We show that the solution in which the plasma remains at the initial equilibrium position is improved when bringing the side walls closer to the plasma. When approximating the ITER first wall as a perfect conductor, the plasma remains stable at the initial equilibrium position far beyond the value predicted by the flat-plates wall limit. When considering an opposite limit in which only the inner shell of the ITER vacuum vessel acts as a perfect conductor, the plasma is displaced during the CQ, but the edge safety factor remains above 2 longer in the current decay compared to the flat-plates wall limit. In all the simulated cases, the VDE is found to be strongly dependent on the plasma current, in agreement with a similar finding in the flat-plates wall limit, showing an important difference with hot VDEs in which the CQ is not a necessary condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transient flow modeling in fractured media using graphs

In this work, we describe a method to simulate transient fluid flows in fractured media using an approach based on graph theory. Our approach builds on past work where the graph-based approach was successfully used to simulate steady-state fluid flows in fractured media. We find a mean computational speedup of the order of 1400 from an ensemble of a 100 discrete fracture networks in contrast to the O(10 4 ) speedup that was obtained for steady-state flows earlier. However, the transient flows considered here involve an additional degree of complexity that was not present in the steady-state flows considered previously with a graph-based approach, that of time marching and solution of the flow equations within a time-stepping scheme. We verify our method with an analytical test case and demonstrate its use on a practical problem related to fluid flows in hydraulically fractured reservoirs. By enabling the study of transient flows, we create an opportunity for a wide set of possibilities where a steady-state approximation is not sufficient, such as the example motivated by hydraulic fracturing that we present here. This work validates the concept that graphs are able to reliably capture the topological properties of the fracture network and serve as effective surrogates in an uncertainty-quantification framework.

58 GEOSCIENCES↗
DOE OSTI2023

Results for the September Bimonthly Calendar Year 2022 Tank 50 Salt Solution Sample

In this Technical Report, the chemical and radionuclide contaminant results from the 2022 September bimonthly sample of Tank 50 salt solution are presented in tabulated form. The information from this characterization will be used by Savannah River Mission Completion (SRMC) for the transfer of aqueous waste from Tank 50 to the Saltstone Production Facility (SPF), where the waste will be treated and disposed in the Saltstone Disposal Facility. This Technical Report compares results, where applicable, to SPF Waste Acceptance Criteria (WAC) LIMITS and TARGETS that were established at the time the Tank 50 sample was obtained. The chemical and radionuclide contaminant results from the characterization of the 2022 September bimonthly sampling of Tank 50 were requested by SRMC personnel via a Task Technical Request (TTR). Details of the testing are presented in the Savannah River National Laboratory (SRNL) Task Technical and Quality Assurance Plan (TTQAP). This Technical Report is the TTR deliverable relating Salt Solution Analysis from the SRMC request. The following facts pertaining to the WAC are drawn from the analytical results, including analytical uncertainty, provided in this report. 1) WAC LIMITS and TARGETS were met for all analyzed chemical and radioactive contaminants for which the detection limits are below the WAC LIMITS and TARGETS; 2) Measured average concentrations of nitrate, nitrite and total mercury are approximately 24%, 9% and 12% of the WAC LIMITS, respectively; 3) Measured average concentrations of Tc-99 and I-129 are approximately 19% and 20% of the WAC LIMITS, respectively; 4) All other radionuclide average concentrations are at 4% or less of the WAC LIMITS and TARGETS.

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

A new method for solving the linearized 1D Vlasov–Poisson system yielding a new class of solutions

We describe a new method for solving the linearized 1D Vlasov–Poisson system by using properties of Cauchy-type integrals. Our method remedies critical flaws of the two standard methods, reveals a previously unrecognized Gaussian-in-time-like decay, and can also account for an externally applied electric field. The Landau approximation involves deforming the Bromwich contour around the poles closest to the real axis due to the analytically continued dielectric function, finding the long-time behavior for a stable system: Landau damping. Jackson's generalization encircles all poles while sending the contour to infinity, assuming its contribution vanishes, which is not true in general. This gives incorrect solutions for physically reasonable configurations and can exhibit pathological behavior, of which we show examples. The van Kampen method expresses the solution for a stable equilibrium as a continuous superposition of waves, resulting in an opaque integral. Case's generalization includes unstable systems and predicts a decaying discrete mode for each growing discrete mode, an apparent contradiction to both the Jackson solution and ours. We show, without imposing additional constraints, that the decaying modes are never present in the time evolution due to an exact cancellation with part of the continuum. Our solution is free of integral expressions, is obtained using algebra and Laurent series expansions, does not rely on analytic continuations, and results in a correct asymptotically convergent form in the case of infinite sums. The analysis used can be readily applied in higher-dimensional, electromagnetic systems and also provides a new technique for evaluating certain inverse Laplace transforms.

Physics↗
DOE OSTI2025

ThO 2 and Th 1– x U x O 2 Nanoscale Materials and Thin Films for Nuclear Science Applications

This study investigates the dynamics and mechanisms of solution combustion synthesis (SCS) for the preparation of nanoscale ThO 2 and Th 1–x U x O 2 materials, utilizing metal nitrates (Th(NO 3 ) 4 and UO 2 (NO 3 ) 2 ) and acetylacetone (C 5 H 8 O 2 ) as reactants dissolved in a 2-methoxyethanol (C 3 H 8 O 2 ) solvent. By combining thermodynamic calculations, dynamic time–temperature profile measurements with differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA), this research reveals how variations in acetylacetone concentration and uranium content influence the structural parameters of the synthesized oxides. The time–temperature measurements show that the heating rate and maximum combustion temperatures are sensitive to acetylacetone concentration. DSC-TGA results indicate shifts in exothermic peak temperatures as the uranium content changes. The complexation between thorium and acetylacetone emerges as a critical factor, impacting combustion parameters and the structural characteristics of the final products. The uniform distribution of Th and U in the Th 1–x U x O 2 solid solution and the formation of nanoscale particles with strained crystallites are considered essential for the low-temperature densification of these materials for nuclear fuel pellet applications. Additionally, high-quality ThO 2 and Th 1–x U x O 2 thin (100–150 nm) films are successfully synthesized via electrospray deposition of combustible solutions followed by a brief period of heat treatment. Furthermore, these films exhibit excellent structural and morphological uniformity, making them ideal candidates for nuclear measurements, irradiation damage studies, and investigations into the physical properties of both pure and mixed oxides.

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

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

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

Development of an Attenuated Total Reflectance–Ultraviolet–Visible Probe for the Online Monitoring of Dark Solutions

Optical spectroscopy is a valuable tool for on-line monitoring of a variety of processes. Ultraviolet-visible (UV-vis) spectroscopy in particular, can monitor the concentration of analytes as well as identify speciation and oxidation state. However, it can be difficult to impossible to employ UV-vis based sensors on chemical systems that are very dark (i.e., high optical density) as exceedingly short pathlengths are required (for transmission approaches) or effective means of backscattering are needed (for reflectance approaches). Examples of processes that would benefit significantly from the use of optical sensors and encounter these challenges include used nuclear fuel recycling and molten salts with high concentrations of dissolved uranium. Utilizing an attenuated total reflectance (ATR) UV-vis approach can overcome these challenges and allow for the measurement of solutions orders of magnitude more concentrated than transmission UV-vis. However, determining ideal sensor specifications for varied processes can be time consuming and expensive. Here, in this study, we evaluate the ability for a novel ATR-UV-vis probe to measure very concentrated solutions of Co(II) and Ni(II) nitrate as well as organic dyes (methylene blue, acid red 1, and crystal violet). This sensor design provides a modular method for exploring possible “pathlengths” by altering the exposed ATR fiber length. Also studied were approaches to loading and measuring the sensor cell. These results are compared to a traditional 1 cm cuvette measured by transmission UV-vis. It was found that the ATR-UV-vis probe was capable of measuring solutions 600 times more concentrated than the 1 cm cuvette. Advanced data analysis in the form of multivariate curve resolution (MCR) was used to analyze the speciation of methylene blue over a large concentration range. The application of this novel ATR-UV-vis probe to the interrogation of dark solutions is a promising avenue for use in on-line monitoring of nuclear processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2021

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗