Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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

Phase field-volumetric lattice Boltzmann model of ion uptake in porous nuclear waste form materials under continuous flow

The flow field within the mesopores of sorbent particles plays a crucial role in radionuclide diffusion and ion uptake kinetics, thus, impacting the overall performance of porous nuclear waste form materials. To fundamentally understand the influence of microstructures and material properties on the radionuclide absorption and retention processes requires a coupled multi-physics model that considers the advection and diffusion within the flow field, the reaction at liquid-solid interfaces, and finally, the solid-state diffusion within a complex nanoporous medium. Here, this study employs the volumetric lattice Boltzmann method (VLBM) to accurately and efficiently calculate the steady state velocity field inside the mesopores of sorbent particles. The obtained velocity field is then utilized to calculate the advection of ions in the steady flow. A phase field (PF) model of ion uptake is used to describe the reaction occurring at the solid-liquid interface and diffusion inside the porous medium. The integrated PF-VLBM model is verified in terms of the mass conservation and numerical efficiency and validated qualitatively with experimental observation data. Then, it is applied to study the influence of thermodynamic and kinetic properties, as well as flow field conditions on the ion uptake kinetics. The numerical results demonstrate that the ion uptake kinetics in porous particles has three distinct stages, which is in agreement with the observations in continuous flow experiments. In the first stage, the kinetics is predominantly controlled by the flow field and ion diffusivity in the liquid phase. The kinetics in the second stage is primarily governed by ion diffusivity in the solid phase. In the third stage the system reaches a dynamic equilibrium with a net zero uptake flux at the interface. It is also found that porous structures significantly affect the efficiency and capacity of ion uptake. The simulation results can help to understand the physics behind the observed ion uptake kinetics in experiments and to facilitate the development of constitutive equations that can account for heterogeneous microstructures in engineering performance codes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Accommodation of Silicon in an Interconnected Copper Network for Robust Li-Ion Storage

Silicon (Si)-based materials are one of the most promising anodes to be applied in rechargeable lithium ion batteries. However, the active Si/electrolyte interface causes continuous side reactions and poor conductivity, which significantly decreases the cycling stability. Cu is the only metallic current collector that has been known to promote electron conduction and lithium-ion transfer without alloying reaction occurrence. However, to the best current knowledge, scalable interface engineering incorporating Cu has not been reported. In this study, this conductive Cu interface (CCI) is constructed through a self-assembly carbothermic reduction method to achieve efficient protection of Si/electrolyte interfaces while allowing for fast Li + diffusion. The energy barrier of lithium-ion diffusion through Cu is calculated to be 0.1965 eV, which is much lower than that through Au, Fe, and Ni films. Benefiting from the enhanced interfacial protection and kinetics of Si with CCI, a fading rate of only 0.068% is maintained for 1000 cycles and an aerial capacity of 4.78 mAh cm -2 is achieved after 280 cycles, which is comparable to the industry standards required for practical application.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tensile behavior of diffusion bonded AA6061 - AA6061 with variation in cooling method

Hot isostatically pressed AA6061 cladding is an important structural component of the high performance, Zr-laminated U-10Mo monolithic fuel system for the application in research and test reactors. In this study, the mechanical behavior of two diffusion bonded aluminum alloy, AA6061, was examined using tensile testing. Solid-to-solid diffusion bonding between two pieces of AA6061 was performed by isothermal annealing at 560 °C for 1.5 h, and diffusion couples were subsequently cooled via three different cooling methods: furnace cooling, air cooling, and water quenching. Dog-bone shaped tensile specimens, with 10 mm in gauge length (with diffusion bonded interface in the middle), and 1.5 × 1.5 mm 2 gauge cross-sections, were fabricated from the diffusion bonded AA6061 by electro-discharge machining. Yield strength (% EL at failure) of furnace cooled, air cooled and water quenched tensile specimens determined was 82 ~ 89 MPa (10 ~ 30%), 112 ~ 116 MPa (10 ~ 14%), and 149 ~ 164 MPa (10 ~ 17%), respectively. This variation in mechanical behavior was examined with cooling-rate dependent, concentrated precipitation of Mg 2 Si at the diffusion bonded interface, with due respect for mechanical properties of the AA6061 alloy that inherently vary as a function of cooling rate from 560 °C. Finite element analysis using ABAQUS was employed to augment experimental findings with the appropriate microstructural constituents and alloy properties. Results suggest that the strength is dominated by matrix/bulk properties of AA6061, while ductility is strongly influenced by the cooling method dependent presence of Mg 2 Si precipitates at the interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1985

Cellular interface morphologies in directional solidification. III - The effects of heat transfer and solid diffusivity

