Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffused interface model”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗
DOE OSTI2025

TEM characterization of two variants of fuel cladding chemical interaction in a HT-9 Clad U-10Zr Fuel. Variant 1: FCCI with a Zr Rind

Here, this study investigated the fuel cladding chemical interaction (FCCI), a key factor that limits operational temperature and burnup, in an HT-9 clad U-10Zr nuclear fuel sample irradiated to a high burnup of 13.1 at.% at a time-averaged peak inner cladding temperature (PICT) of 530 °C. Previous results showed this fuel sample exhibited two distinct levels of FCCI at d. This paper analyzed the FCCI at an azimuthal position showing an interdiffusion layer of <10 µm using transmission electron microscopy to examine chemical and crystallographic nature of phases at the fuel-cladding interface at the nanoscale level. A ZrC layer and a Zr 3 Si phase were identified at the interface; these, along with the relatively low local temperature, potentially contributed to limit interdiffusion, behaving as inhibitors for deleterious interactions. Lanthanides (Ln) partially consumed the ZrC layer and interacted with Fe, forming a Zr-Ln compound and a (Zr,Ce)Fe 2+x phase while also infiltrating up to 4 µm into the cladding. Neither U nor Zr were observed in the cladding, whereas Fe diffused up to 3–5 µm in the fuel. Fe infiltration formed a ternary U-Zr-Fe ε-phase and likely promoted the precipitation of a Cr-rich α’ phase on the cladding interface. Additionally, a Cr-rich χ-phase, likely formed by the dissociation of pre-existing M 23 C 6 carbide precipitates, was identified about 2–5 µm from the fuel-cladding interface. Irradiation-induced nano-voids were also observed in the HT-9 bulk. These findings provide critical insights into FCCI mechanisms at representative irradiation conditions, essential for developing models simulating in-pile metallic fuel behaviors for next-generation reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Nanoscale Turing patterns in a bismuth monolayer

Turing’s reaction–diffusion theory of morphogenesis has been very successful for understanding macroscopic patterns within complex objects ranging from biological systems to sand dunes. However, Turing patterns on microscopic length scales are extremely rare. Here we show that a strained atomic bismuth monolayer assembled on the surface of NbSe 2 —and subject to interatomic interactions and kinetics—displays Turing patterns. Our reaction–diffusion model produces stripe patterns with a period of five atoms (approximately 2 nm) and domain walls with Y-shaped junctions that bear a striking resemblance to what has been experimentally observed. Our work establishes that Turing patterns can occur at the atomic scale in a hard condensed-matter setting.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling the Impact of Grain Size on Corrosion Behavior of Ni-Based Alloys in Molten Chloride Salt via Cellular Automata

Molten chloride salts hold significant promise as both thermal transfer and storage media for next-generation concentrated solar power (CSP) systems. However, molten chlorides pose a considerable corrosion risk to structural materials, particularly Ni-based alloys. One approach to enhancing corrosion resistance is through the optimization of grain structure; however, it remains uncertain whether increasing or decreasing grain size enhances corrosion resistance. A cellular automata (CA) program was developed to evaluate the interplay between grain size and corrosion in Ni-based alloy. Our CA program tracks alloy composition, surface roughness, and thickness loss via a graphical user interface, displaying corrosion and diffusion status, and multiple user input cards for tuning the simulation. CA simulations of Inconel 625 indicate enhanced corrosion resistance with increased grain size, with passivating oxides offering limited protection. Additionally, the temporal evolution of alloy surface roughness demonstrates notable fluctuations, with abrupt increases attributed to corrosion along vertical grain boundaries and sudden decreases to grain detachment from the protective film.

Materials Science↗
DOE OSTI2023

Structure–Transport Relationships of Water–Organic Solvent Co-transport in Carbon Molecular Sieve (CMS) Membranes

