Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIFFUSER”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Modified Closures in Monte Carlo Algorithms for Diffusive Binary Stochastic Media Transport Problems

In a stochastic medium, the material properties at a given spatial location are known only statistically. The most common approach to solving particle transport problems involving binary stochastic media (BSM) is to use the atomic mix approximation in which the transport problem is solved using homogenized (volume-averaged) material properties. A common deterministic model developed for solving BSM particle transport problems is the Levermore-Pomraning (LP) model. Zimmerman and Adams proposed a Monte Carlo algorithm (Algorithm A) that solves the LP equations and another Monte Carlo algorithm (Algorithm B) that locally preserves the sampled material realization; we refer to these Monte Carlo algorithms as the LP and LPLRP (LP local realization preserving) algorithms, respectively. One-dimensional (1D) planar geometry benchmark studies have shown that the LPLRP algorithm is often significantly more accurate than the LP algorithm for problems with an incident angular flux as well as for problems with an interior source. The LPLRP algorithm implementation in these benchmark comparisons made explicit use of the one-dimensional nature of the problem. The LP model is derived by assuming an 'upwind' closure in the coupling term relating the two materials in a binary stochastic medium. Su and Pomraning developed a modified form of this closure by considering the small correlation length limit and requiring the modified closure to produce the correct exponential decay for a source-free half-line albedo problem in rod geometry. They concluded that the modified closure is generally not inferior to the LP closure and in some cases is significantly better. Brantley investigated the use of the Su-Pomraning (SP) closure in the Monte Carlo LP algorithm (LP-SP) for the suite of benchmark problems and concluded that 1) the LP-SP algorithm was somewhat more accurate overall than the LP algorithm and somewhat less accurate overall than the LPLRP algorithm for an incident angular flux benchmark suite, and 2) the LP-SP algorithm was generally the least accurate of the algorithms for an interior source benchmark suite. Larsen et al. performed an asymptotic analysis of the transport equation in 1D planar geometry for the situation in which the physical system is 1) a random binary stochastic medium with material macroscopic total cross sections and mean slab width values of O(1) or smaller and 2) optically thick with weak absorption and sources at each spatial point and therefore globally diffusive. The asymptotic analysis demonstrates that, under these assumptions, the transport equation limits to the conventional diffusion equation with atomically-mixed (volume-averaged) material properties. Larsen et al. also demonstrate that, under these same assumptions, the LP equations asymptotically limit to an atomically-mixed diffusion equation with a diffusion coefficient that is too large, leading to an unphysical flattening of the ensemble-averaged scalar flux distribution. Vasques and Yadav recently performed an asymptotic analysis of an adjusted Levermore-Pomraning closure in which the Markovian transition functions are rescaled such that the equations asymptotically limit to the correct atomically-mixed diffusion equation; we refer to this model as LP-VY. In this paper, we demonstrate that the adjusted LP closure proposed by Vasques and Yadav is a special case of the SP closure obtained by assuming that the absorption cross sections in the two materials are equal. We further demonstrate that the Su-Pomraning closure has the correct asymptotic limit for the diffusive physical system and is therefore an appropriate closure for these problems. We also describe how to incorporate these various closures into the Monte Carlo LP and LPLRP particle transport algorithms, and we present numerical results comparing the accuracy of the Monte Carlo LP, LP-VY, LP-SP, LPLRP, and LPLRP-SP algorithms for the set of diffusive benchmark problems. We demonstrated that the SP closure has the correct asymptotic limit for the diffusive physical system under investigation in this paper and that the VY closure is a special case of the SP closure obtained by assuming equal absorption in the two materials. Both of these closures are readily implemented in the Monte Carlo LP and LPLRP algorithms for particle transport in BSM. Through numerical comparisons to a diffusive benchmark suite, we confirmed that 1) the LP and LPLRP algorithms do not limit to the appropriate diffusion solution and 2) the LP-VY, LP-SP, and LPLRP-SP algorithms limit to the appropriate diffusion solution and are of comparable accuracy. The LPLRP algorithm has been shown to be generally more accurate than the LP algorithm for non-diffusive problems. However, the numerical results in this paper indicate that the LPLRP algorithm does not limit to the correct diffusion equation for the asymptotic limit under consideration. The work in this paper demonstrates that the SP algorithm does limit to the correct diffusion equation for the asymptotic limit considered. Future work will investigate the possibility of whether the LPLRP-SP algorithm is accurate for both diffusive and non-diffusive problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Investigation of Molecular Diffusion at Block Copolymer Thin Films Using Maximum Entropy Method-Based Fluorescence Correlation Spectroscopy and Single Molecule Tracking

