SEARCH · Engineering Papers
Results for “effective diffusivity”
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.
Effects of irradiation on the interface between U-Mo and zirconium diffusion barrier
Irradiated fuel plates were characterized by microscopy that focused on the interface between U–Mo and Zr. Before irradiation, there were three major sub-layers identified in the U–Mo/Zr interface, namely, UZr 2 , Mo 2 Zr, and U with low Mo. The typical total thickness of this U–Mo/Zr interaction is 2–3 microns. The UZr 2 sub-layer formed during fuel fabrication remains stable after irradiation, without large bubbles/porosity accumulation. However, this sub-layer becomes increasingly discontinuous as burnup increases. The low-Mo sub-layer exhibits numerous sub-micron bubbles/porosity at low burnup. Larger, interconnected porosity in this sub-layer was observed in a medium-burnup fuel specimen. However, at higher burnup, regions with the extra-large bubbles/porosity (i.e., larger than 5 ?m) were observed at least 5 microns away from the original location of this sub-layer. The mechanism for the formation of the extra-large bubbles/porosity is still unclear at this time. In general, the U–Mo/Zr interface in monolithic U–Mo fuels is relatively stable after irradiation. As a result, no large detrimental defects, such as large interfacial bubbles or cracks/delamination, were observed in the fuel plates characterized.
Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)
This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The DSA scheme has been implemented in PROTEUS to speed up SN transport calculations. The existing DSA scheme employs a consistent finite element formulation on the same fine mesh structure as the SN solution While the use of consistent discretizations makes the DSA method effective, it does not necessarily make it efficient because the numerical solution of a discretized elliptic diffusion problem itself can be costly, particularly when the problem size becomes large. We developed a new discontinuous Galerkin (DG) discretization of the diffusion equation, called DG-DSA, which can effectively and efficiently accelerate the SN transport iterations. As compared with the previous work, the novelty of our method is that the diffusion equation is solved on a coarse-mesh grid using the DG methods, and the DG diffusion discretization incorporates local hp adaptation, i.e., local adaptation of mesh size and/or polynomial degree, based on local total cross section (or optical thickness). Therefore, the resulting number of degrees of freedom (DOF) of the DG discretization is much less than the conventional consistent DSA discretizations, and thus DG-DSA can achieve significant improvement in computational efficiency. We implemented this scheme in PROTEUS-SN. In addition, we developed a new stabilization scheme named linear prolongation CMFD (lpCMFD). A novel feature of this scheme is that the conventional flat flux ratio–based scaling approach is replaced with a linear interpolation of the scalar flux differences at the coarse-mesh cell edges between the neutron transport and CMFD calculations. Fourier convergence analysis and numerical results show that lpCMFD is unconditionally stable and more effective than other CMFD based on acceleration schemes such as pCMFD, odCMFD. We also developed a new nonlinear diffusion acceleration scheme for solving neutron transport equations. This scheme, called LR-NDA, employs a local refinement approach on the framework of CMFD by solving a local boundary value problem of the scalar flux on the coarse-mesh structure to replace the piecewise constant scalar flux obtained by CMFD. The refined flux is then used to update the scalar flux in the neutron transport source iteration. It has been demonstrated that that LR-NDA is much more effective and stable than CMFD for a wide range of optical thicknesses. LR-NDA is a local adaptive method, which means LR-NDA does not necessarily require local refinement for all the coarse-mesh cells on the problem domain, i.e., it can be used only for relatively optically thick regions where the standard CMFD scheme would encounter the convergence problem.
A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination
Here in this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.
Immobile Pore-Water Storage Enhancement and Retardation of Gas Transport in Fractured Rock
The effect of immobile pore-water on gas transport in fractured rock has implications for numerical modeling of soil vapor extraction, methane leakage, gaseous CO 2 leakage from sequestration operations, radionuclide transport from underground nuclear explosions, and nuclear waste disposal. While the ability for immobile pore-water storage to effect gas transport has been recognized in the past, the details and specific scenarios leading to enhanced, retarded, or unaffected gas transport have not been explored. In this paper, we performed numerical investigations into the enhancement and retardation of gas transport due to immobile pore-water storage in order to identify implications for gas transport applications. To do this, we developed a numerical approach to model gas transport with a single-phase flow solution coupled to the advection–dispersion equation modified to account for immobile pore-water storage. Other than the immobility of pore water, the formulation contains all other physics included in two-phase formulations (advective and diffusive gas transport in fractures and rock matrix and dissolution in immobile pore water). The assumption of immobile pore water is valid here since for many applications involving transport of soluble gases in fractured rock, the rate of aqueous transport is insignificant compared to gas transport. We verify our modeling approach with analytical solutions of: (1) 1D gas diffusion, (2) 1D gas advection, (3) barometric pumping of a fracture, and (4) gas transport with uniform fracture flow and diffusion into the matrix. We account for pore-water storage in our model by implementing a kinetic formulation of gas dissolution wherein the dissolved (aqueous) phase is considered an immobile constituent. Using this formulation, we model the effect of dissolution rate and saturation on the retardation of gas transport during pure diffusion and pure advection. We also demonstrate that although it is commonly believed that pore-water storage will always enhance gas transport in fractures during oscillatory flow (e.g., during reversing pressure gradients such as barometric pumping cycles), our simulations indicate that this may not always be the case. Our numerical investigations indicate that scenarios with lower effective diffusion coefficients ( ⪅<!-- <≈ --> 10 −<!-- - --> 5  m 2 / s ) and lower dissolution coefficients (i.e., the dissolution diffusion coefficient) ( ⪅<!-- <≈ --> 10 −<!-- - --> 11  m 2 / s ) will result in enhanced gas transport. Other combinations of gaseous diffusion and dissolution coefficients result in delayed gas transport or insignificant effect on gas transport. Based on these results, and given the fact that the free-air diffusion coefficients for many gases are slightly above 10 −<!-- - --> 5  m 2 / s (near the boundary of enhanced/non-enhanced gas transport), tortuosity would have to be significant for gas transport enhancement to occur. Similarly, if we consider the free-water diffusion coefficient of gases (around 10 −<!-- - --> 9  m 2 / s for many gases of interest) to be a maximum bound for the effective dissolution coefficients (i.e., the rate of diffusion into the water away from the air/water interface will limit the dissolution rate), the effective dissolution coefficient would have to be at least 2 orders of magnitude less for enhanced gas transport to occur. Otherwise, pore-water storage will delay or have negligible effect on the gas transport. The results also indicate that the assumption of instantaneous equilibrium, often invoked in numerical codes for dissolution processes (i.e., dissolution coefficient is effectively infinite), will fail to accurately capture important details of soluble gas transport. Finally, the implication of our results is that gas breakthrough times in fractured rock during oscillatory flow may be highly dependent on the effective gaseous diffusion and dissolution coefficients.
Short-range thermal magnon diffusion in magnetic garnet
Here, using the spin Seebeck effect (SSE), we study the propagation distance of thermally induced spin currents inside a magnetic insulator thin film in the short-range regime. We disambiguate spin currents driven by temperature and chemical potential gradients by comparing the SSE signal before and after adding a heat-sinking capping layer on the same device. We report that the measured spin decay behavior near the heat source is well accounted for by a diffusion model where the magnon diffusion length is in submicron range, in other words, two orders of magnitude smaller than previous estimates inferred from the long-range behavior. Our results highlight the caveat in applying a diffusive theory to describe thermally generated magnon transport, where a single decay length may not capture the behavior on all length scales.
Size dependent behavior of Fe 3 O 4 crystals during electrochemical (de)lithiation: an in situ X-ray diffraction, ex situ X-ray absorption spectroscopy, transmission electron microscopy and theoretical investigation
Here, the iron oxide magnetite, Fe 3 O 4 , is a promising conversion type lithium ion battery anode material due to its high natural abundance, low cost and high theoretical capacity. While the close packing of ions in the inverse spinel structure of Fe 3 O 4 enables high energy density, it also limits the kinetics of lithium ion diffusion in the material. Nanosizing of Fe 3 O 4 to reduce the diffusion path length is an effective strategy for overcoming this issue and results in improved rate capability. However, the impact of nanosizing on the multiple structural transformations that occur during the electrochemical (de)lithiation reaction in Fe 3 O 4 is poorly understood. In this study, the influence of crystallite size on the lithiation-conversion mechanisms in Fe 3 O 4 is investigated using complementary X-ray techniques along with transmission electron microscopy (TEM) and continuum level simulations on electrodes of two different Fe 3 O 4 crystallite sizes. In situ X-ray diffraction (XRD) measurements were utilized to track the changes to the crystalline phases during (de)lithiation. X-ray absorption spectroscopy (XAS) measurements at multiple points during the (de)lithiation processes provided local electronic and atomic structural information. Tracking the crystalline and nanocrystalline phases during the first (de)lithiation provides experimental evidence that (1) the lithiation mechanism is non-uniform and dependent on crystallite size, where increased Li + diffusion length in larger crystals results in conversion to Fe 0 metal while insertion of Li + into spinel-Fe 3 O 4 is still occurring, and (2) the disorder and size of the Fe metal domains formed when either material is fully lithiated impacts the homogeneity of the FeO phase formed during the subsequent delithiation.
Electric current induced flow of liquid metals: Mechanism and substrate-surface effects
Long range, continuous flow of liquid metals occurs upon application of an electric current. Here, we report experimental results elucidating the mechanism of current-induced liquid metal flow, and its dependence on substrate surface condition. It is shown that the observed flow is diffusion-controlled, with the flow-rate depending linearly on applied current density, indicating that it is driven by electromigration. The effective charge number for liquid electromigration, Z*, of several pure metals, such as Al, Bi, Ga, Sn, and Pb, were deduced from the experimental results and were found to be close to the elemental valency. With the exception of liquid Pb, Z* for all liquid metals tested in this study were positive, indicating that: (i) electron wind contributes much less to Z* in liquid metals than in solids, and (ii) with a few exceptions, liquid metals generally flow in the direction of the electric current. On smooth substrates which are wetted well by the liquid metal, flow occurs in a thin, continuous stream. On rough surfaces which are poorly wetted, on the other hand, discrete beads of liquid form, with mass transport between adjacent beads occurring by surface diffusion on the substrate. A rationale for the role of substrate roughness in fostering this observed transition in flow mechanism is presented.
DISEQUILIBRIUM CARBON, OXYGEN, AND NITROGEN CHEMISTRY IN THE ATMOSPHERES OF HD 189733b AND HD 209458b
We have developed a one-dimensional photochemical and thermochemical kinetics and diffusion model to study the effects of disequilibrium chemistry on the atmospheric composition of 'hot-Jupiter' exoplanets. Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species on HD 189733b and HD 209458b and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on the cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than an mbar due to transport-induced quenching, but CH{sub 4} and NH{sub 3} are photochemically removed at higher altitudes. Disequilibrium chemistry also enhances atomic species, unsaturated hydrocarbons (particularly C{sub 2}H{sub 2}), some nitriles (particularly HCN), and radicals like OH, CH{sub 3}, and NH{sub 2}. In contrast, CO, H{sub 2}O, N{sub 2}, and CO{sub 2} more closely follow their equilibrium profiles, except at pressures {approx}<1 {mu}bar, where CO, H{sub 2}O, and N{sub 2} are photochemically destroyed and CO{sub 2} is produced before its eventual high-altitude destruction. The enhanced abundances of CH{sub 4}, NH{sub 3}, and HCN are expected to affect the spectral signatures and thermal profiles of HD 189733b and other relatively cool, transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficients and discuss further observational consequences of these models.
OBJECT KINETIC MONTE CARLO SIMULATIONS OF RADIATION DAMAGE IN BULK TUNGSTEN
We used our recently developed lattice based OKMC code; KSOME [1] to carryout simulations of radiation damage in bulk W. We study the effect of dimensionality of self interstitial atom (SIA) diffusion i.e. 1D versus 3D on the defect accumulation during irradiation with a primary knock-on atom (PKA) energy of 100 keV at 300 K for the dose rates of 10-5 and 10-6 dpa/s. As expected 3D SIA diffusion significantly reduces damage accumulation due to increased probability of recombination events. In addition, dose rate, over the limited range examined here, appears to have no effect in both cases of SIA diffusion.
Adsorption and Diffusion of Fructose in Zeolite HZSM-5: Selection of Models and Methods for Computational Studies
Explore the source record for details and available documents.
Transition state and product diffusion control by polymer–nanocrystal hybrid catalysts
Effective catalysts stabilize specific transition states and control the transport of species to and from catalytically active sites. Enzymes show these traits thanks to their diverse amino acid functional groups encapsulating metal centres, but are limited in the reaction conditions in which they can operate. Realizing a catalyst with this kinetic and transport control that can be used under demanding industrial conditions is challenging. Here, we show a modular approach for the systematic synthesis of polymer–nanocrystal hybrids, where palladium nanocrystals are encapsulated within tunable microporous polymer layers. The polymer chemistry and morphology control the catalytic performance of the metal sites, affecting the transition state for CO oxidation and controlling the transport of CO 2 away from the active site. In conclusion, this approach can be applied to other polymer–nanocrystal compositions and catalytic applications, and is therefore expected to have an impact in many areas of catalysis.
Interaction transfer of silicon atoms forming Co silicide for Co/√(3)×√(3)R30°-Ag/Si(111) and related magnetic properties
Combined scanning tunneling microscopy, Auger electron spectroscopy, and surface magneto-optic Kerr effect studies were employed to study the microscopic structures and magnetic properties for ultrathin Co/√(3)×√(3)R30°-Ag/Si(111). As the annealing temperature increases, the upward diffusion of Si atoms and formation of Co silicides occurs at temperature above 400 K. Below 600 K, the √(3)×√(3)R30°-Ag/Si(111) surface structure persists. We propose an interaction transferring mechanism of Si atoms across the √(3)×√(3)R30°-Ag layer. The upward transferred Si atoms react with Co atoms to form Co silicide. The step height across the edge of the island, a separation of 0.75 nm from the analysis of the 2 × 2 structure, and the calculations of the normalized Auger signal serve as strong evidences for the formation of CoSi{sub 2} at the interface. The interaction transferring mechanism for Si atoms enhances the possibility of interactions between Co and Si atoms. The smoothness of the surface is advantage for that the easy axis of magnetization for Co/√(3)×√(3)R30°-Ag/Si(111) is in the surface plane. This provides a possible way of growing flat magnetic layers on silicon substrate with controllable silicide formation and shows potential applications in spintronics devices.
Electrochemical Electron Transfer and Proton-Coupled Electron Transfer: Effects of Double Layer and Ionic Environment on Solvent Reorganization Energies
Electron transfer and proton coupled electron transfer (PCET) reactions at electrochemical interfaces play an essential role in a broad range of energy conversion processes. The reorganization energy, which is a measure of the free energy change associated with solute and solvent rearrangements, is a key quantity for calculating rate constants for these reactions. We present a computational method for including the effects of the double layer and ionic environment of the diffuse layer in calculations of electrochemical solvent reorganization energies. This approach incorporates an accurate electronic charge distribution of the solute within a molecular-shaped cavity in conjunction with a dielectric continuum treatment of the solvent, ions, and electrode using the integral equations formalism polarizable continuum model. The molecule-solvent boundary is treated explicitly, but the effects of the electrode-double layer and double layer-diffuse layer boundaries, as well as the effects of the ionic strength of the solvent, are included through an external Green’s function. The calculated total reorganization energies agree well with experimentally measured values for a series of electrochemical systems, and the effects of including both the double layer and ionic environment are found to be very small. This general approach was also extended to electrochemical PCET and produced total reorganization energies in close agreement with experimental values for two experimentally studied PCET systems. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center, funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences.
Rutile solubility in NaF–NaCl–KCl-bearing aqueous fluids at 0.5–2.79GPa and 250–650°C
Abstract not provided
Influence of Hydrogen Bonding on the Surface Diffusion of Molecular Glasses: Comparison of Three Triazines
Surface grating decay measurements have been performed on three closely related molecular glasses to study the effect of intermolecular hydrogen bonds on surface diffusion. The three molecules are derivatives of bis(3,5-dimethyl-phenylamino)-1,3,5-triazine and differ only in the functional group R at the 2-position, with R being C 2 H 5 , OCH 3 , and NHCH 3 , and referred to as “Et”, “OMe”, and “NHMe”, respectively. Of the three molecules, NHMe forms more extensive intermolecular hydrogen bonds than Et and OMe and was found to have slower surface diffusion. For Et and OMe, surface diffusion is so fast that it replaces viscous flow as the mechanism of surface grating decay as temperature is lowered. In contrast, no such transition was observed for NHMe under the same conditions, indicating significantly slower surface diffusion. This result is consistent with the previous finding that extensive intermolecular hydrogen bonds slow down surface diffusion in molecular glasses and is attributed to the persistence of hydrogen bonds even in the surface environment. Here, this result is also consistent with the lower stability of the vapor-deposited glass of NHMe relative to those of Et and OMe and supports the view that surface mobility controls the stability of vapor-deposited glasses.
First-principles studies on vacancy-modified interstitial diffusion mechanism of oxygen in nickel, associated with large-scale atomic simulation techniques
This paper is concerned with the prediction of oxygen diffusivities in fcc nickel from first-principles calculations and large-scale atomic simulations. Considering only the interstitial octahedral to tetrahedral to octahedral minimum energy pathway for oxygen diffusion in fcc lattice, greatly underestimates the migration barrier and overestimates the diffusivities by several orders of magnitude. The results indicate that vacancies in the Ni-lattice significantly impact the migration barrier of oxygen in nickel. Incorporation of the effect of vacancies results in predicted diffusivities consistent with available experimental data. First-principles calculations show that at high temperatures the vacancy concentration is comparable to the oxygen solubility, and there is a strong binding energy and a redistribution of charge density between the oxygen atom and vacancy. Consequently, there is a strong attraction between the oxygen and vacancy in the Ni lattice, which impacts diffusion.
Outcomes of Diffusion Tensor Tractography-Integrated Stereotactic Radiosurgery
Purpose: To analyze the effect of use of tractography of the critical brain white matter fibers created from diffusion tensor magnetic resonance imaging on reduction of morbidity associated with radiosurgery. Methods and Materials: Tractography of the pyramidal tract has been integrated since February 2004 if lesions are adjacent to it, the optic radiation since May 2006, and the arcuate fasciculus since October 2007. By visually confirming the precise location of these fibers, the dose to these fiber tracts was optimized. One hundred forty-four consecutive patients with cerebral arteriovenous malformations who underwent radiosurgery with this technique between February 2004 and December 2009 were analyzed. Results: Tractography was prospectively integrated in 71 of 155 treatments for 144 patients. The pyramidal tract was visualized in 45, the optic radiation in 22, and the arcuate fasciculus in 13 (two tracts in 9). During the follow-up period of 3 to 72 months (median, 23 months) after the procedure, 1 patient showed permanent worsening of pre-existing dysesthesia, and another patient exhibited mild transient hemiparesis 12 months later but fully recovered after oral administration of corticosteroid agents. Two patients had transient speech disturbance before starting integration of the arcuate fasciculus tractography, but no patient thereafter. Conclusion: Integrating tractography helped prevent morbidity of radiosurgery in patients with brain arteriovenous malformations.