Engineering Papers⌕ Search

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.

At least 253 records · Page 14

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗

Linear Stability Analysis of Gravitational Effects on a Low-Density Gas Jet Injected into a High-Density Medium

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The Briggs-Bers criterion was combined with the spatio-temporal stability analysis to determine the nature of the absolute instability of the jet whether absolutely or convectively unstable. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the absolute instability of the jet were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be.

Lawson, Anthony L.↗

Investigation of ion irradiation effects on mineral analogues of concrete aggregates

Irradiation can cause prominent damage to reactor concrete aggregates leading to amorphization, strength and modulus decrease, radiation induced volume expansion (RIVE) and micro-cracking, which limits their long-term performance. Here, to develop an improved understanding of irradiation effects in concrete, three mineral analogues of concrete aggregates (limestone, marble and quartzite) were irradiated by 5.5 MeV He ions and 13 MeV Ni ions to surface doses of 0.011 displacements per atom (dpa) and 0.23 dpa, respectively, at room temperature. The two different ion species allow irradiation spectrum effects (ionizing and displacive) to be examined. Irradiation induced cracks were observed in He irradiated limestone and marble, and Ni irradiated quartzite. Full amorphization was observed in Ni irradiated quartzite with 14.3 % RIVE, and ∼25 % hardness and modulus decrease, while almost no change was observed in He irradiated quartzite except 4.35 % RIVE, revealing a possible ionization enhanced diffusion effect for high energy light ions. Furthermore, partial amorphization was observed in Ni irradiated marble and limestone matrix with a 12 % hardness decrease in marble while no amorphization was observed for He irradiation with a 20 % hardness increase in limestone matrix. The role of knock-on damage and irradiation spectrum on amorphization, volumetric expansion and mechanical property changes are discussed. Moreover, the onset and critical doses for amorphization and RIVE in quartz are obtained for ion irradiations at room temperature. The dose dependence of RIVE exhibits a delay compared to the amorphization behavior. The superior irradiation resistance of calcite phase compared to quartz phase implies there could be advantages to using calcareous aggregates and lowering the usage of siliceous aggregates for concrete in nuclear power plants for extended operation beyond 60 years. However, other effects such as corrosion, aging and reactions during severe accidents should also be considered, and further investigations are needed.

Amorphous↗

An ensemble formulation of PBL fluxes in a GCM

An ensemble approach is applied to Planetary Boundary Layer (PBL) calculations with the bulk Richardson number identified as the key parameter. An ensemble averaging calculation was carried out to rederive the bulk friction and heat transport coefficients for the mean condition. Two simulations are carried out and compared. Significant differences in PBL fluxes low level cloudiness, land surface roughness heights, and surface evaporation are noted between the modified and unmodified simulations. Modifications to the model were: (1) the relationship between actual and potential Effective Temperature (ET) to accord with Sud and Fennessy (1982); (2) maximum permissible instantaneous ET at any time is 1.5 mm per hr; (3) moisture distribution in low level cumulus convection to be consistent with no precipitation; (4) appearance of supersaturation clouds to be consistent with supersaturation condition at that level; (5) invoking a simple function for stomatal diffusion effect in the ET calculation.

Sud, Y. C.↗

Electromagnetic fields from quantum sources in heavy-ion collisions

We compute the electromagnetic field created by an ultrarelativistic charged particle in vacuum at distances comparable to the particle Compton wavelength. The wave function of the particle is governed by the Klein-Gordon equation,for a scalar particle, or the Dirac equation, for a spin-half particle. The produced electromagnetic field is essentially different in magnitude and direction from the Coulomb field, induced by a classical point charge, due to the quantum diffusion effect. Therefore, a realistic computation of the electromagnetic field produced in heavy-ion collisions must be based upon the full quantum treatment of the valence quarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An Integrated Tool for the Coupled Thermal and Mechanical Analysis of Pyrolyzing Heatshield Materials

