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 559 records · Page 31

Reentrant glass transition leading to ultrastable metallic glass

Polyamorphs are often observed in amorphous matters, and a representative example is the reentrant glass transition in colloid systems. For metallic amorphous alloys, however, the cases reported so far are limited to metallic glasses (MGs) that undergo electronic transitions under gigapascal applied pressure, or the presence of two liquids at the same composition. In this paper, we report the first observation of a reentrant glass transition in MGs. This unusual reentrant glass transition transforms an MG from its as-quenched state (Glass I) to an ultrastable state (Glass II), mediated by the supercooled liquid of Glass I. Specifically, upon heating to above its glass transition temperature under ambient pressure, Glass I first transitions into its supercooled liquid, which then transforms into a new Glass II, accompanied by an exothermic peak in calorimetric scan, together with a precipitous drop in volume, electrical resistance and specific heat, as well as clear evidence of local structural ordering on the short-to-medium-range scale revealed via in-situ synchrotron X-ray scattering. Atomistic simulations indicate enhanced ordering of locally favored motifs to establish correlations in the medium range that resemble those in equilibrium crystalline compounds. The resulting lower-energy Glass II has its own glass transition temperature higher than that of Glass I by as much as 50 degrees. This route thus delivers a thermodynamically and kinetically ultrastable MG that can be easily retained to ambient conditions.

36 MATERIALS SCIENCE↗

Topology-enhanced mechanical stability of swelling nanoporous electrodes

Abstract Materials like silicon and germanium offer a 10-fold improvement in charge capacity over conventional graphite anodes in lithium-ion batteries but experience a roughly threefold volume increase during lithiation, which challenges ensuring battery integrity. Nanoporous silicon, created by liquid-metal-dealloying, is a potentially attractive anode design to mitigate this challenge, exhibiting both higher capacity and extended cycle lifetimes. However, how nanoporous structures accommodate the large volume change is unknown. Here, we address this question by using phase-field modeling to produce nanoporous particles and to investigate their elastoplastic swelling behavior and fracture. Our simulations show that enhanced mechanical stability results from the network topology consisting of ligaments connected by bulbous, sphere-like nodes. The ligaments forcefully resist elongation while the nodes, behaving like isolated spherical particles, experience large stresses driving fracture. However, being smaller compared to a sphere of the same volume as the entire nanoporous particle, the nodes are more protected against fracture.

36 MATERIALS SCIENCE↗

Carbon diffusion in molten uranium: an ab initio molecular dynamics study

In this work we used ab initio molecular dynamics (AIMD) within the framework of density functional theory (DFT) and the projector-augmented wave (PAW) method to study carbon diffusion in liquid uranium at temperatures above 1600 K. The electronic interactions of carbon and uranium were described using the local density approximation (LDA). The self-diffusion of uranium based on this approach is compared with literature computational and experimental results for liquid uranium. The temperature dependence of carbon and uranium diffusion in the melt was evaluated by fitting the resulting diffusion coefficients to an Arrhenius relationship. We found that the LDA calculated activation energy for carbon was nearly twice that of uranium: 0.55±0.03 eV for carbon compared to 0.32±0.04 eV for uranium. Structural analysis of the liquid uranium-carbon system is also discussed.

Molecular dynamics computer simulation, Nuclear Pr↗

