SEARCH · Engineering Papers

Results for “Liquid solid interfaces”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Silica and Pyroxene in IVA Irons; Possible Formation of the IVA Magma by Impact Melting and Reduction of L-LL-Chondrite Materials Followed by Crystallization and Cooling

Group IVA is a large magmatic group of iron meteorites. The mean DELTA O-17 (= delta O-17 - 0.52(raised dot) delta O-18) of the silicates is approx. plus or minus 1.2%o, similar to the highest values in L chondrites and the lowest values in LL chondrites; delta O-18 values are also in the L/LL range. This strongly suggests that IVA irons formed by melting L-LL parental material, but the mean Ni content of IVA irons (83 mg/g) is much lower than that of a presumed L-LL parent (approx. 170 mg/g) and the low-Ca pyroxene present in two IVA meteorites is Fs13, much lower than the Fs20-29 values in L and LL chondrites. Thus, formation from L-LL precursors requires extensive addition of metallic Fe, probably produced by reduction of FeS and FeO. Group IVA also has S/Ni, Ga/Ni, and Ge/Ni ratios that are much lower than those in L-LL chondrites or any chondrite group that preserves nebular compositions, implying loss of these volatile elements during asteroidal processing. We suggest that these reduction and loss processes occurred near the surface of the asteroid during impact heating, and resulted partly from reduction by C, and partly from the thermal dissociation of FeS and FeO with loss of O and S. The hot (approx. 1770 K) low-viscosity melt quickly moved through channels in the porous asteroid to form a core. Two members of the IVA group, Sao Joao Nepomuceno (hereafter, SJN) and Steinbach, contain moderate amounts of orthopyroxene and silica, and minor amounts of low-Ca clinopyroxene. Even though SJN formed after approx. 26% crystallization and Steinbach formed after approx. 77% Crystallization of the IVA core, both could have originated within several tens of meters of the core-mantle interface if 99% of the crystallization occurred from the center outwards. Two other members of the group (Gibeon and Bishop Canyon) contain tabular tridymite, which we infer to have initially formed as veins deposited from a cooling SiO-rich vapor. The silicates were clearly introduced into IVA irons after the initial magma crystallized. Because the y-iron crystals in SJN are typically about 5 cm across, an order of magnitude smaller than in IVA irons that do not contain massive silicates, we infer that the metal was in the gamma-iron field when the silicates were injected. The SJN and Steinbach silicate compositions are near the low-Ca-pyroxene/silica eutectic compositions. We suggest that a tectonic event produced a eutectic-like liquid and injected it together with unmelted pyroxene grains into fissures in the solid metal core. Published estimates of IVA metallographic cooling rates range from 20 to 3000 K/Ma, leading to a hypothesized breakup of the core during a major impact followed by scrambling of the core and mantle debris [Haack, H., Scott, E.R.D., Love, S.G., Brearley, A. 1996. Thermal histories of IVA stony-iron and iron meteorites: evidence for asteroid fragmentation and reaccretion. Geochim. Cosmochim. Acta 60, 3103-3113]. This scrambling model is physically implausible and cannot explain the strong correlation of estimated cooling rates with metal composition. Previous workers concluded that the low-Ca clinopyroxene in SJN and Steinbach formed from protopyroxene by quenching at a cooling rate of 10(sup 12) K/Ma, and suggested that this also supported an impact-scrambling model. This implausible spike in cooling rate by a factor of 10(sup 10) can be avoided if the low-Ca clinopyroxene were formed by a late shock event that converted orthopyroxene to clinopyroxene followed by minimal growth in the clinopyroxene field, probably because melt was also produced. We suggest that metallographic cooling-rate estimates (e.g., based on island taenite) giving similar values throughout the metal compositional range are more plausible, and that the IVA parent asteroid can be modeled by monotonic cooling followed by a high-temperature impact event that introduced silicates into the metal and a low-temperature impact event that partially converted thopyroxene into low-Ca clinopyroxene.

Wasson, John T.↗

DOE OSTI2020

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗

DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Probing the thermodynamics and kinetics of ethylene carbonate reduction at the electrode–electrolyte interface with molecular simulations

