Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗

Concentration Dependent Interfacial Chemistry of the NaOH(aq):Gibbsite Interface

Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, gamma-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)—particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the gamma-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5–10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion–surface interactions. The molecular-scale changes to surface speciation and competition between ion–surface vs. ion–ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.

alkaline solutions, IDREAM, gibbsite, molecular dy↗

Growth of Megaspherulites In a Rhyolitic Vitrophyre

Megaspherulites occur in the middle zone of a thick sequence of rhyolitic vitrophyre that occupies a small, late Eocene to early Oligocene volcanic-tectonic basin near Silver Cliff, Custer County, Colorado. Diameters of the megaspherulites range from 0.3 m to over 3.66 m, including a clay envelope. The megaspherulites are compound spherulites. consisting of an extremely large number (3.8 x 10(exp 9) to 9.9 x 10(exp 9)) of individual growth cones averaging 3 mm long by 1.25 mm wide at their termination. They are holocrystalline, very fine- to fine-grained, composed of disordered to ordered sanidine (orthoclase) and quartz, and surrounded by a thin K-feldspar, quartz rich rind, an inner clay layer with mordenite, and an outer clay layer composed wholly of 15 A montmorillonite. Whole rock analyses of the megaspherulites show a restricted composition from their core to their outer edge, with an average analyses of 76.3% SiO2, 0.34% CaO, 2.17% Na2O, 6.92% K2O, 0.83% H2O+ compared to the rhyolitic vitrophyre from which they crystallize with 71.07% SiO2, 0.57% CaO, 4.06% Na2O,4.l0% K2O, and 6.40% H2O+. The remaining oxides of Fe2O3 (total Fe), A12O3, MnO,MgO, TiO2, P2O5, Cr2O3, and trace elements show uniform distribution between the megaspherulites and the rhyolitic vitrophyre. Megaspherulite crystallization began soon after the rhyolitic lava ceased to flow as the result of sparse heterogeneous nucleation, under nonequilibrium conditions, due to a high degree of undercooling, delta T. The crystals grow with a fibrous habit which is favored by a large delta T ranging between 245 C and 295 C, despite lowered viscosity, and enhanced diffusion due to the high H2O content, ranging between 5% and 7%. Therefore, megaspherulite growth proceeded in a diffusion controlled manner, where the diffusion, rate lags behind the crystal growth rate at the crystal-liquid interface, restricting fibril lengths and diameters to the 10 micron to 15 micron and 3 micron and 8 micron ranges respectively. Once diffusion reestablishes itself at the crystallization front, a new nucleation event occurs at the terminated tips of the fibril cones and a new cone begins to develop with a similar orientation (small angle branching) to the earlier cones. During crystallization, these fibril cones impinge upon each other, resulting in fibril cone-free areas. These cone-free areas consist of coarser, fine-grained phases, dominated by quartz, which crystallized from the melt as it accumulated between the crystallizing K-feldspar fibrils of the cones. The anhydrous nature of the disordered to ordered sanidine (orthoclase) and quartz, suggests that water in the vitrophyre moved ahead of the crystallization front, resulting in a water rich fluid being enriched in Si, K, Na, Mg, Ca, Sr, Ba, and Y. The clay layers associated with the megaspherulites are therefore, the result of the deuteric alteration between the fractionated water and the vitrophyre, as indicated by the presence of the minerals mordenite and montmorillonite. This silica rich fluid also resulted in the total silicification of the megaspherulites within the upper 3 m of the vitrophyre.

Smith, Robert K.↗

Anomalous Interlayer Exciton Diffusion in WS 2 /WSe 2 Moiré Heterostructure

