Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface diffusion”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multi-physics Peridynamic Modeling of Damage Processes in Protective Coatings

Protective coatings can revolutionize the aerospace industry by providing a way to build a fuel-efficient and low-emission gas turbine engine operating at elevated temperature. Currently, significant challenges persist in the multi-physics modeling of damage processes in protective coatings exposed to the hot corrosive environment of gas turbine engine. Examples of such detrimental processes include an interface damage between different material layers forming protective coating, a formation of the through-the-thickness vertical cracks, a diffusion of oxygen and moisture through those cracks, an oxidation, and a delamination. In this paper, a feasibility of a multi-physics peridynamic model to provide insight to these detrimental processes has been demonstrated. 3-Dimensional peridynamic equations which are valid everywhere, including dynamically evolving discontinuities, have been presented. In conclusion, validation and demonstration have been shown for conduction heat transfer, reactive oxidation, and delamination of protective coating with a composition borrowed from the open literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Parabolic oxidation kinetics of a plutonium alloy at room temperature

Auger electron spectroscopy (AES) was used to characterize the oxide scale formed on an ~0.5 wt. % Ga plutonium metal alloy. Analysis of changes to the characteristic Auger transitions for Pu and O during sputter depth profiles were used in combination with previously measured sputter yields to determine the thickness of the oxide scale and provide insight into its chemistry. From this data, the parabolic kinetics for the Pu-Ga alloy were determined as a function of relative humidity, RH = 14–92%, and oxygen partial pressure, p(O 2 ) = 31–157 Torr, at 25°C. For each relative humidity investigated, the parabolic rate constant was found to be independent of oxygen partial pressure. Furthermore, power-law fits to the relative humidity effect on the rate constant was found to have a power of n = 1.265, indicating the diffusion species which rate-limits the reaction isn’t molecular water, which subsequent dissociates at the oxide-metal interface, or hydroxyl from dissociation at the gas-surface interface.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2026

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

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

Reactive force fields for aqueous and interfacial magnesium carbonate formation

Here we develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral–water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg 2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite–water and brucite–water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg–CO 3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral–water–CO 2 interface in water-poor conditions.

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

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Interfacial Characterization of Positive Bias Voltage Degradation in PV Modules

Degradation from high system voltage is a prevailing failure mode in fielded photovoltaic modules, and the degradation mechanism is inherently dependent on the bias polarity. Here we report the effects of positive bias. Modules under positive bias demonstrated a significant photocurrent loss caused by two routes. First, delamination and discoloration of the silicon nitride layer, leading to optical loss determined by reflectance measurements. Second, chemical discoloration of the cell gridlines and encapsulant (EVA), which is linked to an electrochemical reaction at the silver electrodes. Chemical compositional analysis using X-ray photoemission spectroscopy demonstrated that the discoloration is attributed to Ag 2 S and/or Ag 2 O. Evidence of Ag ion migration from the cell grid into the encapsulant is observed after shallow depth profiling on the EVA surface. However, Ag was not detected at the EVA/glass interface, inferring limited Ag ion transport through the EVA. The source of sulfur is believed to be ambient air, which diffused into the module through the breathable backsheet.

characterization↗
DOE OSTI2021

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Modeling a ring magnet in ALEGRA

We show here that Sandia's ALEGRA software can be used to model a permanent magnet in 2D and 3D, with accuracy matching that of the open-source commercial software FEMM. This is done by conducting simulations and experimental measurements for a commercial-grade N42 neodymium alloy ring magnet with a measured magnetic field strength of approximately 0.4 T in its immediate vicinity. Transient simulations using ALEGRA and static simulations using FEMM are conducted. Comparisons are made between simulations and measurements, and amongst the simulations, for sample locations in the steady-state magnetic field. The comparisons show that all models capture the data to within 7%. The FEMM and ALEGRA results agree to within approximately 2%. The most accurate solutions in ALEGRA are obtained using quadrilateral or hexahedral elements. In the case where iron shielding disks are included in the magnetized space, ALEGRA simulations are considerably more expensive because of the increased magnetic diffusion time, but FEMM and ALEGRA results are still in agreement. The magnetic field data are portable to other software interfaces using the Exodus file format.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidating the Link Between Alkali Metal Ions and Reaction-Transport Mechanisms in Cathode Electrodes for Alkali-ion Batteries

