SEARCH · Search NASA

Results for “Liquid structure simulation”

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 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2025

New class of tritium breeders for fusion applications: Metal-reinforced composite breeders

Commercial fusion reactors operating on a D-T fuel cycle will require a steady supply of tritium to maintain the burning plasma required for continuous power generation. Since tritium has a short half-life, there is negligible natural abundance which necessitates fusion reactors to produce their own source of tritium. Tritium is most easily produced by surrounding a fusion reactor core with lithium (Li), which reacts under the intense neutron flux leaving the reactor core to form tritium and helium. Due to the hazards and technical challenges associated with surrounding a fusion reactor core with many tons of molten Li, other Li-bearing tritium breeder materials have been pursued. Unfortunately, most of the liquid breeders historically examined are exceedingly corrosive to reactor structural materials while many solid breeders in the form of ceramics are forced to make tradeoffs between Li content and mechanical integrity. In this work, to break the historic limit between Li-density and mechanical integrity of traditional solid breeders, a new class of solid tritium breeders is developed: metal-reinforced composite (MERC) breeders. Specifically, the high Li-density of lithium oxide (Li 2 O) is exploited through the addition of a metal reinforcing phase, producing a composite breeder material exhibiting high splitting tensile strength and quasi-ductility with a Li-density greater than other leading solid breeder candidates, including lithium orthosilicate (Li 4 SiO 4 ) and lithium metatitanate (Li 2 TiO 3 ). Mechanical testing, microstructural characterization, and neutronic simulation results are presented and discussed in light of fusion reactor design considerations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

Effect of Hydrofluoroether Cosolvent Addition on Li Solvation in Acetonitrile-Based Solvate Electrolytes and Its Influence on S Reduction in a Li–S Battery

Li-S batteries are a promising next-generation battery technology. Due to the formation of soluble polysulfides during cell operation, the electrolyte composition of the cell plays an active role in directing the formation and speciation of the soluble lithium polysulfides. Recently, new classes of electrolytes termed "solvates" that contain stoichiometric quantities of salt and solvent and form a liquid at room temperature have been explored due to their sparingly solvating properties with respect to polysulfides. The viscosity of the solvate electrolytes is understandably high limiting their viability; however, hydrofluoroether cosolvents, thought to be inert to the solvate structure itself, can be introduced to reduce viscosity and enhance diffusion. Nazar and co-workers previously reported that addition of 1,1,2,2-tetrafluoroethyl 2,2,3,3-tetrafluoropropyl ether (TTE) to the LiTFSI in acetonitrile solvate, (MeCN) 2 -LiTFSI, results in enhanced capacity retention compared to the neat solvate. Here, we evaluate the effect of TTE addition on both the electrochemical behavior of the Li-S cell and the solvation structure of the (MeCN) 2 -LiTFSI electrolyte. Contrary to previous suggestions, Raman and NMR spectroscopy coupled with ab initio molecular dynamics simulations show that TTE coordinates to Li + at the expense of MeCN coordination, thereby producing a higher content of free MeCN, a good polysulfide solvent, in the electrolyte. Furthermore, the electrolytes containing a higher free MeCN content facilitate faster polysulfide formation kinetics during the electrochemical reduction of S in a Li-S cell likely as a result of the solvation power of the free MeCN.

25 ENERGY STORAGE↗

DOE OSTI2019

Microbial Range Expansions on Liquid Substrates

Despite the importance of fluid flow for transporting and organizing populations, few laboratory systems exist to systematically investigate the impact of advection on their spatial evolutionary dynamics. To address this problem, we study the morphology and genetic spatial structure of microbial colonies growing on the surface of a nutrient-laden fluid 10 4 to 10 5 times more viscous than water in Petri dishes; the extreme but finite viscosity inhibits undesired thermal convection and allows populations to effectively live at the air-liquid interface due to capillary forces. We discover that S. cerevisiae (baker’s yeast) growing on a viscous liquid behave like “active matter”: They metabolically generate fluid flows many times larger than their unperturbed colony expansion speed, and that flow, in turn, can dramatically impact their colony morphology and spatial population genetics. We show that yeast cells generate fluid flows by consuming surrounding nutrients and decreasing the local substrate density, leading to misaligned fluid pressure and density contours, which ultimately generates vorticity via a thresholdless baroclinic instability. Numerical simulations with experimentally measured parameters demonstrate that an intense vortex ring is produced below the colony’s edge. As the viscosity of the substrate is lowered and the self-induced flow intensifies, we observe three distinct morphologies: At the highest viscosity, cell proliferation and movement produces compact circular colonies with, however, a stretched regime of exponential expansion; intermediate viscosities give rise to compact colonies with “fingers” that are usually monoclonal and then break into smaller cell clusters; at the lowest viscosity, the expanding colony fractures into many genetically diverse, mutually repelling, islandlike fragments that can colonize an entire 94-mm-diameter Petri dish within 36 hours. We propose a simple phenomenological model that predicts the early colony dynamics. Our results provide rich opportunities to study the interplay between fluid flow and spatial population genetics for future investigations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