Stacking van der Waals crystals allows for the on-demand creation of a periodic potential landscape to tailor the transport of quasiparticle excitations. We investigate the diffusion of photoexcited electron-hole pairs, or excitons, at the interface of WS 2 /WSe 2 van der Waals heterostructure over a wide range of temperatures. We observe the appearance of distinct interlayer excitons for parallel and antiparallel stacking and track their diffusion through spatially and temporally resolved photoluminescence spectroscopy from 30 to 250 K. While the measured exciton diffusivity decreases with temperature, it surprisingly plateaus below 90 K. Our observations cannot be explained by classical models like hopping in the moiré potential. A combination of ab initio theory and molecular dynamics simulations suggests that low-energy phonons arising from the mismatched lattices of moiré heterostructures, also known as phasons, play a key role in describing and understanding this anomalous behavior of exciton diffusion. Our observations indicate that the moiré potential landscape is dynamic down to very low temperatures and that the phason modes can enable efficient transport of energy in the form of excitons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Designing, Modeling, Constructing, and Testing a Flat Panel Speaker and Sound Diffuser for a Simulator

The goal of this project was to design, model, build, and test a flat panel speaker and frame for a spherical dome structure being made into a simulator. The simulator will be a test bed for evaluating an immersive environment for human interfaces. This project focused on the loud speakers and a sound diffuser for the dome. The rest of the team worked on an Ambisonics 3D sound system, video projection system, and multi-direction treadmill to create the most realistic scene possible. The main programs utilized in this project, were Pro-E and COMSOL. Pro-E was used for creating detailed figures for the fabrication of a frame that held a flat panel loud speaker. The loud speaker was made from a thin sheet of Plexiglas and 4 acoustic exciters. COMSOL, a multiphysics finite analysis simulator, was used to model and evaluate all stages of the loud speaker, frame, and sound diffuser. Acoustical testing measurements were utilized to create polar plots from the working prototype which were then compared to the COMSOL simulations to select the optimal design for the dome. The final goal of the project was to install the flat panel loud speaker design in addition to a sound diffuser on to the wall of the dome. After running tests in COMSOL on various speaker configurations, including a warped Plexiglas version, the optimal speaker design included a flat piece of Plexiglas with a rounded frame to match the curvature of the dome. Eight of these loud speakers will be mounted into an inch and a half of high performance acoustic insulation, or Thinsulate, that will cover the inside of the dome. The following technical paper discusses these projects and explains the engineering processes used, knowledge gained, and the projected future goals of this project

Dillon, Christina↗

Layer-resolved ultrafast extreme ultraviolet measurement of hole transport in a Ni-TiO 2 -Si photoanode

Metal oxide semiconductor junctions are central to most electronic and optoelectronic devices, but ultrafast measurements of carrier transport have been limited to device-average measurements. Here, charge transport and recombination kinetics in each layer of a Ni-TiO 2 -Si junction is measured using the element specificity of broadband extreme ultraviolet (XUV) ultrafast pulses. After silicon photoexcitation, holes are inferred to transport from Si to Ni ballistically in ~100 fs, resulting in characteristic spectral shifts in the XUV edges. Meanwhile, the electrons remain on Si. After picoseconds, the transient hole population on Ni is observed to back-diffuse through the TiO 2 , shifting the Ti spectrum to a higher oxidation state, followed by electron-hole recombination at the Si-TiO 2 interface and in the Si bulk. Electrical properties, such as the hole diffusion constant in TiO 2 and the initial hole mobility in Si, are fit from these transient spectra and match well with values reported previously.

42 ENGINEERING↗

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗

Real-Time KMC Simulation of Vacancy-Mediated Intermixing in Au@Ag Octahedral Core–Cubic Shell Nanocrystals with Ab Initio-Guided Kinetics

Utilization of core–shell rather than monometallic nanocrystals (NCs) facilitates fine-tuning of NC properties for applications. However, compositional evolution via intermixing can degrade these properties prompting recent experimental studies. We develop an atomistic-level stochastic model for vacancy-mediated intermixing exploiting a formalism which allows incorporation at an ab initio density functional theory level of not just the thermodynamics of vacancy formation, but also relevant diffusion barriers for a vast number of possible local environments (in the core and in the shell, at the interface, and in the intermixed phase). This facilitates a predictive treatment and comprehensive understanding of intermixing on the relevant time scale (e.g., 10 1 –10 3 s). In contrast, previous modeling at the atomistic level utilized only unrealistic generic prescriptions of barriers or employed simplified continuum treatments. For Au@Ag octahedral core–cubic shell NCs, our modeling not only captures the experimentally observed rate or time scale for intermixing of ~100 s at 450 °C for 60 nm NCs, but also elucidates the underlying rate controlling processes and the effective intermixing barrier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding coarsening of a post-corrosion microstructure in a molten salt by combining phase-field modeling and in situ tomography

