Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse Scattering”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Negligible magnetic losses at low temperatures in liquid phase epitaxy grown Y 3 Fe 5 O 12 films

Yttrium iron garnet (Y 3 Fe 5 O 12 ; YIG) has a unique combination of low magnetic damping, high spin-wave conductivity, and insulating properties that make it a highly attractive material for a variety of applications in the fields of magnetics and spintronics. While the room-temperature magnetization dynamics of YIG have been extensively studied, there are limited reports correlating the low-temperature magnetization dynamics to the material structure or growth method. Here, in this study, we investigate liquid phase epitaxy grown YIG films and their magnetization dynamics at temperatures down to 10 K. We show there is a negligible increase in the ferromagnetic resonance linewidth down to 10 K, which is unique when compared with YIG films grown by other deposition methods. From the broadband ferromagnetic resonance measurements, polarized neutron reflectivity, and scanning transmission electron microscopy, we conclude that these liquid phase epitaxy grown films have negligible rare-earth impurities present, specifically the suppression of Gd diffusion from the Gd 3 Ga 5 O 12 (GGG) substrate into the Y 3 Fe 5 O 12 film, and therefore negligible magnetic losses attributed to the slow-relaxation mechanism. Overall, liquid phase epitaxy YIG films have a YIG/GGG interface that is five times sharper and have ten times lower ferromagnetic resonance linewidths below 50 K than comparable YIG films by other deposition methods. Thus, liquid phase epitaxy grown YIG films are ideal for low-temperature experiments/applications that require low magnetic losses, such as quantum transduction and manipulation via magnon coupling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Characteristics of Non-Fullerene Acceptor-Based Organic Photovoltaic Active Layers Using X-ray Scattering and Solid-State NMR

With the recent development of non-fullerene acceptors, power conversion efficiencies of bulk-heterojunction organic solar cells have exceeded 18%. The morphology of the BHJ active layer, including packing, ordering, orientation, and phase behavior of the donor(s) and acceptor(s), plays critical roles in determining the device performance. We characterized the morphology of active layers consisting of mixtures of a PTB7/PTB7-Th (donor) and ITIC (acceptor) using grazing incidence wide-angle X-ray scattering, resonant soft X-ray scattering, and solid-state nuclear magnetic resonance (ssNMR) to correlate the morphology with device performance. PTB7-Th/ITIC shows a better device performance than that of PTB7/ITIC due to a smaller π-πstacking distance and smaller domain size of the phase separated morphology. One of PTB7/ITIC samples, processed from mixed solvents (chlorobenzene/benzene = 1:1 by volume), shows a possible partial miscibility and smaller domain size as revealed by ssNMR T 1ρ relaxation times, which is detrimental to device performance. ssNMR results showed that ITIC could crystallize into different forms depending on processing conditions, which may have implications on the manufacturing of devices using it as an ingredient, as well as its long-term stability in service.

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

The Circum-Galactic Medium of MASsive Spirals. I. An overview and a case study of NGC 5908

