Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional diffusion”

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

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

At least 361 records · Page 20

Sagittarius A* as an origin of the Galactic PeV cosmic rays?

Supernova remnants (SNRs) have commonly been considered as a source of the observed PeV cosmic rays (CRs) or a Galactic PeV particle accelerator ('Pevatron'). In this work, we study Sagittarius A* (Sgr A*), which is the low-luminosity active galactic nucleus of the Milky Way Galaxy, as another possible canditate of the Pevatron, because it sometimes became very active in the past. We assume that a large number of PeV CRs were injected by Sgr A* at the outburst about 10{sup 7} yr ago when the Fermi bubbles were created. We constrain the diffusion coefficient for the CRs in the Galactic halo on the condition that the CRs have arrived on the Earth by now, while a fairly large fraction of them have escaped from the halo. Based on a diffusion-halo model, we solve a diffusion equation for the CRs and compare the results with the CR spectrum on the Earth. The observed small anisotropy of the arrival directions of CRs may be explained if the diffusion coefficient in the Galactic disk is smaller than that in the halo. Our model predicts that a boron-to-carbon ratio should be energy-independent around the knee, where the CRs from Sgr A* become dominant. It is unlikely that the spectrum of the CRs accelerated at the outburst is represented by a power-law similar to the one for those responsible for the gamma-ray emission from the central molecular zone (CMZ) around the Galactic center.

79 ASTRONOMY AND ASTROPHYSICS↗

Mixtures of octanol and an ionic liquid: Structure and transport

Ionic liquids (ILs) with long alkyl substituents are amphiphilic, which leads to a bicontinuous liquid structure. The strongly interacting anionic and cationic head groups form a long range charge network, with the hydrocarbon tails forming a nonpolar domain. Such nonpolar domains have been shown to dissolve a variety of neutral organic solvents. In mixtures of ILs with solvents the neutral organic molecules residing in the nonpolar domains experience different environments and friction from the charged cations and anions. Thus, the neutral molecules diffuse much faster than predicted by hydrodynamic scaling using the average viscosity of the mixture. Here, we report studies on the structure and transport properties of mixtures of 1-octanol with the IL trihexyltetradecylphosphonium bis(trifluoromethylsulfonyl)imide (P 6,6,6,14 + /NTf 2 – ). The majority of the atom fraction in the P 6,6,6,14 + cation comprises four hydrocarbon substituents. The unique amphiphilic nature of ILs with the P 6,6,6,14 + cation makes 1-octanol fully miscible with the IL at ambient temperatures. X-ray scattering experiments show that the IL structure persists in the mixtures for 1-octanol mole fractions as large as x oct = 0.90. The self-diffusion coefficients of the three molecular species in the mixtures were measured by NMR experiments. The self-diffusion of the P 6,6,6,14 + cation is well described by the Stokes–Einstein equation, while the diffusivity of the NTf 2 – anion is slightly lower than the hydrodynamic prediction. The measured diffusivities of octanol in these mixtures are 1.3–4 times higher than the hydrodynamic predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Mg impurity in nitride alloys

Although several magnetic resonance studies address the Mg acceptor in GaN, there are few reports on Mg doping in the alloys, where hole production depends strongly on the Al or In content. Our electron paramagnetic resonance (EPR) measurements of the p-type alloys suggest that the Mg impurity retains the axial symmetry, characteristic of a p-type dopant in both alloys; however, In and Al produce additional, different characteristics of the acceptor. In InGaN, the behavior is consistent with a lowering of the acceptor level and increasing hole density as In concentration increases. For AlGaN, the amount of neutral Mg decreases with increasing Al content, which is attributed to different kinetics of hydrogen diffusion thought to occur in samples with higher Al mole fraction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

WE-G-18C-09: Separating Perfusion and Diffusion Components From Diffusion Weighted MRI of Rectum Tumors Based On Intravoxel Incoherent Motion (IVIM) Analysis