Understanding the formation of the solid-electrolyte interphase (SEI) in lithium-ion batteries (LIBs) is an ongoing area of research due to its high degree of complexity and the difficulties encountered by experimental studies. Herein, we investigate the initial stage of SEI growth, the reduction reaction of ethylene carbonate (EC), from both a thermodynamic and kinetic approach with theory and molecular simulation. We employed both the potential distribution theorem (PDT) and the SMD implicit solvent model to EC solvation for the estimation of reduction potentials of Li$^+$, EC, and Li$^+$-solvating EC (s-EC), as well as reduction rate constants of EC and s-EC. We find that solvation effects greatly influence these quantities of interest, particularly the Li$^+$/Li reference electrode potential in EC solvent. Further, we also compute the inner- and outer-sphere reorganization energies for both EC and s-EC at the interface of liquid EC and a hydroxyl-terminated graphite surface, where total reorganization energies are predicted to be 76.6 and 88.9 kcal/mol, respectively. With the computed reorganization energies, we estimate reduction rate constants across a range of overpotentials and show that EC has a larger electron transfer rate constant than s-EC at equilibrium, despite s-EC being more thermodynamically favorable. Furthermore, this manuscript demonstrates how ion solvation effects largely govern the prediction of reduction potentials and electron transfer rate constants at the electrode-electrolyte interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Utilizing Commercial Submersible Mixer Pumps for Sludge Removal in Savannah River Site's Tank 26 - 20289

The Savannah River Site (SRS) Liquid Waste System (LWS) safely manages, stores, treats, and dispositions liquid radioactive waste. The LWS consists of 51 underground waste storage tanks (eight of which are operationally closed and filled with grout), waste evaporators, treatment facilities, and solidification facilities, known as the Defense Waste Processing Facility (DWPF) and Saltstone Production Facility (SPF). One of the waste storage tanks, Tank 26, was placed into service in 1980 as an F Tank Farm (FTF) Evaporator Feed Tank. From 1980 to 2013, Tank 26 received F Canyon receipts and dilute supernate that was then transferred to the 242-16F FTF Evaporator for volume reduction. In 2013, the steam tube bundle failed in the FTF Evaporator, and the decision was made not to repair/resume evaporator operation. During operation of Tank 26, solids built up to a level of 2.03 m in the tank with a volume of 1062.9 kL. Due to the evaporator failure and subsequent shutdown, the solids in Tank 26 were selected to feed Sludge Batch 10, which is collected and prepped prior to being sent to DWPF for final disposition. To accomplish solids removal, Commercial Submersible Mixer Pumps (CSMPs) were selected to slurry the solids in preparation to be sent to Tank 51. Four CSMPs were installed in Tank 26 with each one installed in a separate quadrant of Tank 26. This work would mark the first deployment of CSMPs in an SRS waste tank. The CSMPs were developed in response to operational issues from previous mixing pumps and budgetary constraints. The CSMPs use the concept of modifying commercially available equipment for nuclear waste applications. The CSMPs consist of a 230-horsepower submersible mixing pump (manufactured by GPM, Inc.) joined to a mast fabricated at SRS. For this application, Savannah River Remediation (SRR) design services was challenged to provide a mast design that required no steel superstructure support system on top of Tank 26. Also, SRR design services provided a simplified Tank 26 riser interface that did not require extensive machine parts to adjust height of the suction screen of the CSMP within the tank. Other design improvements include no requirements for tank top radiation monitors and no requirements for flushing of the CSMPs during startup/shutdown operations. The CSMPs were operated at an initial height of 2.29 m above the Tank 26 bottom for 10 days. After sludge sounding, the CSMPs were lowered to a height of 1.52 m with no issues due to proper work planning and the simplified riser interface. With the CSMPs lowered to a level of 1.52 meters above the tank bottom, the CSMPs were operated for 12 days and another sounding was performed. The sounding level came back matching the disturbance depth results of the first pump run. Chemistry samples were analyzed while the tank was left to settle. The analysis revealed enough weight percent solids to provide good feed to Sludge Batch 10. Also, the results revealed the supernate within the mixture would provide good salt solution feed to Salt Waste Processing Facility (SWPF) Salt Batch 3. So, the plan was modified to let the disturbed solids settle and decant the salt solution for addition to SWPF Salt Batch 3. Then water will be added back to Tank 26, and the CSMPs run in order to wash the solids to decrease settling time prior to being sent to Tank 51 for Sludge Batch 10. Utilizing CSMPs for waste removal in an SRS Tank has provided a cost-effective means for further waste removal efforts. In addition, the CSMPs are easier to operate by utilizing a robust and simplified design. The CSMPs performed quite well with no process shutdown or delays during operation. As a result, CSMPs are integral to the future of removing radioactive waste from storage tanks at SRS. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Analysis of an all-solid state nanobattery using molecular dynamics simulations under an external electric field