Alloys corroding in molten salt have been observed to form bicontinuous, nanoporous microstructures via dealloying, which subsequently undergo coarsening due to facile transport in high-temperature conditions. In this work, we describe a methodology to elucidate the underlying transport mechanisms during coarsening of a bicontinuous microstructure via quantitative comparisons between phase-field simulations and four-dimensional in situ experiments, in this case X-ray nanotomography of the coarsening of a dealloyed 80 wt% Ni-20 wt% Cr microwire in molten KCl-MgCl 2 at 800°C. We conduct phase-field simulations initialized from experimental data to model coarsening via three different transport mechanisms: surface diffusion, solid bulk diffusion, and liquid bulk diffusion. These simulations reproduce key features of the experiment, such as the densification of the outer layer of the dealloyed wire and the reduction in radius over time. We quantitatively compare different microstructural characteristics between the simulations and experiment and extract temporal scaling factors that optimally match the time scales of the simulations to that of the experiment. This allows us to evaluate morphological similarity between the simulations and experiment and relate the experimental coarsening kinetics to fundamental material properties. We find that surface diffusion is most likely to be the dominant coarsening mechanism, and its kinetics imply a surface diffusivity of D S = 8.9 x 10 -20 m 3 /s, which is within the range of reported values for Ni-vacuum interfaces at 800°C. However, the difference between the experiment and the surface diffusion simulation increases substantially at late times, suggesting that other mechanisms, such as the dissolution of residual Cr, may be at play.

36 MATERIALS SCIENCE↗

Domain decomposition algorithms and computational fluid dynamics

Some of the new domain decomposition algorithms are applied to two model problems in computational fluid dynamics: the two-dimensional convection-diffusion problem and the incompressible driven cavity flow problem. First, a brief introduction to the various approaches of domain decomposition is given, and a survey of domain decomposition preconditioners for the operator on the interface separating the subdomains is then presented. For the convection-diffusion problem, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is examined.

Chan, Tony F.↗

Understanding LIB Battery Electrodes Through Classical Electrochemical Interface Theory

Classical models of the electrochemical interface can be applied to complex and dynamic electrodes to understand important mechanisms that are relevant to device-level phenomena like calendar aging or cycle life. Here, we show that negative alloy electrodes with wide electrochemical windows cannot be assumed to have static interfaces throughout all states of charge. Under high states of charge, the interface -including the solid electrolyte interphase - has characteristics that resemble a highly polarized electrode far from the point of zero charge. Under these conditions, the interface can be understood through the Helmholtz model. At lower states of charge, the interface capacitance decreases and the space charge layer length increases. This transformation resembles a metal electrode approaching the point of zero charge which is understood through the Gouy-Chapman-Stern-Grahame model of the electrochemical interface. This transformation happens at both silicon and carbon-coated silicon interfaces. From voltage-limited cycling experiments, in the limit of the diffuse double layer, the silicon electrode impedance rises significantly in the first few cycles. This rise is related to the continuous electrochemical reduction of electrolyte components from an interface that is not electronically screened from the electrolyte. In other words, an electrode interface for batteries should resemble a 'Helmholtz-like' structure.

anode↗

Precipitation-site competition in duplex stainless steels: Cu clusters vs spinodal decomposition interfaces as nucleation sites during thermal aging

Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.

Duplex stainless steel, Thermal aging, Precipitati↗

Sulfur transformations at the hydrogen sulfide/oxygen interface in stratified waters and in cyanobacterial mats