Measurement of the (e,e') cross section for12C, 48Ti, 27Al and 40Ar

Many experiments have found that neutrinos have three different flavors, and they can transform from each other. Each flavor is a mixture of three nonzero mass states. In neutrino oscillation experiments, the uncertainty of the extracted oscillation parameters is dependent on the neutrino-nucleus cross section. In recent years various experiments provided precise neutrino and anti-neutrino cross section measurements and they have shown the inability of current Monte Carlo simulations to reproduce the data. A more accurate nuclear structure model of the target employed in neutrino detector is urgently needed to improve the precision of long baseline neutrino experiments. Electron scattering has been one of the most powerful methods of obtaining information about cross sections in the past. In the upcoming deep underground neutrino experiment (DUNE), Liquid Argon Time Projection Chambers (LArTPCs) will be used as the detector technology and argon will be used as the nuclear target. There are very few data available on argon and there is not an accurate nuclear model that describes them accurately, at ~5% level. We have performed an electron-argon scattering experiment at Jefferson Lab (E12-14-012) to provide accurate data that will help build a reliable nuclear model for describing neutrino-argon scattering. The JLab E12-14-012 experiment has successfully taken data in Hall A at the Thomas Jefferson National Accelerator Facility between February and March of 2017. We collected data for the inclusive (e,e') and exclusive (e,e'p) processes for a variety of targets (argon, titanium, aluminum, carbon) at a wide range of kinematic settings. This thesis will present the results of the inclusive double differential cross sections for carbon, titanium, aluminum and argon at beam energy E = 2:222 GeV and scattering angle theta = 15:541 deg with total uncertainties smaller than 5%. The kinematic settings covered a broad range of energy transfers, including deep inelastic scattering (DIS), delta production, and quasielastic scattering.

Dai, Hongxia↗

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen oxygen partial correlation function at the experimental equilibrium density and 311 +/- 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.

Gaiduk, Alex P.↗

Nonadditive Interactions Unlock Small-Particle Mobility in Binary Colloidal Monolayers

We examine the organization and dynamics of binary colloidal monolayers composed of micron-scale silica particles interspersed with smaller-diameter silica particles that serve as minority component impurities. These binary monolayers are prepared at the surface of ionic liquid droplets over a range of size ratios (σ = 0.16–0.66) and are studied with low-dose minimally perturbative scanning electron microscopy (SEM). The high resolution of SEM imaging provides direct tracking of all particle coordinates over time, enabling a complete description of the microscopic state. In these bidisperse size mixtures, particle interactions are nonadditive because interfacial pinning to the droplet surface causes the equators of differently sized particles to lie in separate planes. By varying the size ratio, we control the extent of nonadditivity in order to achieve phase behavior inaccessible to additive 2D systems. Across the range of size ratios, we tune the system from a mobile small-particle phase (σ < 0.24) to an interstitial solid (0.24 < σ < 0.33) and furthermore to a disordered glass (σ > 0.33). These distinct phase regimes are classified through measurements of hexagonal ordering of the large-particle host lattice and the lattice’s capacity for small-particle transport. Altogether, we explain these structural and dynamic trends by considering the combined influence of interparticle interactions and the colloidal packing geometry. Our measurements are reproduced in molecular dynamics simulations of 2D nonadditive disks, suggesting an efficient method for describing confined systems with reduced dimensionality representations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid–Liquid Equilibrium Prediction in Fast Pyrolysis Bio-Oil Systems: A Framework for Incorporating Bio-Oil Complexity

The study of mixtures of bio-oil, water and organic solvents in different proportions can serve as a cost-effective analysis of its content due to the formation of immiscible phases. This manuscript attempts to replicate experimentally determined partition coefficients (K OW ) of relevant species present in fast pyrolysis bio-oil (FPBO). A commercial flowsheeting simulator with surrogate bio-oil model representation is used. Concurrently, pyrolytic lignins in FPBO (‘pyrolignin’) do not have an agreed-upon structural representation, and the literature is ripe with wide variations of said representations. Thus, during the description of FPBO, this pyrolignin fraction was modeled using 20 possible structures (phenolic dimers to tetramers), with the goal of determining the structures for which the experimental data are best described. Two cases were considered: Case 1 normalized the reported experimental mass balance, while Case 2 included the unreported fraction in the mass balance to the total pyroligin. Please, add here a comment on the prediction of the Water oil equilibrium. The best KOW predictions for levoglucosan (LVG) were obtained when the system was modeled with no pyrolignin, presenting an MRE under 10% for both systems WO and BO. Among the possible structures, D2 (dimer), F1 (trimer) and I1, and I3 (tetramers) presented MRE ≤ 13% for both cases.

09 BIOMASS FUELS↗

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗

Design of enhanced porous organic cage solubility in Type 2 porous liquids

Decreasing cost of technologies for direct air capture of carbon can be achieved through the design of new materials with high CO 2 selectivity that can be incorporated into existing industrial processes. An emerging class of materials for these applications are porous liquids (PLs). PLs are mixtures of porous hosts and solvents with intrinsic porosity due to steric exclusion of solvent from inside the porous host. It is currently unknown how solvent -porous host interactions affect porous host solubility in the bulk solvent. Here, in this study, density functional theory simulations were used to investigate interactions between nine solvents and a CC13 porous organic cage (POC). Calculations identified that solvent molecules were the most stable when placed either inside the CC13 POC or in the pore window compared to interfacial binding sites. Structural changes to the CC13 POC correlated with reported experimental solubilities, including expansion of the CC13 POC with solvent molecule infiltration and expansion or contraction of the pore window. Based on these results, new PL design guidelines should include compositions with (1) high concentrations of POCs with flexible cage structures that can expand when solvated and (2) solvent molecule-POC combinations that contract the pore window during solvent molecule-host binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of Rotational Dynamics of Chiral Liquid Crystal Droplets in an Electric Field

The dynamics of the orientational structure of chiral nematic (CN) droplets in an isotropic medium in dc and ac electric fields is investigated by the polarized light microscopy technique. It is shown theoretically that the dynamics of rotational processes in these kinds of systems is determined by electroconvective processes developing due to the flexoelectric polarization associated with the initial configuration of the director field in droplets. It is established experimentally that the linear and quadratic regions of dependence of the rotational velocity of droplets on the electric field strength are explained by the above-mentioned mechanisms. Numerical simulation on the basis of the approach developed gives good agreement with experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hierarchical TiO 2 :Cu 2 O Nanostructures for Gas/Vapor Sensing and CO 2 Sequestration

We investigate the effect of high-surface-area self-assembled TiO 2 :Cu 2 O nanostructures for CO 2 and relative humidity gravimetric detection using polyethylenimine (PEI), 1-ethyl-3-methylimidazolium (EMIM), and polyacrylamide (PAAm). Introduction of hierarchical TiO 2 :Cu 2 O nanostructures on the surface of quartz crystal microbalance sensors is found to significantly improve sensitivity to CO 2 and to H 2 O vapor. The response of EMIM to CO 2 increases fivefold for 100 nm-thick TiO 2 :Cu 2 O as compared to gold. At ambient CO 2 concentrations, the hierarchical assembly operates as a sensor with excellent reversibility, while at higher pressures, the CO 2 desorption rate decreases, suggesting possible application for CO 2 sequestration under these conditions. The gravimetric response of PEI to CO 2 increases by a factor of 3 upon introduction of a 50 nm TiO 2 :Cu 2 O layer. The PAAm gravimetric response to water vapor also increases by a factor of 3 and displays improved reversibility with the addition of 50 nm TiO 2 :Cu 2 O structures. We found that TiO 2 :Cu 2 O can be used to lower the detection limits for CO 2 sensing with EMIM and PEI and lower the detection limits for H 2 O sensing with PAAm by over a factor of 2. Coarse-grained and all-atom molecular dynamics simulations indicate the dissociative character of ionic liquid assembly on TiO 2 :Cu 2 O interfaces and different distributions of CO 2 and H 2 O molecules on bare and ionic liquid-coated surfaces, confirming experimental observations. Overall, our results show high potential of hierarchical assemblies of TiO 2 :Cu 2 O/room temperature ionic liquid and polymer films for sensors and CO 2 sequestration.

36 MATERIALS SCIENCE↗

Amino Acid Sequence Controls Enhanced Electron Transport in Heme-Binding Peptide Monolayers

Metal-binding proteins have the exceptional ability to facilitate long-range electron transport in nature. Despite recent progress, the sequence-structure–function relationships governing electron transport in heme-binding peptides and protein assemblies are not yet fully understood. In this work, the electronic properties of a series of heme-binding peptides inspired by cytochrome bc1 are studied using a combination of molecular electronics experiments, molecular modeling, and simulation. Self-assembled monolayers (SAMs) are prepared using sequence-defined heme-binding peptides capable of forming helical secondary structures. Following monolayer formation, the structural properties and chemical composition of assembled peptides are determined using atomic force microscopy and X-ray photoelectron spectroscopy, and the electronic properties (current density–voltage response) are characterized using a soft contact liquid metal electrode method based on eutectic gallium–indium alloys (EGaIn). Our results show a substantial 1000-fold increase in current density across SAM junctions upon addition of heme compared to identical peptide sequences in the absence of heme, while maintaining a constant junction thickness. These findings show that amino acid composition and sequence directly control enhancements in electron transport in heme-binding peptides. Overall, this study demonstrates the potential of using sequence-defined synthetic peptides inspired by nature as functional bioelectronic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing phonon focusing, thermomechanical behavior, and moiré patterns in van der Waals architectures using surface acoustic waves

Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS 2 . Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS 2 bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS 2 monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS 2 by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS 2 layers.

36 MATERIALS SCIENCE↗

Nanoscale Structure and Dynamics in Geochemical Systems

Neutron scattering is a powerful tool to elucidate the structure and dynamics of systems that are important to geochemists, including ion association in complex aqueous solutions, solvent-exchange reactions at mineral–water interfaces, and reaction and transport of fluids in nanoporous materials. This article focusses on three techniques: neutron diffraction, which can reveal the atomic-level structure of aqueous solutions and solids; quasi-elastic neutron scattering, which measures the diffusional dynamics at mineral–water interfaces; and small-angle neutron scattering, which can show how properties of nanoporous systems change during gas, liquid, and solute imbibition and reaction. Furthermore, the usefulness and applicability of the experimental results are extended by rigorous comparison to computational simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultra-fast interpretable machine-learning potentials

Abstract All-atom dynamics simulations are an indispensable quantitative tool in physics, chemistry, and materials science, but large systems and long simulation times remain challenging due to the trade-off between computational efficiency and predictive accuracy. To address this challenge, we combine effective two- and three-body potentials in a cubic B-spline basis with regularized linear regression to obtain machine-learning potentials that are physically interpretable, sufficiently accurate for applications, as fast as the fastest traditional empirical potentials, and two to four orders of magnitude faster than state-of-the-art machine-learning potentials. For data from empirical potentials, we demonstrate the exact retrieval of the potential. For data from density functional theory, the predicted energies, forces, and derived properties, including phonon spectra, elastic constants, and melting points, closely match those of the reference method. The introduced potentials might contribute towards accurate all-atom dynamics simulations of large atomistic systems over long-time scales.

36 MATERIALS SCIENCE↗

Viscosity in water from first-principles and deep-neural-network simulations

Abstract We report on an extensive study of the viscosity of liquid water at near-ambient conditions, performed within the Green-Kubo theory of linear response and equilibrium ab initio molecular dynamics (AIMD), based on density-functional theory (DFT). In order to cope with the long simulation times necessary to achieve an acceptable statistical accuracy, our ab initio approach is enhanced with deep-neural-network potentials (NNP). This approach is first validated against AIMD results, obtained by using the Perdew–Burke–Ernzerhof (PBE) exchange-correlation functional and paying careful attention to crucial, yet often overlooked, aspects of the statistical data analysis. Then, we train a second NNP to a dataset generated from the Strongly Constrained and Appropriately Normed (SCAN) functional. Once the error resulting from the imperfect prediction of the melting line is offset by referring the simulated temperature to the theoretical melting one, our SCAN predictions of the shear viscosity of water are in very good agreement with experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