Fluorescence correlation spectroscopy (FCS) has been widely used to investigate molecular diffusion behavior in various samples. The use of the maximum entropy method (MEM) for FCS data analysis provides a unique means to determine multiple distinct diffusion coefficients without a priori assumption of their number. Comparison of the MEM-based FCS method (MEM-FCS) with another method will reveal its utility and advantage as an analytical tool to investigate diffusion dynamics. Herein, we measured diffusion of fluorescent probes doped into nanostructured thin films using MEM-FCS, and validated the results with single molecule tracking (SMT) data. The efficacy of the MEM code employed was first demonstrated by analyzing simulated FCS data for systems incorporating one and two diffusion modes with broadly distributed diffusion coefficients. The MEM analysis accurately afforded the number of distinct diffusion modes and their mean diffusion coefficients. These results contrasted with those obtained by fitting the simulated data to conventional two-component and anomalous diffusion models, which yielded inaccurate estimates of the diffusion coefficients. Subsequently, the MEM analysis was applied to FCS data acquired from hydrophilic dye molecules incorporated into microphase-separated polystyrene-block-poly(ethylene oxide) (PS-b-PEO) thin films characterized under a water-saturated nitrogen atmosphere. The MEM analysis revealed distinct fast and slow diffusion components attributable to molecules diffusing on the film surface and inside the film, respectively. SMT studies of the same materials yielded trajectories for mobile molecules that appear to follow the curved PEO microdomains. Diffusion coefficients obtained from the SMT data were consistent with those obtained for the slow diffusion component detected by MEM-FCS. Furthermore, these results highlight the utility of MEM-FCS and SMT for gaining complementary information on molecular diffusion processes in heterogeneous material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A First-Principles Density Function Theory Study of Tritium Diffusion in Li 2 ZrO 3 : Application for Producing Tritium

The ceramic material Li 2 ZrO 3 has superior thermo-physical and thermo-chemical properties and is highly compatible with other blanket materials used in nuclear reactors. Like LiAlO 2 , it could be used in the form of an annular pellet in tritium-producing burnable absorber rods (TPBARs) to produce tritium( 3 H) upon thermal neutron irradiation of lithium isotope ( 6 Li). The radiation damaged pellet crystal contains vacancies, defects of its constituent elements, and several other trapping sites which hinder the 3 H diffusion process and releasing behavior. In this study, we investigate the diffusion mechanisms of 3 H and O 3 H species in Li 2 ZrO 3 ceramic pellet in order to understand the effects on diffusion barriers and diffusion coefficients due to the presence of interstitial and substitutional Li defects, hydroxide (O– 3 H) vacancy defect, and of the interactions of 3 H with O-vacancies in the radiation damaged Li 2 ZrO 3 crystal. We consider several possible diffusion pathways of interstitial and substitutional 3 H and its species in a defective supercell and calculate the activation energy barrier profiles. We find that the smallest activation energy barrier is 0.3 eV and corresponding diffusion coefficient at 600 K is 1.93 × 10 –9 m 2 /s for 3 H substitutional diffusion. The smallest 3 H interstitial diffusion barrier and diffusion coefficient are found to be 0.34 eV and 3.25 × 10 –9 m 2 /s. By examining several channels for diffusion, our results show the most likely diffusion mechanism of 3 H migration is occurring along the c-direction by exchange of 3 H within and between different Li sites. Here, our calculated results reveal that the smallest energy barrier for O 3 H diffusion is 0.75 eV which is when O 3 H is diffusing to the O–Li vacancy pair and the corresponding diffusion coefficient is 6.31 × 10 –13 m 2 /s. In terms of 3 H diffusion, the obtained results indicate that the performance of Li 2 ZrO 3 could be better than γ-LiAlO 2 , a widely used ceramic material in TPBAR for producing 3 H.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Predicting the influence of long-range molecular interactions on macroscopic-scale diffusion by homogenization of the Smoluchowski equation

