SEARCH · Engineering Papers
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.
This work proposes an analytical model considering the effects of hydrodynamic drag and energy barriers induced by liquid solvation forces to predict the in-plane translational diffusivity of a nanoparticle physically adsorbed on a wetted surface.
Abstract Constructing an artificial solid electrolyte interphase (SEI) on lithium metal electrodes is a promising approach to address the rampant growth of dangerous lithium morphologies (dendritic and dead Li 0 ) and low Coulombic efficiency that plague development of lithium metal batteries, but how Li + transport behavior in the SEI is coupled with mechanical properties remains unknown. We demonstrate here a facile and scalable solution-processed approach to form a Li 3 N-rich SEI with a phase-pure crystalline structure that minimizes the diffusion energy barrier of Li + across the SEI. Compared with a polycrystalline Li 3 N SEI obtained from conventional practice, the phase-pure/single crystalline Li 3 N-rich SEI constitutes an interphase of high mechanical strength and low Li + diffusion barrier. We elucidate the correlation among Li + transference number, diffusion behavior, concentration gradient, and the stability of the lithium metal electrode by integrating phase field simulations with experiments. We demonstrate improved reversibility and charge/discharge cycling behaviors for both symmetric cells and full lithium-metal batteries constructed with this Li 3 N-rich SEI. These studies may cast new insight into the design and engineering of an ideal artificial SEI for stable and high-performance lithium metal batteries.
Orthopyroxene crystals in a number of meteorites exhibit compositional zoning of Fe and Mg, which provide important constraint on their cooling rates. However, attempts to model cooling rate of these crystals from Fe-Mg zoning profiles suffer from the lack of any measured or theoretically well constrained Fe-Mg interdiffusion data in OP(x) It has been assumed that Fe-Mg interdiffusion in OP(x) only slightly slower than that in olivine. The purpose of this paper is to (1) calculate the Fe-Mg fractionation, and (2) provide analytical formulation relating cooling rate to the length of the diffusion zone across the interface of the overgrowth of a mineral on itself with application to Mg diffusion profile across OP(x) growth on OP(x) in certain mesosiderites.
We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.
Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.
Polydicyclopentadiene, p(DCPD), is a high‐performance thermoset valued for its exceptional toughness, strength, and stiffness. When copolymerized with 1,5‐cyclooctadiene (COD), its mechanical properties can be tuned from glassy to rubbery at room temperature. While frontal polymerization enables a rapid and energy‐efficient route to 3D print DCPD‐based materials, challenges such as ink shelf life and gravitational distortion, especially in direct ink writing of soft COD‐rich formulations, must be considered. Here, a complementary chemical strategy is presented, embedded 3D printing, that enables localized in situ polymerization of printed DCPD/COD inks within a reactive support matrix. The matrix provides both physical support and a reservoir of chemical activator, which diffuses into the ink, activates a latent bis(N‐heterocyclic carbene) Ru precatalyst, and initiates ring‐opening metathesis polymerization. Curing begins at the ink–matrix interface and propagates inward via diffusion, stabilizing the interface and preventing capillary‐driven deformation regardless of the matrix yield stress. This approach eliminates the need for cold storage, external curing, or photoinitiation, significantly expanding the processing window. Using this method, diverse thermosetting and elastomeric architectures are fabricated with features as small as 5 µm and aspect ratios of 100, including interlinked chains, shallow spherical shells exhibiting snap‐through buckling, and hair‐like fin arrays inaccessible through traditional techniques.
Schottky diodes of Palladium/SiC are good candidates for hydrogen and hydrocarbon gas sensors at elevated temperature. The detection sensibility of the diodes has been found heavily temperature dependent. In this work, emphasis has been put on the understanding of changes of physical and chemical properties of the Schottky diodes with variation of temperature. Schottky diodes were made by depositing ultra-thin palladium films onto silicon carbide substrates. The electrical and chemical properties of Pd/SiC Schottky contacts were studied by XPS and AES at different annealing temperatures. No significant change in the Schottky barrier height of the Pd/SiC contact was found in the temperature range of RT-400 C. However, both palladium diffused into SiC and silicon migrated into palladium thin film as well as onto surface were observed at room temperature. The formation of palladium compounds at the Pd/SiC interface was also observed. Both diffusion and reaction at the Pd/SiC interface became significant at 300 C and higher temperature. In addition, silicon oxide was found also at the interface of the Pd/SiC contact at high temperature. In this report, the mechanism of diffusion and reaction at the Pd/SiC interface will be discussed along with experimental approaches.
This work focuses on the thermal oxidation of silicon near 1273 K using the double-tracer oxidation method. The results confirm that oxidation occurs by the transport of electrically neutral nonnetwork oxygen through the interstitial space of the vitreous silica (v-SiO2) scale. Simultaneously, self- (or isotopic-) diffusion occurs in the network, resulting in characteristic isotopic fraction distributions near the gas-scale interface. The self-diffusion coefficients calculated from these profiles agree with those reported for tracer diffusion in v-SiO2, and the diffusion coefficient calculated from the scale growth is consistent with reported O2 permeation data. An important parameter that describes the double-oxidation behavior is a ratio relating the scale thickness grown during the second oxidation, the network self-diffusion coefficient for oxygen, and the time of the second oxidation.
Silicon photodiodes have been developed by defect-free phosphorus diffusion having practically no carrier recombination at the SiSiO2 interface or in the front diffused region. The quantum efficiency of these photodiodes was found to be around 120 percent at 100 nm. Unlike the previously tested silicon photodiodes, the developed photodiodes exhibit extremely stable quantum efficiency over extended periods of time. The possibility of using these photodiodes as vacuum ultraviolet detector standards is being currently investigated.
The dealloying corrosion behavior of the FCC Ni20Cr (wt%) in molten LiF-NaF-KF (FLiNaK) salts at 600 °C under varying applied potentials was investigated. Using in-operando electrochemical techniques and a multi-modal suite of characterization methods, we connect electrochemical potential, thermodynamic stability, and electro-dissolution kinetics to the corrosion morphologies. Notably, under certain potential regimes, a micron-scale bicontinuous structure, characterized by a network of interconnected pores and ligaments riched with the composition of the more noble (MN) element, becomes prominent. At other potentials both MN and less noble (LN) elements dealloy but at different rates. The dealloying process consists of lattice and grain boundary diffusion of Cr to the metal/salt interface, interphase Cr oxidation, accompanied by surface diffusion of Ni to form interconnected ligaments. At higher potentials, the bicontinuous porous structure undergoes further surface coarsening. Concurrently, Cr(II), Cr(III), and Ni(II) begin to dissolve, with the dissolution of Ni occurring at a significantly slower rate. When solid-state transport of Cr is exceeded by the interfacial rates, dealloying depths are limited.
Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.
Amphiphilic bottlebrush block copolymers (BBCPs), having a hydrophilic bottlebrush polymer (BP) linked covalently to a hydrophobic BP, were found to segregate to liquid-liquid interfaces to minimize the free energy of the system. The key parameter influencing the outcome of the experiments is the ratio between the degree of polymerization of the backbone (NBB) and that of the side-chain brushes (NSC). Specifically, a spherical, star-like configuration results when NBB < NSC, while a cylindrical, bottlebrush-like shape is preferred when NBB > NSC. Dynamic interfacial tension (γ) and fluorescence recovery after photobleaching (FRAP) measurements show that the BBCP configuration influences the areal density and in-plane diffusion at the fluid interface. The characteristic relaxation times associated with BBCP adsorption (τA) and reorganization (τR) were determined by fitting time-dependent interfacial tension measurements to a sum of two exponential relaxation functions. Both τA and τR initially increased with NBB up to 92 repeat units, due to the larger hydrodynamic radius in solution and slower in-plane diffusivity, attributed to a shorter cross-sectional diameter of the side-chains near the block junction. This trend reversed at NBB = 190, with shorter τA and τR attributed to increased segregation strength and exposure of the bare water/toluene interface due to tilting and/or wiggling of the backbone chains, respectively. The adsorption energy barrier decreased with higher NBB, due to a reduced BBCP packing density at the fluid interface. In conclusion, this study provides fundamental insights into macromolecular assembly at fluid interfaces, as it pertains to unique bottlebrush block architectures.
Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.
Numerical solutions of the governing diffusion equation for two-phase concentration dependent diffusion coefficients are examined. Solutions were also calculated for planar, cylindrical, and spherical geometries to compare the effect of interface geometries with those caused by concentration-dependent diffusion coefficients, and two methods of averaging D were considered to determine the best averaging method for different types of D-variations. The effects of interface-location criteria on mass conservation and convergence of interface location, diffusion coefficient variation in the alpha and beta-phases of a two-phase binary alloy system, effect of D(alpha) variation in a cylindrical couple on beta-phase thickness, and geometry and D-variation effects on the degree of homogenization were determined. It is concluded that typical D(alpha)-variations can have a greater influence on the kinetics of interdiffusion than the geometry.
The ability to fabricate additively manufactured laser waveguides with sharp dopant concentration interfaces is limited by diffusion of the dopants at the temperatures required to fully densify the material. Compositional analysis of bilayer samples, where each layer was either undoped YAG or YAG doped with Yb, Lu, or Nd, were fabricated such that all combinations were available for testing. Samples were fabricated at both 1750°C and 1850°C to determine the diffusion behavior of each dopant alone and also in the presence of a second dopant. It was found that the experimental concentration profiles exhibited both intragranular (bulk) and grain boundary contributions, and thus fitting to a complementary error function equation required the use of two diffusion coefficients. Nd always diffused further along grain boundaries than the other dopants, due in part to its small segregation coefficient in YAG. It is shown that the presence of Nd as a counter dopant inhibits intragranular diffusion of other dopants while enhancing their grain boundary diffusion. All the observed trends were attributed to a combination of intragranular lattice strain due to: dopant ions replacing yttrium substitutionally, the relative driving forces for the segregation of dopant ions to grain boundaries, and the ability of one dopant to affect the diffusion of a different dopant in the other direction to maintain charge neutrality.
Metal biogeochemistry in the sediment, water, and the sediment–water interface (SWI) was studied in a free water surface constructed wetland. Concentrations of labile copper (Cu), zinc (Zn), sulfate, chloride, and dissolved organic carbon (DOC) were measured with the diffusive gradients in thin film (DGT) and peeper. A good agreement between peeper- and DGT-measured metals was observed for Cu (regression r 2 = 0.3, 95% CI of the slopes > 0), but not for Zn (95% CI of the slopes overlapped with 0), which was attributed to the different complexed compounds between Cu and Zn in porewater. The depth profile of labile metals in sediment porewater varied with time and was consistent with the solid-phase metal deposition, showing higher concentrations in the surface layer (3 to –3 cm) than in the bottom layer (–4 to –13 cm). The depth-averaged labile Cu and Zn concentrations measured by DGT were 1.0 and 3.1 µg/L, and labile sulfate, chloride, and DOC concentrations measured by peeper were 1.8, 3.6, and 2.1 mg/L, respectively. A sharp decrease in sulfate occurred in September when sulfate concentrations became the lowest among sampling months. This was caused by the seasonal sulfur cycles in the wetland, where the dominant sulfur reaction is sulfate reduction in warm seasons and sulfide oxidation in cold seasons. Different metal-removal mechanisms were observed in the two wetland cells; sulfur dynamics controlled the removal processes in the cell without frequent disturbance but failed to influence metal removal in the cell with frequent disturbance due to the interruption of anoxic layers. The flux ratios that compare labile element concentrations between the water column and the SWI (R-Cu, R-Zn, R-DOC, R-sulfate, and R-chloride) were generated to determine metal diffusive fluxes at the interface. Labile Zn was mostly transported from the water to the SWI during all seasons (R-Zn < 1 for all months except January). Labile Cu moved from the SWI to the water during the warm months (R-Cu < 1), which was explained by the bioturbation-induced transport of organic matter based on the positive correlations between R-Cu and R-DOC. In general, sediment can serve either as a sink or a source depending on the environmental conditions, metal speciation, and presence of living organisms. Here, metal flux at the SWI is a key component in the biogeochemical cycling of a constructed wetland.
Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.