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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2021

Bulk Glassy GeTe 2 : A Missing Member of the Tetrahedral GeX 2 Family and a Precursor for the Next Generation of Phase-Change Materials

Vitreous germanium disulfide GeS 2 and diselenide GeSe 2 belong to canonical chalcogenide glasses extensively studied over the past half century. Their high-temperature orthorhombic polymorphs are congruently melting compounds, and the tetrahedral crystal and glass structure is largely preserved in the melt. In contrast, the ditelluride counterpart is absent in the Ge-Te phase diagram, which shows only a single compound, monotelluride GeTe. Phase-change materials based on GeTe have become a technologically important class of solids, and their structure and properties are also widely studied. Surprisingly, very scarce information is available for alloys having GeTe 2 stoichiometry. Using a fast quenching procedure in silica capillaries, high-energy X-ray diffraction, and Raman spectroscopy supported by first-principles simulations, we show that bulk glassy GeTe 2 differs substantially from the lighter GeX 2 members, revealing 46% of trigonal germanium, 31% of three-fold coordinated tellurium, and only 20% of edge-sharing tetrahedra or pyramids. The fraction of homopolar Ge-Ge bonds is low; however, the population of dominant Te-Te dimers and Te n oligomers, n <= 10, appears to be significant. The complex structural and chemical topology of g-GeTe 2 is directly related to the thermodynamic metastability of germanium ditelluride, schematically represented by the following reaction: GeTe 2 $\rightleftarrows$ GeTe + Te. Disproportionation is complete above liquidus in the temperature range of semiconductor-metal transition, and the dense metallic GeTe 2 liquid, mostly consisting of five-fold coordinated Ge species, exhibits high fluidity, strong fragility (m = 99 ± 5), and presumably a fast structural transformation rate combined with low atomic mobility in the vicinity of the glass transition temperature, favorable for reliable long-term data retention in nonvolatile memories. The observed and predicted characteristic features make GeTe 2 a promising precursor for the next generation of phase-change materials, especially coupled with additional metal doping, depolymerizing the tetrahedral interconnected glass network and accelerating (sub)nanosecond crystallization.

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

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2020

Review of Emerging Potassium–Sulfur Batteries

Abstract This is the first review on potassium–sulfur (K–S) batteries (KSBs), which are emerging metal battery (MB) systems. Since KSBs are quite new, there are fundamental questions regarding the electrochemistry of S‐based cathode and of K metal anode, as well as the holistic aspects of full‐cell performance. The manuscript begins with a critical discussion regarding the potassium–sulfur electrochemistry and on how it differs from the much better‐known lithium–sulfur. Cathodes are discussed next, focusing on the role of sulfur structure, carbon host chemistry and porosity, and electrolytes in establishing the reversible potassium sulfide K 2 S n phase sequence, the parasitic polysulfide shuttle, pulverization‐driven capacity fade, etc. Following is a discussion of solid‐state electrolytes (SSEs), including of hybrid solid–liquid systems that show much promise. Potassium metal anodes are then critically reviewed, emphasizing electrolyte reactions to form stable versus unstable solid electrolyte interphase (SEI), covering the current understanding of potassium dendrites, and highlighting the deep‐eutectic K–Na alloying approaches for room temperature liquid anodes. The manuscript concludes with K–S batteries, focusing on cell architectures and providing quantitative performance comparisons as master plots. Unanswered scientific/technological questions are identified, emerging research opportunities are discussed, and potential experimental and simulation‐based studies that can unravel these unknowns are proposed.

Ding, Jia↗
DOE OSTI2022

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

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

Critical role of slags in pitting corrosion of additively manufactured stainless steel in simulated seawater

Abstract Pitting corrosion in seawater is one of the most difficult forms of corrosion to identify and control. A workhorse material for marine applications, 316L stainless steel (316L SS) is known to balance resistance to pitting with good mechanical properties. The advent of additive manufacturing (AM), particularly laser powder bed fusion (LPBF), has prompted numerous microstructural and mechanical investigations of LPBF 316L SS; however, the origins of pitting corrosion on as-built surfaces is unknown, despite their utmost importance for certification of LPBF 316L SS prior to fielding. Here, we show that Mn-rich silicate slags are responsible for pitting of the as-built LPBF material in sodium chloride due to their introduction of deleterious defects such as cracks or surface oxide heterogeneities. In addition, we explain how slags are formed in the liquid metal and deposited at the as-built surfaces using high-fidelity melt pool simulations. Our work uncovers how LPBF changes surface oxides due to rapid solidification and high-temperature oxidation, leading to fundamentally different pitting corrosion mechanisms.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Protein Crystallization: Specific Phenomena and General Insights on Crystallization Kinetics

Experimental and simulation studies of the nucleation and growth kinetics of proteins have revealed phenomena that are specific for macromolecular crystallization, and others that provide a more detailed understanding of solution crystallization in general. The more specific phenomena, which include metastable liquid-liquid phase separations and gelation prior to solid nucleation, are due to the small ratio of the intermolecular interaction-range to the size of molecules involved. The apparently more generally applicable mechanisms include the cascade-like formation of macrosteps, as an intrinsic morphological instability that roots in the coupled bulk transport and nonlinear interface kinetics in systems with mixed growth rate control. Analyses of this nonlinear response provide (a) criteria for the choice of bulk transport conditions to minimize structural defect formation, and (b) indications that the "slow" protein crystallization kinetics stems from the mutual retardation of growth steps.