The macroscopic diffusion constant for a charged diffuser is in part dependent on (1) the volume excluded by solute “obstacles” and (2) long-range interactions between those obstacles and the diffuser. Increasing excluded volume reduces transport of the diffuser, while long-range interactions can either increase or decrease diffusivity, depending on the nature of the potential. We previously demonstrated [P. M. Kekenes-Huskey et al., Biophys. J. 105, 2130 (2013)] using homogenization theory that the configuration of molecular-scale obstacles can both hinder diffusion and induce diffusional anisotropy for small ions. As the density of molecular obstacles increases, van der Waals (vdW) and electrostatic interactions between obstacle and a diffuser become significant and can strongly influence the latter's diffusivity, which was neglected in our original model. Here, we extend this methodology to include a fixed (time-independent) potential of mean force, through homogenization of the Smoluchowski equation. We consider the diffusion of ions in crowded, hydrophilic environments at physiological ionic strengths and find that electrostatic and vdW interactions can enhance or depress effective diffusion rates for attractive or repulsive forces, respectively. Additionally, we show that the observed diffusion rate may be reduced independent of non-specific electrostatic and vdW interactions by treating obstacles that exhibit specific binding interactions as “buffers” that absorb free diffusers. Finally, we demonstrate that effective diffusion rates are sensitive to distribution of surface charge on a globular protein, Troponin C, suggesting that the use of molecular structures with atomistic-scale resolution can account for electrostatic influences on substrate transport. This approach offers new insight into the influence of molecular-scale, long-range interactions on transport of charged species, particularly for diffusion-influenced signaling events occurring in crowded cellular environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Charge Carrier Diffusion Length Determination in 4H-SiC Schottky Alpha Detectors

When developing semiconductor radiation detectors detailed characterization of the detector behavior is required. The ability of a semiconductor radiation detector to collect the charge liberated by incident ionizing radiation is a critical characteristic. Several factors influence the charge collection efficiency of a semiconductor radiation detector, including electron-hole trapping, applied electric field strength, and dopant concentrations and dimensions of low doped regions. Applying a reverse bias electric potential to a low doped semiconductor region creates a free charge carrier depleted region and establishes an electric field which can remove subsequently freed charge in the form of an electrical signal. When ionizing radiation is incident upon a semiconductor electron-hole pairs are generated. These charge carriers are created along the path of the ionizing radiation regardless of the doping concentration. However, only the charge freed in the low doped epitaxial region of a detector will be collected as a result of the applied electric field. When the electric field is insufficient to develop a depleted region equal to or larger than the range of the ionizing radiation in the detector material some charge freed by the ionizing radiation outside of the sensitive region will be lost, thus potentially resulting in less than complete charge collection. The incomplete collection of charge complicates the analysis of the incident particle energy spectrum. Charge freed in the highly doped bulk region beyond the low doped epitaxial region may diffuse back into the low doped depleted region and be subsequently included in the integral pulse. The amount of charge collected as a result of this diffusion mechanism is dependent upon the diffusion length of the charge carriers in the material. The diffusion length is subsequently dependent on the material properties of the detector including the number of defects, trapping centers, and temperature. In order to fully characterize the properties of 4H-SiC as a radiation detector the diffusion length was calculated by adjusting the theoretical diffusion length value used in the calculation of the charge collection efficiency to agree with empirically observed charge collection efficiencies. The discrepancy between the observed and theoretical charge collection efficiency at 50 to 150 volts applied bias indicates an unexplained charge collection phenomenon which is contributing to the larger than predicted charge collection efficiency of the detectors. This voltage range corresponds to an intermediate applied bias, indicating that either the theoretical depletion width is too short, the diffusion length varies with the composition of the region that the charge carriers are diffusing through, or some other phenomenon is involved. If the charge carrier diffusion length varies depending on whether or not the diffusing charge is freed in unbiased epitaxy or bulk layer, then a function describing the variation of the diffusion layer could be determined by measuring the diffusion length at the applied voltages used previously, but over a range of epitaxial layer thicknesses. Future work includes the characterization of the effects of high temperature operation on the diffusion length of 4H-SiC in addition to measuring the effects of varying energy alpha emitter electrodeposition on the observed diffusion length. Additional Schottky barrier contacts, namely those with a nickel-platinum contact, will be tested in various sizes which have been developed in order to eliminate certain failure mechanisms and explore the potential effects of Schottky contact size on charge collection efficiency. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2017

