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 37 records · Page 2

Microstructure, mechanical and corrosion behaviors of AlCoCuFeNi-(Cr,Ti) high entropy alloys

The equimolar AlCoCuFeNi-(Cr,Ti) high entropy alloys (HEAs) were synthesized by nonconsumable arc melting to investigate the effects of Cr and Ti on the mechanical and corrosion properties of HEAs. The results showed that as-cast AlCoCuFeNi-(Cr,Ti) HEAs have a multi-phase microstructure, of which the solid-solution face-centered cubic (FCC), body-centered cubic (BCC) phases, and intermetallics can be observed. Ab initio molecular-dynamics (AIMD) simulations exhibit the existence of the preferred short-range ordering of Al-Ni, Co-Cr, Cr-Fe, and Ti-Co pairs in the AlCoCuFeNiCrTi liquid structure. The AIMD simulations are consistent with the experimental observation during solidification. The segregations and the FCC Cu-rich phase appear in the AlCoCuFeNiCrTi alloy, which is in agreement with AIMD calculations. The Cr addition to AlCoCuFeNi facilitates the formation of the BCC phases in the AlCoCuFeNiCr alloy, which can be explained by the larger Ω and smaller δ values. The addition of large Ti atoms facilitates the formation of the FCC phase, which is due to the fact that Ti will easily induce the breakdown of the BCC solid-solution of the AlCoCuFeNi alloy in terms of decreasing the Ω value and increasing the δ value. Finally, the Cr addition improves the corrosion resistance of AlCoCuFeNi alloys.

36 MATERIALS SCIENCE↗

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↗

Effects of Si solute on the glass formation and atomic structure of Pd liquid

Molecular dynamics simulations were performed to study the effects of Si solute on the glass formation and crystallization of Pd liquid. Pure Pd sample prepared by quenching process with cooling rate of 10 13 K/s can be in an amorphous state but the structural analysis indicates there is nearly no glass-forming motif in the sample. However, doping a small amount Si (Si concentration ~ 4%) the sample can be vitrified at a cooling rate of 10 12 K/s. The glass-forming motifs such as Pd-centered Z13, Si-centered Z9-like and Mixed-ICO-Cube clusters with 5-fold local symmetry are found to be the dominant short-range orders in the glassy samples. With the increasing of the Si-doping concentration, these glass-forming motifs tend to aggregate and connect with each other forming a network structure. Our calculated results revealed that Si solutes in liquid Pd can significantly enhance the glass-forming ability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of water confined in lyotropic liquid crystals: Molecular dynamics simulations of the dynamic structure factor

The properties of water under confinement are of practical and fundamental interest. Here in this work we study the properties of water in the self-assembled lyotropic phases of gemini surfactants with a focus on testing the standard analysis of quasi-elastic neutron scattering (QENS) experiments. In QENS experiments the dynamic structure factor is measured and fit to models to extract the translational diffusion constant, D T , and rotational relaxation time, τ R . We test this procedure by using simulation results for the dynamic structure factor, extracting the dynamic parameters from the fit as is typically done in experiments, and comparing the values to those directly measured in the simulations. We find that the decoupling approximation, where the intermediate scattering function is assumed to be a product of translational and rotational contributions, is quite accurate. The jump-diffusion and isotropic rotation models, however, are not accurate when the degree of confinement is high. In particular, the exponential approximations for the intermediate scattering function fail for highly confined water and the values of D T and τ R can differ from the measured value by as much as a factor of two. Other models have more fit parameters, however, and with the range of energies and wave-vectors accessible to QENS, the typical analysis appears to be the best choice. In the most confined lamellar phase, the dynamics are sufficiently slow that QENS does not access a large enough time scale and neutron spin echo measurements would be a valuable technique in addition to QENS.

36 MATERIALS SCIENCE↗

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↗

Structure and lifetimes in ionic liquids and their mixtures

With the aid of molecular dynamics simulations, we study the structure and dynamics of different ionic liquid systems, with focus on hydrogen bond, ion pair and ion cage formation. To do so, we report here radial distribution functions, their number integrals, and various time–correlation functions, from which we extract well-defined lifetimes by means of the reactive flux formalism. We explore the influence of polarizable force fields vs. non-polarizable ones with downscaled charges (±0.8) for the example of 1-butyl-3-methylimidazolium bromide. Furthermore, we use 1-butyl-3-methylimidazolium trifluoromethanesulfonate to investigate the impact of temperature and mixing with water as well as with the chloride ionic liquid. Smaller coordination numbers, larger distances, and tremendously accelerated dynamics are observed when the polarizable force field is applied. The same trends are found with increasing temperature. Adding water decreases the ion–ion coordination numbers whereas the water–ion and water–water coordination is enhanced. A domain analysis reveals that the nonpolar parts of the ions are dispersed and when more water is added the water clusters increase in size. The dynamics accelerate in general upon addition of water. In the ionic liquid mixture, the coordination number around the cation changes between the two anions, but the number integrals of the cation around the anions remain constant and the dynamics slow down with increasing content of the chloride ionic liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