Present Li-ion battery (LIB) technology requires strong improvements in performance, energy capacity, charging-time, and cost to expand their application to e-mobility and grid storage. Li-metal is one of the most promising materials to replace commercial anodes such as graphite because of its 10 times higher specific capacity. However, Li-metal has high reactivity with commercial liquid electrolytes; thus, new solid materials are proposed to replace liquid electrolytes when Li-metal anodes are used. We present a theoretical analysis of the charging process in a full nanobattery, containing a LiCoO 2 cathode, a Li 7 P2S 8 I solid-state electrolyte (SSE), a Li-metal anode as well as Al and Cu collectors for the cathode and anode, respectively. In addition, we added a Li 3 P/Li 2 S film as a solid electrolyte interphase (SEI) layer between the Li-anode and SSE. Thus, we focus this study on the SEI and SSE. We simulated the charging of the nanobattery with an external voltage by applying an electric field. We estimated temperature profiles within the nanobattery and analyzed Li-ion transport through the SSE and SEI. Here, we observed a slight temperature rise at the SEI due to reactions forming $PS_{3}^{–}$ and $P_{2}S_{7}^{4}$$^{–}$ fragments at the interfaces; however, this temperature profile changes due to the charging current under the presence of the external electric field ε = 0.75 V Å –1 . Without the external field, the calculated open-circuit voltage (OCV) was 3.86 V for the battery, which is within the range of values of commercial cobalt-based LIBs. This voltage implies a spontaneous fall of available Li-ions from the anode to the cathode (during discharge). The charge of this nanobattery requires overcoming the OCV plus an additional voltage that determines the charging current. Thus, we applied an external potential able to neutralize the OCV, plus an additional 1.6 V to induce the transport of Li + from the cathode up to the anode. Several interesting details about Li + transport paths through the SSE and SEI are discussed.

25 ENERGY STORAGE↗

NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗

NASA NTRSDate not supplied

Validation of Two-Phase Flow Modeling Techniques in Loci-Stream Using Axial Jet Pressure Control and Autogenous Pressurization Experiments

Improvements to computational fluid dynamics modeling of cryogenic fluid management operations are presented with experimental validations. Modeling improvements includes utilizing improved delayed detaches eddy simulation to more effectively model turbulence compared to the Reynolds-averaged Navier Stokes modeling, modeling temperature dependent tank solid properties to capture large variation in thermodynamic properties at cryogenic conditions, modeling real gas equation of state for hydrogen instead of ideal gas behavior, and an on-going effort to improve interface turbulence transport for the volume of fluid method. Two flight-scale ground tests, five K-site tank axial jet pressure control experiments and an Engineering Design Unit tank autogenous pressurization experiment are validated and show the effectiveness of the improved modeling techniques. The successful validation establishes further trust in these models for flight project support and provides a pathway for maturation of in-development interface turbulence damping model for the volume of fluid method for more generalized applications.

CFD↗

NASA NTRSDate not supplied

Validation of Two-phase Flow Modeling Techniques in Loci-Stream using Axial Jet Pressure Control and Autogenous Pressurization Experiment

CFD↗

DOE OSTI2024

Electrochemical formation of bis(fluorosulfonyl)imide-derived solid-electrolyte interphase at Li-metal potential