Purpose: Pseudodiffusion arises from the microcirculation of blood in the randomly oriented capillary network and contributes to the signal decay acquired using a multi-b value diffusion weighted (DW)-MRI sequence. This effect is more significant at low b-values and should be properly accounted for in apparent diffusion coefficient (ADC) calculations. The purpose of this study was to separate perfusion and diffusion component based on a biexponential and a segmented monoexponential model using IVIM analysis Methods. The signal attenuation is modeled as S(b) = S0[(1−f)exp(−bD) + fexp(−bD*)]. Fitting the biexponetial decay leads to the quantification of D, the true diffusion coefficient, D*, the pseudodiffusion coefficient, and f, the perfusion fraction. A nonlinear least squares fit and two segmented monoexponential models were used to derive the values for D, D*,‘and f. In the segmented approach b = 200 s/mm{sup 2} was used as the cut-off value for calculation of D. DW-MRI's of a rectum cancer patient were acquired before chemotherapy, before radiation therapy (RT), and 4 weeks into RT and were investigated as an example case. Results: Mean ADC for the tumor drawn on the DWI cases was 0.93, 1.0 and 1.13 10{sup −3}×mm{sup 2}/s before chemotherapy, before RT and 4 weeks into RT. The mean (D.10{sup −3} × mm{sup 2}/s, D* 10{sup −3} × mm{sup 2}/s, and f %) based on biexponential fit was (0.67, 18.6, and 27.2%), (0.72, 17.7, and 28.9%) and (0.83,15.1, and 30.7%) at these time points. The mean (D, D* f) based on segmented fit was (0.72, 10.5, and 12.1%), (0.72, 8.2, and 17.4%) and (.82, 8.1, 16.5%) Conclusion: ADC values are typically higher than true diffusion coefficients. For tumors with significant perfusion effect, ADC should be analyzed at higher b-values or separated from the perfusion component. Biexponential fit overestimates the perfusion fraction because of increased sensitivity to noise at low b-values.

60 APPLIED LIFE SCIENCES↗

Experimental investigation of the effect of pore size distribution on nano-particle capture efficiency within ceramic particulate filters

The effect of the pore size distribution on size-resolved filtration efficiency was investigated for two ceramic particulate filters using particulate matter (PM) generated by a spark-ignition direct-injection engine fueled with gasoline. The cordierite filter tested had a porosity of 43%, a median pore diameter (d_50) of 12 µm, and a wide pore size distribution with a lognormal standard deviation (s’) of 0.3. The aluminum titanate filter had very similar porosity, d_50, and thickness, but significantly narrower pore size distribution (s’ = 0.1). The testing of two filters under identical experimental conditions enabled the impact of the pore size distribution on filtration performance to be evaluated. Filtration experiments were performed focusing on just the filter wall, starting from a clean filter until the transition to cake filtration (filtration efficiency > 99%). Time-resolved particle size distribution measurements were used to evaluate the progression of filtration performance and estimate trapped mass within the filter. The aluminum titanate filter, with a narrow pore size distribution, exhibited significantly better diffusion capture efficiency. The negative impact of higher flow velocity on diffusion capture efficiency was more pronounced for a narrower pore size distribution. Flow distributions measured using capillary-flow porometry were used to develop a cylindrical pore flow model to understand the impact of the differences in pore size distribution on observed trends in diffusion capture efficiency within a clean filter. The model helps demonstrate how large pores on the high end of a wider pore size distribution can carry a large fraction of the flow, with less efficient capture of small particles by diffusion. The experimental results and data demonstrate that the bubble point diameter and width of the pore size distribution significantly influence diffusion capture efficiencies for filters with very similar median pore diameter, porosity, and thickness.

aerosol filtration, Aftertreatment, ceramic exhaus↗

Polarimetry of the Interstellar Medium