Comparison of Experimental Methods for Estimating Matrix Diffusion Coefficients for Contaminant Transport Modeling

Diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2018

Comparison of experimental methods for estimating matrix diffusion coefficients for contaminant transport modeling

Here, diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating effective matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of effective matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than effective matrix diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields effective matrix diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Density functional theory calculations of self- and Xe diffusion in U 3 Si 2

Uranium silicide, U 3 Si 2 , has been proposed as an advanced nuclear fuel to be used in light water reactors (LWRs). Development of this alternative to the predominant current fuel, UO 2 , is motivated by enhanced accident tolerance as a result of higher thermal conductivity as well as improved fuel cycle economics through increased uranium density. In order to accurately model the fuel performance of U 3 Si 2 , the diffusion rate of point defects, which is related to self-diffusion, and of fission gas atoms must be determined. DFT calculations are used to predict the U and Si point defect concentrations, the corresponding self-diffusivities, the preferred Xe trap site and the Xe diffusivity. Effects of irradiation are not considered. A low defect formation energy and a high entropy for Si interstitials give rise to Si-rich non-stoichiometry at elevated temperatures. Both U and Si self-diffusion and Xe diffusion are anisotropic as a consequence of the tetragonal crystal structure of U 3 Si 2 . Si diffusion occurs by interstitial mechanisms in both the a-b plane and along the c axis, while the U c axis diffusion rate is controlled by a vacancy mechanism. Interstitial diffusion of U is very fast in the a-b plane of the U 3 Si 2 crystal structure. Xe atoms prefer to occupy U vacancy trap sites. The highest Xe diffusion rate occurs by a vacancy mechanism in both the a-b plane and along the c axis. The diffusion rate is similar in the a-b plane and along the c axis. U and Si self-diffusion and Xe diffusion are all faster in U3Si2 than intrinsic U and Xe diffusion in conventional UO 2 nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
NASA NTRS1974

Boundary-layer analysis of subsonic inlet diffuser geometries for engines nacelles

Theoretical Mach number distributions and boundary-layer parameters are presented for subsonic nacelle inlet diffuser geometries with length to exit diameter ratios ranging from 0.4 to 1.6 and diffuser exit area to throat area ratios ranging from 1.1 to 2.0. The major portion of the study was done with a cubic diffuser contour with the inflection point at the midpoint of the diffuser, a diffuser throat Mach number of 0.6, and a free-stream Mach number of 0.12. Calculations were performed at both model (diffuser exit diameter, 30.5 cm) and full-scale (diffuser exit diameter, 183 cm) sizes. Separation limits were defined by establishing a separation boundary on plots of diffuser area ratio as a function of diffuser length to diameter ratio. The effects of diffuser contour, inlet lip geometry, and throat Mach number on the boundary-layer characteristics are illustrated. The major results of the study indicate that the separation boundary is shifted to greater area ratios by (1) increasing the diffuser length, (2) increasing the scale of the diffuser and, (3) moving the inflection point of the diffuser contour to or ahead of the midpoint of the diffuser.

