Search NASA⌕ Search

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.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗
DOE OSTI2021

Ab initio and force field molecular dynamics study of bulk organophosphorus and organochlorine liquid structures

Here, we performed ab initio molecular dynamics (AIMD) simulations to benchmark bulk liquid structures and to evaluate results from all-atom force field molecular dynamics (FFMD) simulations with the generalized Amber force field (GAFF) for organophosphorus (OP) and organochlorine (OC) compounds. Our work also addresses the current and important topic of force field validation, applied here to a set of nonaqueous organic liquids. Our approach differs from standard treatments, which validate force fields based on thermodynamic data. Utilizing radial distribution functions (RDFs), our results show that GAFF reproduces the AIMD-predicted asymmetric liquid structures moderately well for OP compounds that contain bulky alkyl groups. Among the OCs, RDFs obtained from FFMD overlap well with AIMD results, with some offsets in position and peak structuring. However, re-parameterization of GAFF for some OCs is needed to reproduce fully the liquid structures predicted by AIMD. The offsets between AIMD and FFMD peak positions suggest inconsistencies in the developed force fields, but, in general, GAFF is able to capture short-ranged and long-ranged interactions of OPs and OCs observed in AIMD. Along with the local coordination structure, we also compared enthalpies of vaporization. Overall, calculated bulk properties from FFMD compared reasonably well with experimental values, suggesting that small improvements within the FF should focus on parameters that adjust the bulk liquid structures of these compounds.

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

The Molecular Structures of Liquid and Glassy Nifedipine and Felodipine and Their Incorporation into PVP

Background: Amorphous drug formulations are commonly used to improve the solubility and bioavailability of poorly soluble molecular pharmaceuticals, yet less is known about their molecular conformations and local bonding interactions than their crystalline phases. Methods: High-energy X-ray diffraction structure factor measurements have been made on liquid and glassy nifedipine (NIF), felodipine (FEL), NIF 1:3 polyvinylpyrrolidone (PVP), and FEL 1:3 PVP wt.% mixtures. The corresponding X-ray pair distribution functions have been interpreted using empirical potential structure refinement using different models and density functional theory conformer calculations. Results: In both NIF and FEL, the NH···O inter-molecular hydrogen bonds between the pyridyl nitrogen and ester carbonyls are found to be considerably weaker than those observed in the crystalline polymorphs. For nifedipine, it is proposed that either inter-molecular NH…ON nitro bonds are present and/or a fraction (<20%) of conformational changes, with the aryl ring flipped, occur in the liquid state. For felodipine, the models indicate significant disorder associated with the methyl and ethyl side chains in the liquid state, with the main peak intensity at 3.0 Å arising from intra-molecular Cl-Cl atom pairs. When nifedipine molecules are incorporated into PVP, our models show they possess stronger NH···O bonds to the PVP polymer than felodipine molecules, which have stronger affinity for bonding to the polymer than to other felodipine molecules. Conclusions: The amorphous forms of both NIF and FEL show much weaker hydrogen bonding than found in their crystalline phases. Liquid NIF also exhibits configurations which are not observed in the crystal phases.

Amorphous Solid Dispersion↗
DOE OSTI2019

Coordination changes in liquid tin under shock compression determined using in situ femtosecond x-ray diffraction

Little is known regarding the liquid structure of materials compressed to extreme conditions, and even less is known about liquid structures undergoing rapid compression on nanosecond timescales. Here, we report on liquid structure factor and radial distribution function measurements of tin shock compressed to 84(19) GPa. High-quality, femtosecond x-ray diffraction measurements at the Linac Coherent Light Source were used to extract the liquid diffuse scattering signal. From the radial distribution function, we find that the structural evolution of the liquid with increasing pressure mimics the evolution of the solid phase. With increasing pressure, we find that the liquid structure evolves from a complex structure, with a low coordination number, to a simple liquid structure with a coordination number of ~12. We provide a pathway for future experiments to study liquids at elevated pressures using high-energy lasers to shock compress materials beyond the reach of static diamond anvil cell techniques.

X-ray diffraction↗
DOE OSTI2022

Nature of the bonded-to-atomic transition in liquid silica to TPa pressures

