Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2

A new interatomic potential for mixed Mg-Al-Ga-In spinels

While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. Here, in this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO 4 structure while correctly predicting that neither MgAlInO 4 nor MgGaInO 4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.

36 MATERIALS SCIENCE↗

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↗

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗

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

In this study, we will calculate the 3H2 and 3H2O diffusion mechanism in the Ni coated layer and the dominant 3H species dissolution into the Zircaloy-4 getter. We will optimize the Ni/Zr interface and introduce Sn and O impurities. To achieve the objectives, we will perform the investigation by using spin-polarized density functional theory (DFT) methods to capture the effects of Ni magnetic properties on 3H diffusion. We will introduce NiO/Ni(OH) on the Ni side closer to the interface. Then, the 3H diffusion pathways and barriers from Ni side to the zircalory-4 will be explored. We will investigate how NiO/Ni(OH) affects the 3H diffusion across the interface.

Paudel, Hari P.↗

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

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↗