The talk will review what is known about the composition of ices and organics in the dense and diffuse interstellar media (ISM). Mixed molecular ices make up a significant fraction of the solid materials in dense molecular clouds and it is now known that thermal and radiation processing of these ices results in the production of more complex organic species, some of which may survive transport into forming stellar systems and the diffuse ISM. Molecular species identified in interstellar ices include H2O, CH3OH, CO, CH4, CO2, and somewhat surprisingly, H2. Theoretical and laboratory studies of the processing of interstellar analog ices containing these species indicate that species like HCO, H2CO, CH3, and NH3 are readily made and should also be present. The irradiation of mixed molecular ices containing these species, when followed by warming, leads to the production of a large variety of more complex species, including ethanol (CH3CH2OH), formamide (HC(=O)NH2), acetamide (CH3C(=O)NH2), nitriles or isonitriles (R-CN or R-NC hexamethylenetetramine (HMT; C6H12N4), a number of polymeric species related to polyoxymethylene [POM,(-CH2O-)n], and ketones {R-C(=O)-R'}. Spectral studies of dust in the diffuse ISM indicate the presence of fairly complex organics, some of which may be related to the organics produced in dense molecular clouds. Spectral comparisons indicate that the diffuse ISM organics may be quite similar to meteoritic kerogens, i.e. they may consist largely of aromatic moieties interlinked by short aliphatic bridges. Interestingly, recent evidence indicates that the galactic distribution of this material closely matches that of silicates, but does not correlate directly with visual extinction. This implies that a large fraction of the visual extinction is caused by a material other than these organics and silicates and that this other material has a significantly different distribution within the galaxy.

Sandford, Scott↗

Periodic Excitation for Jet Vectoring and Enhanced Spreading

The effects of periodic excitation on the evolution of a turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit and excitation was introduced at the junction between the jet exit and the diffuser inlet. The introduction of high amplitude periodic excitation at the jet exit enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Vectoring is achieved by applying the excitation over a fraction of the circumference of the circular jet, enhancing its spreading rate on the excited side and its tendency to reattach to that side. Static deflection studies demonstrate that the presence of the wide-angle diffuser increases the effectiveness of the added periodic momentum due to a favorable interaction between the excitation, the jet shear-layer and the diffuser wall. This point was further demonstrated by the evolution of a wave packet that was excited in the jet shear-layer. Strong amplification of the wave packet was measured with a diffuser attached to the jet exit. The turbulent jet responds quickly (10-20 msec) to step changes in the level of the excitation input. The response scales with the jet exit velocity and is independent of the Reynolds number. Jet deflection angles were found to be highly sensitive to the relative direction between the excitation and the jet flow and less sensitive to the excitation frequency. The higher jet deflection angles were obtained for a diffuser length of about two diameters and for diffusers with half-angles greater than 15 degrees.

Pack, LaTunia G.↗

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗

Method and system for gas flow mitigation of molecular contamination of optics

A computer-implemented method for determining an optimized purge gas flow in a semi-conductor inspection metrology or lithography apparatus, comprising receiving a permissible contaminant mole fraction, a contaminant outgassing flow rate associated with a contaminant, a contaminant mass diffusivity, an outgassing surface length, a pressure, a temperature, a channel height, and a molecular weight of a purge gas, calculating a flow factor based on the permissible contaminant mole fraction, the contaminant outgassing flow rate, the channel height, and the outgassing surface length, comparing the flow factor to a predefined maximum flow factor value, calculating a minimum purge gas velocity and a purge gas mass flow rate from the flow factor, the contaminant mass diffusivity, the pressure, the temperature, and the molecular weight of the purge gas, and introducing the purge gas into the semi-conductor inspection metrology or lithography apparatus with the minimum purge gas velocity and the purge gas flow rate.

42 ENGINEERING↗

CO DIFFUSION INTO AMORPHOUS H{sub 2}O ICES

The mobility of atoms, molecules, and radicals in icy grain mantles regulates ice restructuring, desorption, and chemistry in astrophysical environments. Interstellar ices are dominated by H{sub 2}O, and diffusion on external and internal (pore) surfaces of H{sub 2}O-rich ices is therefore a key process to constrain. This study aims to quantify the diffusion kinetics and barrier of the abundant ice constituent CO into H{sub 2}O-dominated ices at low temperatures (15–23 K), by measuring the mixing rate of initially layered H{sub 2}O(:CO{sub 2})/CO ices. The mixed fraction of CO as a function of time is determined by monitoring the shape of the infrared CO stretching band. Mixing is observed at all investigated temperatures on minute timescales and can be ascribed to CO diffusion in H{sub 2}O ice pores. The diffusion coefficient and final mixed fraction depend on ice temperature, porosity, thickness, and composition. The experiments are analyzed by applying Fick’s diffusion equation under the assumption that mixing is due to CO diffusion into an immobile H{sub 2}O ice. The extracted energy barrier for CO diffusion into amorphous H{sub 2}O ice is ∼160 K. This is effectively a surface diffusion barrier. The derived barrier is low compared to current surface diffusion barriers in use in astrochemical models. Its adoption may significantly change the expected timescales for different ice processes in interstellar environments.

74 ATOMIC AND MOLECULAR PHYSICS↗

Limiting Current in Nanostructured Block Copolymer Electrolytes

The next-generation electrolytes for lithium batteries must be able to conduct ions at sufficiently high current densities; yet this regime is rarely studied directly. The limiting current density of an electrolyte quantifies the highest possible rate of ion transport under an applied dc potential. Herein, we report on the limiting current density in twelve nanostructured polystyrene-block-poly(ethylene oxide) (PS-b-PEO, or SEO) copolymer electrolytes. We find that the limiting current at a given salt concentration increases systematically with increasing volume fraction of the PEO block (φ EO ). In contrast, the effective-medium theory, commonly used to analyze conductivity in block copolymer electrolytes, predicts that limiting current is independent of φ EO . To resolve this conundrum, the ionic conductivity, the mutual diffusion coefficient of the salt, and the steady-state current fraction of the block copolymer electrolytes were measured. These measurements enable predictions of limiting current with no adjustable parameters using the concentrated solution theory. We found quantitative agreement between experimentally measured limiting current densities and predictions based on the concentrated solution theory. This work sheds light on how to reliably measure and predict limiting current density in composite electrolytes.

25 ENERGY STORAGE↗

Evolution of Microstructure and Creep Behavior in an Fe-Ni-Cr-Nb-C Alloy during Service in Hydrocarbon Cracker Tubes

Heat resistant, cast austenitic Fe-Ni-Cr-Nb-C alloy tubes are used by the petrochemical industry to transport hydrocarbons at temperatures between 850 and 1050 °C. Exposure to carbon from the process stream and nitrogen from the air causes severe carburization and nitridation in the tube walls, which influences creep resistance and degrades their ductility. A better understanding of these processes would enable their service life to be extended. To this end, we report the results of an experimental and computational study of the changes in microstructure and creep behavior in a representative alloy during service. The microstructure of the fresh material is compared to that of 8- and 10-year-old tubes. The results confirm the presence of M{sub 23}C{sub 6}, M{sub 7}C{sub 3}, and M{sub 2}(CN) phases reported in prior studies, and also reveal the nucleation and growth of the carbo-nitride phase M{sub 6}(CN). We show that the spatial variations of carbide and nitride particle volume fractions in aged tubes can be predicted by a simple reaction–diffusion equation. In addition, the uniaxial compressive creep response of the fresh and aged tubes was obtained as functions of temperature and strain rate. The materials deform by power-law creep, with strain rate versus stress relation as ε-dot =Aσ{sup n}. Aging reduces the stress exponent from n = 8 in fresh tubes to n = 5 in aged tubes, and also reduces the pre-exponential coefficient A. Tests on model alloys with composition identical to that of the matrix combined with finite element simulations of creep in model microstructures reveal that Cr depletion is responsible for the change in stress exponent, while the change in A is caused by the increase in particle volume fraction. By combining the microstructure-scale simulations of creep with the predictions of the reaction–diffusion model of microstructure evolution, the creep behavior of tubes left in service beyond 10 years is estimated.

36 MATERIALS SCIENCE↗

Ion Thermal Decoupling and Species Separation in Shock-Driven Implosions

Anomalous reduction of the fusion yields by 50% and anomalous scaling of the burn-averaged ion temperatures with the ion-species fraction has been observed for the first time in DHe 3 -filled shock-driven inertial confinement fusion implosions. Two ion kinetic mechanisms are used to explain the anomalous observations: thermal decoupling of the D and He 3 populations and diffusive species separation. The observed insensitivity of ion temperature to a varying deuterium fraction is shown to be a signature of ion thermal decoupling in shock-heated plasmas. The burn-averaged deuterium fraction calculated from the experimental data demonstrates a reduction in the average core deuterium density, as predicted by simulations that use a diffusion model. Accounting for each of these effects in simulations reproduces the observed yield trends.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis techniques for tracer studies of oxidation

Analysis techniques to obtain quantitative diffusion data from tracer concentration profiles were developed. Mass balance ideas were applied to determine the mechanism of oxide growth and to separate the fraction of inward and outward growth of oxide scales. The process of inward oxygen diffusion with exchange was theoretically modelled and the effect of lattice diffusivity, grain boundary diffusivity and grain size on the tracer concentration profile was studied. The development of the tracer concentration profile in a growing oxide scale was simulated. The double oxidation technique was applied to a FeCrAl-Zr alloy using 0-18 as a tracer. SIMS was used to obtain the tracer concentration profile. The formation of lacey oxide on the alloy was discussed. Careful consideration was given to the quality of data required to obtain quantitative information.

Basu, S. N.↗

Effects of MOF linker rotation and functionalization on methane uptake and diffusion

The flexible degrees of freedom in metal–organic frameworks (MOFs) can have significant effects on guest molecule behavior. However, in the majority of studies applying molecular simulations to MOFs, the framework is assumed to be rigid in order to minimize computational cost. Here, in this work, we assess the significance of this assumption on a representative example of methane uptake and diffusion in UiO-66. We introduce an open-source code to modify MOFs through functionalization and linker rotation and we perform Grand Canonical Monte Carlo and molecular dynamics simulations of methane in each of the functionalized and linker-rotated derivatives of UiO-66. We find that linker rotation moderately influences methane uptake and significantly influences methane diffusion. Our assessment provides ranges of property values that serve as measures of uncertainty of these two properties associated with linker rotation. We further determine that void volume fraction and minimum pore size are the features that govern methane uptake and diffusion, respectively. These findings illustrate the impact of linker rotation on MOFs and provide design principles to guide future investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of silver, cesium, and strontium release predictions using PARFUME with results from the AGR-1 irradiation experiment

Here, the PARFUME (PARticle FUel ModEl) code was used to predict the release of fission products silver, cesium, and strontium from tristructural isotropic coated fuel particles and compacts during the first irradiation experiment (AGR-1) of the Advanced Gas Reactor Fuel Development and Qualification program. The PARFUME model for the AGR-1 experiment used the fuel compact volume average temperature for each of the 620 days of irradiation to calculate the release of silver, cesium, and strontium from a representative particle for a select number of AGR-1 compacts. Post-irradiation examination measurements provided data on release of these fission products from fuel compacts and fuel particles, and retention of silver in the compacts outside of the silicon carbide (SiC) layer. PARFUME-predicted fractional release of silver, cesium, and strontium was determined and compared to the PIE measurements. For silver, comparisons show a trend of over-prediction at low burnup and under-prediction at high burnup. PARFUME has limitations in the modeling of the temporal and spatial distributions of the temperature and burnup across the compacts, which affects the accuracy of its predictions. Nevertheless, the comparisons on silver release lie in the same order of magnitude. Results show an overall over-prediction of the fractional release of cesium by PARFUME. For particles with failed SiC layers, the over-prediction is by a factor of up to 3, corresponding to a potential over-estimation of the diffusivity in uranium oxycarbide (UCO) by a factor of up to 250. For intact particles, whose release is much lower, the over-prediction is by a factor of up to 100, which could be attributed to an over-estimated diffusivity in SiC by about 40% on average. The release of strontium from intact particles is also over-predicted by PARFUME, which also points towards an over-estimated diffusivity of strontium in either SiC or UCO, or possibly both. The measured strontium fractional release from intact particles varied considerably from compact to compact, making it difficult to assess the effective over-estimation of the diffusivities. Moreover, the release of strontium from particles with failed SiC is difficult to observe experimentally due to the release from intact particles, preventing any conclusions to be made on the accuracy or validity of the PARFUME predictions and the modeled diffusivity of strontium in UCO.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structure of turbulent non-premixed flames modeled with two-step chemistry

Direct numerical simulations of turbulent diffusion flames modeled with finite-rate, two-step chemistry, A + B yields I, A + I yields P, were carried out. A detailed analysis of the turbulent flame structure reveals the complex nature of the penetration of various reactive species across two reaction zones in mixture fraction space. Due to this two zone structure, these flames were found to be robust, resisting extinction over the parameter ranges investigated. As in single-step computations, mixture fraction dissipation rate and the mixture fraction were found to be statistically correlated. Simulations involving unequal molecular diffusivities suggest that the small scale mixing process and, hence, the turbulent flame structure is sensitive to the Schmidt number.

Chen, J. H.↗