Rosenberger, F.↗
DOE OSTI2023

A numerical study of an atomizing jet in a resonant acoustic field

By applying state-of-the-art, high-fidelity methods for compressible multiphase flows, this work features a parametric study of a turbulent liquid jet atomizing under the influence of standing acoustic waves with five different amplitude values. We perform these simulations in a semi-periodic domain to isolate interactions between acoustic and atomization processes and to emulate a pure liquid fuel spray placed at the pressure node of a resonant acoustic field. Characterization of the flow includes measurements of the mean dimensions, droplet number, and surface area of the liquid distribution, alongside qualitative depictions of the interface evolution, breakup events, and hydrodynamic fields. As acoustic forcing is provided to the flow, the added energy provokes atomization when shear would otherwise be insufficient for breakup. The magnitude of the forcing determines the rate of droplet production, which accelerates as the evolving liquid surface modifies the induced motion of the gas phase. At the nodal plane, the acoustic field produces an organizing effect on the liquid structures.

42 ENGINEERING↗
DOE OSTI2018

Structure, Polarization, and Sum Frequency Generation Spectrum of Interfacial Water on Anatase TiO 2

The photocatalytic activity of TiO 2 for water splitting has been known for decades, yet the adsorption structure and hydrogen bonding of water at the interface with TiO 2 have remained controversial. We investigate the prototypical aqueous interface with anatase TiO 2 (101) using ab initio molecular dynamics (AIMD) with the strongly constrained and appropriately normed (SCAN) density functional, recently shown to provide an excellent description of the properties of bulk liquid water. We find that water forms a stable bilayer of intact molecules with ice-like dynamics and enhanced dipole moment and polarizability on the anatase surface. Here, the orientational order and H-bond environment of interfacial water are reflected in the computed sum frequency generation (SFG) spectrum, which agrees well with recent measurements in the OH stretching frequency range (3000–3600 cm –1 ). Additional AIMD simulations for a model interface with 66% of dissociated water in the contact layer show that surface hydroxyls disrupt the order in the bilayer and lead to a much faster orientational dynamics of interfacial water. Nonetheless, the computed SFG spectrum for the hydroxylated surface also agrees with experiment, suggesting that SFG measurements in a wider frequency range would be necessary to unambiguously identify the character of interfacial water on anatase.

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

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

ZIF-8 Vibrational Spectra: Peak Assignments and Defect Signals

In this study, Zeolitic Imidazolate Framework-8 (ZIF-8) is a promising material for gas separation applications. It also serves as a prototype for numerous ZIFs, including amorphous ones, with a broader range of possible applications, including sensors, catalysis, and lithography. It consists of zinc coordinated with 2-methylimidazolate (2mIm) and has been synthesized with methods ranging from liquid-phase to solvent-free synthesis, which aim to control its crystal size and shape, film thickness and microstructure, and incorporation into nanocomposites. Depending on the synthesis method and post-synthesis treatments, ZIF-8 materials may deviate from the nominal defect-free ZIF-8 crystal structure due to defects like missing 2mIm, missing zinc, and physically adsorbed 2mIm trapped in the ZIF-8 pores, which may alter its performance and stability. Infra-red (IR) spectroscopy has been used to assess the presence of defects in ZIF-8 and related materials. However, conflicting interpretations by various authors persist in the literature. Here, we systematically investigate ZIF-8 vibrational spectra by combining experimental IR spectroscopy and first-principles molecular dynamics simulations, focusing on assigning peaks and elucidating the spectroscopic signals of putative defects present in ZIF-8 material. We attempt to resolve conflicting assignments from the literature and to provide a comprehensive understanding of the vibrational spectra of ZIF-8 and its defect-induced variations, aiming towards more precise quality control and design of ZIF-8-based materials for emerging applications.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
DOE OSTI2022

Multiplexed Anodic Stripping Voltammetry Detection of Heavy Metals in Water Using Nanocomposites Modified Screen-Printed Electrodes Integrated With a 3D-Printed Flow Cell