Materials that pyrolyze at elevated temperature have been commonly used as thermal protection materials in hypersonic flight, and advanced pyrolyzing materials for this purpose continue to be developed. Because of the large temperature gradients that can arise in thermal protection materials, significant thermal stresses can develop. Advanced applications of pyrolytic materials are calling for more complex heatshield configurations, making accurate thermal stress analysis more important, and more challenging. For non-pyrolyzing materials, many finite element codes are available and capable of performing coupled thermal-mechanical analyses. These codes do not, however, have a built-in capability to perform analyses that include pyrolysis effects. When a pyrolyzing material is heated, one or more components of the original virgin material pyrolyze and create a gas. This gas flows away from the pyrolysis zone to the surface, resulting in a reduction in surface heating. A porous residue, referred to as char, remains in place of the virgin material. While the processes involved can be complex, it has been found that a simple physical model in which virgin material reacts to form char and pyrolysis gas, will yield satisfactory analytical results. Specifically, the effects that must be modeled include: (1) Variation of thermal properties (density, specific heat, thermal conductivity) as the material composition changes; (2) Energy released or absorbed by the pyrolysis reactions; (3) Energy convected by the flow of pyrolysis gas from the interior to the surface; (4) The reduction in surface heating due to surface blowing; and (5) Chemical and mass diffusion effects at the surface between the pyrolysis gas and edge gas Computational tools for the one-dimensional thermal analysis these materials exist and have proven to be reliable design tools. The objective of the present work is to extend the analysis capabilities of pyrolyzing materials to axisymmetric configurations, and to couple thermal and mechanical analyses so that thermal stresses may be efficiently and accurately calculated.

Pronchick, Stephen W.↗

Non-equilibrium effects on the chemistry of nebular condensates - Implications for the planets and asteroids

Kinetic effects, for example nucleation constraints and slow reactions, should have been important in nebular condensation. Consideration of these effects leads to the prediction of pressure-dependent compositions and physical properties of nebular condensates which is consistent with (1) the differences between different classes of chondritic meteorites, (2) some of the differences between planets, and (3) the presence of oxidized iron on the moon and in the eucrite parent body (presumably an asteroid) despite the low abundance of volatiles. Diffusion effects appear to be important for understanding oxygen isotope anomalies in refractory inclusions in Allende. The consideration of kinetic effects leads to more information concerning nebular processes than if equilibrium is assumed.

Blander, M.↗

Effect of finite-rate chemistry and unequal Schmidt numbers on turbulent non-premixed flames modeled with single-step chemistry

The interaction between a quasi-laminar flame and a turbulent flowfield is investigated through direct numerical simulations (DNS) of reacting flow in two- and three-dimensional domains. Effects due to finite-rate chemistry are studied using a single step global reaction A (fuel) + B (oxidizer) yields P (product), and by varying a global Damkoehler number, as a result of which the turbulence-chemistry interaction in the flame is found to generate a wide variety of conditions, ranging from near-equilibrium to near-extinction. Differential diffusion effects are studied by changing the Schmidt number of one reactive species to one-half. It is observed that laminar flamelet response is followed within the turbulent flowfield, except in regions where transient effects seem to dominate.

Chen, J. H.↗

Instability Analysis of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Lawson, Anthony Layiwola↗

An Experimental Method to Measure Gaseous Diffusivity in Tight and Partially Saturated Porous Media via Continuously Monitored Mass Spectrometry

Detection of radioxenon and radioargon produced by underground nuclear explosions is one of the primary methods by which the Comprehensive Nuclear-Test–Ban Treaty (CTBT) monitors for nuclear activities. However, transport of these noble gases to the surface via barometric pumping is a complex process relying on advective and diffusive processes in a fractured porous medium to bring detectable levels to the surface. To better understand this process, experimental measurements of noble gas and chemical surrogate diffusivity in relevant lithologies are necessary. However, measurement of noble gas diffusivity in tight or partially saturated porous media is challenging due to the transparent nature of noble gases, the lengthy diffusion times, and difficulty maintaining consistent water saturation. Here, the quasi-steady-state Ney–Armistead method is modified to accommodate continuous gas sampling via effusive flow to a mass spectrometer. An analytical solution accounting for the cumulative sampling losses and induced advective flow is then derived. Experimental results appear in good agreement with the proposed theory, suggesting the presence of retained groundwater reduces the effective diffusivity of the gas tracers by 10–1000 times. Finally, by using a mass spectrometer, the method described herein is applicable to a broad range of gas species and porous media.

47 OTHER INSTRUMENTATION↗

An analysis of curvature effects for the control of wall-bounded shear flows

