Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse interface”

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

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2022

Spatiotemporal Visualization and Chemical Identification of the Metal Diffusion Layer at the Electrochemical Interface

The diffusion layer created by transition metal dissolution is ubiquitous at the electrochemical solid-liquid interface and plays a key role in determining electrochemical performance. Tracking the spatiotemporal dynamics of the diffusion layer has remained an unresolved challenge. With spatially resolved synchrotron X-ray fluorescence microscopy and micro-X-ray absorption spectroscopy, we demonstrate the in situ visualization and chemical identification of the dynamic diffusion layer near the electrode surface under electrochemical operating conditions. Finally, our method allows for direct mapping of the reactive electrochemical interface and provides insights into engineering the diffusion layer for improving electrochemical performance.

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

On the importance of species immiscibility in mixing-layer dynamics at supercritical pressures

The mixing dynamics of injected propellants is a key factor in determining the ignition performance and combustion-instability response of rocket engines, internal combustion engines, and gas turbines. A key source of uncertainty in the prediction of phase-exchange dynamics is the immiscibility of the injected propellants at supercritical pressure conditions, which induces liquid/vapor phase separation and surface-tension dynamics. While experimental observations indicate the presence of liquid/vapor interfacial structures, this interfacial dynamics is typically neglected in numerical analyses. To address this issue, the objective of the present study is to systematically evaluate the importance of species immiscibility on the phase-exchange dynamics of cryogenic LOX/GH2 mixing layers at typical rocket engine injection conditions. This is accomplished by comparing simulations of (i) a recently developed interface-capturing Regularized-Interface Method (RIM) formulation, and (ii) the commonly employed Diffuse-Interface Method (DIM) formulation. Analysis shows that the interfacial dynamics significantly impact the atomization and mixing of the propellants in the near-injector region, which is not captured by the DIM formulation. The findings of this study extend to other immiscible injection systems, such as LOX/kerosene and hydrocarbon-fuel/air, thereby highlighting the importance of resolving species immiscibility in simulating high-pressure combustion engines.

diffused-interface method↗
NASA NTRS1992

Effect of Index of Refraction on Radiation Characteristics in a Heated Absorbing, Emitting, and Scattering Layer

The index of refraction can considerably influence the temperature distribution and radiative heat flow in semitransparent materials such as some ceramics. For external radiant heating, the refractive index influences the amount of energy transmitted into the interior of the material. Emission within a material depends on the square of its refractive index, and hence this emission can be many times that for a biackbody radiating into a vacuum. Since radiation exiting through an interface into a vacuum cannot exceed that of a blackbody, there is extensive reflection at the internal surface of an interface, mostly by total internal reflection. This redistributes energy within the layer and tends to make its temperature distribution more uniform. The purpose of the present analysis is to show that, for radiative equilibrium in a gray layer with diffuse interfaces, the temperature distribution and radiative heat flux for any index of refraction can be obtained very simply from the results for an index of refraction of unity. For the situation studied here, the layer is subjected to external radiative heating incident on each of its surfaces. The material emits, absorbs, and isotropically scatters radiation. For simplicity the index of refraction is unity in the medium surrounding the layer. The surfaces of the layer are assumed diffuse. This is probably a reasonable approximation for a ceramic layer that has not been polished. When transmitted radiation or radiation emitted from the interior reaches the inner surface of an interface, the radiation is diffused and some of it thereby placed into angular directions for which there is total internal reflection. This provides a trapping effect for retaining energy within the layer and tends to equalize its temperature distribution. An analysis of temperature distributions in absorbing-emitting layers, including index of refraction effects, was developed by Gardon (1958) to predict cooling and heat treating of glass plates. The interfaces were optically smooth; the resulting specular reflections were computed from the Fresnel reflection laws. This provides a somewhat different behavior than for diffuse interfaces. A similar application was for heating that occurs in a window of a re-entry vehicle (Fowle et al., 1969). A number of recent papers (Rokhsaz and Dougherty, 1989; Ping and Lallemand, 1989; Crosbie and Shieh, 1990) further examined the effects of Fresnel boundary reflections and nonunity refractive index. Other examples of analyses of both steady and transient heat transfer to single or multiple plane layers (Amlin and Korpela, 1979; Tarshis et al., 1969) have used diffuse assumptions at the interfaces as in the present study

Siegel, R.↗
DOE OSTI2020

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗
DOE OSTI2026

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗
DOE OSTI2023

Analysis of ducted fuel injection at high-pressure transcritical conditions using large-eddy simulations

Ducted fuel injection (DFI) is a proposed fuel injection concept for achieving substantial reductions in emissions. In this concept, the fuel is injected through a coannular duct, resulting in increased fuel-air mixing and minimized formation of soot and other unwanted combustion products. Apart from comprehensive experimental investigations on DFI, so far computational studies have been limited to single-point Reynolds-averaged Navier Stokes simulations. Therefore, the objective of this work is to complement these studies by performing large-eddy simulations using a diffuse-interface method to examine the physical mechanisms and combustion processes of DFI, specifically focusing on the mixing process and the effect of fuel-ducting on combustion and pollutant emissions. To this end, finite-rate chemistry simulations are performed of the DFI configuration corresponding to the Engine Combustion Network Spray A injector at transcritical conditions (n-dodecane fuel, 60 bar pressure and 1000 K temperature chamber conditions). A two-equation soot model is employed for the qualitative analysis of soot emissions. Direct comparisons of averaged and instantaneous flow field results with the Spray A configuration are performed to assess the effect of DFI on the first- and second-stage ignition and soot formation. Compared to the free-spray condition, the results show that the DFI case exhibits a combination of (i) increased mass flow rate and entrained air, (ii) larger pressure drop magnitude and flow velocity, and (iii) a closer-to-stoichiometric mixture composition (both globally and locally), each of which is conjectured to contribute toward reduced soot production.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1987