The shape and stability of two-dimensional finite-amplitude cellular interfaces arising during directional solidification are compared for several solidification models that account differently for latent heat released at the interface, unequal thermal conductivities of melt and solid, and solute diffusivity in the solid. Finite-element analysis and computer-implemented perturbation methods are used to analyze the families of steadily growing cellular forms that evolve from the planar state. In all models a secondary bifurcation between different families of finite-amplitude cells exists that halves the spatial wavelength of the stable interface. The quantitative location of this transition is very dependent on the details of the model. Large amounts of solute diffusion in the solid retard the growth of large-amplitude cells.

Ungar, Lyle H.↗
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 OSTI2024

Probing elemental diffusion and radiation tolerance of perovskite solar cells via non-destructive Rutherford backscattering spectrometry

Mixed organic–inorganic halide perovskite-based solar cells have attracted interest in recent years due to their potential for both terrestrial and space applications. Analysis of interfaces is critical to predicting device behavior and optimizing device architectures. Most advanced tools to study buried interfaces are destructive in nature and can induce further degradation. Ion beam techniques, such as Rutherford backscattering spectrometry (RBS), is a useful non-destructive method to probe an elemental depth profile of multilayered perovskite solar cells (PSCs) as well as to study the inter-diffusion of various elemental species across interfaces. Additionally, PSCs are becoming viable candidates for space photovoltaic applications, and it is critical to investigate their radiation-induced degradation. RBS can be simultaneously utilized to analyze the radiation effects induced by He+ beam on the device, given their presence in space orbits. In the present work, a 2 MeV He+ beam was used to probe the evidence of elemental diffusion across PSC interfaces with architecture glass/ITO/SnO2/Cs0.05(MA0.17FA0.83)0.95Pb(I0.83Br0.17)3/spiro-OMeTAD/MoO3/Au. During the analysis, the device active area was exposed to an irradiation equivalent of up to 1.62 × 1015 He+/cm2, and yet, no measurable evidence (with a depth resolution ∼1 nm) of beam-induced ion migration was observed, implying high radiation tolerance of PSCs. On the other hand, aged PSCs exhibited indications of the movement of diverse elemental species, such as Au, Pb, In, Sn, Br, and I, in the active area of the device, which was quantified with the help of RBS.

14 SOLAR ENERGY↗
DOE OSTI2020

Electric current–assisted direct joining of silicon carbide

Conventional direct joining technologies are difficult to use with silicon carbide (SiC) materials, especially for fiber composite forms of SiC, because of the harsh conditions required. Therefore, to reduce the temperature and/or process time required for the direct joining process, an electric current–assisted joining (ECAJ) method was studied. Joining of low–resistivity grade, nitrogen doped β-SiC was demonstrated at a relatively low nominal temperature of 1750 °C with a 10 min hold by enhancing the passage of current through the material. The joining mechanism is discussed in terms of localized overheating and accelerated self-diffusion at the interface. In the case of joining at 2160 °C for 1 min, rapid crystal growth of textured SiC was found at the interface. This study indicates that rapid ECAJ-based direct joining is a practical and appropriate method for joining SiC-based materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enable High-Energy LiNi0.5Co0.2Mn0.3O2 by Ultra-Thin Coating through Wet Impregnation

A high cut-off voltage is required for nickel-rich layered oxide LiNixCoyMnzO2 (NCM) to meet the high energy density requirement of lithium-ion batteries in electric vehicles. However, such a high voltage application leads to an unstable interface between NCM and liquid electrolytes. To stabilize the interface, the facile wet impregnation method has been developed to apply an ultra-thin Al2O3 coating layer on the NCM particles. This coating layer was found to have a strong interaction with the NCM and resulted in Al-doped NCM at the surface structure of NCM. The change of surface structure can not only reduce the surface resistance of lithium diffusion of LiNi0.5Co0.2Mn0.3O2 (NCM523), but also stabilize the solid electrolyte interface between NCM523 and the electrolyte with the cut-off voltage of 4.5 V vs. Li/Li+. Compared to other coating methods, wet impregnation coating can provide an ultra-thin and uniform coating with surface doping on NCM particles. Furthermore, this scalable coating method can be applied to various electrode materials without adding much additional cost.

25 ENERGY STORAGE↗
DOE OSTI2024

Kinetics and Mechanism of Ce(IV) Phase Transfer by the Neutral Extractants Tributyl Phosphate and N,N -Di(2-ethylhexyl)isobutyramide