Albers, J. A.↗
NASA NTRS1993

Irradiance calibration with solar diffuser

The sun's energy is used in combination of movable and fixed diffuser plates, windows and apertures which are positioned in a series of test sequences (modes) for reflectance monitoring and calibration without the use of man-made sources. There are three embodiments, or implementations, of the invention--one embodiment uses two diffusers--a working diffuser and a secondary diffuser--the second embodiment uses three diffusers, a working diffuser, a secondary diffuser and a reference diffuser--and the third embodiment uses two diffusers--a working diffuser and a secondary diffuser, the latter also functioning as a cover for the working diffuser. The movable diffusers are mounted on rotatable cones and, in all embodiments, the sun is blocked from reaching the diffusers when not in use. Thus, the sun is used as a stable source for calibration and monitoring and the sun/diffuser combination is used in such a way that the response of all elements of the optical subsystem of the TOMS can be unambiguously and efficiently characterized with high accuracy and precision.

Haring, Robert E.↗
NASA NTRS2014

Spring 2014 Internship Diffuser Data Analysis

J-2X engine testing on the A-2 test stand at the NASA John C. Stennis Space Center (SSC) has recently concluded. As part of that test campaign, the engine was operated at lower power levels in support of expanding the use of J-2X to other missions. However, the A-2 diffuser was not designed for engine testing at the proposed low power levels. To evaluate the risk of damage to the diffuser, computer simulations were created of the rocket engine exhaust plume inside the 50ft long, water-cooled, altitude-simulating diffuser. The simulations predicted that low power level testing would cause the plume to oscillate in the lower sections of the diffuser. This can possibly cause excessive vibrations, stress, and heat transfer from the plume to the diffuser walls. To understand and assess the performance of the diffuser during low power level engine testing, nine accelerometers and four strain gages were installed around the outer surface of the diffuser. The added instrumentation also allowed for the verification of the rocket exhaust plume computational model. Prior to engine hot-fire testing, a diffuser water-flow test was conducted to verify the proper operation of the newly installed instrumentation. Subsequently, two J-2X engine hot-fire tests were completed. Hot-Fire Test 1 was 11.5 seconds in duration, and accelerometer and strain data verified that the rocket engine plume oscillated in the lower sections of the diffuser. The accelerometers showed very different results dependent upon location. The diffuser consists of four sections, with Section 1 being closest to the engine nozzle and Section 4 being farthest from the engine nozzle. Section 1 accelerometers showed increased amplitudes at startup and shutdown, but low amplitudes while the diffuser was started. Section 3 accelerometers showed the opposite results with near zero G amplitudes prior to and after diffuser start and peak amplitudes to +/- 100G while the diffuser was started. Hot-Fire Test 1 strain gages showed different data dependent on section. Section 1 strains were small, and were in the range of 50 to 150 microstrain, which would result in stresses from 1.45 to 4.35 ksi. The yield stress of the material, A-285 Grade C Steel, is 29.7 ksi. Section 4 strain gages showed much higher values with strains peaking at 1600 microstrain. This strain corresponds to a stress of 46.41 ksi, which is in excess of the yield stress, but below the ultimate stress of 55 to 75 ksi. The decreased accelerations and strain in Section 1, and the increased accelerations and strain in Sections 3 and 4 verified the computer simulation prediction of increased plume oscillations in the lower sections of the diffuser. Hot-Fire Test 2 ran for a duration of 125 seconds. The engine operated at a slightly higher power level than Hot-Fire Test 1 for the initial 35 seconds of the test. After 35 seconds the power level was lowered to Hot-Fire Test 1 levels. The acceleration and strain data for Hot-Fire Test 2 was similar during the initial part of the test. However, just prior to the engine being lowered to the Hot-Fire Test 1 power level, the strain gage data in Section 4 showed a large decrease to strains near zero microstrain from their peak at 1500 microstrain. Future work includes further strain and acceleration data analysis and evaluation.

