Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “liquid dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phonon dynamics in the generalized Kitaev spin liquid

Candidate materials for the Kitaev spin liquid generically have residual interactions beyond the Kitaev coupling. It therefore becomes necessary to understand how signatures of the quantum spin liquid, e.g., fractionalization of the spin excitations, are affected by the presence of these interactions. Recently it was shown that phonon dynamics is an indirect but effective probe to study fractionalized excitations in the Kitaev spin liquid. Ultrasound experiments can measure sound attenuation, which should show characteristic temperature and angular dependence of the sound-attenuation coefficient if the scattering of phonons happens predominantly on Majorana fermions. So far the computation of the sound attenuation was only done in the pure spin-phonon coupled Kitaev model, without taking into account residual interactions. In order to understand experimental signatures, here we present a mean-field study of the sound attenuation in the generalized J-K_Γ model, which is relevant to candidate materials. Further, our findings show that as long as the system is in the spin liquid phase, characteristic features of the sound attenuation remain observable even in the presence of residual interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A Method for Obtaining Liquid–Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2 O

Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas–solid adsorption in both theory and experiment, such rates are still not well studied for liquid–solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid–solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid–solid adsorption systems.

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

Role of Surface Features on the Initial Dissolution of CH 3 NH 3 PbI 3 Perovskite in Liquid Water: An Ab Initio Molecular Dynamics Study

In this study, the degradation of CH 3 NH 3 PbI 3 (MAPbI 3 ) hybrid organic inorganic perovskite (HOIP) by water has been the major issue hampering its use in commercial perovskites solar cells (PSCs) as MAPbI 3 HOIP have been known to easily degrade in the presence of water molecules. However, even though there have been numerous studies investigating this phenomenon, there is still no consensus on the mechanisms of initial stages of dissolution. Here, we attempt to consolidate differing mechanistic interpretations previously reported in the literature through the use of the first-principles constrained ab-initio molecular dynamics (AIMD) to study the mechanisms, kinetics, and thermodynamics that accompany the degradation of MAPbI 3 HOIP in liquid water. We consider not only the dissolution of the species found on the pristine MAPbI 3 HOIP surfaces, but also the dissolution from defect sites to imitate the successive degradation steps, and propose a sequence of events in the initial phase of MAPbI 3 HOIP dissolution. By comparing the dissolution free energy barrier between surface species of different surficial types, we find that the dominant dissolution mechanisms of surface species varies widely based on the specific surface features adjacent to the dissolving ion. The high sensitivity of dissolution mechanism to surface features has contributed to the many dissolution mechanisms proposed in the literature. In contrast, the dissolution free energy barriers are mainly determined by the dissolving species rather than the type of surfaces, and the type of surfaces the ions are dissolving from are inconsequential toward the dissolution free energy barrier. However, the presence of surface defects such as vacancy sites are found to significantly lower the dissolution free energy barriers. Based on the estimated dissolution free energy barriers from different types of surfaces that we investigated in this study, we proposed that the dissolution of MAPbI 3 HOIP in liquid water originates from surface defect sites that propagate laterally along the surface layer of the MAPbI 3 HOIP crystal.

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

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

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.

Coarse-grained simulations↗
DOE OSTI2020

Tuning the dynamics of imidazolium-based ionic liquids via hydrogen bonding. I. The viscous regime