An effective air–liquid contactor for CO 2 direct air capture using aqueous solvents

Here, the development of a cost-effective, corrosion resistant, high-flux direct air capture (HiDAC) contactor for a solvent-based direct air capture (DAC) process is reported. Literature technoeconomic analyses suggest that the air–liquid contactor can cost over 20 % of the overall DAC plant’s annualized capital costs. To bring down the overall cost of DAC, it is imperative that an effective contactor is developed. A hybrid contactor consisting of a commercial polyvinylchloride structured packing enhanced with stainless-steel 410 random packing has been developed to provide a high surface area for air–liquid contact. The contactor geometry, wettability, corrosion resistance, pressure drop, along with its CO 2 uptake efficiency, CO 2 uptake rates, and extended loading potential using potassium sarcosinate solutions are investigated. Results show that the HiDAC contactor has a relatively high specific surface area (885 m 2 /m 3 ), which allows for CO 2 uptake efficiencies of up to 75 % and capture rates of up to 550 g of CO 2 per day for a 0.3 × 0.25 × 0.3 m contactor. The contactor also exhibits high levels of wettability and corrosion resistance with amino acid-based DAC solvents. A CO 2 uptake model was developed, and modeling results are compared to experimental data to simulate and predict the performance of the contactor in a DAC process using potassium sarcosinate solvent. The system’s overall mass transfer coefficient and theoretical pressure drop were also calculated, and these results were compared to DAC data found in literature. The results presented indicate that the HiDAC contactor is well suited for DAC due to its high specific surface area, resistance to corrosion, and high degree of wettability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗

DOE OSTI2016

Obtaining the solid-liquid interfacial free energy via multi-scheme thermodynamic integration: Ag-ethylene glycol interfaces

The solid-liquid interfacial free energy γ sl is an important quantity in wetting, nucleation, and crystal growth. Although various methods have been developed to calculate γ sl with atomic-scale simulations, such calculations still remain challenging for multi-component interfaces between molecular fluids and solids. In this paper, we present a multi-scheme thermodynamic integration method that is inspired by the “cleaving-wall” method and aimed at obtaining γ sl for such systems using open-source simulation packages. This method advances two aspects of its predecessor methods. First, we incorporate separate schemes to resolve difficulties when manipulating periodic boundary conditions of the supercell using open-source simulation packages. Second, we introduce a numerical approximation to obtain thermodynamic integrands for complex force fields when an analytical differentiation is not readily available. To demonstrate this method, we obtain γ sl for interfaces between Ag(100) and Ag(111) and ethylene glycol (EG). These interfacial free energies mirror interfacial potential energies for each facet. We also estimate entropies of interface formation and these are consistent with theoretical predictions in signs and trends. For the Ag-EG systems, we find that the largest contribution to γ sl is the free energy to create the bare metal surfaces. The second-largest contribution to γ sl is from the liquid-solid interaction. Finally, this user-friendly method will accelerate investigation in a broad range of research topics, such as the thermodynamic effect of structure-directing agents in solution-phase shape-controlled nanocrystal syntheses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗

NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank

Alok Majumdar↗

DOE OSTI2023

Phase coexistence implications of violating Newton’s third law

Newton’s third law, action = reaction, is a foundational statement of classical mechanics. However, in natural and living systems, this law appears to be routinely violated for constituents interacting in a nonequilibrium environment. Here, we use computer simulations to explore the macroscopic phase behavior implications of breaking microscopic interaction reciprocity for a simple model system. We consider a binary mixture of attractive particles and introduce a parameter that is a continuous measure of the degree to which interaction reciprocity is broken. In the reciprocal limit, the species are indistinguishable, and the system phase separates into domains with distinct densities and identical compositions. Furthermore, increasing nonreciprocity is found to drive the system to explore a rich assortment of phases, including phases with strong composition asymmetries and three-phase coexistence. Many of the states induced by these forces, including traveling crystals and liquids, have no equilibrium analogs. By mapping the complete phase diagram for this model system and characterizing these unique phases, our findings offer a concrete path forward toward understanding how nonreciprocity shapes the structures found in living systems and how this might be leveraged in the design of synthetic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Microstructure development during rapid alloy solidification