We explore the effects of the carbon molecular sieve (CMS) microstructure on the separation performance and transport mechanism of water–organic mixtures. Specifically, we utilize PIM-1 dense films and integrally skinned asymmetric hollow fiber membranes as polymer precursors for the CMS materials. The PIM-1 membranes were pyrolyzed under several different pyrolysis atmospheres (argon, carbon dioxide, and diluted hydrogen gas) and at multiple pyrolysis temperatures. Detailed gas physisorption measurements reveal that membranes pyrolyzed under 4% H 2 and CO 2 had broadened ultramicropore distributions (pore diameter <7 Å) compared to Ar pyrolysis, and pyrolysis under CO 2 increased ultramicropore volume and broadened micropore distributions at increased pyrolysis temperatures. Gravimetric water and p-xylene sorption and diffusion measurements reveal that the PIM-1-derived CMS materials are more hydrophilic than other CMS materials that have been previously studied, which leads to sorption-diffusion estimations showing water-selective permeation. Water permeation in the vapor phase, pervaporation, and liquid-phase hydraulic permeation reveal that the isobaric permeation modes (vapor permeation and pervaporation) are reasonably well predicted by the sorption-diffusion model, whereas the hydraulic permeation mode is significantly underpredicted (>250×). Conversely, the permeation of p-xylene is well predicted by the sorption-diffusion model in all cases. The collection of pore size analysis, vapor sorption and diffusion, and permeation in different modalities creates a picture of a combined transport mechanism in which water–under high transmembrane pressures–permeates via a Poiseuille-style mechanism, whereas p-xylene solutes in the mixture permeate via sorption-diffusion.

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

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. While simulations have been able to capture key aspects of these processes to reproduce film morphology, these models typically ignore any subsurface diffusion that occurs. The results of these diffusion and intermixing processes are often observed experimentally, but cannot be explained by invoking bulk lattice diffusion at the deposition temperature. In this work, we present quantitative experimental data that reveals anion and cation intermixing over long length scales during the deposition of epitaxial Fe 2 O 3 and Cr 2 O 3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18 O and/or 57 Fe, and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. These reveal that adatoms on the film surface act to “pull up” subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. A simple model incorporating these events qualitatively agrees with the experimental intermixing results. In addition to film deposition, these intermixing mechanisms may to be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Spontaneous Lithiation of Binary Oxides during Epitaxial Growth on LiCoO 2

Epitaxial growth is a powerful tool for synthesizing heterostructures and integrating multiple functionalities. Interfacial mixing can readily occur at temperatures required for complete film crystallization and can significantly modify the properties of layered structures, particularly for those containing energy storage materials with smaller cations. Here we show a two-step sequence involving the growth of an epitaxial LiCoO 2 cathode layer followed by the deposition of a binary transition metal oxide (WO 3 , TiO 2 , and others) in which controlled lithiation of the binary oxide occurs. Orientation-controlled epitaxial synthesis of the model solid-state-electrolyte Li 2 WO 4 and model anode material Li 4 Ti 5 O 12 occurs as WO 3 and TiO 2 nucleate their respective host lattices and attract Li ions from the underlying cathode. The cathode layer readily provides a tunable amount of Li ions for the controlled lithiation of the subsequent layer. We demonstrate that this approach can be used for energy materials discovery and exploring different combinations of epitaxial interfaces that can serve as well-defined model systems for mechanistic studies of energy storage and conversion processes.

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

Efficiency and Forward Voltage of Blue and Green Lateral LEDs with V-shaped Defects and Random Alloy Fluctuation in Quantum Wells