First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO 2 characterize the bonded-to-atomic transition at 0.1–1.6 TPa and 10 4 –10 5 K (1–9 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. Furthermore, these results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and “rocky” (oxide) constituents of large planets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Characterization of three phases of liquid carbon by tight-binding molecular dynamics simulations

In this work, we have performed systematic molecular dynamics simulations to study the structures of liquid carbon at 5000 K with the weight density ranging from 1.4 to 3.5 g cm –3 , using a three-center tight-binding potential of carbon. The simulation results show that the bonding characteristics of the liquid changes predominately from twofold to threefold, and then to fourfold coordination as the density increases. Signals of two structural changes at the densities of about 1.9 and 3.0 g cm –3 respectively are revealed by the slope changes in the density dependence of structural, electronic and dynamical properties. Our simulation results suggest that there are three distinct liquid carbon phases in this density range. However, further thermodynamics calculations, e.g. , free energy calculations, would be required to clarify the possible liquid–liquid transitions.

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

Incorporating bubble growth volume feedback to improve simulation of the response of a structure containing liquid and gas to sudden energy input

The SNS target module is a stainless-steel structure that contains and directs mercury. The mercury is struck with short, intense proton pulses to create neutrons. The pulses deposit energy and cause loads on the target structure and mercury cavitation. Helium bubbles are introduced into the mercury to reduce loads and cavitation. This gas complicates the prediction of the target module’s physical response. The current method used to simulate the structural response to the proton pulse incorporates a simple but effective approach to approximate mercury cavitation effects. This method cannot account for the changes from injected non-condensable gas. The ability to model the target vessel response is a prerequisite for estimating its life. Known approaches for simulating a bubbly fluid mixture’s response are computationally expensive and impractical for detailed engineering models. A constitutive model is proposed for target mercury with injected gas bubbles, with the intent of providing more accurate mercury–vessel response behavior without the computing costs required to track individual bubbles. The model’s assumptions and a computational model embodiment are described. The bubbly mercury computational model is subjected to test problems to assess its potential utility. Here, the model was able to better predict the effect of gas on wave propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

First Principles Calculation of the Entropy of Liquid Aluminum

The information required to specify a liquid structure equals, in suitable units, its thermodynamic entropy. Hence, an expansion of the entropy in terms of multi-particle correlation functions can be interpreted as a hierarchy of information measures. Utilizing first principles molecular dynamics simulations, we simulate the structure of liquid aluminum to obtain its density, pair and triplet correlation functions, allowing us to approximate the experimentally measured entropy and relate the excess entropy to the information content of the correlation functions. We discuss the accuracy and convergence of the method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Liquid crystal torons in Poiseuille-like flows

Three-dimensional (3D) simulations of the structure of liquid crystal (LC) torons, topologically protected distortions of the LC director field, under material flows are rare but essential in microfluidic applications. Here, we show that torons adopt a steady-state configuration at low flow velocity before disintegrating at higher velocities, in line with experimental results. Furthermore, we show that under partial slip conditions at the boundaries, the flow induces a reversible elongation of the torons, also consistent with the experimental observations. These results are in contrast with previous simulation results for 2D skyrmions under similar flow conditions, highlighting the need for a 3D description of this LC soliton in relation to its coupling to the material flow. These findings pave the way for future studies of other topological solitons, like hopfions and heliknotons, in flowing soft matter systems.

Science & Technology - Other Topics↗
DOE OSTI2023

Water, Ion, and Proton Dynamics and Interactions in Nanoconfined Systems

The research was directed at understanding the role that nanoconfinement has on the dynamics and structure of liquids confined in nanometer size structures. Experiments were conducted using ultrafast infrared methods. The methods were two dimensional infrared (2D IR) spectroscopy and polarization pump-probe (PSPP) experiments. In addition FT-IR spectroscopy was used. To augment the experimental work, high level molecular dynamics simulations of the systems were conducted. The system used for confinement were nanoporous silica with essentially monodispersed pores. Pores of different diameters were studies, i.e., 2.5 nm to 8 nm. Water and organic liquids were investigated. The experiments were performed on a vibrational probe molecule, i.e., the CN stretch of the anion SeCN-. The 2D IR experiments measured spectral diffusion, which directly reported on the structural dynamics of the confined liquids. 2D IR was also used to perform chemical exchange measurements, which observed the formation and dissociation of molecular complexes with surface silanols. The PSPP experiments measured orientational relaxation of the probe molecule as well as it vibrational lifetime. All of the observables depended on the diameter of the pores and the type of liquid. The simulations of the observables provided detailed structural information on the structure of the confined liquids and there interactions with the surfaces. In addition to the study of the confinement in nanopores, water in hydrogels was also studied. The same experimental methods were used along with MD simulations. It was found the influence of hydrogels on water did not come from confinement of the water in the hydrogel polymer network, but rather from the interactions of the water with the individual monomers that compose the polymer. In addition, the temperature dependence of a polymer/water system was tracked from below to above the lower critical solution temperature to observe the nature of the phase separated polymer and how water responded to the phase separation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Effect of an external electric field on the dynamics and intramolecular structures of ions in an ionic liquid

Simulations of the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide in an external electric field have been performed using a Drude particle polarizable force field. The structure of the ions has been analyzed, with close attention paid to the configurations of the ions. The “charge arm” concept is used to explain certain changes of these configurations that would be difficult to rationalize otherwise, e.g., trans → cis isomerization of the bis(trifluoromethylsulfonyl)imide anion and extension of the alkyl chain of the cation. It has also been shown that the ions orient themselves so that their charge arms align with and stretch out along the field, and these effects occur at lower external electric field strengths than cause a change in the inherent diffusion of the ions. Here, the dynamics of the system parallel and perpendicular to the field were analyzed, and it was found that the applied field affected the diffusion normal to the field. This is explained as a secondary effect of a change in the ion cage formed by the surrounding counterions of a given ion in the ionic liquid. The breakdown of the ion cages was rationalized by correlating changes to the inherent diffusion of the ions with other changes to the diffusion and bulk structure of the liquid, as well as considering the average forces on the ions compared to the force the ions would be expected to experience in an electric field. Parallel to the field, a drift was observed at every electric field studied. In electric fields with no changes to the ion cage structure, the relationship between the drift and electric field was found to be linear, becoming nonlinear as the ion cage structure breaks down.

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

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Tricoupled hybrid lattice Boltzmann model for nonisothermal drying of colloidal suspensions in micropore structures

A tricoupled hybrid lattice Boltzmann model (LBM) is developed to simulate colloidal liquid evaporation and colloidal particle deposition during the nonisothermal drying of colloidal suspensions in micropore structures. An entropic multiple-relaxation-time multirange pseudopotential two-phase LBM for isothermal interfacial flow is first coupled to an extended temperature equation for simulating nonisothermal liquid drying. Then the coupled model is further coupled with a modified convection diffusion equation to consider the nonisothermal drying of colloidal suspensions. Two drying examples are considered. First, drying of colloidal suspensions in a two-pillar micropore structure is simulated in two dimensions (2D), and the final configuration of colloidal particles is compared with the experimental one. Good agreement is observed. Second, at the temperature of 343.15 K ( 70&deg;C ), drying of colloidal suspensions in a complex spiral-shaped micropore structure containing 220 pillars is simulated (also in 2D). The drying pattern follows the designed spiral shape due to capillary pumping, i.e., transport of the liquid from larger pores to smaller ones by capillary pressure difference. Because the colloidal particles are passively carried with liquid, they accumulate at the small menisci as drying proceeds. As liquid evaporates at the small menisci, colloidal particles are deposited, eventually forming solid structures between the pillars (primarily), and at the base of the pillars (secondarily). As a result, the particle deposition conforms to the spiral route. Qualitatively, the simulated liquid and particle configurations agree well with the experimental ones during the entire drying process. Quantitatively, the model demonstrates that the evaporation rate and the particle accumulation rate slowly decrease during drying, similar to what is seen in the experimental results, which is due to the reduction of the liquid-vapor interfacial area. In conclusion, the hybrid model shows the capability and accuracy for simulating nonisothermal drying of colloidal suspensions in a complex micropore structure both qualitatively and quantitatively, as it includes all the required physics and captures all the complex features observed experimentally. Such a tricoupled LBM has a high potential to become an efficient numerical tool for further investigation of real and complex engineering problems incorporating drying of colloidal suspensions in porous media.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Computational characterization of the secondary droplets formed during the impingement of a train of ethanol drops

This paper uniquely characterizes the secondary droplets formed during the impingement of a train of ethanol drops, using three-dimensional direct numerical simulations performed under conditions studied experimentally by Yarin and Weiss. Our numerical results have been previously validated against experimental data demonstrating the ability to accurately capture the splashing dynamics. In this work, the predictive ability of the model is leveraged to gain further insight into secondary droplet formation. We present a robust post-processing algorithm, which scrutinizes the liquid volume fraction field in the volume-of-fluid method and quantifies the number, volume and velocity of secondary droplets. The high-resolution computational simulations enable secondary droplet characterization within close proximity of the impingement point at small length and time scales, which is extremely challenging to achieve experimentally. By studying the temporal evolution of secondary droplet formation, direct connections are made between liquid structures seen in the simulation and the instantaneous distribution of secondary droplets, leading to detailed insight into the instability-driven breakup process of lamellae. Time-averaged secondary droplet characteristics are also studied to describe the global distribution of secondary droplets. Such analysis is vital to understanding fuel drop impingement in direct injection engines, facilitating the development of highly accurate spray–wall interaction models for use in Lagrangian solvers.

42 ENGINEERING↗
DOE OSTI2019

EXPLORING THE STRUCTURE OF GLASS-FORMING LIQUIDS USING HIGH ENERGY X-RAY DIFFRACTION, CONTAINERLESS METHODOLOGY AND STATE-OF-THE-ART MOLECULAR DYNAMICS SIMULATION

High energy X-ray diffraction can be combined with containerless techniques to provide information on the atomic arrangements in glass-forming liquids in stable and metastable regimes. The high incident energies provide bulk diffraction data to high values of scattering vector which enables significantly more robust analysis of the local and medium-range order that influences important physical properties such as viscosity and crystal nucleation. These combined techniques have been applied to a range of oxide liquids. In this contribution we illustrate addition of further dimensions to phase space by controlling the partial pressure of oxygen that permits the study liquids containing iron. The advantages of rapid data acquisition are also demonstrated in a study of tellurite glass-forming systems where a transition from ergodic to non-ergodic regimes in the deeply supercooled liquid is shown. Finally we demonstrate how descriptions of the liquid structure can be developed by combining HEXRD with molecular dynamics simulations.

containerless techniques↗
DOE OSTI2024

Detecting underscreening and generalized Kirkwood transitions in aqueous electrolytes

We establish the connection between the measured small angle x-ray scattering signal and the charge–charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a 0 , and the inverse periodicity length, Q 0 . Above the KTs, we find a universal scaling of a 0 ∝ c – $\sqrt{ζ/3}$ and Q 0 ∝ c 1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