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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

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

Numerical determination of the transmissibility characteristics of a squeeze film damped forced vibration system

Numerical solutions of the governing equations of motion of a liquid squeeze film damped forced vibration system were carried out to examine the feasibility of using a liquid squeeze film to cushion and protect large structures, such as buildings, located in areas of high seismic activity. The mathematical model used was that for a single degree of freedom squeeze film damped spring mass system. The input disturbance was simulated by curve fitting actual seismic data with an eleventh order Lagranging polynomial technique. Only the normal component of the seismic input was considered. The nonlinear, nonhomogeneous governing differential equation of motion was solved numerically to determine the transmissibility over a wide range of physical parameters using a fourth-order Runge-Kutta technique. It is determined that a liquid squeeze film used as a damping agent in a spring-mass system can significantly reduce the response amplitude for a seismic input disturbance.

Sutton, M. A.↗
DOE OSTI2025

Structure of Novel Phosphonium-Based Ionic Liquids with S and O Substitutions from Experiments and a Mixed Quantum-Classical Approach

This article presents experimental characterization information and synchrotron X-ray scattering measurements on a set of novel O- and S-substituted phosphonium-based ionic liquids (ILs) all coupled with the bis(fluorosulfonyl)imide (FSI-) anion. The ILs include the ethoxyethyltriethylphosphonium (P 222(2O2) + ) and triethyl[2-(ethylthio)ethyl]phosphonium (P 222(2S2) + ) cations, and we contrast results on these with those for unsubstituted triethylpentylphosphonium (P 2225 + ). The article also introduces a physics-based protocol that combines classical force field studies on larger simulation boxes with classical and first-principles studies on smaller boxes. The method produces significantly improved S(q) functions in the regime which in prior publications we have associated with inter- and intraionic adjacency correlations. By understanding which shorter-range structural changes improve S(q) in the q-regime of interest, we are also able to pinpoint specific deficiencies in the classical force field model. The approach we take should be quite general and could help study other complex liquids on different length scales.

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

Corrosion on Mars: An Investigation of Corrosion Mechanisms Under Relevant Simulated Martian Environments

This one-year project was selected by NASA's Science Innovation Fund in FY17 to address Corrosion on Mars which is a problem that has not been addressed before. Corrosion resistance is one of the most important properties in selecting materials for landed spacecraft and structures that will support surface operations for the human exploration of Mars. Currently, the selection of materials is done by assuming that the corrosion behavior of a material on Mars will be the same as that on Earth. This is understandable given that there is no data regarding the corrosion resistance of materials in the Mars environment. However, given that corrosion is defined as the degradation of a metal that results from its chemical interaction with the environment, it cannot be assumed that corrosion is going to be the same in both environments since they are significantly different. The goal of this research is to develop a systematic approach to understand corrosion of spacecraft materials on Mars by conducting a literature search of available data, relevant to corrosion in the Mars environment, and by performing preliminary laboratory experiments under relevant simulated Martian conditions. This project was motivated by the newly found evidence for the presence of transient liquid brines on Mars that coincided with the suggestion, by a team of researchers, that some of the structural degradation observed on Curiosity's wheels may be caused by corrosive interactions with the brines, while the most significant damage was attributed to rock scratching. An extensive literature search on data relevant to Mars corrosion confirmed the need for further investigation of the interaction between materials used for spacecraft and structures designed to support long-term surface operations on Mars. Simple preliminary experiments, designed to look at the interaction between an aerospace aluminum alloy (AA7075-T73) and the gases present in the Mars atmosphere, at 20degC and a pressure of 700 Pa, showed that there is an interaction between the small amount of oxygen present in the Mars gas and the alloy when there is a scratch that removes the protective aluminum oxide film. Further studies are needed to consider many other important components of the Mars environment that can affect this interaction such as: the effect of oxidants, the effect of radiation on their oxidizing properties and the possible catalytic effects of the clays present in the Martian regolith. The results of this one-year project provide strong justification for further investigation of the corrosion mechanism of materials relevant to long-term surface operations in support of future human exploration missions on Mars.

Perchlorates↗
DOE OSTI2025

Thermal mechanical assessment of a SiC-SiC-composite clad fuel pin concept in a light water reactor environment

Accident Tolerant Fuels (ATFs) are designed to increase coping time following an accident scenario while preserving or improving current steady state reactor operational performance. A potential ATF concept is SiC-SiC composite claddings. Fuel performance simulations were conducted on a SiC-SiC based cladding concept utilizing a multilayered approach for improved performance. This cladding concept referred to as the Duplex concept is a duplex structure composed of a monolithic SiC layer placed on the outside of a SiC-SiC composite. A liquid metal is added to fuel-cladding gap for improved heat dissipation from the fuel. The monolithic SiC layer is used to improve the coolant corrosion characteristics and protect the SiC-SiC composite layer from exposure to the coolant. The fuel performance code BISON was used to conduct fuel performance simulations on the cladding concepts. Comparisons are made with a current prototypic fuel rod design (UO 2 fuel enclosed in Zircaloy-4 cladding). Representative steady-state cases were considered for normal power and two cycle power histories. Additionally, a PCI ramp case was simulated to analyze potential anticipated operational occurrences. Transient response during a Loss of Coolant Accident and a Reactivity Initiated Accident were also simulated. This computational study demonstrated that for normal operating conditions, the SiC concept cladding performed as well as the baseline for the standard power cases evaluated. The ramping evaluations indicate potential fracturing of the SiC-SiC composite of the composite cladding compared to the Zircaloy-4 cladding due to the temperature gradient and the subsequent differential thermal conductivity degradation and swelling across the composite thickness. In conclusion, the rod fails early at low enthalpy for RIA but survives a LOCA with minimal material loss due to high temperature steam corrosion.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Characterization of an Optically Accessible RP/GOx Rotating Detonation Engine