The Circum-Galactic Medium of MASsive Spirals is a project studying the overall content, physical and chemical properties, and spatial distributions of the multi-phase circumgalactic medium around a small sample of the most massive (M{sub ∗}≳2×10{sup 11} M{sub ⊙}, v{sub rot}>300 km s{sup −1}) isolated spiral galaxies in the local universe. We introduce the sample and present a detailed case study of the XMM-Newton observation of the hot gas halo of NGC 5908. After data calibration, point source removal, and background analysis, we find that the diffuse soft X-ray emission of NGC 5908 is significantly more extended than the stellar light in the vertical direction. The 0.5–1.25 keV radial intensity profile tracing hot gas emission can be detected above the background out to about 2{sup ′}, or 30 kpc from the nucleus. The unresolved soft X-ray emission can be characterized using a β-model with a slope of β≈0.68. The unresolved 0.5–2 keV luminosity within 1{sup ′} is 6.8×10{sup 39} erg s{sup −1}, but young stellar sources may contribute significantly to it. Assuming a metallicity of 0.2 solar, an upper limit (without subtracting the very uncertain young stellar contribution) to the mass of hot gas within this radius is 2.3×10{sup 9} M{sub ⊙}. The cooling radius is r{sub cool}≈25 kpc or ≈0.06r{sub 200}, within which the hot gas could cool radiatively in less than 10 Gyr, and the cooling of hot gas could significantly contribute to replenishing the gas consumed in star formation. The hot gas accounts for ≈1.9% of the baryons detected within the cooling radius. By comparing NGC 5908 to other galaxies, we find that its X-ray luminosity per stellar mass is consistent with lower-mass non-starburst field spiral galaxies. However, a large scatter in hot gas soft X-ray emissivity is indicated for spiral galaxies with M{sub ∗}≳2×10{sup 11} M{sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SIMULATING ASTRO-H OBSERVATIONS OF SLOSHING GAS MOTIONS IN THE CORES OF GALAXY CLUSTERS

Astro-H will be the first X-ray observatory to employ a high-resolution microcalorimeter, capable of measuring the shift and width of individual spectral lines to the precision necessary for estimating the velocity of the diffuse plasma in galaxy clusters. This new capability is expected to bring significant progress in understanding the dynamics, and therefore the physics, of the intracluster medium. However, because this plasma is optically thin, projection effects will be an important complicating factor in interpreting future Astro-H measurements. To study these effects in detail, we performed an analysis of the velocity field from simulations of a galaxy cluster experiencing gas sloshing and generated synthetic X-ray spectra, convolved with model Astro-H Soft X-ray Spectrometer (SXS) responses. We find that the sloshing motions produce velocity signatures that will be observable by Astro-H in nearby clusters: the shifting of the line centroid produced by the fast-moving cold gas underneath the front surface, and line broadening produced by the smooth variation of this motion along the line of sight. The line shapes arising from inviscid or strongly viscous simulations are very similar, indicating that placing constraints on the gas viscosity from these measurements will be difficult. Our spectroscopic analysis demonstrates that, for adequate exposures, Astro-H will be able to recover the first two moments of the velocity distribution of these motions accurately, and in some cases multiple velocity components may be discerned. The simulations also confirm the importance of accurate treatment of point-spread function scattering in the interpretation of Astro-H/SXS spectra of cluster plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Electrical Properties of Barium and Zirconium Modified NBT Ferroelectric Ceramics

Recently a new wave of interest has risen on relaxor ferroelectrics with complex perovskite structure due to their wide use in fabrication of multilayer ceramic capacitors, electrostrictive actuators, and electromechanical transducers. The polycrystalline 0.93(Bi{sub 0.5}Na{sub 0.5})Ba{sub 0.07}Z{sub 0.04}T{sub 0.96}O{sub 3}(0.07BNBZT) ceramic material, which is in the vicinity of the morphotropic phase boundary (MPB) has been prepared by using high temperature solid state reaction method. The tolerance factor has been estimated and found to be 0.815. XRD analysis revealed a rhombohedral perovskite type structure. SEM micrographs showed highly dense grains with rectangular shape. The average grain size is found to be 1.51{mu}m. Dielectric studies in the material ha indicated relaxor behaviour with diffuse phase transition. High value of {epsilon}{sub m}>1958 is found at 1kHz, Tm (phase transition temperature) 335 deg. C, The diffuseness parameter was established to be 1.60 revealing the relaxor behaviour. Further, to confirm the relaxor behaviour in the material, Vogel-Fulcher (V-F) relation has been used. Estimated V-F parameters are found to be T{sub f} = 138 deg. C, E{sub a} = 0.080 eV and {nu}{sub o} = 2.32x10{sup 8} Hz. Cole-Cole analysis has shown a non-Debye type relaxation in the system. Conductivity studies in the material obeyed the Jonscher's power law in frequency range of (45Hz-5MHz) and temperature range of (35 deg. C - 600 deg. C). The electric conduction in the system may be due to hopping/mobility/ transportation of charge carriers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-resolution Near-infrared Spectroscopy of Diffuse Sources around MWC 1080

To reveal the origins of diffuse Hα emissions observed around the Herbig star MWC 1080, we have performed a high-resolution near-infrared (NIR) spectroscopic observation using the Immersion Grating Infrared Spectrograph. In the NIR H and K bands, we detected various emission lines (six hydrogen Brackett lines, seven H{sub 2} lines, and an [Fe ii] line) and compared their spatial locations with the optical (Hα and [S ii]) and radio ({sup 13}CO and CS) line maps. The shock-induced H{sub 2} and [Fe ii] lines indicate the presence of multiple outflows, consisting of at least three associated young stars in this region. The kinematics of H{sub 2} and [Fe ii] near the northeast (NE) cavity edge supports the idea that the NE main outflow from MWC 1080A is the blueshifted one with a low inclination angle. The H{sub 2} and [Fe ii] lines near the southeast molecular region newly reveal that additional highly blueshifted outflows originate from other young stars. The fluorescent H{sub 2} lines were found to trace photodissociation regions formed on the cylindrical surfaces of the main outflow cavity, which are expanding outward with a velocity of about 10–15 km s{sup −1}. For the Hα emission, we identify its components associated with two stellar outflows and two young stars in addition to the dominant component of MWC 1080A scattered by dust. We also report a few faint Hα features located ∼0.4 pc away in the southwest direction from MWC 1080A, which lie near the axes of the NE main outflow and one of the newly identified outflows.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Characterizing the surface compositions of supported bimetallic PtSn clusters: Effects of cluster-support interactions and surface adsorbates

PtSn bimetallic clusters on TiO2(110) and highly oriented pyrolytic graphite (HOPG) surfaces have been characterized by scanning tunneling microscopy, low energy ion scattering (LEIS), Xray photoelectron spectroscopy, and temperature programmed desorption (TPD); density functional theory (DFT) calculations have also been performed to better understand adsorption of CO and D2 on the PtSn surfaces. On TiO2 at coverages of 2 ML of Pt and 2 ML of Sn, exclusively bimetallic clusters are formed for both orders of deposition because clusters of the first metal completely cover the surface such that all atoms of the second metal are incorporated into the existing clusters. In contrast, on HOPG, the high mobility and weak cluster-support interactions on HOPG result in much larger 2 ML monometallic clusters (~30 Å high) that do not completely cover the surface, and deposition of the second metal produces larger clusters as well as smaller ones. Despite the difference in cluster morphologies for the different orders of deposition and supports, the LEIS experiments demonstrate that in all cases, the PtSn clusters are rich in Sn at the surface, as expected based on the lower surface free energy for Sn compared to Pt. Furthermore, the +0.2 eV shift in the Sn(3d5/2) binding energy observed on all surfaces in the presence of Pt is consistent with PtSn alloy formation. Deposition of 2 ML of Sn on TiO2 produces two-dimensional clusters with oxidation of Sn and reduction of titania at the clustersupport interface, but addition of Pt to the Sn clusters causes Sn to diffuse away from this interface, leaving Sn in the metallic state. TPD experiments on 2 ML Pt/TiO2 with increasing coverages of Sn show that the number of adsorption sites for D2 sharply decreases to nearly zero at 0.5 ML, while CO adsorption decreases to zero only at much higher Sn coverages of 2 ML. DFT studies for Sn modified Pt surfaces and bulk structures demonstrate that for CO adsorption at low Sn coverages (<0.25 ML), the strong Pt-CO interactions induce diffusion of Pt to the cluster surface and the formation of a bulk Pt3Sn alloy, whereas D2 adsorption does not lead to interactions with the Pt surface that are strong enough to induce alloy formation. A single Sn adatom prevents D2 adsorption on four neighboring Pt atoms via site-blocking and the donation of electron density to Pt.

Li, Fangliang↗
DOE OSTI2017

The Diversity of Diffuse Ly α Nebulae around Star-forming Galaxies at High Redshift

We report the detection of diffuse Ly α emission, or Ly α halos (LAHs), around star-forming galaxies at z ≈ 3.78 and 2.66 in the NOAO Deep Wide-Field Survey Boötes field. Our samples consist of a total of ∼1400 galaxies, within two separate regions containing spectroscopically confirmed galaxy overdensities. They provide a unique opportunity to investigate how the LAH characteristics vary with host galaxy large-scale environment and physical properties. We stack Ly α images of different samples defined by these properties and measure their median LAH sizes by decomposing the stacked Ly α radial profile into a compact galaxy-like and an extended halo-like component. We find that the exponential scale-length of LAHs depends on UV continuum and Ly α luminosities, but not on Ly α equivalent widths or galaxy overdensity parameters. The full samples, which are dominated by low UV-continuum luminosity Ly α emitters ( M {sub UV} ≳ −21), exhibit LAH sizes of 5–6 kpc. However, the most UV- or Ly α- luminous galaxies have more extended halos with scale-lengths of 7–9 kpc. The stacked Ly α radial profiles decline more steeply than recent theoretical predictions that include the contributions from gravitational cooling of infalling gas and from low-level star formation in satellites. However, the LAH extent matches what one would expect for photons produced in the galaxy and then resonantly scattered by gas in an outflowing envelope. The observed trends of LAH sizes with host galaxy properties suggest that the physical conditions of the circumgalactic medium (covering fraction, H i column density, and outflow velocity) change with halo mass and/or star formation rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

An Optical Spectrum of the Diffuse Galactic Light from BOSS and IRIS

We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 Å, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-2 increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with the features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemispheres. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions, but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 Å, as evidence for luminescence in the form of extended red emission. The observed strength of the 4000 Å break indicates that at most ~7% of the dust-correlated light at 4000 Å can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence, independent of measurements of extinction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Characterization and Modeling of Grain Boundary Chemistry Evolution in Ferritic Steels under Irradiation

Ferritic/martensitic (FM) steels such as HT-9, T-91 and NF12 with chromium concentrations in the range of 9-12 at.% Cr and high Cr ferritic steels (oxide dispersion strengthened steels with 12-18% Cr) are receiving increasing attention for advanced nuclear applications, e.g. cladding and duct materials for sodium fast reactors, pressure vessels in Generation IV reactors and first wall structures in fusion reactors, thanks to their advantages over austenitic alloys. Predicting the behavior of these alloys under radiation is an essential step towards the use of these alloys. Several radiation-induced phenomena need to be taken into account, including phase separation, solute clustering, and radiation-induced segregation or depletion (RIS) to point defect sinks. RIS at grain boundaries has raised significant interest because of its role in irradiation assisted stress corrosion cracking (IASCC) and corrosion of structural materials. Numerous observations of RIS have been reported on austenitic stainless steels where it is generally found that Cr depletes at grain boundaries, consistently with Cr atoms being oversized in the fcc Fe matrix. While FM and ferritic steels are also subject to RIS at grain boundaries, unlike austenitic steels, the behavior of Cr is less clear with significant scatter and no clear dependency on irradiation condition or alloy type. In addition to the lack of conclusive experimental evidence regarding RIS in F-M alloys, there have been relatively few efforts at modeling RIS behavior in these alloys. The need for predictability of materials behavior and mitigation routes for IASCC requires elucidating the origin of the variable Cr behavior. A systematic detailed high-resolution structural and chemical characterization approach was applied to ion-implanted and neutron-irradiated model Fe-Cr alloys containing from 3 to 18 at.% Cr. Atom probe tomography analyses of the microstructures revealed slight Cr clustering and segregation to dislocations and grain boundaries in the ion-irradiated alloys. More significant segregation was observed in the neutron irradiated alloys. For the more concentrated alloys, irradiation did not affect existing Cr segregation to grain boundaries and segregation to dislocation loops was not observed perhaps due to a change in the dislocation loop structure with increasing Cr concentration. Precipitation of α’ was observed in the neutron irradiated alloys containing over 9 at.% Cr. However ion irradiation appears to suppress the precipitation process. Initial low dose ion irradiation experiments strongly suggest a cascade recoil effect. The systematic analysis of grain boundary orientation on Cr segregation was significantly challenged by carbon contamination during ion irradiation or by existing carbon and therefore carbide formation at grain boundaries (sensitization). The combination of the proposed systematic experimental approach with atomistic modeling of diffusion processes directly addresses the programmatic need for developing and benchmarking predictive models for material degradation taking into account atomistic kinetics parameters

36 MATERIALS SCIENCE↗
DOE OSTI2020

The energetics and dynamics of confinement in flexible frameworks and molecular confinement

Porous frameworks form the chemical and structural basis for critical technologies in separations, catalysis, nuclear waste containment and biomedical applications. Hundreds of zeolites and metal organic frameworks (MOFs) have been synthesized, and their ability to separate and store hydrogen, methane and carbon dioxide has been investigated both experimentally and theoretically. Nevertheless, a fundamental and systematic molecular-level understanding of the thermodynamic and structural factors governing the stability and guest-host interactions in these materials lags behind focused studies of specific systems. Because the guest molecules interact with each other and with the host framework, molecular confinement is a finely balanced and complex phenomenon. The ability of the guest molecules to bind and diffuse through the pores is determined by the nature of the host framework which, in turn, responds to the nature and concentration of guest molecules and to pressure and temperature. The work explores how framework flexibility, tailored by structure, composition, temperature and pressure, is a general phenomenon, similar in nature but variable in extent, in both zeolites and MOFs and is part of a free energy landscape in which framework-guest interactions, pressure, and temperature result in changes in framework geometry and, in some cases, phase transitions. These subtle and/or pronounced changes in lattice geometry, energetics, and dynamics can play a decisive role in confinement and in differentiating the binding of molecules of similar size. The free energy landscape created by these structural changes links polymorphism, amorphization, “breathing,” “gate opening” and confinement. Specifically, the generality of such behavior arises from commonalities in lattice dynamics and energetics of frameworks containing a combination of strong rigid bonds and weaker more flexible deformation modes. Identifying and describing these common and collective phenomena is the focus of the research on a selected group of zeolites and MOFs. Structural studies using X-ray and neutron diffraction explore the mechanical functionality of these important materials, specifically how framework materials respond to changes in temperature, pressure and guest loading. These structural studies are combined with calorimetric measurements, using techniques uniquely developed in the participating laboratories, of heats of formation, heat capacities and entropies, and guest-host interactions. The experimental thermodynamic studies are complemented by inelastic neutron scattering studies of the lattice dynamics related both to framework vibrations and to guest-host interactions. The report below summarizes the work done at UC Davis, which is complemented by work done at BYU (Woodfield) and Virginia Tech (Ross).

36 MATERIALS SCIENCE↗
DOE OSTI2022

The energetics and dynamics of confinement in flexible frameworks and molecular confinement

Porous frameworks form the chemical and structural basis for critical technologies in separations, catalysis, nuclear waste containment and biomedical applications. Hundreds of zeolites and metal organic frameworks (MOFs) have been synthesized, and their ability to separate and store hydrogen, methane and carbon dioxide has been investigated both experimentally and theoretically. Nevertheless, a fundamental and systematic molecular-level understanding of the thermodynamic and structural factors governing the stability and guest-host interactions in these materials lags behind focused studies of specific systems. Because the guest molecules interact with each other and with the host framework, molecular confinement is a finely balanced and complex phenomenon. The ability of the guest molecules to bind and diffuse through the pores is determined by the nature of the host framework which, in turn, responds to the nature and concentration of guest molecules and to pressure and temperature. The work explores how framework flexibility, tailored by structure, composition, temperature and pressure, is a general phenomenon, similar in nature but variable in extent, in both zeolites and MOFs and is part of a free energy landscape in which framework-guest interactions, pressure, and temperature result in changes in framework geometry and, in some cases, phase transitions. These subtle and/or pronounced changes in lattice geometry, energetics, and dynamics can play a decisive role in confinement and in differentiating the binding of molecules of similar size. The free energy landscape created by these structural changes links polymorphism, amorphization, “breathing,” “gate opening” and confinement. Specifically, the generality of such behavior arises from commonalities in lattice dynamics and energetics of frameworks containing a combination of strong rigid bonds and weaker more flexible deformation modes. Identifying and describing these common and collective phenomena is the focus of the research on a selected group of zeolites and MOFs. Structural studies using X-ray and neutron diffraction explore the mechanical functionality of these important materials, specifically how framework materials respond to changes in temperature, pressure and guest loading. These structural studies are combined with calorimetric measurements, using techniques uniquely developed in the participating laboratories, of heats of formation, heat capacities and entropies, and guest-host interactions. The experimental thermodynamic studies are complemented by inelastic neutron scattering studies of the lattice dynamics related both to framework vibrations and to guest-host interactions. This research was a collaboration between UC Davis (A. Navrotsky) BYU (B. Woodfield) and Virginia Tech (N. Ross).

36 MATERIALS SCIENCE↗
DOE OSTI2012

EARLY-TYPE GALAXIES AT z = 1.3. I. THE LYNX SUPERCLUSTER: CLUSTER AND GROUPS AT z = 1.3. MORPHOLOGY AND COLOR-MAGNITUDE RELATION

We confirm the detection of three groups in the Lynx supercluster, at z Almost-Equal-To 1.3, through spectroscopic follow-up and X-ray imaging, and we give estimates for their redshifts and masses. We study the properties of the group galaxies compared to the two central clusters, RX J0849+4452 and RX J0848+4453. Using spectroscopic follow-up and multi-wavelength photometric redshifts, we select 89 galaxies in the clusters, of which 41 are spectroscopically confirmed, and 74 galaxies in the groups, of which 25 are spectroscopically confirmed. We morphologically classify galaxies by visual inspection, noting that our early-type galaxy (ETG) sample would have been contaminated at the 30%-40% level by simple automated classification methods (e.g., based on Sersic index). In luminosity-selected samples, both clusters and groups show high fractions of bulge-dominated galaxies with a diffuse component that we visually identified as a disk and which we classified as bulge-dominated spirals, e.g., Sas. The ETG fractions never rise above Almost-Equal-To 50% in the clusters, which is low compared to the fractions observed in other massive clusters at z Almost-Equal-To 1. In the groups, ETG fractions never exceed Almost-Equal-To 25%. However, overall bulge-dominated galaxy fractions (ETG plus Sas) are similar to those observed for ETGs in clusters at z {approx} 1. Bulge-dominated galaxies visually classified as spirals might also be ETGs with tidal features or merger remnants. They are mainly red and passive, and span a large range in luminosity. Their star formation seems to have been quenched before experiencing a morphological transformation. Because their fraction is smaller at lower redshifts, they might be the spiral population that evolves into ETGs. For mass-selected samples of galaxies with masses M > 10{sup 10.6} M{sub Sun} within {Sigma} > 500 Mpc{sup -2}, the ETG and overall bulge-dominated galaxy fractions show no significant evolution with respect to local clusters, suggesting that morphological transformations might occur at lower masses and densities. The ETG mass-size relation shows evolution toward smaller sizes at higher redshift in both clusters and groups, while the late-type mass-size relation matches that observed locally. When compared to the clusters, the group ETG red sequence shows lower zero points (at {approx}2{sigma}) and larger scatters, both expected to be an indication of a younger galaxy population. However, we show that any allowed difference between the age in groups and clusters would be small when compared to the differences in age in galaxies of different masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗
DOE OSTI2019

Oxygen point defect accumulation in single-phase U O 2 + x

UO 2.07 was characterized using neutron total scattering in order to elucidate defect morphology in the low oxygen-to-metal regime (x < 0.125 for UO 2+x ). Data were collected at temperatures (600 and 1000 °C) coinciding with the single-phase UO 2+x region of the established phase diagram, and results were compared with data of stoichiometric UO 2 collected at near-identical temperatures. Experimental data were modeled and interpreted using a holistic approach employing a combination of analyses that characterized multiple spatial length scales. Preferential modeling of long-range atomic arrangements with Rietveld refinement suggests the existence of primarily monointerstitials in UO 2.07 , whereas preferential modeling of short-range atomic structures with small-box pair distribution function (PDF) refinement indicates the presence of defect clusters in UO 2.07 . Simultaneous modeling of multiple length scales using complementary reverse Monte Carlo (RMC) and molecular dynamics (MD) methods confirms that excess oxygen atoms in UO 2.07 exist as small defects, such as monointerstitials and di-interstitials. RMC and MD results agree with diffraction analysis but differ significantly from small-box PDF refinements, which may be related to a lack of intermediate- and long-range structural information gained from the small-box PDF refinement procedure. Employing a combination of analysis methods with varying length-scale sensitivities enabled more accurate assessment of the UO 2+x defect structure. Furthermore, our findings provide experimental support for previously predicted di-interstitial defect morphologies in UO 2+x that highly influence the accurate prediction of bulk physiochemical properties of UO 2+x , such as oxygen diffusivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Electron transport in gaseous detectors with a Python-based Monte Carlo simulation code

Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (‘swarm’) properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Classical transport equations for burning gas-metal plasmas

Thermonuclear inertial confinement fusion plasmas confined by a heavy metal shell may be subject to the mixing of metal into the gas with a resulting degradation of fusion yield. Classical plasma diffusion driven by a number of gradients can provide a physical mechanism to produce atomic mix, possibly in concert with complex hydrodynamic structures and/or turbulence. This paper gives a derivation of the complete dissipative plasma hydrodynamics equations from kinetic theory, for a binary ionic mixture plasma consisting of electrons, e, a light (hydrogenic gas) ion species, i, and a heavy, high ZI plasma metal species, I. A single mean ionization state for the heavy metal, ZI, is assumed to be provided by some independent thermodynamic model of the heavy metal ZI=ZI(ni,nI,Te). The kinetic equations are solved by a generalized Chapman-Enskog expansion that assumes small Knudsen numbers for all species: NKe≡λe/L≪1,NKi≡λi/L≪1. The small electron to ion mass ratio, me/mi≪1, is utilized to account for electron-ion temperature separation, Te≠Ti, and to decouple the electron and ion transport coefficient calculations. This produces a well ordered perturbation theory for the electrons, resulting in the well known “Spitzer” problem of Spitzer and collaborators and solved independently by Braginskii. The formulation in this paper makes clear the inherent symmetry of the transport and gives an analytic solution for all values of the effective charge Zeff, including Zeff<1. The electron problem also determines the ambipolar electric field and the “thermal forces” on both ion species that are needed for the ion kinetic solution. The ion transport problem makes use of the small mass ratio between ion species, mi/mI≪1, to identify an “ion Spitzer problem” that is mathematically identical to that for the electrons but with different thermodynamic forces. The ionic scattering parameter, ΔI≡nIZI2/ni, replaces the Zeff of the electron problem, but has an extended domain, 0≤ΔI<∞, to cover all mixture fractions from the pure gas to the pure metal plasma. The extension of the Spitzer problem to include this extended domain is given in this work. The resulting transport equations for the binary gas-metal plasma mixture are complete and accurate through second order. All transport coefficients are provided in analytic form.

Molvig, Kim↗