Stratified water bodies allow the development of several microbial plates along the water column. The microbial plates develop in relation to nutrient availability, light penetration, and the distribution of oxygen and sulfide. Sulfide is initially produced in the sediment by sulfate-reducing bacteria. It diffuses along the water column creating a zone of hydrogen sulfide/oxygen interface. In the chemocline of Solar Lake oxygen and sulfide coexist in a 0 to 10 cm layer that moves up and down during a diurnal cycle. The microbial plate at the chemocline is exposed to oxygen and hydrogen sulfide, alternating on a diurnal basis. The cyanobacteria occupying the interface switch from anoxygenic photosynthesis in the morning to oxygenic photosynthesis during the rest of the day which results in a temporal build up of elemental sulfur during the day and disappears at night due to both oxidation to thiosulfate and sulfate by thiobacilli, and reduction to hydrogen sulfide by Desulfuromonas sp. and anaerobically respiring cyanobacteria. Sulfate reduction was enhanced in the light at the surface of the cyanobacterial mats. Microsulfate reduction measurements showed enhanced activity of sulfate reduction even under high oxygen concentrations of 300 to 800 micrometer. Apparent aerobic SO sub 4 reduction activity is explained by the co-occurrence of H sub 2. The physiology of this apparent sulfate reduction activity is studied.

Cohen, Y.↗

Directional solidification in an imperfect furnace

The upward directional solidification of a binary alloy in which the rejected solute is more dense than the solvent is considered. The system is stably stratified both thermally and solutally. If heat is lost through the sidewalls, horizontal thermal gradients drive a steady double-diffusive convection. Solidification cells with small aspect ratios (width/length) A are considered, and asymptotic solutions to the steady double-diffusive system are obtained as A tends to 0 under the assumption heat losses are also small. When the thermal Grashof number is large, buoyancy boundary layers are at the sidewalls. The shape of the crystal-melt interface is determined by the coupled effects of local thermodynamic equilibrium, diffusion and convection.

Brattkus, K.↗

Solidification and crystal growth of solid solution semiconducting alloys

Problems associated with the solidification and crytal growth of solid-solution semiconducting alloy crystals in a terrestrial environment are described. A detailed description is given of the results for the growth of mercury cadmium telluride (HgCdTe) alloy crystals by directional solidification, because of their considerable technological importance. A series of HgCdTe alloy crystals are grown from pseudobinary melts by a vertical Bridgman method using a wide range of growth rates and thermal conditions. Precision measurements are performed to establish compositional profiles for the crystals. The compositional variations are related to compositional variations in the melts that can result from two-dimensional diffusion or density gradient driven flow effects ahead of the growth interface. These effects are discussed in terms of the alloy phase equilibrium properties, the recent high temperature thermophysical data for the alloys and the highly unusual heat transfer characteristics of the alloy/ampule/furnace system that may readily lead to double diffusive convective flows in a gravitational environment.

Lehoczky, S. L.↗

Thermodynamics and kinetics of core-shell versus appendage co-precipitation morphologies

What determines precipitate morphologies in co-precipitating alloy systems? We focus on alloys of two precipitating phases, with the fast-precipitating phase acting as heterogeneous nucleation sites for a second phase manifesting slower kinetics. Kinetic lattice Monte Carlo simulations show that the interplay between interfacial and ordering energies, plus active diffusion paths, strongly affect the selection of core-shell verses appendage morphologies. We study a FeCuMnNiSi alloy using the combination of atom probe tomography and simulations, and show that the ordering energy reduction of the MnNiSi phase heterogeneously nucleated on a pre-existing copper-rich precipitate exceeds the energy penalty of a predominantly Fe/Cu interface, leading to initial appendage, rather than core-shell, formation. Diffusion of Mn, Ni and Si around and through the Cu core towards the ordered phase results in subsequent appendage growth. We further show that in cases with higher primary precipitate interface energies and/or suppressed ordering, the coreshell morphology is favored.

Shu, Shipeng↗

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