Laigaie, Robert T.↗
DOE OSTI2022

Effect of charge and anisotropic diffusivity on ion exchange kinetics in nuclear waste form materials

Motifs of radionuclide absorption particles (absorbers or fillers) for the nuclear waste treatment usually have crystal structures with nano-sized tunnels for fast ion diffusion. On the other hand, the accumulation of ions generates an electric field which affects the ion diffusion. Therefore, the strong anisotropic kinetic properties and the electric field affect ion exchange kinetics, especially in an assembly of these particles which may block the fast diffusion tunnels for each other. In this work, we developed a mesoscale model to describe the effect of anisotropic kinetic properties and electric field on ion exchange kinetics. With the model we simulated the effect of anisotropic diffusivity, electrochemical potential, particle morphology, size and aggregation on ion exchange kinetics during batch experiments. Here, tt is found that 1) ion exchange results in charge accumulation inside the particle because different ions have different diffusivities and different diffusion driving forces. The electric neutrality is not preserved; 2) the electric field speeds up the slower diffusion ions (either uptake or release) and slows down the faster diffusion ions to reduce charge accumulation and maintains charge neutrality; 3) the ion exchange becomes very slow at the later stages because a) diffusion driving force (electrochemical gradient, Cs + concentration in the solution and Na + inside the particle) continuously decreases, as observed in batch experiments; and b) the diffusivity decreases due to the Cs + concentration dependence of ion diffusivity. This can explain the observed experimental phenomena that the radionuclide capacity usually cannot be reached; 4) for needle-like particle and their clusters, both surface area and diffusion distance along the direction with the fastest diffusivities have significant impact on the ion exchange kinetics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Perturbative thermal diffusivity from partial sawtooth crashes in Alcator C-Mod

Perturbative thermal diffusivity has been measured on Alcator C-Mod via the use of the extended-time-to-peak method on heat pulses generated by partial sawtooth crashes. Perturbative thermal diffusivity governs the propagation of heat pulses through a plasma. It differs from power balance thermal diffusivity, which governs steady state thermal transport. Heat pulses generated by sawtooth crashes have been used extensively in the past to study heat pulse thermal diffusivity, but the details of the sawtooth event typically lead to non-diffusive 'ballistic' transport, making them an unreliable measure of perturbative diffusivity on many tokamaks. Partial sawteeth are common on numerous tokamaks, and generate a heat pulse without the 'ballistic' transport that often accompanies full sawteeth. This is the first application of the extended-time-to-peak method of diffusivity calculation to partial sawtooth crashes. Furthermore, this analysis was applied to over 50 C-Mod shots containing both L- and I-Mode. Results indicate correlations between perturbative diffusivity and confinement regime (L- versus I-mode), as well as correlations with local temperature, density, the associated gradients, and gradient scale lengths (a/L Te and a/L n ). In addition, diffusivities calculated from partial sawteeth are compared to perturbative diffusivities calculated with the nonlinear gyrokinetic code GYRO. We find that standard ion-scale simulations (ITG/TEM turbulence) under-predict the perturbative thermal diffusivity, but new multi-scale (ITG/TEM coupled with ETG) simulations can match the experimental perturbative diffusivity within error bars for an Alcator C-Mod L-mode plasma. Perturbative diffusivities extracted from heat pulses due to partial sawteeth provide a new constraint that can be used to validate gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Analytical and experimental studies of a short compact subsonic diffuser for a two-dimensional supersonic inlet