Abstract Solidification processing of structural alloys can take place over an extremely wide range of solid–liquid interface velocities spanning six orders of magnitude, from the low-velocity constitutional supercooling limit of microns/s to the high-velocity absolute stability limit of m/s. In between these two limits, the solid–liquid interface is morphologically unstable and typically forms cellular-dendritic microstructures, but also other microstructures that remain elusive. Rapid developments in additive manufacturing have renewed the interest in modeling the high-velocity range, where approximate analytical theories provide limited predictions. In this article, we discuss recent advances in phase-field modeling of rapid solidification of metallic alloys, including a brief description of state-of-the-art experiments used for model validation. We describe how phase-field models can cope with the dual challenge of carrying out simulations on experimentally relevant length- and time scales and incorporating nonequilibrium effects at the solid–liquid interface that become dominant at rapid rates. We present selected results, illustrating how phase-field simulations have yielded unprecedented insights into high-velocity interface dynamics, shedding new light on both the absolute stability limit and the formation of banded microstructures that are a hallmark of rapid alloy solidification near this limit. We also discuss state-of-the-art experiments used to validate those insights. Graphical abstract

36 MATERIALS SCIENCE↗

NASA NTRS1996

The Structural Properties of Vapor Deposited Water Ice and Astrophysical Implications

Films of vapor deposited water ice at low temperature (T<30 K) show a number of interesting structural changes during a gradual warmup. We would like to talk about the structure of the low temperature high density amorphous form of water ice, the process of crystallization, and some recent work on the morphological changes of water ice films at high temperature. The studies of the high density amorphous form are from in-situ electron microscopy as well as numerical simulations of molecular dynamics and have lead to new insights into the physical distinction between this high density amorphous form and the low density amorphous form. For the process of crystallization, we propose a model that describes the crystallization of water ice from the amorphous phase to cubic ice in terms of the nucleation of small domains in the ice. This model agrees well with the behavior of water ice in our electron microscopy studies and finds that pure water above the glass transition is a strong liquid. In more recent work, we have concentrated on temperatures above the crystallization temperature and we find interesting morphological changes related to the decrease in viscosity of the amorphous component in the cubic crystalline regime. Given enough time, we would like to put these results in an astrophysical context and discuss some observed features of the frost on interstellar grains and the bulk ice in comets.

Jenniskens, P.↗

DOE OSTI2017

Direct experimental observation of mesoscopic fluorous domains in fluorinated room temperature ionic liquids

Fluorinated room temperature ionic liquids (FRTILs) represent a class of solvent media that are attracting great attention due to their IL-specific properties as well as features stemming from their fluorous nature. Medium-to-long fluorous tails constitute a well-defined apolar moiety in the otherwise polar environment. Similarly to the case of alkyl tails, such chains are expected to result in the formation of self-assembled fluorous domains. So far, however, no direct experimental observation has been made of the existence of such structural heterogeneities on the nm scale. We report here the first experimental evidence of the existence of mesoscopic spatial segregation of fluorinated domains, on the basis of highly complementary X-ray and neutron scattering data sets (highlighting the importance of the latter probe) and NMR spectroscopy. Data are interpreted using atomistic molecular dynamics simulations, emphasizing the existence of a self-assembly mechanism that delivers segregated fluorous domains, where preferential solubilisation of fluorinated compounds can occur, thus paving the way for several smart applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Unsupervised learning of representative local atomic arrangements in molecular dynamics data

Molecular dynamics (MD) simulations present a data-mining challenge, given that they can generate a considerable amount of data but often rely on limited or biased human interpretation to examine their information content. By not asking the right questions of MD data we may miss critical information hidden within it. Here we combine dimensionality reduction (UMAP) and unsupervised hierarchical clustering (HDBSCAN) to quantitatively characterize prevalent coordination environments of chemical species within MD data. By focusing on local coordination, we significantly reduce the amount of data to be analyzed by extracting all distinct molecular formulas within a given coordination sphere. We then efficiently combine UMAP and HDBSCAN with alignment or shape-matching algorithms to partition these formulas into structural isomer families indicating their relative populations. The method was employed to reveal details of cation coordination in electrolytes based on molecular liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Accelerated path integral methods for atomistic simulations at ultra-low temperatures