Combining results from impedance spectroscopy and oscillatory shear rheology, this study focuses on the relation between the mass and charge flows and on how these are affected by the H-bonding in viscous ionic liquids (ILs). In particular, we compare the relaxational behaviors of the paradigmatic IL 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) and its OH-functionalized counterpart 1-(2-hydroxyethyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (OHEMIM-TFSI). Our results and their analysis demonstrate that the presence of cationic OH-groups bears a strong impact on the overall dynamics of OHEMIM-TFSI, although no signatures of suprastructural relaxation modes could be identified in their dielectric and mechanical responses. To check whether at the origin of this strong variation is the H-bonding or merely the difference between the corresponding cation sizes (controlling both the hydrodynamic volume and the inter-charge distance), the present study includes 1-propyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (PMIM-TFSI), mixtures of EMIM-TFSI and PMIM-TFSI with lithium bis(trifluoromethylsulfonyl)imide (Li-TFSI), and mixtures of OHEMIM-TFSI with PMIM-TFSI. Their investigation clearly reveals that the dynamical changes induced by H-bonding are significantly larger than those that can be attributed to the change in the ion size. Moreover, in the mixtures of OHEMIM-TFSI with PMIM-TFSI, a dilution of the OH-groups leads to strong deviations from ideal mixing behavior, thus highlighting the common phenomenological ground of hydroxy-functionalized ILs and other H-bonded liquids.

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

Influence of Non-Equilibrium Conditions on Liquid Hydrogen Storage Tank Behavior [Slides]

Understanding liquid hydrogen tank fluid dynamics is key for modeling liquid hydrogen systems. The tank is the source for nearly all liquid hydrogen systems. Accurate flow modeling out of the tank is needed to predict flows through downstream components. Tank contains liquid and gas that may not be at equilibrium. Questions to be addressed are: Does heat and mass transfer between liquid and vapor affect the flow rate? Is boiling an important consideration? For what conditions is a pressure relief valve (PRV) sufficient to relieve pressure and when is the burst disc needed?

08 HYDROGEN↗
DOE OSTI2023

Electrical breakdown dynamics in an argon bubble submerged in conductive liquid for nanosecond pulsed discharges

Abstract This study delves into the dynamics of cold atmospheric plasma and their interaction within conductive solutions under the unique conditions of nanosecond pulsed discharges (22 kV peak voltage, 10 ns FWHM, 4.5 kV ns −1 rate-of-rise). The research focuses on the electrical response, breakdown, and discharge propagation in an argon bubble, submerged in a NaCl solution of varying conductivity. Full or partial discharges were observed at conductivities of 1.5 µ S cm −1 (deionized water) to 1.6 mS cm −1 , but no breakdown was observed at 11.0 mS cm −1 when reducing the electrode gap. It is demonstrated that at higher conductivity electric breakdown is observed only when the gas bubble comes into direct contact with the electrode and multiple emission nodes were observed at different timescales. These nodes expanded in the central region of the bubble over timescales longer than the initial high-voltage pulse. This work offers a temporal resolution of 2 ns exposure times over the first 30 ns of the initial voltage pulse, and insight into plasma formation over decaying reflected voltage oscillations over 200 ns.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamic Motions of Topological Defects in Nematic Liquid Crystals under Spatial Confinement

Abstract Topological defects (TDs) have become a sensational topic due to their significant influence on the unusual optical and physiochemical characteristics of the material. To facilitate their application across a wide range of disciplines it is desirable to analyze and gain fundamental understanding of TDs in both equilibrium and nonequilibrium systems. Liquid crystals (LCs) are considered an ideal system for study given the direct visualization of TDs and a straightforward analyzation process. In addition to the equilibrium morphology of LC TDs, it is also of great interest to track and control the formation and annihilation of defects during thermodynamic processes. However, controlling the dynamic behavior of formed defects remains a challenge. Here, nematic LCs confined in a cell exhibiting surfaces with periodic anchoring conditions containing surface topography are relied on. The effects of patterned surface characteristics such as width, periodicity, and degree of curvature on defects dynamic motion, stabilization, and annihilation are explored. The computational experiments recapitulate the TDs transition path and provide free energy‐based predictions of critical distances for defect annihilation. Taken together, this simple approach offers a promising opportunity to control the dynamics of TDs in LCs through chemical patterned surfaces with topography.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Performance of a dynamic single bubbler in single and two-phase immiscible liquids

Ensuring nonproliferation and safeguards of special nuclear materials (SNM) is a critical aspect of advancing the nuclear fuel cycle. Traditional bubbler systems used to estimate liquid levels and densities in nuclear recycling processes have limitations, particularly in harsh environments where dip-tube corrosion and buildup necessitate frequent maintenance and recalibration. This study explores the Dynamic Single Bubbler (DSB) method, which utilizes a single dip-tube attached to a linear actuator to estimate liquid properties dynamically. This approach is extended to estimate liquid-liquid interfaces in immiscible liquids and employs a linear regression method to reduce uncertainties and improve accuracy. The DSB method achieved density estimate uncertainties of less than 0.5% and surface level estimate uncertainties typically under 0.5%, across various fluids including water, acetone, methanol, mineral oil, glycerol, and aqueous salt solutions. Results indicate that the DSB method provides accurate and robust estimates of liquid density and surface levels with minimal maintenance and without the need for calibration. Additionally, the method's applicability to immiscible liquids and various dip-tube geometries was demonstrated, showing promise for widespread use in nuclear and other industrial applications.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗
DOE OSTI2022

Regimes of evaporation and mixing behaviors of nanodroplets at transcritical conditions

The objective of this paper is to examine the fundamental mechanisms responsible for the transition between subcritical evaporation and supercritical dense-fluid-mixing in the absence of convection effects, specifically focusing on the liquid–vapor interfacial dynamics. To isolate the dynamics of this transition process, we characterize the different physical behaviors exhibited by an -dodecane nanoscale droplet placed in different nitrogen ambient conditions across the fuel’s critical point. We employ a continuum-based interface-resolving diffuse-interface method to explore the underlying phase-exchange mechanisms that bring about such distinct dynamics. Following the comparison against molecular dynamics simulations and experiments of evaporating droplets and experimental data for vapor–liquid equilibria, a parametric study at various ambient conditions and droplet sizing is performed to identify four regimes of evaporation/mixing behaviors: sub- and supercritical droplet evaporation, and sub- and supercritical dense-fluid-mixing. It is shown that the distinction in the phase-exchange mechanisms in these four regimes are brought about by the different thermodynamic phases the droplet center can exhibit during the evaporation/mixing process: subcritical liquid, supercritical liquid-like, subcritical gaseous, and supercritical gas-like, respectively. It is shown that the subcritical dense-fluid-mixing behavior is a direct result of nanoconfinement of the liquid–vapor interfacial structure and thus is not present for large droplet sizes. Finally, the present study also shows that the supercritical phase-exchange dynamics can follow two different pathways: supercritical droplet-like evaporation and supercritical dense-fluid-mixing. Furthermore, promoting the early transition to supercritical dense-fluid-mixing can significantly expedite the phase-exchange process through the disintegration of the liquid-like droplet core.

42 ENGINEERING↗
DOE OSTI2020

Isotope effects in molecular structures and electronic properties of liquid water via deep potential molecular dynamics based on the SCAN functional

Here, Feynman path-integral deep potential molecular dynamics (PI-DPMD) calculations have been employed to study both light (H 2 O) and heavy water (D 2 O) within the isothermalisobaric ensemble. In particular, the deep neural network is trained based on ab initio data obtained from the strongly constrained and appropriately normed (SCAN) exchange-correlation functional. Because of the lighter mass of hydrogen than deuteron, the properties of light water are more influenced by nuclear quantum effect than those of heavy water. Clear isotope effects are observed and analyzed in terms of hydrogen-bond structure and electronic properties of water that are closely associated with experimental observables. The molecular structures of both liquid H 2 O and D 2 O agree well with the data extracted from scattering experiments. The delicate isotope effects on radial distribution functions and angular distribution functions are well reproduced as well. Our approach demonstrates that deep neural network combined with SCAN functional based ab initio molecular dynamics provides an accurate theoretical tool for modeling water and its isotope effects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Role of Fast Dynamics in Conductivity of Polymerized Ionic Liquids

Polymerized ionic liquids (PolyILs) are promising candidates for a broad range of technologies. However, the relatively low conductivity of PolyILs at room temperature has strongly limited their applications. In this work, we provide new insights into the roles of various microscopic parameters controlling ion transport in these polymers, which are crucial for their rational design and practical applications. Using broadband dielectric spectroscopy and neutron and light scattering techniques, we found a clear connection between the activation energy for conductivity, fast dynamics, and high-frequency shear modulus in PolyILs at their glass transition temperature (T g ). In particular, our analysis reveals a correlation between conductivity and the amplitude of fast picosecond fluctuations at T g , suggesting the possible involvement of fast dynamics in lowering the energy barrier for ion conductivity. We also demonstrate that both the activation energy for ion transport and the amplitude of the fast fluctuations depend on the high-frequency shear moduli of PolyILs, thus identifying a practically important parameter for tuning conductivity. Finally, the parameters recognized in this work and their connection to the ionic conductivity of PolyILs set the stage for a deeper understanding of the mechanism of ion transport in PolyILs in the glassy state.

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

Investigation of γ - ( U , Z r ) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics

In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ–(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at. % and 100 at. %. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at. % which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ–(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ–(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. Finally, this work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks

The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.

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

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

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

Maximizing ion dynamics and electrochemical performance of ionic liquid-acetonitrile electrolyte in Ti 3 C 2 T x MXene

Modification of the structure and morphology of MXene electrodes and the formulation of the electrolytes used in their supercapacitor configurations are significant factors affecting the performance of electrochemical devices. In this study, we investigated the electrochemical performance and ion dynamics of 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [EmimTFSI], ionic liquid in the presence of acetonitrile (ACN) at different concentrations in Ti 3 C 2 T x MXene supercapacitor. We found an optimum concentration of ACN, at which more cations from the ionic liquid attach to the MXene electrode surface, providing higher electrochemical performance. This higher capacitance is also associated with increased microscopic dynamics of the cation away from the pore wall. These findings give a guideline to optimize the performance of MXene-based supercapacitors using organic solvents-ionic liquid-based electrolyte systems.

36 MATERIALS SCIENCE↗