Our long-term goal is to improve the reliability of electrode materials and their ability to transport and store various metal ions for electrochemical energy storage applications. The main objective of this work was to investigate the intrinsic relationship between the role of alkali metal ions and electrochemically driven mechanical stability and kinetic properties of battery materials. The overall question was “What is the role of alkali metal ions on the electrochemical and mechanical behavior of cathode electrodes? Our guiding hypothesis was that intercalation of larger alkali metal ions (Na and K) inevitably alters the coupled transport-reaction processes during battery operation in organic electrolytes, leading to more intensive chemo-mechanical instabilities in cathode electrodes, resulting in rapid capacity fade. To validate the hypothesis, we experimentally characterized the reaction-transport processes and governing forces driving the instability of electrode materials in different alkali metal-ion environments. The project had three main tasks. The first one was to investigate intercalation-induced strains and associated stress generation, and their impact on structural deformations in composite cathode electrodes. The second task focused on identifying potential-dependent dynamic changes in the electrode-electrolyte interface in alkali metal ion batteries. The last task was focused on determining how larger alkali metal ions with slower diffusivity affect the transport-mechanics coupling at faster scan rates, compared to smaller ions with faster diffusivity in electrodes. We shortly provided the outcome of each task in the accomplishment section. This project produced 10 peer-reviewed publications (9 research papers and one review manuscript) and supported two Ph.D. students, who graduated from Oklahoma State University.

25 ENERGY STORAGE↗
DOE OSTI2020

Investigating the Linear Thermal Expansion of Additively Manufactured Multi-Material Joining between Invar and Steel

This work investigated the linear thermal expansion properties of a multi-material specimen fabricated with Invar M93 and A36 steel. A sequence of tests was performed to investigate the viability of additively manufactured Invar M93 for lowering the coefficient of thermal expansion (CTE) in multi-material part tooling. Invar beads were additively manufactured on a steel base plate using a fiber laser system, and samples were taken from the steel, Invar, and the interface between the two materials. The CTE of the samples was measured between 40 °C and 150 °C using a thermomechanical analyzer, and the elemental composition was studied with energy dispersive X-ray spectroscopy. The CTE of samples taken from the steel and the interface remained comparable to that of A36 steel; however, deviations between the thermal expansion values were prevalent due to element diffusion in and around the heat-affected zone. The CTEs measured from the Invar bead were lower than those from the other sections with the largest and smallest thermal expansion values being 10.40 μm/m-K and 2.09 μm/m-K. In each of the sections, the largest CTE was measured from samples taken from the end of the weld beads. An additional test was performed to measure the aggregate expansion of multi-material tools. Invar beads were welded on an A36 steel plate. The invar was machined, and the sample was heated in an oven from 40 °C and 160 °C. Strain gauges were placed on the surface of the part and were used to analyze how the combined thermal expansions of the invar and steel would affect the thermal expansion on the surface of a tool. There were small deviations between the expansion values measured by gauges placed in different orientations, and the elongation of the sample was greatest along the dimension containing a larger percentage of steel. On average, the expansion of the machined Invar surface was 42% less than the expansion of the steel surface. The results of this work demonstrate that additively manufactured Invar can be utilized to decrease the CTE for multi-material part tooling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Turbulent mixing by microscopic self-assembled spinners

A system for mixing particles that includes a liquid comprising inert particles and defining a liquid and air interface; magnetic microparticles suspended at the liquid and air interface; and a magnetic source configured to apply a uniaxial alternating magnetic field parallel to the liquid and air interface, wherein the uniaxial alternating magnetic field promotes a turbulent motion of the magnetic microparticles, which in turn promotes a diffusive motion of the inert particles.

Snezhko, Oleksiy↗
DOE OSTI2025

Dissolved CO 2 Modulates the Electrochemical Capacitance on Gold Electrodes

