SEARCH · Engineering Papers
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.
This article presents a perspective on what we think are key topics related to the structure and structural dynamics of ILs and to some extent high-temperature molten salts.
This research advances the field of geothermal energy by developing and evaluating novel fluid systems designed explicitly for high-temperature subsurface environments to enhance fracture flow distribution. The project integrates chemical synthesis, environmental assessment, rheological testing, microfluidic analysis, and reservoir-scale laboratory testing and analysis to identify optimal ionic and nonionic fluids for Enhanced Geothermal Systems (EGS).
Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.
We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.
The local atomic structure of NaF-ZrF 4 (53–47 mol%) molten system and its evolution with temperature are examined with x-ray scattering measurements which are then used to validate the quality of ab initio and neural network-based molecular dynamics (NNMD) calculations in the temperature range 515–700°C. The machine-learning enhanced NNMD calculations offer improved efficiency while maintaining accuracy at higher distances compared to ab initio calculations. Looking at the evolution of the pair distribution function with increasing temperature, a fundamental change in the liquid structure within the selected temperature range, accompanied by a slight decrease in overall correlation is revealed. NNMD calculations indicate the coexistence of three different fluorozirconate complexes: [ZrF 6 ] 2– , [ZrF 7 ] 3– , and [ZrF 8 ] 4– , with a shift in the dominant coordination state from the 7-coordinated Zr cation toward a 6-coordinated cation with increasing temperature. The study also highlights the metastability of different local coordination structures, with frequent interconversions between the 6- and 7-coordinate states. Analysis of the Zr-F-Zr angular distribution function reveals the presence of both “edge-sharing” and “corner-sharing” fluorozirconate complexes with specific bond angles and distances in accord with previous studies, while the next-nearest-neighbor cation-cation correlations demonstrate a clear preference for unlike cations as nearest-neighbor pairs, emphasizing nonrandom arrangement. Finally, these findings contribute to a comprehensive understanding of the complex local structure of the molten salt, providing insights into temperature-dependent preferences and correlations within the molten system.
Effective monitoring of internal corrosion risk is crucial to ensuring the safety and longevity of natural gas pipeline infrastructure. While electrochemical sensors are commonly used to assess corrosion rates and corrosion indicators in aqueous fluids, they are rarely used in gas pipelines as these fluids lack the ionic conductivity needed for electrochemical measurements. The inclusion of ion-conductive membranes into electrochemical sensors can extend their functionality into humidified gas streams, providing critical information about emerging corrosion events that are common during withdrawal season in pipeline systems downstream from underground storage facilities. In parallel, new protective films, like those obtained through cold spray coating, are being developed to protect oil and gas pipelines and recover losses in structural integrity due to corrosion damage. Herein, we demonstrate how membrane-based electrochemical sensors (MBES) can be used to monitor fluid corrosivity by examining their response to changes in water content for a wide range of fluid compositions. It was found that MBES readings were highly sensitive to water content changes with membrane conductivity measurements varying from 10 –6 to 10 –1 S cm -1 , and corrosion rate measurements which varied from 10 –7 to 1 mm y -1 . Electron microscopy confirmed that the self-healing characteristics of metal coating films were still active despite their inclusion into an MBES probe. In conclusion, these findings indicate that membrane-based corrosion monitoring can be expanded to monitor coated-pipeline materials and provide early detection of emerging corrosion upsets relevant to underground gas storage facilities.
The directed self-assembly of colloidal nanoparticles (NPs) using external fields guides the formation of sophisticated hierarchical materials but becomes less effective with decreasing particle size. As an alternative, electron-beam-driven assembly offers a potential avenue for targeted nanoscale manipulation, yet remains poorly controlled due to the variety and complexity of beam interaction mechanisms. Here, we investigate the beam–particle interaction of silica NPs pinned to the fluid–vacuum interface of ionic liquid droplets. In these experiments, scanning electron microscopy of the droplet surface resolves NP trajectories over space and time while simultaneously driving their reorganization. With this platform, we demonstrate the ability to direct particle transport and create transient, reversible colloidal patterns on the droplet surface. By tuning the beam voltage, we achieve precise control over both the strength and sign of the beam–particle interaction, with low voltages repelling particles and high voltages attracting them. This response stems from the formation of well-defined solvent flow fields generated from trace radiolysis of the ionic liquid, as determined through statistical analysis of single-particle trajectories under varying solvent composition. Altogether, electron-beam-guided assembly introduces a versatile strategy for nanoscale colloidal manipulation, offering new possibilities for the design of dynamic, reconfigurable systems with applications in adaptive photonics and catalysis.
Electrically induced bearing damage (EIBD) is a growing concern in electric drivetrains, where electrical discharge currents generate localized pitting that accelerates wear and premature failure. While ionic liquids (ILs) have shown promise as lubricant additives due to their natural physical adsorption, strong capabilities of tribofilm formation, and high molecular tunability, their role in mitigating EIBD is little known. This gap is especially evident because lubricant viscosity and film thickness affect discharge behavior, and additive performance in these regimes has not been well explored. In this study, tribological tests were performed using a ball-on-disk pure rolling-contact system under an applied voltage on polyalphaolefin oils of three viscosities (4, 10, and 150 cSt measured at 100 °C), with and without the addition of a phosphonium-phosphate IL. Surface damage was analyzed using scanning electron microscopy (SEM) and stylus profilometry. Damage quantification was based on pit morphology analysis of SEM images. The results provide new insights into how oil viscosity governs the discharge behavior and how an IL additive could influence the pit formation. While viscosity seems to play a significant role, an IL additive shows potential to reduce EIBD and deserves further study.
Waste CO 2 can be repurposed as a carbon feedstock for synthesizing valuable chemicals via CO 2 electrolysis. Porous solid electrolyte (PSE) CO 2 electrolysis has been demonstrated as an economically viable method to produce high purity products. This work applies an interfacial engineering approach to determine key factors to improve performance in PSE CO 2 electrolyzers. We standardize the assembly by binding the ionic resin into an ionomer wafer and utilize Computational Fluid Dynamics (CFD) to design gaskets for uniform fluid flow. Here, we employ the distribution of relaxation times (DRT) method to determine that anionic-conducting interfaces are the primary contributor to energy losses. To address this, we demonstrate that enhancing the contact between the cathode and the anion exchange membrane (AEM) and the AEM-ionic resin interface allows for low overpotential in deionized water operation.
We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature nonuniform distributions of charge nonuniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator that is powered by enzyme-coated patches that catalyze reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion, which disrupts a nonuniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.
Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.
A tube-in-tube heat exchanger utilizes a selectively permeable tube having a selective permeable layer to allow the refrigerant to transfer into an ionic liquid to generate heating or cooling. The ionic liquid then provides heating or cooling to a heat transfer fluid through a non-permeable layer or tube. The system may be configured as a shell and tube design, with the third fluid free to flow on the outside of the shell, or as a shell and tube-in-tube, with a central tube containing a first liquid, a second tube containing a second liquid, and an outer shell containing the third liquid. The selectively permeable tube may include an anion or cation selectively permeable layer and this layer may be supported by a support layer or tube.
Here, we evaluate the exponentially rare fluctuations of the ionic current for a dilute electrolyte by means of macroscopic fluctuation theory. We consider the fluctuating hydrodynamics of a fluid electrolyte described by a stochastic Poisson–Nernst–Planck equation. We derive the Euler–Lagrange equations that dictate the optimal concentration profiles of ions conditioned on exhibiting a given current, whose form determines the likelihood of that current in the long-time limit. For a symmetric electrolyte under small applied voltages, number density fluctuations are small, and ionic current fluctuations are Gaussian with a variance determined by the Nernst–Einstein conductivity. Under large applied potentials, the ionic current distribution is generically non-Gaussian. Its structure is constrained thermodynamically by Gallavotti–Cohen symmetry and the thermodynamic uncertainty principle.
Molten salts are ionic liquids that are used for the electrolytic pyroprocessing of metals and as heat transfer fluids in very high temperature processes. Recently, halide-based molten salt reactors (MSR) have gained momentum for high density and environmentally responsible electricity generation. The function of these reactors and their fuel cycle depend on a knowledge of the halide salt’s thermodynamic properties [1]. Therefore, high accuracy phase equilibria of MSR relevant base halide and actinide containing salts are needed. Modern thermal analysis of the phase transitions of halide salts is usually done during heating at relatively high scan rates with commercial devices. Normally such measurements are adequate for pure or pseudo-binary salts. However, as the number of components in the mixture increases, accurately resolving the liquidus becomes increasingly difficult. Reversing the scanning mode greatly increases the sensitivity of phase transition measurements but can decrease accuracy due to undercooling [2] — a common occurrence in molten halide systems [3]. This work presents the development of a differential thermal analysis (DTA) cell intended for use in radiological gloveboxes. Non-contact ultrasonic agitation of the halide salt is implemented to limit kinetic limitations on crystallization during cooling to minimize undercooling [4]. Measurements on halide salts are also presented to elucidate the effect of non-contact mixing on undercooling. Reference [1] S. Boyd and C. Taylor, “3 - Chemical fundamentals and applications of molten salts,” in Molten Salt Reactors and Thorium Energy, T. J. Dolan, Ed., Woodhead Publishing, 2017, pp. 29–91. doi: 10.1016/B978-0-08-101126-3.00003-8. [2] K. Nitsch, A. Cihlár, and M. Rodová, “Molten state and supercooling of lead halides,” J. Cryst. Growth, vol. 264, no. 1, pp. 492–498, Mar. 2004, doi: 10.1016/j.jcrysgro.2004.01.011. [3] L. Rycerz, “Practical remarks concerning phase diagrams determination on the basis of differential scanning calorimetry measurements,” J. Therm. Anal. Calorim., vol. 113, no. 1, pp. 231–238, Jul. 2013, doi: 10.1007/s10973-013-3097-0. [4] Md. H. Zahir, S. A. Mohamed, R. Saidur, and F. A. Al-Sulaiman, “Supercooling of phase-change materials and the techniques used to mitigate the phenomenon,” Appl. Energy, vol. 240, pp. 793–817, Apr. 2019, doi: 10.1016/j.apenergy.2019.02.045.
Not Available
Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.
Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.
The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.