The accurate characterization of the dynamic response of liquid jets to transient detonation wave passage is crucial for optimizing and modeling liquid-fueled rotating detonation combustors. To address this need, we have developed an optically accessible rocket RDE (RP-RDE), which utilizes liquid RP-1 and gaseous oxygen as propellants and operates at elevated pressures. The experimental setup features optical access from the oxidizer injection plenum to the combustor exit, allowing for detailed characterization of injection, mixing, and combustion processes and their impact on the detonation wave structure. To design and develop the liquid injection scheme for the RP-RDE test article, we employed a combined experimental and numerical approach that included cold-flow simulations. Our aim was to optimize the liquid injection process and ensure that the resulting fuel-oxidizer mixture is homogeneous and stable. To this end, high-speed optical imaging will be employed to provide qualitative visualization of the liquid and gaseous phases of the propellant as well as the local heat release during detonation wave propagation through the non-premixed medium. This work will provide critical insights into the behavior of the transient detonation waves in more relevant rocket conditions, which is essential for the optimization and modeling of rotating detonation combustors.

Propulsion↗
NASA NTRS1974

Immiscible materials and alloys

Only immiscible materials which involve a fluid phase have a potential advantage when considered for space processing. Processing of such systems at low gravity can lead to a material which is very homogeneous and contains a finely dispersed mixture of phases. Materials with such structures may exhibit potentially useful properties for such applications as superconductors, dispersion-strengthened materials, superplastic materials, permanent magnets, etc. These conclusions have resulted from review of past efforts in the field of low-g processing of materials containing a liquid phase miscibility gap and from present efforts dealing with such systems. These studies have concentrated on precipitation of liquid droplets in a host liquid during cooling through the miscibility gap. This effort has been to study the agglomeration of the droplets both experimentally and through computer simulation in order to understand the mechansism involved and to deduce the effect of gravity on the droplet distribution and resulting structure of the solidified material. The agglomeration mechanisms treated in detail are those due to diffusional growth and collision processes caused by Stokes on convection current-induced droplet migration.

Markworth, A. J.↗
DOE OSTI2023

Ion Solvation and Hydrogen Bonding in Liquid Electrolytes

The primary goal of the proposed research was to explore the use of far-infrared (FIR) spectroscopy for providing new informal on about the solvation structure and dynamics of small ions such as Li + , Na + , and Mg 2+ in electrolyte solutions commonly used in energy-related applications. Classical and ab initio molecular dynamics simulations were to be carried out together with spectroscopic measurements to aid interpretation of the spectra generated.

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

Finite-size effects in the static structure factor S ( k ) and S ( 0 ) for a two-dimensional Yukawa liquid

Finite-size effects in the static structure factor S(k) are analyzed for an amorphous substance. As the number of particles is reduced, S(0) increases greatly, up to an order of magnitude. Meanwhile, there is a decrease in the height of the first peak S peak. These finite-size effects are modeled accurately by the Binder formula for S(0) and our empirical formula for S peak. Procedures are suggested to correct for finite-size effects in S(k) data and in the hyperuniformity index H ≡ S(0)/S peak. These principles generally apply to S(k) obtained from particle positions in noncrystalline substances. As a result, the amorphous substance we simulate is a two-dimensional liquid, with a soft Yukawa interaction modeling a dusty plasma experiment.

2-dimensional systems↗
DOE OSTI2023

Direct Extraction of Sodium Hydroxide by Calix[4]pyrrole-Based Ion-Pair Receptors

Here, described in this work, are calix[4]pyrrole-based ion-pair receptors, cis/trans-1 and cis/trans-2, designed for the extraction of sodium hydroxide. An X-ray diffraction analysis of a single crystal of the cis-1·NaOH isomer isolated from a mixture of cis/trans-1 revealed a unique dimeric supramolecular structure. An average dimer in toluene-d 8 solution was inferred on the basis of diffusion-ordered spectroscopy (DOSY). Support for the proposed stoichiometry came from density functional theory (DFT) calculations. The structural stability of the dimeric cis-1·NaOH complex in toluene solution was further confirmed by ab initio molecular dynamics (AIMD) simulation with explicit representation of solvent. Under conditions of liquid–liquid extraction (LLE), purified receptors cis- and trans-2 were both found to remove NaOH from a pH 11.01 aqueous source phase into toluene with extraction efficiencies (E%) of 50–60% when used equimolar to NaOH. However, in all cases, precipitation was observed. Complexities associated with precipitation could be avoided by immobilization of the receptors onto a chemically inert poly(styrene) resin by means of solvent impregnation. The use of solvent-impregnated resins (SIRs) eliminated precipitation in solution while retaining the extraction efficiency toward NaOH. This allowed both the pH and salinity of the alkaline source phase to be lowered.

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

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational Analysis of the Effect of Structured Packing Design on Absorption Column Hydrodynamics for Post-Combustion Carbon Capture Applications