In this study, we present multiplexed anodic stripping voltammetry (ASV) detection of heavy metal ions (HMIs)—As(III), Cd(II), and Pb(II)—using a homemade electrochemical cell consisting of dual working, reference and counter screen-printed electrodes (SPE) on polyimide substrate integrated with a 3D-printed flow cell. Working and counter electrodes were fabricated by the screen-printing of graphite paste while the Ag/AgCl paste was screen-printed as a reference electrode (Ag/AgCl quasi-reference electrode). The working electrodes were modified with (BiO) 2 CO 3 -reduced graphene oxide (rGO)-Nafion [(BiO) 2 CO 3 -rGO-Nafion] and Fe 3 O 4 magnetic nanoparticles (Fe 3 O 4 MNPs) decorated Au nanoparticles (AuNPs)-ionic liquid (IL) (Fe 3 O 4 -Au-IL) nanocomposites separately to enhance HMIs sensing. Electrochemical detection was achieved using square wave ASV technique. The desired structure of the flow electrochemical cell was optimized by the computational fluid dynamic (CFD). Different experimental parameters for stripping analysis of HMIs were optimized including deposition time, deposition potential and flow rate. The linear range of calibration curves with the sensing nanocomposites modified SPE for the three metal ions was from 0–50 μg/L. The limits of detection (S/N = 3) were estimated to be 2.4 μg/L for As(III), 1.2 μg/L for Pb(II) and 0.8 μg/L for Cd(II). Furthermore, the homemade flow anodic stripping sensor platform was used to detect HMIs in simulated river water with a 95–101% recovery, indicating high selectivity and accuracy and great potential for applicability even in complex matrices.

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

Creating the Thermal Environment for Safely Testing the James Webb Space Telescope at the Johnson Space Center's Chamber A

Chamber A is the largest thermal vacuum chamber at the Johnson Space Center and is one of the largest space environment chambers in the world. The chamber is 19.8 m (65 ft) in diameter and 36.6 m (120 ft) tall and is equipped with cryogenic liquid nitrogen panels (shrouds) and gaseous helium shrouds to create a simulated space environment. The chamber was originally built to support testing of the Apollo Service and Command Module for lunar missions, but underwent major modifications to be able to test the James Webb Space Telescope in a simulated deep space environment. To date seven tests have been performed in preparation of testing the flight optics for the James Webb Space Telescope (JWST). Each test has had a uniquie thermal profile and set of thermal requirements for cooling down and warming up, controlling contamination, and releasing condensed air. These range from temperatures from 335K to 15K, with tight uniformity and controllability for maintining thermal stability and pressure control. One unique requirement for two test was structurally proof loading hardware by creating thermal gradients at specific temperatures. This paper will discuss the thermal requirements and goals of the tests, the original requirements of the chamber thermal systems for planned operation, and how the new requirements were met by the team using the hardware, system flexiblilty, and engineering creativity. It will also discuss the mistakes and successes to meet the unique goals, especially when meeting the thermal proof load.

Homan, Jonathan L.↗
DOE OSTI2017

Energy landscape-driven non-equilibrium evolution of inherent structure in disordered material

Complex states in glasses can be neatly expressed by the potential energy landscape (PEL). But, because PEL is highly multi-dimensional it is difficult to describe how the system moves around in PEL. We demonstrate that it is possible to predict the evolution of macroscopic state in a metallic glass, such as ageing and rejuvenation, through a set of simple equations describing excitations in the PEL. The key to this simplification is the realization that the step of activation from the initial state to the saddle point in PEL and the following step of relaxation to the final state are essentially decoupled. Furthermore, the model shows that the interplay between activation and relaxation in PEL is the key driving force that simultaneously explains both the equilibrium of supercooled liquid and the thermal hysteresis observed in experiments. It further predicts anomalous peaks in truncated thermal scanning, validated by independent molecular dynamics simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Revealing the Size-Dependent d–d Excitations of Cobalt Nanoparticles Using Soft X-ray Spectroscopy

Cobalt-based catalysts are widely used to produce liquid fuels through the Fischer-Tropsch (FT) reaction. However, the cobalt nanocatalysts can exhibit intriguing size-dependent activity whose origin remains heavily debated. To shed light on this issue, the electronic structures of cobalt nanoparticles with size ranging from 4 to 10 nm are studied using soft X-ray absorption (XAS) and resonant inelastic X-ray scattering (RIXS) spectroscopies. The RIXS measurements reveal the significant size-dependent d-d excitations, from which we determine that the crystal-field splitting energy 10Dq changes from 0.6 to 0.9 eV when the particle size is reduced from 10 to 4 nm. The finding that larger Co nanoparticles have smaller 10Dq value is further confirmed by the Co L-edge RIXS simulations with atomic multiplet code. Furthermore, our RIXS results demonstrate a stronger Co-O bond in smaller Co nanoparticles, which brings in further insight into their size-dependent catalytic performance.

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

Fracture Mechanics Analysis of LH2 Feed Line Flow Liners

Inspections of the Space Shuttle Main Engine revealed fatigue cracks growing from slots in the flow liner of the liquid hydrogen (LH2) feed lines. During flight, the flow liners experience complex loading induced by flow of LH2 and the resonance characteristics of the structure. The flow liners are made of Inconel 718 and had previously not been considered a fracture critical component. However, fatigue failure of a flow liner could have catastrophic effect on the Shuttle engines. A fracture mechanics study was performed to determine if a damage tolerance approach to life management was possible and to determine the sensitivity to the load spectra, material properties, and crack size. The load spectra were derived separately from ground tests and material properties were obtained from coupon tests. The stress-intensity factors for the fatigue cracks were determined from a shell-dynamics approach that simulated the dominant resonant frequencies. Life predictions were obtained using the NASGRO life prediction code. The results indicated that adequate life could not be demonstrated for initial crack lengths of the size that could be detected by traditional NDE techniques.

James, Mark A.↗