Path integral methods provide a rigorous and systematically convergent framework to include the quantum mechanical nature of atomic nuclei in the evaluation of the equilibrium properties of molecules, liquids, or solids at finite temperature. Such nuclear quantum effects are often significant for light nuclei already at room temperature, but become crucial at cryogenic temperatures such as those provided by superfluid helium as a solvent. Unfortunately, the cost of converged path integral simulations increases significantly upon lowering the temperature so that the computational burden of simulating matter at the typical superfluid helium temperatures becomes prohibitive. Here we investigate how accelerated path integral techniques based on colored noise generalized Langevin equations, in particular the so-called path integral generalized Langevin equation thermostat (PIGLET) variant, perform in this extreme quantum regime using as an example the quasi-rigid methane molecule and its highly fluxional protonated cousin, CH{sub 5}{sup +}. We show that the PIGLET technique gives a speedup of two orders of magnitude in the evaluation of structural observables and quantum kinetic energy at ultralow temperatures. Moreover, we computed the spatial spread of the quantum nuclei in CH{sub 4} to illustrate the limits of using such colored noise thermostats close to the many body quantum ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Functionalization of Cathode–Electrolyte Interface with Ionic Liquids for High-Performance Quasi-Solid-State Lithium–Sulfur Batteries: A Low-Sulfur Loading Study

We introduce a quasi-solid-state electrolyte lithium-sulfur (Li–S) battery (QSSEB) based on a novel Li-argyrodite solid-state electrolyte (SSE), Super P–Sulfur cathode, and Li-anode. The cathode was prepared using a water-based carboxymethyl cellulose (CMC) solution and styrene butadiene rubber (SBR) as the binder while Li6PS5F0.5Cl0.5 SSE was synthesized using a solvent-based process, via the introduction of LiF into the argyrodite crystal structure, which enhances both the ionic conductivity and interface-stabilizing properties of the SSE. Ionic liquids (IL) were prepared using lithium bis(trifluoromethyl sulfonyl)imide (LiTFSI) as the salt, with pre-mixed pyrrolidinium bis(trifluoromethyl sulfonyl)imide (PYR) as solvent and 1,3-dioxolane (DOL) as diluent, and they were used to wet the SSE–electrode interfaces. The effect of IL dilution, the co-solvent amount, the LiTFSI concentration, the C rate at which the batteries are tested and the effect of the introduction of SSE in the cathode, were systematically studied and optimized to develop a QSSEB with higher capacity retention and cyclability. Interfacial reactions occurring at the cathode–SSE interface during cycling were also investigated using electrochemical impedance spectroscopy, cyclic voltammetry, and X-ray photoelectron spectroscopy supported by ab initio molecular dynamics simulations. This work offers a new insight into the intimate interfacial contacts between the SSE and carbon–sulfur cathodes, which are critical for improving the electrochemical performance of quasi-solid-state lithium–sulfur batteries.

25 ENERGY STORAGE↗

NASA NTRS1994

Passive microwave remote sensing of rainfall with SSM/I: Algorithm development and implementation

A physically based algorithm sensitive to emission and scattering is used to estimate rainfall using the Special Sensor Microwave/Imager (SSM/I). The algorithm is derived from radiative transfer calculations through an atmospheric cloud model specifying vertical distributions of ice and liquid hydrometeors as a function of rain rate. The algorithm is structured in two parts: SSM/I brightness temperatures are screened to detect rainfall and are then used in rain-rate calculation. The screening process distinguishes between nonraining background conditions and emission and scattering associated with hydrometeors. Thermometric temperature and polarization thresholds determined from the radiative transfer calculations are used to detect rain, whereas the rain-rate calculation is based on a linear function fit to a linear combination of channels. Separate calculations for ocean and land account for different background conditions. The rain-rate calculation is constructed to respond to both emission and scattering, reduce extraneous atmospheric and surface effects, and to correct for beam filling. The resulting SSM/I rain-rate estimates are compared to three precipitation radars as well as to a dynamically simulated rainfall event. Global estimates from the SSM/I algorithm are also compared to continental and shipboard measurements over a 4-month period. The algorithm is found to accurately describe both localized instantaneous rainfall events and global monthly patterns over both land and ovean. Over land the 4-month mean difference between SSM/I and the Global Precipitation Climatology Center continental rain gauge database is less than 10%. Over the ocean, the mean difference between SSM/I and the Legates and Willmott global shipboard rain gauge climatology is less than 20%.

Ferriday, James G.↗

DOE OSTI2021

Phase behaviours of superionic water at planetary conditions

Most water in the Universe may be superionic, and its thermodynamic and transport properties are crucial for planetary science but difficult to probe experimentally or theoretically. We use machine learning and free-energy methods to overcome the limitations of quantum mechanical simulations and characterize hydrogen diffusion, superionic transitions and phase behaviours of water at extreme conditions. Furthermore, we predict that close-packed superionic phases, which have a fraction of mixed stacking for finite systems, are stable over a wide temperature and pressure range, whereas a body-centred cubic superionic phase is only thermodynamically stable in a small window but is kinetically favoured. Our phase boundaries, which are consistent with existing—albeit scarce—experimental observations, help resolve the fractions of insulating ice, different superionic phases and liquid water inside ice giants.

36 MATERIALS SCIENCE↗