The presence of CO 2 at an electrified interface between an aqueous electrolyte and a metal electrode is the prerequisite for many electrochemical CO 2 capture technologies. To understand the behavior of dissolved CO 2 at an aqueous electrified interface, we characterized the electrochemical interface of planar gold electrodes with cyclic voltammetry, electrochemical impedance spectroscopy (EIS), electrochemical surface plasmon resonance (EC–SPR), and attenuated total reflectance surface-enhanced infrared spectroscopy (ATR–SEIRAS). Under all investigated conditions, we observed a decrease in the electrochemical capacitance upon saturation of the electrolyte with CO 2 , as compared to an electrolyte saturated with Ar. EIS and EC–SPR showed that this capacitance reduction was also potential dependent: it reached a minimum near the point of zero charge and became more significant as the applied potential moved further away from the point of zero charge. Hybrid quantum–classical simulations of the gold/aqueous electrolyte interface indicate that bicarbonate decreases the capacitance and modifies the composition of the electric double layer. In addition to the binding of bicarbonate under positive bias, we propose that molecular CO 2 can be induced by applied potential to concentrate in the diffuse layer of the electric double layer, leading to a reduction in the electrochemical capacitance under both negative and positive bias. Furthermore, this work advances the understanding of non-Faradaic effects of dissolved CO 2 at aqueous electrified interfaces of relevance for electrochemical CO 2 capture.

25 ENERGY STORAGE↗
DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Microstructure characteristics of LPBF&HIP fabricated graded composite transition joint between ferritic steel and austenitic stainless steel

Graded composite transition joints (GCTJs) offer a promising alternative to conventional dissimilar metal welds (DMWs) by enabling smooth compositional and microstructural transitions. However, GCTJs fabricated solely through additive manufacturing (AM) face challenges such as heat accumulation, complex parameter control, and elemental segregation. In this study, we propose a novel approach that relies on AM to design a spatially graded structure in one alloy and then employs hot isostatic pressing (HIP) as a diffusion bonding method to join it with a second alloy. Here, this method combines the flexibility of AM with the powder net-shaping advantage of HIP. Specifically, a series of closely packed austenitic stainless steel 304 conical structures were printed using laser powder bed fusion (LPBF) and then combined with ferritic steel P91 powder via HIP. By using electron backscatter diffraction (EBSD), electron probe microanalysis (EPMA), and transmission electron microscopy (TEM) techniques, the microstructure characteristics of the GCTJ of 304&P91, especially the interdiffusion zone (IDZ), have been systematically investigated. The microstructure at the interface transitions from austenite-ferrite (A+F) to austenite-martensite-ferrite (A+M+F), and finally to martensite-ferrite (M+F) due to diffusion. Additionally, the diffusion width between 304 and P91 increases with the volume fraction of P91. This unique design also ensures a gradual transition in both hardness and thermal expansion coefficient from 304 to P91, thereby enabling a smooth gradient in functional properties. Overall, this study proposes a novel approach for fabricating GCTJs and contributes to advancing design concepts in the field of dissimilar metal joining.

Additively manufacturing (AM)↗
DOE OSTI2023

Fundamental study of gas species transport in the oxygen electrode of solid oxide fuel and electrolysis cells

A fundamental analysis of multicomponent gas transport models was performed in application to the oxygen electrodes of solid oxide cells. It is common practice to neglect the effect of pressure gradients within oxygen electrodes, even though a net molar flux at the electrolyte surface implies that a pressure gradient must exist. The influence of both Darcy velocity and Knudsen flux are considered in the context of ordinary (Fickian) diffusion, the dusty gas model, and the binary friction model. Comparisons between the models and different sets of assumptions are made via parametric studies on operating load, oxygen partial pressure, microstructural properties, and electrode thickness. Results show that the pressure gradient will have a significant impact on the oxygen concentration distribution and therefore the concentration overpotential. In electrolysis mode, pressure increases up to 1 atm are predicted, indicating that pressure at the electrode/electrolyte interface could contribute to electrode delamination. Additionally, it is found that Darcy's law is insufficient for calculating the pressure distribution without accounting for the flux due to Knudsen diffusion. Additionally, it is found that for the range of properties typical of oxygen electrodes, there is a negligibly small difference between the dusty gas model and binary friction model from a practical standpoint.

08 HYDROGEN↗