An experimental study of a two-dimensional supersonic inlet with a short compact subsonic diffuser, length to exit diameter (dl/d) ratio of 1.25, was conducted to investigate the impact of the short diffuser on inlet performance at low speeds and to assess the diffuser subsonic performance for a simulated diffuser flow corresponding to high-speed inlet conditions near the design flight Mach number of 2.2. For the low-speed testing, a drooped lip was employed to improve the inlet performance at a high angle of attack. For the simulated high-speed testing, air was blown through slots or discrete nozzles as an active boundary-layer control. The results from the low-speed performance test were compared with the results from a previous test program on the same inlet with a long subsonic diffuser (dl/d = 4.5). The comparison indicates that inlet recovery was not affected by the use of the short diffuser for either the baseline (no droop) or the drooped cowl lip configuration. However, the inlet baseline distortion for the short diffuser configuration was substantially higher than for the long diffuser. A comparison of the two configurations with a 70 deg drooped lip showed no significant difference in distortion. For the portion of the experimental program in which diffuser conditions for high-speed flight were simulated, diffuser-induced flow separation occurred. This separation was predicted from an analytical study that used the Hess potential flow panel method and the Herring two-dimensional boundary-layer analysis computer codes. The flow separated mainly on the diffuser ramp. Subsequent tests in which boundary-control systems were utilized showed that blowing with either slots or discrete nozzles could suppress the flow separation in the short subsonic diffuser, thereby substantially improving the diffuser performance.

Iek, Chanthy↗
NASA NTRS2000

Laminar Diffusion Flame Studies (Ground- and Space-Based Studies)

Laminar diffusion flames are of interest because they provide model flame systems that are far more tractable for analysis and experiments than more practical turbulent diffusion flames. Certainly, understanding flame processes within laminar diffusion flames must precede understanding these processes in more complex turbulent diffusion flames. In addition, many properties of laminar diffusion flames are directly relevant to turbulent diffusion flames using laminar flamelet concepts. Laminar jet diffusion flame shapes (luminous flame boundaries) have been of particular interest since the classical study of Burke and Schumann because they are a simple nonintrusive measurement that is convenient for evaluating flame structure predictions. Thus, consideration of laminar flame shapes is undertaken in the following, emphasizing conditions where effects of gravity are small, due to the importance of such conditions to practical applications. Another class of interesting properties of laminar diffusion flames are their laminar soot and smoke point properties (i.e., the flame length, fuel flow rate, characteristic residence time, etc., at the onset of soot appearance in the flame (the soot point) and the onset of soot emissions from the flame (the smoke point)). These are useful observable soot properties of nonpremixed flames because they provide a convenient means to rate several aspects of flame sooting properties: the relative propensity of various fuels to produce soot in flames; the relative effects of fuel structure, fuel dilution, flame temperature and ambient pressure on the soot appearance and emission properties of flames; the relative levels of continuum radiation from soot in flames; and effects of the intrusion of gravity (or buoyant motion) on emissions of soot from flames. An important motivation to define conditions for soot emissions is that observations of laminar jet diffusion flames in critical environments, e.g., space shuttle and space station facilities, cannot involve soot emitting flames in order to ensure that test chamber windows used for experimental observations are not blocked by soot deposits, thereby compromising unusually valuable experimental results. Another important motivation to define conditions where soot is present in diffusion flames is that flame chemistry, transport and radiation properties are vastly simplified when soot is absent, making such flames far more tractable for detailed numerical simulations than corresponding soot-containing flames. Motivated by these observations, the objectives of this phase of the investigation were as follows: (1) Observe flame-sheet shapes (the location of the reaction zone near phi=1) of nonluminous (soot free) laminar jet diffusion flames in both still and coflowing air and use these results to develop simplified models of flame-sheet shapes for these conditions; (2) Observe luminous flame boundaries of luminous (soot-containing) laminar jet diffusion flames in both still and coflowing air and use these results to develop simplified models of luminous flame boundaries for these conditions. In order to fix ideas here, maximum luminous flame boundaries at the laminar smoke point conditions were sought, i.e., luminous flame boundaries at the laminar smoke point; (3) Observe effects of coflow on laminar soot- and smoke-point conditions because coflow has been proposed as a means to control soot emissions and minimize the presence of soot in diffusion flames.

Dai, Z.↗