Solvent based post-combustion carbon capture technologies have a potential for reducing carbon emissions from fossil-fuel-fired power plants and industrial sources where CO2 emissions are inherently harder to mitigate, such as steel or cement industries. A prominent technology to achieve this is by retrofitting absorption columns to the existing infrastructure. While these systems have been among the less costly alternatives for carbon capture, they still impose a considerable energy penalty to the operation of power plants or industrial facilities. The optimization of CO2 capture rate in solvent-based absorption process is complex as it depends on several factors including CO2 solubility, solvent reaction kinetics and temperature effects on the solubility, reaction rates, surface tension, and thermophysical properties of the solvent and the flue gas. The overall heat and mass transfer also depends on the hydrodynamics, which in turn, is affected by the packing geometry. In the current work, we systematically quantify the effects the design of the structured packing has on the column hydrodynamics, by performing detailed CFD simulations for different geometrical configurations and operating conditions. We then obtain relationships between the key hydrodynamic metrics, such as liquid holdup, interfacial area, wetted area, and pressure drop to the parameters defining the packing geometries and identify new more effective packing designs for the given operating conditions.

Shah, Yash Girish↗
DOE OSTI2018

Phase Field Simulations of Microstructure Evolution in IN718 using a Surrogate Ni–Fe–Nb Alloy during Laser Powder Bed Fusion

The solidification microstructure in IN718 during additive manufacturing was modeled using phase field simulations. The novelty of the research includes the use of a surrogate Ni–Fe–Nb alloy that has the same equilibrium solidification range as IN718 as the model system for phase field simulations, the integration of the model alloy thermodynamics with the phase field simulations, and the use of high-performance computing tools to perform the simulations with a high enough spatial resolution for realistically capturing the dendrite morphology and the level of microsegregation seen under additive manufacturing conditions. Heat transfer and fluid flow models were used to compute the steady state temperature gradient and an average value of the solid-liquid (s-l) interface velocity that were used as input for the phase field simulations. The simulations show that the solidification morphology is sensitive to the spacing between the columnar structures. Spacing narrower than a critical value results in continued growth of a columnar microstructure, while above a critical value the columnar structure evolves into a columnar dendritic structure through the formation of secondary arms. These results are discussed in terms of the existing columnar to dendritic transition (CDT) theories. The measured interdendritic Nb concentration, the primary and secondary arm spacing is in reasonable agreement with experimental measurements performed on the nickel-base superalloy IN718.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Stability Analysis of Intertank Formed Skin/Stringer Compression Panel with Simulated Damage

The External Tank (ET) is a component of the Space Shuttle launch vehicle that contains fuel and oxidizer. During launch, the ET supplies the space shuttle main engines with liquid hydrogen and liquid oxygen. In addition to supplying fuel and oxidizer, it is the backbone structural component of the Space Shuttle. It is comprised of a liquid hydrogen (LH2) tank and a liquid oxygen (LOX) tank, which are separated by an Intertank. The Intertank is a stringer-stiffened cylindrical structure with hat-section stringers that are roll formed from aluminum-lithium alloy Al-2090. Cracks in the Intertank stringers of the STS-133 ET were noticed after a November 5, 2010 launch attempt. The cracks were approximately nine inches long and occurred on the forward end of the Intertank (near the LOX tank), along the fastener line, and were believed to have occurred while loading the ET with the cryogenic propellants. These cracks generated questions about the structural integrity of the Intertank. In order to determine the structural capability of the Intertank with varying degrees of damage, a finite element model (FEM) simulating a 1995 compression panel test was analyzed and correlated to test data. Varying degrees of damage were simulated in the FEM, and non-linear stability analyses were performed. The high degree of similarity between the compression panel and the Intertank provided confidence that the ET Intertank would have similar capabilities.

Harper, David W.↗
DOE OSTI2019

Effects of dopants on the glass forming ability in Al-based metallic alloy

The effect of dopants on the metallic glass forming ability is usually considered based on the analysis of changes in the liquid structure or thermodynamics. What is missing in such considerations is an analysis of how a dopant changes the properties of the crystal phases which can form instead of the glass. In order to illuminate this aspect we performed molecular dynamics simulations to study the effects of Mg and Sm dopants on the crystal nucleation in Al. The simulation data were found to be consistent with the experimental observations that addition of Mg to Al does not lead to vitrification but addition of only 8% Sm does. The significant effect of Sm doping was related to the intolerance of Al to this dopant. This leads to increase in the solid-liquid interfacial free energy, and therefore, to increase in the nucleation barrier and to a dramatic decrease in the nucleation rate. Here, the intolerance mechanism also significantly affects the growth kinetics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

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

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