Tributyl phosphate and N,N-di(2-ethylhexyl)isobutyramide (DEHiBA) are uranium-selective extractants that could be used to recover low-enriched uranium from commercial used nuclear fuel. The economic viability of liquid-liquid extraction processes for chemical separations depends on both thermodynamic and kinetic considerations. Here we used microfluidic droplet extraction to measure interfacial phase transfer rate constants for the extraction and stripping of Ce(IV), a non-radioactive actinide surrogate. The Damkohler numbers for each system were calculated and the phase transfer processes were determined to be diffusion-limited under all conditions, contrary to prior studies using constant interface stirred cells and single drop methods. A well-extracted 2:1 DEHiBA/Ce(IV) complex in the organic phase was determined from batch distribution studies, suggesting a different extraction mechanism for Ce(IV) compared with poorly extracted Pu(IV). Finally, phase disengagement rates were found to be rapid in both systems, in agreement with other studies.

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

An extended classical solution of the droplet growth problem

Problems of applying the classical kinetic theory to the growth of small droplets from vapor are examined. A solution for the droplet growth equation is derived which is based on the assumption of a diffusive field extending to the drop surface. The method accounts for partial thermal and mass accommodation at the interface and the kinetic limit to the mass and heat fluxes, and it avoids introducing the artifact of a discontinuity in the thermal and vapor field near the droplet. Consideration of the environmental fields in spherical geometry utilizing directional fluxes yields boundary values in terms of known parameters and a new Laplace transform integral.

Anderson, B. J.↗
DOE OSTI2021

An efficient reconstruction algorithm for diffusion on triangular grids using the nodal discontinuous Galerkin method

High-energy-density (HED) hydrodynamics studies such as those relevant to inertial confinement fusion and astrophysics require highly disparate densities, temperatures, viscosities, and other diffusion parameters over relatively short spatial scales. This presents a challenge for high-order accurate methods to effectively resolve the hydrodynamics at these scales, particularly in the presence of highly disparate diffusion. A significant volume of engineering and physics applications use an unstructured discontinuous Galerkin (DG) method developed based on the finite element mesh generation and algorithmic framework. This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection–diffusion equations. However, to the best of the authors’ knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. As a result, the algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.

Computational efficiency↗
DOE OSTI2022

Method for fabricating embedded nanostructures with arbitrary shape

A layered heterostructure, comprising alternating layers of different semiconductors, wherein one of the atom species of one of the semiconductors has a faster diffusion rate along an oxidizing interface than an atom species of the other semiconductor at an oxidizing temperature, can be used to fabricate embedded nanostructures with arbitrary shape. The result of the oxidation will be an embedded nanostructure comprising the semiconductor having slower diffusing atom species surrounded by the semiconductor having the higher diffusing atom species. The method enables the fabrication of low- and multi-dimensional quantum-scale embedded nanostructures, such as quantum dots (QDs), toroids, and ellipsoids.

Wang, George T.↗
DOE OSTI2022

Multifunctional transitional metal-based phosphide nanoparticles towards improved polysulfide confinement and redox kinetics for highly stable lithium-sulfur batteries

The shuttle effect and the sluggish redox kinetics of lithium polysulfides (LiPSs) are the major issues impeding the practical applications of lithium-sulfur batteries (LSBs). Herein, a highly-efficient Ni 2 P electrocatalyst supported on N, P co-doped graphene (Ni 2 P@NPG) is developed via a simple “recrystallization-self-assembly” method to address the above issues. Here, the ultrafine Ni 2 P nanoparticles ensure abundant adsorption-diffusion-conversion interfaces for accelerating LiPSs transformation and Li 2 S deposition, which extremely decreases the accumulation of LiPSs in the electrolyte and therefore prevents the migration of LiPSs. Their superior catalytic performance is demonstrated by reduced Gibbs free energy changes of rate-limiting step based on the systematic theoretical calculations and the reduced shuttle effect is tested by the three-dimensional reconstructions of Raman depth profiles. Benefiting from these synergistic effects, the LSBs with Ni 2 P@NPG modified separators present a superior cycling performance with an average capacity decay rate of 0.048 % per cycle at 1C around the 400 cycles and a high-rate capacity of 731 mAh/g at 2C. Even with a high-sulfur loading of 3.53 mg cm –2 , the cell can still contain a reversible capacity of 809 mAh/g at 0.2C with a remarkable columbic efficiency of 98.4 %.

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

Brazing method produces solid-solution bond between refractory metals

Brazing two refractory metals by diffusion bonding minimizes distortion and avoids excessive grain growth in the metals. This method requires the selection of an interface metal that forms intermediate low-melting eutectics or solid solutions with the metals to be brazed.

Source record↗
DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

42 ENGINEERING↗
DOE OSTI2023

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