The Reynolds stress transport equations are used to predict the effects of simultaneous and sequential combinations of distortions on turbulent boundary layers. The equations are written in general orthogonal curvilinear coordinates, with the curvature terms expressed in terms of the principal radii of curvature of the respective coordinate surfaces. Results are obtained for the cases of two-dimensional and three-dimensional flows in the limit where production and pressure-strain redistribution dominate over diffusion effects.

Gatski, T. B.↗

Evaluation of Lightning Induced Effects in a Graphite Composite Fairing Structure

Defining the electromagnetic environment inside a graphite composite fairing due to lightning is of interest to spacecraft developers. This paper is the first in a two part series and studies the shielding effectiveness of a graphite composite model fairing using derived equivalent properties. A frequency domain Method of Moments (MoM) model is developed and comparisons are made with shielding test results obtained using a vehicle-like composite fairing. The comparison results show that the analytical models can adequately predict the test results. Both measured and model data indicate that graphite composite fairings provide significant attenuation to magnetic fields as frequency increases. Diffusion effects are also discussed. Part 2 examines the time domain based effects through the development of a loop based induced field testing and a Transmission-Line-Matrix (TLM) model is developed in the time domain to study how the composite fairing affects lightning induced magnetic fields. Comparisons are made with shielding test results obtained using a vehicle-like composite fairing in the time domain. The comparison results show that the analytical models can adequately predict the test and industry results.

Trout, Dawn H.↗

Theoretical aspects of tidal and planetary wave propagation at thermospheric heights

