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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sub‐5 Ångstrom Porosity Tuning in Calixarene‐Derived Porous Liquids via Supramolecular Complexation Construction

Abstract Sub‐Ångstrom‐level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub‐5 Ångstrom scale is created via sophisticated porosity tuning in calixarene‐derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy‐, X‐ray‐, and neutron‐scattering techniques. The presence of permanent porosity in calixarene‐derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub‐5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. The approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

Li, Errui [Department of Chemistry University of T↗
DOE OSTI2024

Numerical Simulation of Irradiation Induced Swelling for STAR 4.1 Blanket

Here, this article presents the preliminary results of FEM implementation of irradiation-induced swelling, hardening, and creep effect using ANSYS user programmable features (UPFs) applied to material properties. The first liquid metal breeding blanket model for STAR 4.1 tokamak as a part of the virtual prototyping system is under development at Princeton Plasma Physics Laboratory (PPPL). Intense neutron irradiation produces significant changes in the physical and mechanical properties of Fe-(8%–9%)Cr-based reduced activation ferritic martensitic (RAFM) steels. An ANSYS model of the breeding blanket was built to involve these irradiation-induced material property changes with neutron fluence distribution mapped from the MCNP model. Simulation results show that significant structure deformation forms from nuclear swelling, stress, and plastic strain arising from DPA gradient through wall thickness. The material degradation effect is not negligible. This modeling feature can either assist in analyzing the structural behaviors with the influence of nuclear swelling or provide guidance to design the structure to withstand irradiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Design and engineering of dispersion and loss in photonic crystal fiber 1 × 4 power splitter (PCFPS) based on hole size alteration and optofluidic infiltration

We have presented a technique based on optofluidic infiltration and air-holes diameter variation together to design a 1 × 4 photonic crystal fiber power splitter (PCFPS) which have very low dispersion (D): (0 ≤ D ≤ 2.5 (ps/nm/km) and very low loss (L): 0 ≤ L ≤ 0.025(dB/cm) in a wide range of wavelengths (1100–1700 nm). This approach allows us to control the dispersion of the fundamental mode in a PCF beam splitter by choosing appropriate refractive indices for liquids and suitable diameters for air-holes in PCF power splitter. In fact, the techniques, used in this paper are complementary of each other and give us more excellent results which are better than other reported results in researchers’ works so far. In this work, a new design of 1 × 4 photonic crystal fiber power splitter is proposed by using beam propagation method. An optical Gaussian signal at a wavelength of third communication window range (1550 nm) is inserted into the central core and equally is divided into four core (25% of the total input power interred to each core). In addition, the physical behavior of coupling characteristics is obtained by using coupled mode analysis. Numerical simulations show that input optical signal can be equally divided in photonic crystal fiber structure with low dispersion and low loss. The total size of proposed PCFPS is 30 µm × 30 µm × 1.2 mm, too.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Study of the Effect of Mold Corner Shape on the Initial Solidification Behavior of Molten Steel Using Mold Simulator

The chamfered mold with a typical corner shape (angle between the chamfered face and hot face is 45 deg) was applied to the mold simulator study in this paper, and the results were compared with the previous results from a well-developed right-angle mold simulator system. The results suggested that the designed chamfered structure would increase the thermal resistance and weaken the two-dimensional heat transfer around the mold corner, causing the homogeneity of the mold surface temperatures and heat fluxes. In addition, the chamfered structure can decrease the fluctuation of the steel level and the liquid slag flow around the meniscus at mold corner. The cooling intensities at different longitudinal sections of shell are close to each other due to the similar time-average solidification factors, which are 2.392 mm/s{sup 1/2} (section A-A: chamfered center), 2.372 mm/s{sup 1/2} (section B-B: 135 deg corner), and 2.380 mm/s{sup 1/2} (section D-D: face), respectively. For the same oscillation mark (OM), the heights of OM roots at different positions (profile L1 (face), profile L2 (135 deg corner), and profile L3 (chamfered center)) are very close to each other. The average value of height difference (H{sub D}) between two OMs roots for L1 and L2 is 0.22 mm, and for L2 and L3 is 0.38 mm. Finally, with the help of metallographic examination, the shapes of different hooks were also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗
DOE OSTI2026

Probing the atomic dynamics of ultrafast melting with femtosecond electron diffraction

Melting is an every-day phase transition that is determined by thermodynamic parameters like temperature and pressure. In contrast, ultra-fast melting is governed by the microscopic response to a rapid energy input and, thus, can reveal the strength and dynamics of atomic bonds as well as the energy flow rate to the lattice. Accurately describing these processes remains challenging and requires detailed insights into transient states encountered. Here, we present data from femtosecond electron diffraction measurements that capture the structural evolution of copper during the ultrafast solid-to-liquid phase transformations. At absorbed energy densities 2-4 times the melting threshold, melting begins at the surface slightly below the nominal melting point followed by rapid homogeneous melting throughout the volume. Molecular dynamics simulations reproduce these observations and reveal a weak electron-lattice energy transfer rate for the given experimental conditions. Both simulations and experiments show no indications of rapid lattice collapse when its temperature surpasses proposed limits of superheating, providing evidence that the inherent dynamics limits the speed of disordering in ultrafast melting of metals.

FOS: Physical sciences↗
DOE OSTI2023

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chirality and dislocation effects in single nanostructures probed by whispering gallery modes

Nanostructures such as nanoribbons and -wires are of interest as components for building integrated photonic systems, especially if their basic functionality as dielectric waveguides can be extended by chiroptical phenomena or modifications of their optoelectronic properties by extended defects, such as dislocations. However, conventional optical measurements typically require monodisperse (and chiral) ensembles, and identifying emerging chiral optical activity or dislocation effects in single nanostructures has remained an unmet challenge. Here we show that whispering gallery modes can probe chirality and dislocation effects in single nanowires. Wires of the van der Waals semiconductor germanium(II) sulfide (GeS), obtained by vapor–liquid–solid growth, invariably form as growth spirals around a single screw dislocation that gives rise to a chiral structure and can modify the electronic properties. Cathodoluminescence spectroscopy on single tapered GeS nanowires containing joined dislocated and defect-free segments, augmented by numerical simulations and ab-initio calculations, identifies chiral whispering gallery modes as well as a pronounced modulation of the electronic structure attributed to the screw dislocation. Further, our results establish chiral light-matter interactions and dislocation-induced electronic modifications in single nanostructures, paving the way for their application in multifunctional photonic architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Shock growth of ice crystal near equilibrium melting pressure under dynamic compression

Crystal growth is governed by an interplay between macroscopic driving force and microscopic interface kinetics at the crystal–liquid interface. Unlike the local equilibrium growth condition, the interplay becomes blurred under local nonequilibrium, which raises many questions about the nature of diverse crystal growth and morphological transitions. Here, we systematically control the growth condition from local equilibrium to local nonequilibrium by using an advanced dynamic diamond anvil cell (dDAC) and generate anomalously fast growth of ice VI phase with a morphological transition from three- to two-dimension (3D to 2D), which is called a shock crystal growth. Unlike expected, the shock growth occurs from the edges of 3D crystal along the (112) crystal plane rather than its corners, which implies that the fast compression yields effectively large overpressure at the crystal–liquid interface, manifesting the local nonequilibrium condition. Molecular dynamics (MD) simulation reproduces the faster growth of the (112) plane than other planes upon applying large overpressure. Moreover, the MD study reveals that the 2D shock crystal growth originates from the similarity of the interface structure between water and the (112) crystal plane under the large overpressure. As a result, this study provides insight into crystal growth under dynamic compressions, which makes a bridge for the unknown behaviors of crystal growth between under static and dynamic pressure conditions.

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

Classical and quantum spin dynamics of the honeycomb Γ model

Quantum-to-classical crossover is a fundamental question in dynamics of quantum many-body systems. In frustrated magnets, for example, it is highly nontrivial to describe the crossover from the classical spin liquid with a macroscopically degenerate ground-state manifold to the quantum spin liquid phase with fractionalized excitations. This is an important issue, as we often encounter the demand for a sharp distinction between the classical and quantum spin liquid behaviors in real materials. Here we take the example of the classical spin liquid in a frustrated magnet with novel bond-dependent interactions to investigate the classical dynamics, and critically compare it with quantum dynamics in the same system. In particular, we focus on signatures in the dynamical spin structure factor. Combining Landau-Lifshitz dynamics simulations and the analytical Martin-Siggia-Rose approach, we show that the low-energy spectra are described by relaxational dynamics and highly constrained by the zero mode structure of the underlying degenerate classical manifold. Further, the higher energy spectra can be explained by precessional dynamics. Surprisingly, many of these features can also be seen in the dynamical structure factor in the quantum model studied by finite-temperature exact diagonalization. Here, we discuss the implications of these results and their connection to recent experiments on frustrated magnets with strong spin-orbit coupling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Thermo-Fluid Modeling Framework for Supercomputer Digital Twins: Part 1, Demonstration at Exascale

A thermo-fluid modeling framework is being developed for ExaDigiT---an open-source framework for developing comprehensive digital twins of liquid-cooled supercomputers. The work is being conducted in two parts, and discussion is divided into two companion papers. The work documented in this paper focuses on the development of a cooling system library in Dymola for the Frontier supercomputer at Oak Ridge National Laboratory. The second part, outlined in a companion paper, focuses on a templating structure called Auto-CSM for easily creating model-agnostic, physics-based thermo-fluid cooling system models for liquid-cooled supercomputers using a text-based schema. The cooling model is being developed using primarily the open-source Transient Simulation Framework of Reconfigurable Models (TRANSFORM) library. The library follows the templating architecture developed within the TRANSFORM library for modeling subsystems. A full-system validation was performed to validate a very simple model that is integrated with the system controls, and the results are presented herein.

Kumar, Vineet↗
DOE OSTI2025

Selective and Stable Ethanol Synthesis via Electrochemical CO 2 Reduction in a Solid Electrolyte Reactor

Electrochemical CO 2 reduction to ethanol faces challenges such as low selectivity, a product mixture with liquid electrolyte, and poor catalyst/reactor stability. Here, we developed a grain-rich zinc-doped Cu 2 O precatalyst that presented a high ethanol Faradaic efficiency of over 40% under a current density of 350 mA·cm –2 . Our density functional theory (DFT) simulation suggested that Zn atoms inside the structure have a greater carbophilicity than the Cu atoms to help facilitate *CHCHO formation, a key reaction intermediate toward ethanol instead of other C 2 products. Here, a high Faradaic efficiency ratio between ethanol and ethylene (FE EtOH /FE C2H4 ) reached 2.34 in the zinc-doped Cu 2 O precatalyst, representing an over 4-fold improvement compared to bare Cu 2 O precatalyst. By integrating this Cu-based catalyst into a porous solid electrolyte (PSE) reactor with a salt-managing design, we achieved stable ethanol production for over 180 h under a current density of 250 mA·cm –2 while maintaining ethanol selectivity at ~30%.

09 BIOMASS FUELS↗
DOE OSTI2021

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗
DOE OSTI2018

Experimental and Computational Investigation of Subcritical Near-Nozzle Spray Structure and Primary Atomization in the Engine Combustion Network Spray D

In order to improve understanding of the primary atomization process for diesel-like sprays, a collaborative experimental and computational study was focused on the near-nozzle spray structure for the Engine Combustion Network (ECN) Spray D single-hole injector. These results were presented at the 5th Workshop of the ECN in Detroit, Michigan. Application of x-ray diagnostics to the Spray D standard cold condition enabled quantification of distributions of mass, phase interfacial area, and droplet size in the near-nozzle region from 0.1 to 14 mm from the nozzle exit. Using these data, several modeling frameworks, from Lagrangian-Eulerian to Eulerian-Eulerian and from Reynolds-Averaged Navier-Stokes (RANS) to Direct Numerical Simulation (DNS), were assessed in their ability to capture and explain experimentally observed spray details. Due to its computational efficiency, the Lagrangian-Eulerian approach was able to provide spray predictions across a broad range of conditions. In general, this “engineering-level” simulation was able to reproduce the details of the droplet size distribution throughout the spray after calibration of the spray breakup model constants against the experimental data. Complementary to this approach, higher-fidelity modeling techniques were able to provide detailed insight into the experimental trends. For example, interface-capturing multiphase simulations were able to capture the experimentally observed bimodal behavior in the transverse interfacial area distributions in the near-nozzle region. Further analysis of the spray predictions suggests that peaks in the interfacial area distribution may coincide with regions of finely atomized droplets, whereas local minima may coincide with regions of continuous liquid structures. Lastly, the results from this study highlight the potential of x-ray diagnostics to reveal salient details of the near-nozzle spray structure and to guide improvements to existing primary atomization modeling approaches.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

The FastEddy® Resident-GPU Accelerated Large-Eddy Simulation Framework: Moist Dynamics Extension, Validation and Sensitivities of Modeling Non-Precipitating Shallow Cumulus Clouds

Herein we describe the moist dynamics formulation implemented within the graphics processing unit-resident large-eddy simulation FastEddy® model, which includes a simple saturation adjustment scheme for condensation and evaporation processes. Two LES model intercomparison exercises for non-precipitating shallow cumulus clouds are simulated in order to validate this model extension, including a static forcing and a time-dependent forcing case. Overall, we find our dynamical, thermodynamical and microphysical quantities, along with turbulence variability and fluxes, to be commensurate with the corresponding model intercomparison results. In addition, sensitivities to specific model settings are investigated. Among these settings, it is shown that boundary layer and cloud layer structure and characteristics are sensitive to use of higher-order advection schemes impacting the vertical distribution of cloud content and associated turbulence statistics. Increasing the timescale of the saturation scheme leads to enhanced liquid water presence and decreases vertical velocity variance within the cloud deck. In some cases, these sensitivities agree with the model-to-model variability reported in the intercomparison exercises, highlighting the important role of specific model implementation choices in the context of shallow cumulus convection simulations. These analyses and findings also provide the basis for future extensions and applications of FastEddy® for modeling moist convection and precipitation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sub-5 Ångstrom Porosity Tuning in Calixarene-Derived Porous Liquids via Supramolecular Complexation Construction

Sub-Ångstrom-level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub-5 Ångstrom scale is created via sophisticated porosity tuning in calixarene-derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy-, X-ray-, and neutron-scattering techniques. The presence of permanent porosity in calixarene-derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub-5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. Further, the approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

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

Lithium Solvation and Mobility in Ionic Liquid Electrolytes with Asymmetric Sulfonyl-Cyano Anion

The solvation structure and transport properties of Li + in ionic liquid (IL) electrolytes based on n-methyl-n-butylpyrrolidinium cyano(trifluoromethanesulfonyl)imide [PYR 14 ][CTFSI] and [Li][CTFSI] (0 ≤ x Li ≤ 0.7) were studied by Raman and Nuclear Magnetic Resonance (NMR) diffusometry, and molecular dynamics (MD) simulations. At x Li < 0.3, Li + coordination is dominated by the cyano group. As x Li is increased, free cyano-sites become limited, resulting in increased coordination via the sulfonyl group. Here, the 1:1 mixture of the symmetric anions bis(trifluoromethanesulfonyl)imide ([TFSI]) and dicyanamide ([DCA]) results in similar physical properties as the IL with [CTFSI]. However, anion asymmetry is shown to increase Li-salt solubility and promote Li+ transference. The lifetimes of Li + -cyano coordination for [CTFSI] are calculated to be shorter than those for [DCA], indicating that the competition from the sulfonyl group weakens its solvation with Li + . This resulted in higher Li + transference for the electrolyte with [CTFSI]. In relation to the utility of these electrolytes in energy storage, the Li–LiFePO 4 half cells assembled with IL electrolyte (x Li = 0.3, 0.5, and 0.7) demonstrated a nominal capacity of 140 mAh/g at 0.1C rate and 90 °C where the cell with x Li = 0.7 IL electrolyte demonstrated 61% capacity retention after 100 cycles and superior rate capability owing to increased electrochemical stability.

25 ENERGY STORAGE↗