Liquid Structure of CO 2 –Reactive Aprotic Heterocyclic Anion Ionic Liquids from X-ray Scattering and Molecular Dynamics

A combination of X-ray scattering experiments and molecular dynamics simulations were conducted to investigate the structure of ionic liquids (ILs) which chemically bind CO 2 . The structure functions were measured and computed for four different ILs consisting of two different phosphonium cations, triethyloctylphosphonium ([P 2228 ] + ) and trihexyltetradecylphosphonium ([P 66614 ] + ), paired with two different aprotic heterocyclic anions which chemically react with CO 2 , 2-cyanopyrrolide, and 1,2,4-triazolide. Simulations were able to reproduce the experimental structure functions, and by deconstructing the simulated structure functions, further information on the liquid structure was obtained. All structure functions of the ILs studied had three primary features which have been seen before in other ILs: a prepeak near 0.3–0.4 Å–1 corresponding to polar/nonpolar domain alternation, a charge alternation peak near 0.8 Å–1, and a peak near 1.5 Å–1 due to interactions of adjacent molecules. The liquid structure functions were only mildly sensitive to the specific anion and whether or not they were reacted with CO 2 . Upon reacting with CO 2 , small changes were observed in the structure functions of the [P 2228 ] + ILs, whereas virtually no change was observed upon reacting with CO 2 in the corresponding [P 66614 ] + ILs. When the [P 2228 ] + cation was replaced with the [P 66614 ] + cation, there was a significant increase in the intensities of the prepeak and adjacency interaction peak. While many of the liquid structure functions are similar, the actual liquid structures differ as demonstrated by computed spatial distribution functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Low-Dimensional Network Formation in Molten Sodium Carbonate

Molten carbonates are highly inviscid liquids characterized by low melting points and high solubility of rare earth elements and volatile molecules. An understanding of the structure and related properties of these intriguing liquids has been limited to date. We report the results of a study of molten sodium carbonate (Na2CO3) which combines high energy X-ray diffraction, containerless techniques and computer simulation to provide insight into the liquid structure. Total structure factors (F-x(Q)) are collected on the laser-heated carbonate spheres suspended in flowing gases of varying composition in an aerodynamic levitation furnace. The respective partial structure factor contributions to Fx(Q) are obtained by performing molecular dynamics simulations treating the carbonate anions as flexible entities. The carbonate liquid structure is found to be heavily temperature-dependent. At low temperatures a low-dimensional carbonate chain network forms, at T = 1100 K for example similar to 55% of the C atoms form part of a chain. The mean chain lengths decrease as temperature is increased and as the chains become shorter the rotation of the carbonate anions becomes more rapid enhancing the diffusion of Na+ ions.

Wilding, Martin C.↗

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↗

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°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↗

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↗

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

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↗

Impact of water dilution and cation tail length on ionic liquid characteristics: Interplay between polar and non-polar interactions

The recalcitrance of lignocellulosic biomass poses a major challenge that hinders the economical utilization of biomass for the production of biofuel, plastics, and chemicals. Ionic liquids have become a promising solvent that addresses many issues in both the pretreatment process and the hydrolysis of the glycosidic bond for the deconstruction of cellulosic materials. However, to make the use of ionic liquids economically viable, either the cost of ionic liquids must be reduced, or a less expensive solvent (e.g., water) may be added to reduce the overall amount of ionic liquid used in addition to reducing the viscosity of the binary liquid mixture. In this work, we employ atomistic molecular dynamics simulations to investigate the impact of water dilution on the overall liquid structure and properties of three imidazolium based ionic liquids. It is found that ionic liquid-water mixtures exhibit characteristics that can be grouped into two distinct regions, which are a function of the ionic liquid concentration. The trends observed in each region are found to correlate with the ordering in the local structure of the ionic liquid that arises from the dynamic interactions between the ion pairs. Simulation results suggest that there is a high level of local ordering in the molecular structure at high concentrations of ionic liquids that is driven by the aggregation of the cationic tails and the anion-water interactions. It is found that as the concentration of ionic liquids in the binary mixture is decreased, there is a point at which the competing self and cross interaction energies between the ionic liquid and water shifts away from a cation-anion dominated regime, which results in a significant change in the mixture properties. This break point, which occurs around 75% w/w ionic liquids, corresponds to the point at which water molecules percolate into the ionic liquid network disrupting the ionic liquids’ nanostructure. It is observed that as the cationic alkyl tail length increases, the changes in the binary mixtures’ properties become more pronounced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