A simple semiquantitative model is presented which allows analytic solutions of tidal and planetary wave propagation at thermospheric heights. This model is based on perturbation approximation and mode separation. The effects of viscosity and heat conduction are parameterized by Rayleigh friction and Newtonian cooling. Because of this simplicity, one gains a clear physical insight into basic features of atmospheric wave propagation. In particular, we discuss the meridional structures of pressure and horizontal wind (the solutions of Laplace's equation) and their modification due to dissipative effects at thermospheric heights. Furthermore, we solve the equations governing the height structure of the wave modes and arrive at a very simple asymptotic solution valid in the upper part of the thermosphere. That 'system transfer function' of the thermosphere allows one to estimate immediately the reaction of the thermospheric wave mode parameters such as pressure, temperature, and winds to an external heat source of arbitrary temporal and spatial distribution. Finally, the diffusion effects of the minor constituents due to the global wind circulation are discussed, and some results of numerical calculations are presented.

Volland, H.↗

High-Subsonic Performance Characteristics and Boundary-Layer Investigations of a 12 10-Inch-Inlet-Diameter Conical Diffuser

Performance and boundary-layer data were taken in a 12 degree 10-inch inlet-diameter conical diffuser of 2:1 exit- to inlet-area ratio. These data were taken for two inlet-boundary-layer conditions. The first condition was that of a thinner inlet boundary later (boundary-layer displacement thickness, delta* approximately equal to 0.034) produced by an inlet section approximately 1 inlet diameter in length between the entrance bell and the diffuser. The second condition was a thicker inlet boundary layer (delta* approximately equal to 0.120) produced by an additional inlet section length of approximately 6 diameters. Longitudinal static-pressure distributions were measured fro wall static orifices. Transverse total- and static-pressure surveys were made at the inlet and exit stations. Boundary-layer velocity distributions were measured at seven stations between the inlet and exit. These data were obtained for a Reynolds number (based on inlet diameter) range of 1 x 10(exp 6) to 3.9 x 10(exp 6). The corresponding Mach number range was from M = 0.2 to choking. At the maximum-power-available condition supersonic flow was obtained as far as 4.5 inches downstream from the diffuser inlet with a maximum Mach number of M approximately equal to 1.5. The total-pressure loss through the diffuser in percentage of inlet dynamic pressure was approximately 2.5 percent for the thinner inlet boundary later and 5.5 percent for the thicker inlet boundary later over the lower subsonic range. These valued increased with increasing flow rate- the values for the thicker inlet boundary later more than those for the thinner inlet boundary layer. The diffuser effectiveness, expressed as the ratio of the actual static-pressure rise to the ideal static-pressure rise, was about 85 percent for the thinner inlet boundary layer and about 67 percent for the thicker inlet boundary later in the lower subsonic range. These values decrease with increasing flow rate. Separated flow was observed for both inlet-boundary-layer conditions in the region of adverse pressure gradient just downstream of the transition curvature from inlet section to diffuser. The flow for the thinner-inlet-boundary-layer condition did not fully re-establish itself along the diffuser walls. The thicker inlet-boundary-layer flow, while not completely re-establishing the normal flow pattern downstream of the separated region, did re-establish more successfully than the thinner inlet boundary layer.

Little, B. H., Jr.↗

The effects of kinetic instabilities on the electron cyclotron emission from runaway electrons

In this paper we show that the kinetic instabilities associated with runaway electron beams play an essential role for the production of high-level non-thermal electron–cyclotron-emission (ECE) radiation. Most of the non-thermal ECE comes from runaway electrons in the low-energy regime with large pitch angle, which are strongly scattered by the excited whistler waves. The power of ECE from runaway electrons is obtained using a synthetic diagnostic model based on the reciprocity method. The electron distribution function is calculated using a kinetic simulation model including the whistler wave instabilities and the quasilinear diffusion effects. Simulations based on DIII-D low-density discharge reproduces the rapid growth of the ECE signals observed in DIII-D experiments. Unlike the thermal ECE where radiation for a certain frequency is strongly localized inside the resonance region, the non-thermal ECE radiation from runaway electrons is nonlocal, and the emission-absorption ratio is higher than that of thermal electrons. The runaway electron tail is more significant for ECE with higher frequencies, and the ECE spectrum becomes flatter as RE population grows. The nonlinear behavior of the kinetic instabilities is illustrated in the oscillations of the ECE waves. In conclusion, the good agreement with the DIII-D experimental observations after including the kinetic instabilities clearly illustrate the significance of the scattering effects from wave-particle interactions, which can also be important for runaway electrons produced in disruptions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Water Saturation Relations and Their Diffusion-Limited Equilibration in Gas Shale: Implications for Gas Flow in Unconventional Reservoirs

Large volumes of water are used for hydraulic fracturing of low permeability shale reservoirs to stimulate gas production, with most of the water remaining unrecovered and distributed in a poorly understood manner within stimulated regions. Because water partitioning into shale pores controls gas release, we measured the water saturation dependence on relative humidity (rh) and capillary pressure (P c ) for imbibition (adsorption) as well as drainage (desorption) on samples of Woodford Shale. Experiments and modeling of water vapor adsorption into shale laminae at rh = 0.31 demonstrated that long times are needed to characterize equilibrium in larger (5 mm thick) pieces of shales, and yielded effective diffusion coefficients from 9 × 10 -9 to 3 × 10 -8 m 2 s -1 , similar in magnitude to the literature values for typical low porosity and low permeability rocks. Most of the experiments, conducted at 50°C on crushed shale grains in order to facilitate rapid equilibration, showed significant saturation hysteresis, and that very large P c (~1 MPa) are required to drain the shales. These results quantify the severity of the water blocking problem, and suggest that gas production from unconventional reservoirs is largely associated with stimulated regions that have had little or no exposure to injected water. Finally, gravity drainage of water from fractures residing above horizontal wells reconciles gas production in the presence of largely unrecovered injected water, and is discussed in the broader context of unsaturated flow in fractures.

58 GEOSCIENCES↗

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗

Loss of solar He and Ne from IDPS in subducting sediment: Diffusion and the effect of phase changes

The results of the diffusion experiment for solar He and Ne in IDP's in a magnetic separate from Pacific Ocean sediment suggest that solar He and Ne would be easily released from IDP grains and hence lost from subducting slabs at shallow depths. However, since the diffusion experiments was conducted under high vacuum, there may be a possibility that magnetite grains, which are supposedly the main constituent of the magnetic fraction, might be partly reduced to form a metal phase due to low oxygen fugacity in the experimental condition. If this is the case, such a phase change might affect the gas release and hence the results of the diffusion coefficients. In order to examine whether or not such a phase change really occurred in the condition of the diffusion experiment, I conducted a heating experiment for a magnetic separate from Pacific Ocean sediment. In the same condition as in the diffusion experiment, and the run products were examined with an x ray diffraction method. Three samples were prepared: they were wrapped with platinum foil, put in a vacuum line, and heated in a molybdenum crucible for two hours at 500 C, 800 C,and 950 C, respectively. After cooling the furnace, the samples were taken out from the crucible and analyzed with an x ray diffraction method.

Hiyagon, H.↗