Lithium bis(fluorosulfonyl)imide-based liquid electrolytes are promising for realizing high Coulombic efficiency and long cycle life for next-generation Li-metal batteries. However, the role of anions in the formation of solid-electrolyte interphase remains unclear. Herein, we combine electrochemical analyses and X-ray photoelectron spectroscopy measurements both with and without sample washing, together with computational simulations, to propose the reaction pathways of electrolyte decomposition and correlate the interphase component solubility with the efficacy of passivation. Additionally, we discovered that not all the products derived from interphase-forming reactions were incorporated into the resulting passivation layer, with a notable portion present in the liquid electrolyte. Further, we found that the high-performance electrolytes can afford a sufficiently passivating interphase with minimized electrolyte decomposition, by incorporating more anion-decomposition products. Overall, this work presents a systematic approach of coupling electrochemical and surface analyses to paint a comprehensive picture of solid-electrolyte interphase formation, while identifying key attributes of high-performance electrolytes for guiding future designs.

Yu, Weilai [Stanford University, CA (United States↗

NASA NTRSDate not supplied

Overview of the SLS Core Stage Thrust Vector Control System Design

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) consists of four independent hydraulic systems. The SLS CS TVC system is comprised of 8 mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. Each hydraulic system nominally provides hydraulic power to one RS-25 engine and two actuators. Additionally, each system provides redundant control capability to one actuator on each of its neighboring systems. The RS-25 uses hydraulic power to control propellant valves, and the TVC actuators are used to move the engine in the pitch and yaw gimbal planes. The TVC system design leverages hardware from the Space Shuttle program as well as new hardware designed specifically for the Core Stage. The Space Shuttle heritage hardware directly reused on SLS includes the Orbiter TVC hydraulic servo-actuators (with two slight design modifications), the Orbiter hydraulic circulation pumps, the Orbiter gimbal block/bearing, and the Solid Rocket Booster hydraulic pumps. The Solid Rocket Booster APU turbines are powered by hot gas produced by a catalyzed hydrazine decomposition. The SLS Core Auxiliary Power Unit (CAPU) is derived from the Space Shuttle Orbiter Auxiliary Power Unit (APU); on the SLS Core Stage, the CAPU turbine is spun using cold gas tapped-off from the RS-25 to CS liquid hydrogen autogenous pressurization line. The remaining hardware in the TVC system (hydraulic Filter Manifold (FM), hydraulic Supply Accumulator (SA), hydraulic Return Accumulator (RA), Hydraulic Reservoir, Exhaust Gas Heat Exchanger (EGHE)) as well as the avionics providing control and telemetry (TVC Actuator Controller (TAC) and CAPU Controller (CAPUC) are new components developed for SLS. This paper is the first installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, the overall design architecture of the CS TVC is presented, with a focus on the interfaces between the TVC actuators, the engines, their hydraulic power systems, and the avionics that provide commands from the SLS Vehicle Management (VM) software to effect stable and robust flight control for the integrated SLS launch vehicle.

Thrust Vector Control↗

NASA NTRSDate not supplied

Overview of the SLS Core Stage Thrust Vector Control System Design

Thrust Vector Control↗

DOE OSTI2020

The organic redox transistor for neuromorphic computing

Inspired by the in-memory computing architectures of biological systems, neuromorphic computing using crossbar arrays of artificial synapses based on non-volatile memory (NVM) devices with variable conductances has emerged as a new paradigm to enable massively parallel and ultra-low power computing hardware for data centric applications. Although inference has been demonstrated successfully using crossbars based on a variety of NMV technologies, efficient learning and scaling to large arrays (>10 6 elements) remains a challenge due to the synaptic elements' non-ideal electrical characteristics which degrades ANN accuracy. A further challenge is that in the conductive state memristors draw large currents >μA resulting in significant voltage drops in the interconnect wires and increased probability of failure in scaled arrays. We suggest the organic polymer redox transistor (RT) is an alternate approach that could solve many of these challenges, enabling both inference and parallel outer product updates, as recently demonstrated by Fuller et al. An RT consists of redox-active channel and gate electrodes in contact with a liquid or solid electrolyte. lon insertion through the electrolyte controls the channel electronic conductivity, while electron transfer through an external circuit maintains overall charge neutrality. Unlike a rechargeable battery, in the RT the voltage built-up across the electrolyte is kept to a minimum (typically <100 mV) by using the same material for the gate and channel. Elimination of the voltage offset simplifies integration of the RT into programmable arrays by enabling the use of various selectors. RTs based on inorganic and organic materials have been recently demonstrated with conductance tuning occurring at potentials of just a few mV and hundreds to thousands of linearly and symmetrically programmable conductance states, enabling near ideal accuracy in neural network simulations. Introduced in the 1980's, redox transistors with metallic gate electrodes and organic channel materials, also known as organic electrochemical transistors (OECTs), have been explored for a variety of applications such as chem- and bio-sensing, neural interfaces, and low cost printed circuits. A typical channel material for OECTs is the conducting polymer poly(3,4-ethylenedioxythiophene) doped with poly(styrene sulfonate) (PEDOT:PSS). PEDOT is a p-type semiconducting polymer with mobile positively charged polarons that hop chain-to-chain.

Talin, Albert Alec↗

DOE OSTI2020

Development of Operando Nuclear Magnetic Resonance Capabilities for Solid-State Battery Research

Mechanistic information regarding solid-solid interfaces in solid-state batteries remains elusive as its observation requires a non-destructive operando technique. Such a technique does not currently exist. Solid-state nuclear magnetic resonance (NMR) is uniquely suited to address this challenge as it is sensitive to chemistry, structure, and dynamics and can measure buried interfaces. However solid-state NMR requires high speed spinning of the sample at the ‘magic angle’ (MAS) to obtain adequate resolution. The metallic current collectors in a battery are detrimental to spinning in magnetic fields as they generate eddy currents which prevent spinning and dampen the NMR signal. Furthermore, the spinning nature of the MAS NMR experiment requires a wireless means of charging for operando capabilities. These challenges were addressed in this research and operando MAS NMR enabled by optical wireless charging was demonstrated. By limiting the size of the metallic components, we were able to spin a sample holder with electrical leads and a solid-state battery up to a few kilohertz, fast enough for high-resolution NMR. We achieved wireless charging by attaching a photovoltaic cell to the head of the NMR rotor and used a laser to transmit power to the photovoltaic cell. Demonstration of the optical wireless charging was performed by polarizing a liquid electrolyte during an operando MAS NMR measurement. This novel capability is a transformative technique for investigating the mechanistic electrochemical processes in solid-state samples.

25 ENERGY STORAGE↗

DOE OSTI2022

Wet chemical synthesis and properties of argyrodite sulfide solid electrolytes for solid state lithium batteries

The commercialization of the lithium-ion battery (LIB) in 1991 was responsible for the explosion in portable electronic technologies that has been seen over the past 30 years. With the advent of electric vehicles and other high-powered technologies, there is tremendous demand for LIBs with higher energy density and high safety. To achieve this, new electrode materials must be explored. The obvious choice of anode material would be pure metal lithium, which has a theoretical specific capacity of 3860 mAh g-1 . Unfortunately, metal lithium anodes have not been widely commercialized due to their tendency to react violently with the flammable liquid electrolytes used in today’s batteries. Battery safety can best be achieved by adopting solid electrolytes in place of liquid electrolytes. Solid electrolytes are nonvolatile and nonflammable, safely allowing for the combination of high-capacity cathode materials with a Li metal anode. Argyrodite sulfide solid electrolytes such as halogen-doped Li6PS5X (X = Cl, Br, I) are noted for their high ionic conductivity. But before sulfides can be commercially adopted, they possess several disadvantages which must be addressed, including time- and energy-consuming synthesis processes, poor electrochemical stability, and intrinsically poor air stability. This dissertation seeks to address each of these challenges through materials design an synthesis strategies. In this work, we pioneer a solvent-based approach for the synthesis of argyrodite solid electrolytes Li7PS6 and Li6PS5Xinstead of a stringent solid-state synthesis. Nontoxic ethanol is employed as the solvent, enabling a rapid synthetic approach to produce argyrodite solid electrolytes with high phase purity and compositional flexibility. Compared with Li7PS6, halogen doping (i.e. X = F, Cl, Br, I) not only increases the ionic conductivity, but also enhances the electrochemical stability at the interface towards Li metal. Specifically, F-doped argyrodites produce a robust SEI layer containing LiF, contributing to enhanced interfacial stability. Finally, to address the air instability challenge, argyrodite-incorporated composite solid electrolytes (CSEs) are designed and prepared to produce stable and flexible membranes that are demonstrated in solid-state Li metal batteries. These advances push argyrodite sulfide solid electrolyte research further and pave the way for the proliferation of next generation lithium metal batteries.

25 ENERGY STORAGE↗

DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗

DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