For nitride-based blue and green light-emitting diodes (LEDs), the forward voltage V for is larger than expected, especially for green LEDs. This is mainly due to the large barriers to vertical carrier transport caused by the total polarization discontinuity at multiple quantum well and quantum barrier interfaces. The natural random alloy fluctuation in quantum wells has proven to be an important factor reducing V for . However, this does not suffice in the case of green LEDs because of their larger polarization-induced barrier. V-shaped defects (V-defects) have been proposed as another key factor in reducing V for to allow lateral injection into multiple quantum wells, thus bypassing the multiple energy barriers incurred by vertical transport. In this paper, to model carrier transport in the whole LED, we consider both random-alloy and V-defect effects. A fully two-dimensional drift-diffusion charge-control solver is used to model both effects. The results indicate that the turn-on voltages for blue and green LEDs are both affected by random alloy fluctuations and the V-defect density. For green LEDs, V for decreases more due to V-defects, where the smaller polarization barrier at the V-defect sidewall is the major path for lateral carrier injection. Then, we discuss how the V-defect density and size affects the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Black and gray box learning of amplitude equations: Application to phase field systems

Here, we present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

42 ENGINEERING↗
DOE OSTI2020

Converting 2D Nanofiber Membranes to 3D Hierarchical Assemblies with Structural and Compositional Gradients Regulates Cell Behavior

Abstract New methods are described for converting 2D electrospun nanofiber membranes to 3D hierarchical assemblies with structural and compositional gradients. Pore‐size gradients are generated by tuning the expansion of 2D membranes in different layers with incorporation of various amounts of a surfactant during the gas‐foaming process. The gradient in fiber organizations is formed by expanding 2D nanofiber membranes composed of multiple regions collected by varying rotating speeds of mandrel. A compositional gradient on 3D assemblies consisting of radially aligned nanofibers is prepared by dripping, diffusion, and crosslinking. Bone mesenchymal stem cells (BMSCs) on the 3D nanofiber assemblies with smaller pore size show significantly higher expression of hypoxia‐related markers and enhanced chondrogenic differentiation compared to BMSCs cultured on the assemblies with larger pore size. The basic fibroblast growth factor gradient can accelerate fibroblast migration from the surrounding area to the center in an in vitro wound healing model. Taken together, 3D nanofiber assemblies with gradients in pore sizes, fiber organizations, and contents of signaling molecules can be used to engineer tissue constructs for tissue repair and build biomimetic disease models for studying disease biology and screening drugs, in particular, for interface tissue engineering and modeling.

Chen, Shixuan↗
DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

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

Modeling the aqueous interface of amorphous TiO 2 using deep potential molecular dynamics

Amorphous titanium dioxide (a-TiO 2 ) is widely used as a coating material in applications such as electrochemistry and self-cleaning surfaces where its interface with water has a central role. However, little is known about the structures of the a-TiO 2 surface and aqueous interface, particularly at the microscopic level. In this work, we construct a model of the a-TiO 2 surface via a cut-melt-and-quench procedure based on molecular dynamics simulations with deep neural network potentials (DPs) trained on density functional theory data. After interfacing the a-TiO 2 surface with water, we investigate the structure and dynamics of the resulting system using a combination of DP-based molecular dynamics (DPMD) and ab initio molecular dynamics (AIMD) simulations. Both AIMD and DPMD simulations reveal that the distribution of water on the a-TiO 2 surface lacks distinct layers normally found at the aqueous interface of crystalline TiO 2 , leading to an ~10 times faster diffusion of water at the interface. Bridging hydroxyls (Ti 2 –O b H) resulting from water dissociation decay several times more slowly than terminal hydroxyls (Ti–O w H) due to fast Ti–O w H 2 → Ti–O w H proton exchange events. These results provide a basis for a detailed understanding of the properties of a-TiO 2 in electrochemical environments. Moreover, the procedure of generating the a-TiO 2 -interface employed here is generally applicable to studying the aqueous interfaces of amorphous metal oxides.

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

Topographic effects on reflected acoustic waves from the OSIRIS-REx reentry observed from stratospheric balloons

During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.