Predicting diffusion paths and interface motion in gamma/gamma + beta, Ni-Cr-Al diffusion couples

A simplified model has been developed to predict Beta recession and diffusion paths in ternary gamma/gamma + beta diffusion couples (gamma:fcc, beta: NiAl structure). The model was tested by predicting beta recession and diffusion paths for four gamma/gamma + beta, Ni-Cr-Al couples annealed for 100 hours at 1200 C. The model predicted beta recession within 20 percent of that measured for each of the couples. The model also predicted shifts in the concentration of the gamma phase at the gamma/gamma + beta interface within 2 at. pct Al and 6 at. pct Cr of that measured in each of the couples. A qualitative explanation based on simple kinetic and mass balance arguments has been given which demonstrates the necessity for diffusion in the two-phase region of certain gamma/gamma + beta, Ni-Cr-Al couples.

Nesbitt, J. A.↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Inelastic phonon transport across atomically sharp metal/semiconductor interfaces

Understanding thermal transport across metal/semiconductor interfaces is crucial for the heat dissipation of electronics. The dominant heat carriers in non-metals, phonons, are thought to transport elastically across most interfaces, except for a few extreme cases where the two materials that formed the interface are highly dissimilar with a large difference in Debye temperature. In this work, we show that even for two materials with similar Debye temperatures (Al/Si, Al/GaN), a substantial portion of phonons will transport inelastically across their interfaces at high temperatures, significantly enhancing interface thermal conductance. Moreover, we find that interface sharpness strongly affects phonon transport process. For atomically sharp interfaces, phonons are allowed to transport inelastically and interface thermal conductance linearly increases at high temperatures. With a diffuse interface, inelastic phonon transport diminishes. Our results provide new insights on phonon transport across interfaces and open up opportunities for engineering interface thermal conductance specifically for materials of relevance to microelectronics.

42 ENGINEERING↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-Principles Studies of Tritium Species Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

α-Zr and its alloys are known as best 3 H getters due to their excellent corrosion resistance under chemically corrosive environment, low thermal neutron absorption cross-section and robust mechanical strength at high temperature. In this work, we conducted a systematic study on constructing and optimizing Ni and Zr surfaces, interfacing the stable Ni and Zr surfaces to create an optimal Ni-Zr interface, and an understanding the diffusion mechanisms for 3 H across Ni(111)-Zr(0001) interface with and without vacancies and impurities such as Sn and O.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

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

Notes on Interface Growth Kinetics 50 Years after Burton, Cabrera and Frank

This overview is devoted to some unresolved basic problems in crystal growth kinetics. The density wave approach to propagation of a spatially diffuse interface between a growing crystal and its simple (e.g., metallic) melt is discussed is Section 2. This approach allows for the calculation of kinetic coefficients and is an alternative to the localized interface concept in which each atom belongs to either a solid or a liquid. Sections 3 and 4 deals mainly with layer growth from solution. Mutual retardation of the growth steps via their bulk and d a c e diffusion fields is the major subject. The influence of solution flow on step bunching (Section 4) suggests the essential influence of bulk diffusion on the surface morphology. The flow within the solution boundary layer enhances step-step interaction, influences the step bunching process and the resulting step pattern morphology on the growing surface. Recent experiments on the rates at which strongly polygonized steps on protein and small molecule crystals propagate during growth from solution are analyzed in Section 5 . We have shown that the step segments may be "singular" and that "one-dimensional nucleation" may be the rate limiting stage for the segments that are shorter or comparable in length to the thermodynamically equilibrium interlink distance. In this case, the reciprocal dependence of the segment propagation rate on the segment length that follow from the Gibbs-Thompson law, should be replaced by an abrupt switch from zero to a finite constant velocity. Until recently, the Kossel crystal remained the only model used in crystal growth theory. In such Kossel Gibbs-Thomson law, should be replaced by an abrupt switch &om zero to a finite constant velocity. crystals, all kinks at the steps are identical and the kink rate is a linear function of the supersaturation. In the non-Kossel crystals, there may be several kink configurations characterized by different geometries and energies. These configurations should appear in a specific sequence when each new lattice unit cell is filled. As a result of such a cooperative interaction within the unit cell, a non-linear dependence of the kink rate on the vapor pressure or solution concentration in excess over the equilibrium value should be expected.

Chernov, A. A.↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2022

Emergent interface vibrational structure of oxide superlattices

As the length scales of materials decrease, the heterogeneities associated with interfaces become almost as important as the surrounding materials. This has led to extensive studies of emergent electronic and magnetic interface properties in superlattices. However, the interfacial vibrations that affect the phonon-mediated properties, such as thermal conductivity, are measured using macroscopic techniques that lack spatial resolution. Although it is accepted that intrinsic phonons change near boundaries, the physical mechanisms and length scales through which interfacial effects influence materials remain unclear. Here we demonstrate the localized vibrational response of interfaces in strontium titanate–calcium titanate superlattices by combining advanced scanning transmission electron microscopy imaging and spectroscopy, density functional theory calculations and ultrafast optical spectroscopy. Structurally diffuse interfaces that bridge the bounding materials are observed and this local structure creates phonon modes that determine the global response of the superlattice once the spacing of the interfaces approaches the phonon spatial extent. Our results provide direct visualization of the progression of the local atomic structure and interface vibrations as they come to determine the vibrational response of an entire superlattice. Direct observation of such local atomic and vibrational phenomena demonstrates that their spatial extent needs to be quantified to understand macroscopic behaviour. Tailoring interfaces, and knowing their local vibrational response, provides a means of pursuing designer solids with emergent infrared and thermal responses.

36 MATERIALS SCIENCE↗