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.
Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.
Not Available
The understanding of the sorption/desorption kinetics is essential for practical applications of moisture-controlled CO 2 sorption. We introduce an analytic model of the kinetics of moisture-controlled CO 2 sorption and its interpretation in two limiting cases. In one case, chemical reaction kinetics on pore surfaces dominates, in the other case, diffusive transport through the sorbent defines the kinetics. Here, we show that reaction kinetics, which is dominant in the first case, can be expressed as a linear combination of 1st and 2nd order kinetics in agreement with the static isotherm equation derived and validated in a previous paper. The interior transport kinetics can be described by non-linear diffusion equations. By combining all carbon species into a single equation, we can eliminate — in certain limits — the source terms associated with chemical reactions. In this case, the governing equation is ∂θ/∂t = –∇ · (–D eff ∇ θ ). For a sorbent in a form of a flat sheet or a membrane, one can maintain the same functional form of a diffusion equation by introducing a generalized effective diffusivity D M that combines contributions from both surface chemical reaction kinetics and interior diffusive transport kinetics. Experimental data of transient CO 2 flux in a preconditioned commercial anion exchange membrane fit well to the 1st order model as long as very dry states are avoided, validating the theory. The observed DM for a preconditioned commercial anion exchange membrane ranges from 6.6× 10 -14 to 7.1× 10 -14 m 2 s -1 at 35°C. These small values compared to typical ionic diffusivities imply a very slow kinetics, which will be the largest issue that needs to be addressed for practical application. The collected transient CO 2 flux data are used to predict the magnitude of a continuous CO 2 pumping flux in an active membrane that transports CO 2 against a CO 2 concentration gradient. The pumped CO 2 flux is supported by water flux due to a water concentration gradient.
Sodium diffusion has been studied in p-type 4H-SiC. Heat treatments have been performed from 1200 °C to 1800 °C for 1 min to 4 h. Secondary ion mass spectrometry has been used to measure the sodium distribution. We show that sodium has a considerable mobility at 1200 °C in p-type 4H-SiC. On the other hand for sodium atoms trapped at suitable sites the mobility is limited up to 1800 °C. Trap limited diffusion kinetics is suggested and an effective diffusivity has been extracted with an activation energy of 4 eV for sodium diffusion in p-type 4H-SiC.
Standard second-order self-adjoint forms of the transport equation, such as the even-parity, odd-parity, and self-adjoint angular flux equation, cannot be used in voids. Perhaps more important, they experience numerical convergence difficulties in near-voids. Here we present a new form of a second-order self-adjoint transport equation that has an advantage relative to standard forms in that it can be used in voids or near-voids. Our equation is closely related to the standard least-squares form of the transport equation with both equations being applicable in a void and having a nonconservative analytic form. However, unlike the standard least-squares form of the transport equation, our least-squares equation is compatible with source iteration. It has been found that the standard least-squares form of the transport equation with a linear-continuous finite-element spatial discretization has difficulty in the thick diffusion limit. Here we extensively test the 1D slab-geometry version of our scheme with respect to void solutions, spatial convergence rate, and the intermediate and thick diffusion limits. We also define an effective diffusion synthetic acceleration scheme for our discretization. Our conclusion is that our least-squares S n formulation represents an excellent alternative to existing second-order S n transport formulations
The TRISO particle for HTRs has been developed to an advanced state where the coating withstands internal gas pressures and retains fission products during irradiation and under postulated accidents. However, one exception is Ag that has been found to be released from high quality TRISO coated particles when irradiated and can also during high temperature accident heating tests. Although out- of- pile laboratory tests have never hither to been able to demonstrate a diffusion process of Ag in SiC, effective diffusion coefficients have been derived to successfully reproduce measured Ag-110m releases from irradiated HTR fuel elements, compacts and TRISO particles It was found that silver transport through SiC does not proceed via bulk volume diffusion. Presently grain boundary diffusion that may be irradiation enhanced either by neutron bombardment or by the presence of fission products such as Pd, are being investigated. Recent studies of irradiated AGR-1 TRISO fuel using scanning transmission electron microscopy (STEM), transmission kukuchi diffraction (TKD) patterns and high resolution transmission electron microscopy (HRTEM) have been used to further the understanding of Ag transport through TRISO particles. No silver was observed in SiC grains, but Ag was identified at triple-points and grain boundaries of the SiC layer in the TRISO particle. Cadmium was also found in some of the very same triple junctions, but this could be related to silver behavior as Ag-110m decays to Cd-110. Palladium was identified as the main constituent of micron-sized precipitates present at the SiC grain boundaries and in most SiC grain boundaries and the potential role of Pd in the transport of Ag will be discussed.
Gas sorption and non-Darcy flow are two important issues for shale gas reservoirs. The sorption consists of dissolution and adsorption. Dissolved gas and adsorbed gas are different. The former is dissolved in the shale matrix, while the latter is concentrated near the solid walls of pores. In this paper, the Langmuir equation is used to describe adsorption and Henry’s law is used to describe dissolution. The K coefficient in Henry’s law of 0.052 mmol/(MPa g TOC) is obtained by matching experimental data. The amount of dissolved gas increases linearly when pressure increases. Using only the Langmuir equation without considering dissolution can lead to a significant underestimation of the amount of sorbed gas in shales. For non-Darcy gas flow, the apparent permeability model for free gas is established by combining slip flow and Knudsen flow. For adsorbed gas, the surface diffusion effect is also considered in this model. Here, the surface diffusion coefficient is suggested to be of the same scale as the gas self-diffusion coefficient, and the corresponding effective permeability is derived. When
Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. This study aims at calculating diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO 2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This paper estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 10 2 to 10 4 times as porosity decreases from 0.1 to 0.03.
Evaluation of the long-term retention mechanisms and potential release rates for the primary constituents of potential concern (COPCs) (i.e., Tc, I, Se, and nitrate) is necessary to determine if Cast Stone, a radioactive waste form, can meet performance objectives under near-surface disposal scenarios. In this work, a mineral and parameter set accounting for the solubility of I and Se in Cast Stone was developed based on pH-dependent and monolithic diffusion leaching test results, to extend a geochemical speciation model previously developed. The impact of oxidation and carbonation as environmental aging processes on the retention properties of Cast Stone for primary COPCs was systematically estimated. Physically, the effective diffusion coefficients of 4 COPCs in Cast Stone were increased after carbonation and/or oxidation, reflecting an increase in permeability to diffusion. Chemically, i) pH & pe conditions in the original Cast Stone were favorable for the stabilization of Tc, but not for I, Se, and N; ii) oxidation (with/without carbonation) of Cast Stone changed the pe & pH conditions to be detrimental for Tc stabilization; and iii) carbonation (with/without oxidation) of Cast Stone modified the pH & pe conditions to be beneficial for the stabilization of I (in system with Ag added) and Se.
Diffusive gradient in thin-film (DGT) passive samplers are frequently used to monitor the concentrations of metals such as mercury and zinc in sediments and other aquatic environments. The application of these samplers generally presumes that they quantify only the dissolved fraction and not particle-bound metal species that are too large to migrate into the sampler. However, metals associated with very small nanoparticles (smaller than the pore size of DGT samplers) can be abundant in certain environments, yet the implications of these nanoparticles for DGT measurements are unclear. The objective of this study was to determine how the performance of the DGT sampler is affected by the presence of nanoparticulate species of Hg and Zn. DGT samplers were exposed to solutions containing known amounts of dissolved Hg(II) and nanoparticulate HgS (or dissolved Zn(II) and nanoparticulate ZnS). The amounts of Hg and Zn accumulated onto the DGT samplers were quantified over hours to days, and the rates of diffusion of the dissolved metal (i.e., the effective diffusion coefficient D) into the sampler’s diffusion layer were calculated and compared for solutions containing varying concentrations of nanoparticles. The results suggested that the nanoparticles deposited on the surface of the samplers might have acted as sorbents, slowing the migration of the dissolved species into the samplers. The consequence was that the DGT sampler data underestimated the dissolved metal concentration in the solution. In addition, X-ray absorption spectroscopy was employed to determine the speciation of the Hg accumulated on the sampler binding layer, and the results indicated that HgS nanoparticles did not appear to directly contribute to the DGT measurement. Overall, our findings suggest that the deployment of DGT samplers in settings where nanoparticles are relevant (e.g., sediments) may result in DGT data that incorrectly estimated the dissolved metal concentrations. Models for metal uptake into the sampler may need to be reconsidered.
Iron-chromium-aluminum (FeCrAl) alloys such as APMT (advanced powder metallurgy tubing) are candidate materials to replace zirconium alloys for the light water reactor fuel cladding. This alloy meets the requirements to be a material more tolerant of high-temperature accidents than the current zirconium alloys. One concern is that the use of FeCrAl may result in an increase in tritium presence in the coolant compared to the current design. The aim of the current research was to obtain effective diffusion coefficients (D eff ) for hydrogen through APMT using the Devanathan-Stachurski cell. Results showed that at 30°C the D eff value was 2.8 × 10 -8 cm 2 /s. Results also showed a linear relationship between the permeated hydrogen flux and the inverse of the test specimen thickness, which demonstrated the validity of the permeation measurements.
In this paper, a one-dimensional computational framework is developed that can solve for the evolution of voltage and current in a lithium-ion battery electrode under different operating conditions. A reduced order model is specifically constructed to predict the growth of mechanical degradation within the active particles of the carbon anode as a function of particle size and C-rate. Using an effective diffusivity relation, the impact of microcracks on the diffusivity of the active particles has been captured. Reduction in capacity due to formation of microcracks within the negative electrode under different operating conditions (constant current discharge and constant current constant voltage charge) has been investigated. At the beginning of constant current discharge, mechanical damage to electrode particles predominantly occurs near the separator. As the reaction front shifts, mechanical damage spreads across the thickness of the negative electrode and becomes relatively uniform under multiple discharge/charge cycles. Mechanical degradation under different drive cycle conditions has been explored. It is observed that electrodes with larger particle sizes are prone to capacity fade due to microcrack formation. Finally, under drive cycle conditions, small particles close to the separator and large particles close to the current collector can help in reducing the capacity fade due to mechanical degradation.
This report summarizes the development of new fission gas diffusion models from lower length scale simulations and assessment of these models in terms of annealing experiments and fission gas release simulations using the BISON fuel performance code. Based on the mechanisms established from density functional theory (DFT) and empirical potential calculations, continuum models for diffusion of xenon (Xe) in UO 2 were derived for both intrinsic conditions and under irradiation. The importance of the large X eU 3 O cluster (a Xe atom in a uranium + oxygen vacancy trap site with two bound uranium vacancies) is emphasized, which is a consequence of its high mobility and stability. These models were implemented in the MARMOT phase field code, which is used to calculate effective Xe diffusivities for various irradiation conditions. The effective diffusivities were used in BISON to calculate fission gas release for a number of test cases. The results are assessed against experimental data and future directions for research are outlined based on the conclusions.
Knowledge of the exchange of carbon isotopes between dissolved inorganic carbon and calcite minerals is of long-standing importance for the interpretation of sedimentary paleoclimate records and 14 C transport in the geosphere. To assess the mechanism and rates of carbon isotope exchange, we equilibrated 12 C-pure synthetic calcite particles in water, first in a glovebox and then in contact with atmospheric P CO 2 and 13 C/ 12 C ratio, at two different temperatures. Cavity ring-down infrared spectroscopy δ 13 C measurements of the solid revealed sustained 13C incorporation for over a period of 500 days (21 °C) and 125 days (50 °C). We developed a quantitative model for recrystallization and isotope exchange, assuming that the interfacial free energy provides a thermodynamic driving force for the growth of larger particles at the expense of smaller ones. Furthermore, this Ostwald ripening model did not reproduce the kinetics of 13C uptake and required greater coarsening than observed. Rather, the data were best explained by a mechanism involving surface exchange and solid-phase diffusion of 13C into the particles with an inferred effective diffusion constant at 21 °C of about 10 –25 m 2 /s. Although this work cannot rule out the possibility that structural or chemical aspects of the synthetic particles enabled faster 13 C uptake than could be observed in natural systems, this study adds to the body of the recent work, suggesting that fast exchange processes are possible, likely through grain boundaries and other defects.
In this paper, we describe our work on spin polarization in high-energy electron storage rings which we base on the Full Bloch equation (FBE) for the polarization density and which aims towards the [Formula: see text] option of the proposed Future Circular Collider (FCC-ee) and the proposed Circular Electron Positron Collider (CEPC). The FBE takes into account non spin-flip and spin-flip effects due to synchrotron radiation including the spin-diffusion effects and the Sokolov–Ternov effect with its Baier–Katkov generalization as well as the kinetic-polarization effect. This mathematical model is an alternative to the standard mathematical model based on the Derbenev–Kondratenko formulas. For our numerical and analytical studies of the FBE, we develop an approximation to the latter to obtain an effective FBE. This is accomplished by finding a third mathematical model based on a system of stochastic differential equations (SDEs) underlying the FBE and by approximating that system via the method of averaging from perturbative ODE theory. We also give an overview of our algorithm for numerically integrating the effective FBE. This discretizes the phase space using spectral methods and discretizes time via the additive Runge–Kutta (ARK) method which is a high-order semi-implicit method. We also discuss the relevance of the third mathematical model for spin tracking.
Generation IV reactors and future fusion reactor designs have led to more demanding materials performance requirements due to their increased operating temperatures, corrosive coolants, and increased radiation doses compared to the current light-water reactor fleet. Among the innovative nuclear technologies under development, molten salt reactors stand out for their potential to offer superior fuel utilization, intrinsic safety characteristics, and economic viability. Of the proposed Generation IV designs, the gas fast reactor operates at 450 to 850°C and the molten salt reactor operates at 565 to 850°C, with the molten salt reactor design needing molten salt corrosion resistant materials [1, 2]. These increased temperatures and more extreme corrosion environments necessitate higher material performance, such as creep strength, radiation-tolerant microstructures, corrosion resistance, and high-temperature tensile properties. Hastelloy-N, a nickel-based alloy with additions of molybdenum and chromium, has been successfully employed to contain molten fluoride salt at temperatures up to 705°C. However, Hastelloy-N becomes embrittled upon neutron irradiation, primarily due to the accumulation of helium produced by (n,a) transmutation reactions. Furthermore, the corrosive nature of molten fluoride and chloride salts presents a formidable challenge, as these salts can react with and dissolve alloying elements such as Cr, Mo, and Fe, leading to selective leaching, loss of protective oxide layers, and accelerated degradation. High entropy alloys (HEAs) and refractory high entropy alloys (RHEAs) have emerged as a prominent area of interest, due to their ability to achieve tailored chemical compositions for specific applications. Unlike conventional alloys, HEAs are characterized by having multiple principal elements in equimolar or near equimolar ratios, leading to an unconventional alloying strategy [3]. This alloying strategy is believed to promote unique properties, such as single-phase stabilization of chemically compatible elements, lattice distortion effects due to atomic radius differences, and proposed sluggish diffusion effects. For extreme-environment applications, RHEAs have garnered much research interest because of the possibility of creating relatively ductile materials that can operate in extremely high-temperature environments, beyond the operating temperatures where other Ni-based superalloys begin to lose strength [4-6]. Idaho National Laboratory (INL) initiated a joint international effort with the Czech Republic to explore the feasibility of manufacturing HEAs for high-temperature nuclear applications using advanced manufacturing. This effort was funded at INL by the United States Department of Energy's Office of Nuclear Energy under the Advanced Reactor Technologies and Advanced Materials and Manufacturing Technologies (AMMT) Program. The HEAs were specifically designed for the corrosive and irradiation environments experienced in gas-cooled fast reactors, molten salt reactors, and fusion power. These alloys have been manufactured by multiple processes to determine the impact of manufacturing processes on the performance of the alloys in corrosive and irradiation environments. Preliminary molten salt corrosion testing showed that equimolar MoNbTiV and MoNbTi alloys exhibit exceptional performance, with arc-melted variants demonstrating only minimal degradation after 1000 hours of exposure to molten chloride salt at 700°C. Conversely, Nb2TiVZr2 showed significant molten salt corrosion susceptibility and microstructural instability during high-temperature molten salt exposures, and was, therefore deemed unfit for molten salt reactor applications. The MoNbTiV, MoNbTi, and Nb2TiVZr2 alloys were further evaluated through ion irradiation experiments conducted at the Michigan Ion Beam Laboratory at the University of Michigan. The microstructural stability and the evolution of irradiation-induced defects were characterized to assess the irradiation resistance of each of these alloys.
We investigate wave propagation in rotationally symmetric tubes with a periodic spatial modulation of cross section. Using an asymptotic perturbation analysis, the governing quasi-two-dimensional reaction-diffusion equation can be reduced into a one-dimensional reaction-diffusion-advection equation. Assuming a weak perturbation by the advection term and using projection method, in a second step, an equation of motion for traveling waves within such tubes can be derived. Both methods predict properly the nonlinear dependence of the propagation velocity on the ratio of the modulation period of the geometry to the intrinsic width of the front, or pulse. As a main feature, we observe finite intervals of propagation failure of waves induced by the tube’s modulation and derive an analytically tractable condition for their occurrence. For the highly diffusive limit, using the Fick-Jacobs approach, we show that wave velocities within modulated tubes are governed by an effective diffusion coefficient. Furthermore, we discuss the effects of a single bottleneck on the period of pulse trains. We observe period changes by integer fractions dependent on the bottleneck width and the period of the entering pulse train.
We propose a mathematical model for uptake and diffusion of air plasma-generated reactive nitrogen species (RNS) in a model keratinized membrane, such as a thin slice of bovine hoof or human nail material. An experimental system consisting of a surface microdischarge (SMD) in air was designed for the purpose of developing and validating a reaction-diffusion model to describe this system. Key variables such as membrane effective diffusivity, surface reaction rate coefficients, and other parameters are determined through comparison between the model predictions and experimental measurements. The model results yield spatial and temporal concentration profiles of RNS inside the keratinized membrane, leading to an improved understanding of transport and reaction of plasma generated RNS through the membrane. This work offers insights into possible mechanisms underlying plasma treatment of toenail fungus.