Lees, Jonathan M. [University of North Carolina, C↗
DOE OSTI2021

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
DOE OSTI2023

Transport and Energetics of Carbon Dioxide in Ionic Liquids at Aqueous Interfaces

A major hurdle in utilizing carbon dioxide (CO 2 ) lies in separating it from industrial flue gas mixtures and finding suitable storage methods that enable its application in various industries. To address this issue, we utilized a combination of molecular dynamics simulations and experiments to investigate the behavior of CO 2 in common room-temperature ionic liquids (RTIL) when in contact with aqueous interfaces. Here, our investigation of RTILs, [EMIM][TFSI] and [OMIM][TFSI], and their interaction with a pure water layer mimics the environment of a previously developed ultrathin enzymatic liquid membrane for CO 2 separation. We analyzed diffusion constants and viscosity, which reveals that CO 2 molecules exhibit faster mobility within the selected ILs compared to what would be predicted solely based on the viscosity of the liquids using the standard Einstein–Stokes relation. Moreover, we calculated the free energy of translocation for various species across the aqueous–IL interface, including CO 2 and HCO 3 – . Free energy profiles demonstrate that CO 2 exhibits a more favorable partitioning behavior in the RTILs compared to that in pure water, while a significant barrier hinders the movement of HCO 3 – from the aqueous layer. Experimental measurement of the CO 2 transport in the RTILs corroborates the model. These findings strongly suggest that hydrophobic RTILs could serve as a promising option for selectively transporting CO 2 from aqueous media and concentrating it as a preliminary step toward storage.

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

Exchange of Ions across the TiN/TaO x Interface during Electroformation of TaO x -Based Resistive Switching Devices

The valence change model describes the resistive switching in metal oxide-based devices as due to electroreduction of the oxide and subsequent electromigration of oxygen vacancies. Here, we present cross-sectional X-ray energy-dispersive spectroscopy elemental maps of Ta, O, N, and Ti in electroformed TiN/TaO 2.0 /TiN structures. O, N, and Ti were exchanged between the anode and the functional oxide in devices formed at high power (~1 mW), but the exchange was below the detection limit at low power (<0.5 mW). All structures exhibit a similar Ta-enriched and O-depleted filament formed by the elemental segregation in the functional oxide by the temperature gradient. The elemental interchange is interpreted as due to Fick’s diffusion caused by high temperatures in the gap of the filament and is not an essential part of electroformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Combinations of Single Chain Variable Fragments From HIV Broadly Neutralizing Antibodies Demonstrate High Potency and Breadth

Broadly neutralizing antibodies (bNAbs) are currently being assessed in clinical trials for their ability to prevent HIV infection. Single chain variable fragments (scFv) of bNAbs have advantages over full antibodies as their smaller size permits improved diffusion into mucosal tissues and facilitates vector-driven gene expression. We have previously shown that scFv of bNAbs individually retain significant breadth and potency. Here we tested combinations of five scFv derived from bNAbs CAP256-VRC26.25 (V2-apex), PGT121 (N332-supersite), 3BNC117 (CD4bs), 8ANC195 (gp120-gp41 interface) and 10E8v4 (MPER). Either two or three scFv were combined in equimolar amounts and tested in the TZM-bl neutralization assay against a multiclade panel of 17 viruses. Experimental IC 50 and IC 80 data were compared to predicted neutralization titers based on single scFv titers using the Loewe additive and the Bliss-Hill model. Like full-sized antibodies, combinations of scFv showed significantly improved potency and breadth compared to single scFv. Combinations of two or three scFv generally followed an independent action model for breadth and potency with no significant synergy or antagonism observed overall although some exceptions were noted. The Loewe model underestimated potency for some dual and triple combinations while the Bliss-Hill model was better at predicting IC 80 titers of triple combinations. Given this, we used the Bliss-Hill model to predict the coverage of scFv against a 45-virus panel at concentrations that correlated with protection in the AMP trials. Using IC 80 titers and concentrations of 1μg/mL, there was 93% coverage for one dual scFv combination (3BNC117+10E8v4), and 96% coverage for two of the triple combinations (CAP256.25+3BNC117+10E8v4 and PGT121+3BNC117+10E8v4). Combinations of scFv, therefore, show significantly improved breadth and potency over individual scFv and given their size advantage, have potential for use in passive immunization.

60 APPLIED LIFE SCIENCES↗