Results for “constrained variational analysis” — Page 17

DOE OSTI2025

A Review of Quantum Computing Technologies in Power System Optimization

As modern power grids increasingly integrate variable renewable generation, distributed energy resources, and energy storage systems, classical optimization techniques are facing unprecedented challenges. This review examines the emerging application of quantum computing to overcome these challenges in power system optimization, including optimal power flow (OPF), unit commitment (UC), economic dispatch (ED), and intelligent switching and topology optimization (IS-TO). Recent research has introduced various quantum methodologies—such as gate-based, annealing-based, variational algorithms, and quantum-inspired algorithms—to address the combinatorial complexity inherent in grid reconfiguration and energy management. The review summaries the quantum algorithms, quantum devices and the power system test cases, highlighting hybrid quantum–classical strategies that leverage the complementary strengths of both paradigms. Some quantum advantages have been observed, including theoretical speedup, accurate simulation results, scalable qubit usage, efficient QUBO mapping. In particular, the review emphasizes the importance of integrating quantum optimization techniques with classical control frameworks, these hybrid approaches demonstrate the potential to improve real-time grid management and operational reliability. A significant portion of the analysis is devoted to the practical limitations of current quantum devices. Present-day quantum hardware, operating in the noisy intermediate-scale quantum (NISQ) era, remains highly sensitive to noise and limited in qubit connectivity, which constrains the scale and accuracy of implemented algorithms. The review delves into specific challenges such as the need for qubit-efficient encoding techniques and error mitigation strategies that are critical for handling real-world grid optimization problems. In addition, the work draws attention to the performance discrepancies between theoretical quantum speedups and experimental validations, underscoring the importance of rigorous benchmark studies using representative power grid test cases. In summary, this review highlights both the promise and limitations of quantum computing for power system optimization. It provides a comprehensive overview of the state-of-the-art technologies, categorizes recent advancements in algorithm design, and discusses practical considerations for implementation, and serves as an informative resource on current research. Future research directions include developing robust hybrid frameworks, advancing qubit-efficient formulations, and scaling up experimental demonstrations to confirm the theoretical advantages of quantum methods in large-scale power system operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Visualisation and outlier detection for probability density function ensembles

Abstract Exploratory data analysis (EDA) for functional data—data objects where observations are entire functions—is a difficult problem that has seen significant attention in recent literature. This surge in interest is motivated by the ubiquitous nature of functional data, which are prevalent in applications across fields such as meteorology, biology, medicine and engineering. Empirical probability density functions (PDFs) can be viewed as constrained functional data objects that must integrate to one and be nonnegative. They show up in contexts such as yearly income distributions, zooplankton size structure in oceanography and in connectivity patterns in the brain, among others. While PDF data are certainly common in modern research, little attention has been given to EDA specifically for PDFs. In this paper, we extend several methods for EDA on functional data for PDFs and compare them on simulated data that exhibit different types of variation, designed to mimic that seen in real‐world applications. We then use our new methods to perform EDA on the breakthrough curves observed in gas transport simulations for underground fracture networks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Weak Lensing by Galaxy Clusters: from Pixels to Cosmology

The story of the origin and evolution of our Universe is told, equivalently, by space-time itself and by the structures that grow inside of it. Clusters of galaxies are the frontier of bottom-up structure formation. They are the most massive objects to have collapsed at the present epoch. By that virtue, their abundance and structural parameters are highly sensitive to the composition and evolution of the Universe. The most common probe of cluster cosmology, abundance, uses samples of clusters selected by some observable. Applying a mass-observable relation (MOR), cosmological parameters can be constrained by comparing the sample to predicted cluster abundances as a function of observable and redshift. Arguably, however, cluster probes have not yet entered the era of per cent level precision cosmology. The primary reason for this is our imperfect understanding of the MORs. The overall normalization, the slope of mass vs. observable, the redshift evolution, and the degree and correlation of intrinsic scatters of observables at fixed mass have to be constrained for interpreting abundances correctly. Mass measurement of clusters by means of the differential deflection of light from background sources in their gravitational field, i.e. weak lensing, is a powerful approach for achieving this. This thesis presents new methods for and scientific results of weak lensing measurements of clusters of galaxies. The former include, on the data reduction side, (i) the correction of CCD images for non-linear effects due to the electric fields of accumulated charges and (ii) a method for masking artifact features in sets of overlapping images of the sky by comparison to the median image. Also, (iii) I develop a method for the selection of background galaxy samples based on their color and apparent magnitude that includes a new correction for contamination with cluster member galaxies. The main scientific results are the following. (i) For the Hubble Frontier Field cluster RXC J2248.7--4431 our lensing analysis constrains mass and concentration of the cluster halo and we confirm the large mass predicted by X-ray and Sunyaev-Zel’dovich (SZ) observations. The study of cluster members shows the relation of galaxy morphology to luminosity and environment. (ii) Our lensing mass measurements for 12 clusters are consistent with X-ray masses derived under the assumption of hydrostatic equilibrium of the intra-cluster gas. We confirm the MORs derived by the South Pole Telescope collaboration for the detection significance of the cluster SZ signal in their survey. We find discrepancies, however, with the Planck SZ MOR. We hypothesize that these are related either to a shallower slope of the MOR or a size-, redshift- or noise-dependent bias in SZ signal extraction. (iii) Finally, using a combination of simulations and theoretical models for the variation of cluster profiles at fixed mass, we find that the latter is a significant contribution to the uncertainty of cluster lensing mass measurements. A cosmic variance model, such as the one we develop, is necessary for MOR constraints to be accurate at the level required for future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1989

High-resolution microwave images of Saturn

An analysis of high-resolution microwave images of Saturn and Saturn's individual rings is presented. Radio interferometric observations of Saturn taken at the Very Large Array in New Mexico at wavelengths of 2 and 6 centimeters reveal interesting new features in both the atmosphere and rings. The resulting maps show an increase in brightness temperature of about 3 K from equator to pole at both wavelengths, while the 6-centimeter map shows a bright band at northern midlatitudes. The data are consistent with a radiative transfer model of the atmosphere that constrains the well-mixed, fully saturated, NH3 mixing ratio to be 0.00012 in a region just below the NH3 clouds, while the observed bright band indicates a 25 percent relative decrease of NH3 in northern midlatitudes. Brightness temperatures for the classical rings are presented. Ring brightness shows a variation with azimuth and is linearly polarized at an average value of about 5 percent. The variations in ring polarization suggest that at least 20 percent of the ring brightness is the result of a single scattering process.

Grossman, A. W.↗

DOE OSTI2025

Aerosol hygroscopicity over the South-East Atlantic Ocean during the biomass burning season – Part 2: Influence of sea salt and burning conditions on CCN hygroscopicity

Biomass burning (BB) significantly influences cloud condensation nuclei (CCN) concentrations over the southeastern Atlantic; however, aerosol hygroscopicity (κ) – a key factor for CCN activation – remains poorly constrained during the BB season. This study investigates κ variability using in situ measurements from Ascension Island during the 2016 and 2017 BB seasons. Results show substantial monthly variability, with κ values lowest in August and increasing through October. On average, κ was significantly higher in 2017 (∼ 0.55) than in 2016 (∼ 0.33), suggesting that the aerosols in 2017 were more hygroscopic and more easily activated as CCN. Sulfate and sea salt were the two dominant contributors to κ and the primary drivers of its interannual variability. During the 2017 BB season, sulfate – the major inorganic component – accounted for ∼ 34 % of the submicron aerosol mass, while sea salt, estimated via κ-closure analysis, contributed ∼ 17 %. The higher κ in 2017 was largely attributed to increased sea salt, likely driven by stronger marine winds. Approximately 67 % of sulfate was linked to BB emissions. Variations in BB combustion efficiency, modulated by regional meteorology, influenced sulfate fraction and thus κ values. Specifically, higher relative humidity and lower wind speeds over BB source regions in 2017 favored smoldering combustion, explaining the higher sulfate fraction. Overall, the observed interannual differences in aerosol hygroscopicity reflect the combined impacts of BB combustion characteristics and sea salt emissions, underscoring the critical roles of both BB and marine aerosol sources in regulating aerosol-cloud interactions over the southeastern Atlantic.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Analysis of coseismic surface displacement gradients using radar interferometry: New insights into the Landers earthquake

The map of the coseismic displacement field generated by interferometric processing of synthetic aperture radar (SAR) images taken before and after the June 28, 1992, Landers earthquake sequence brings new insights into the nature of deformation caused by these earthquakes. We use the interferometric map generated by Massonnet et al. (1993) to analyze the surface displacement field in the vicinity of the fault trace. Complexities in the fringe pattern near the fault reflect short-wavelength variations of the surface rupture and slip distribution, and attest to large displacement gradients. Along two sections of the fault, characteristic fringe patterns can be recognized, contrasting in density and direction with patterns observed away from the rupture. In order to understand the observed fringe patterns, we compute synthetic interferograms in three simple cases: (1) rigid-body rotations about a vertical axis, (2) about a horizontal axis (tilt), and (3) distributed, simple shear. The orientation and spatial separation of interferometric fringes predicted by these models help constrain near-field deformation and rupture parameters. Where the Kickapoo fault connects with the Homestead Valley fault, the interferogram shows a clear pattern of parallel N20 deg W fringes separated by about 160 m. This pattern and vertical offsets measured along the Kickapoo fault suggest that the block between this fault and the Johnson Valley fault may have been tilted, down to the west. A 5-km block lifted by 1 m on one side would be tilted by an angle of 0.01 deg (190 microrad), producing fringes separated by about 160 m, parallel to the tilt axis. Such a tilt, parallel to a N20 deg W direction, would account for the gradual, northward increase of the vertical slip component observed along the Kickapoo fault. This tilt may also explain the 1 m of reverse slip observed along the 'slip gap' section of the Homestead Valley break. Between the southern end of the Johnson Valley fault and the Eureka Peak fault, where no surface rupture has been mapped, the dense pattern of fringes implies distributed shear, probably resulting from fault slip at depth. The density and direction of the fringes in the gap are consistent with a right-lateral slip of 1.2-3.8 m on a blind fault locked above the depth of 1.5-2 km. Such observations of small wavelength features in the SAR interferogram bring new insights into the near-field displacement gradient and thus on response of the uppermost crust to seismic rupture.

Peltzer, Gilles↗

DOE OSTI2013

COSMIC SHEAR RESULTS FROM THE DEEP LENS SURVEY. I. JOINT CONSTRAINTS ON {Omega} {sub M} AND {sigma}{sub 8} WITH A TWO-DIMENSIONAL ANALYSIS

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 deg{sup 2} fields with two National Optical Astronomy Observatory (NOAO) 4 m telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the {Omega} {sub M}-{sigma}{sub 8} likelihood contour tightly, providing a joint constraint of {Omega} {sub M} = 0.262 {+-} 0.051 and {sigma}{sub 8} = 0.868 {+-} 0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the {Omega} {sub M}-{sigma}{sub 8} degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) likelihood data, we obtain {Omega} {sub M} = 0.278 {+-} 0.018 and {sigma}{sub 8} = 0.815 {+-} 0.020.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2022

Inferring iron-oxide species content in atmospheric mineral dust from DSCOVR EPIC observations

The iron-oxide content of dust in the atmosphere and most notably its apportionment between hematite (α-Fe2O3) and goethite (α-FeOOH) are key determinants in quantifying dust’s light absorption, its top of atmosphere UV radiances used for dust monitoring, and ultimately shortwave dust direct radiative effects (DRE). Hematite and goethite column mass concentrations and iron-oxide mass fractions of total dust mass concentration were retrieved from the DeepSpace Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) measurements in the ultraviolet–visible (UV–Vis) channels. The retrievals were performed for dust-identified aerosol plumes over land using aerosol optical depth (AOD) and spectral imaginary refractive index provided by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm over six continental regions (North America, North Africa, West Asia, Central Asia, EastAsia, and Australia). The dust particles are represented as an internal mixture of non-absorbing host and absorbing hematite and goethite. We use the Maxwell–Garnett effective medium approximation with carefully selected complex refractive indices of hematite and goethite that produce mass fractions of iron oxides species consistent with in situ values found in the literature to derive the hematite and goethite volumetric/mass concentrations from MAIAC EPIC products. We compared the retrieved hematite and goethite concentrations with in situ dust aerosol mineralogical content measurements, as well as with published data. Our data display variations within the published range of hematite, goethite, and iron-oxide mass fractions for pure mineral dust cases. A specific analysis is presented for 15 sites over the main dust source regions. Sites in the central Sahara, Sahel, and Middle East exhibit a greater temporal variability of iron oxides relative to other sites. Niger site(13.52°N, 2.63°E) is dominated by goethite over Harmattan season with median of ~2 weight percentage (wt.%) of iron-oxide. Saudi Arabia site (27.49°N, 41.98°E) over Middle East also exhibited surge of goethite content with the beginning of Shamal season. The Sahel dust is richer in iron-oxide than Saharan and northern China dust except in Summer. The Bodélé Depression area shows a distinctively lower iron-oxide concentration (~1 wt.%) throughout the year. Finally, we show that EPIC data allow to constrain the hematite refractive index. Specifically, we select 5 out of 13 different number of hematite refractive indices widely variable in published laboratory studies by constraining the iron-oxide mass ratio to the known measured values. Provided climatology of hematite and goethite mass fractions across main dust regions of the Earth will be useful for dust shortwave DRE studies and climate modeling.

Sujung Go↗

NASA NTRS1989

Spaceborne radar observations: A guide for Magellan radar-image analysis

Geologic analyses of spaceborne radar images of Earth are reviewed and summarized with respect to detecting, mapping, and interpreting impact craters, volcanic landforms, eolian and subsurface features, and tectonic landforms. Interpretations are illustrated mostly with Seasat synthetic aperture radar and shuttle-imaging-radar images. Analogies are drawn for the potential interpretation of radar images of Venus, with emphasis on the effects of variation in Magellan look angle with Venusian latitude. In each landform category, differences in feature perception and interpretive capability are related to variations in imaging geometry, spatial resolution, and wavelength of the imaging radar systems. Impact craters and other radially symmetrical features may show apparent bilateral symmetry parallel to the illumination vector at low look angles. The styles of eruption and the emplacement of major and minor volcanic constructs can be interpreted from morphological features observed in images. Radar responses that are governed by small-scale surface roughness may serve to distinguish flow types, but do not provide unambiguous information. Imaging of sand dunes is rigorously constrained by specific angular relations between the illumination vector and the orientation and angle of repose of the dune faces, but is independent of radar wavelength. With a single look angle, conditions that enable shallow subsurface imaging to occur do not provide the information necessary to determine whether the radar has recorded surface or subsurface features. The topographic linearity of many tectonic landforms is enhanced on images at regional and local scales, but the detection of structural detail is a strong function of illumination direction. Nontopographic tectonic lineaments may appear in response to contrasts in small-surface roughness or dielectric constant. The breakpoint for rough surfaces will vary by about 25 percent through the Magellan viewing geometries from low to high Venusian latitudes. Examples of anomalies and system artifacts that can affect image interpretation are described.

Ford, J. P.↗

DOE OSTI2016

The Southern Ocean as a constraint to reduce uncertainty in future ocean carbon sinks

Earth system model (ESM) simulations exhibit large biases compares to observation-based estimates of the present ocean CO 2 sink. The inter-model spread in projections increases nearly 2-fold by the end of the 21st century and therefore contributes significantly to the uncertainty of future climate projections. In this study, the Southern Ocean (SO) is shown to be one of the hot-spot regions for future uptake of anthropogenic CO 2 , characterized by both the solubility pump and biologically mediated carbon drawdown in the spring and summer. Here, we show, by analyzing a suite of fully interactive ESMs simulations from the Coupled Model Intercomparison Project phase 5 (CMIP5) over the 21st century under the high-CO 2 Representative Concentration Pathway (RCP) 8.5 scenario, that the SO is the only region where the atmospheric CO 2 uptake rate continues to increase toward the end of the 21st century. Furthermore, our study discovers a strong inter-model link between the contemporary CO 2 uptake in the Southern Ocean and the projected global cumulated uptake over the 21st century. This strong correlation suggests that models with low (high) carbon uptake rate in the contemporary SO tend to simulate low (high) uptake rate in the future. None the less, our analysis also shows that none of the models fully capture the observed biophysical mechanisms governing the CO 2 fluxes in the SO. The inter-model spread for the contemporary CO 2 uptake in the Southern Ocean is attributed to the variations in the simulated seasonal cycle of surface p CO 2 . Two groups of model behavior have been identified. The first one simulates anomalously strong SO carbon uptake, generally due to both too strong a net primary production and too low a surface p CO 2 in December–January. The second group simulates an opposite CO 2 flux seasonal phase, which is driven mainly by the bias in the sea surface temperature variability. Furthermore, we show that these biases are persistent throughout the 21st century, which highlights the urgent need for a sustained and comprehensive biogeochemical monitoring system in the Southern Ocean to better constrain key processes represented in current model systems.

58 GEOSCIENCES↗

DOE OSTI2023

Third-order Hall effect in the surface states of a topological insulator

Time reversal and inversion symmetric materials fail to yield linear and nonlinear responses since they possess net zero Berry curvature. However, higher-order Hall response can be generated in these systems upon constraining the crystalline symmetries. Motivated by the recently discovered third-order Hall (TOH) response mediated by Berry connection polarizability, namely, the variation of the Berry connection with respect to an applied electric field, here, we investigate the existence of such a Hall effect in the surface states of hexagonal warped topological insulators (e.g., Bi 2 Te 3 ) under the application of only electric field. Using the semiclassical Boltzmann formalism, we investigate the effect of tilt and hexagonal warping on the Berry connection polarizability tensor and, consequently, the TOH effect, provided the Dirac cone remains gapless. We find that the magnitude of the response increases significantly with increasing tilt strength and warping, and therefore, they can provide the tunability of this effect. In addition, we also explore the effect of chemical doping on the TOH response in this system. Interestingly, we show, based on the symmetry analysis, that the TOH can be the leading-order response in this system, which can directly be verified in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Highly constrained ferroelectric [BaTiO{sub 3}]{sub (1−x)Λ}/[BaZrO{sub 3}]{sub xΛ} superlattices: X-ray diffraction and Raman spectroscopy

We report an x-ray diffraction (XRD) and a Raman-scattering investigation of ferroelectric/paraelectric superlattices [BaTiO{sub 3}] {sub (1−x)Λ}/[BaZrO{sub 3}]{sub xΛ} for which the composition varied, 0.15 ≤ x ≤ 0.85, while the superlattice (SL) modulation period Λ was kept constant at about 100 Å. The samples were epitaxially grown by pulsed laser deposition on MgO substrates buffered with La{sub 0.5}Sr{sub 0.5}CoO{sub 3}. Based on the XRD analysis and on polarized Raman spectra, we have showed that the large strain in SLs induced ferroelectricity in BaZrO{sub 3} (BZ) for all SLs, a material that is paraelectric in the bulk form at any temperature and in the single film. The induced polar axis in BZ layers is perpendicular to the plane of substrate while BaTiO{sub 3} (BT) layers exhibit in-plane polar orientation. Raman spectroscopy revealed a lattice ordering in SLs due to the misfit strain generated by the large lattice mismatch between the alternating BZ and BT layers. This strain induced a huge upward frequency of the lowest E(1TO) soft mode from 60 cm{sup −1} in the BT single film to 215 cm{sup −1} in the SL with x = 0.85. These results show that in spite of relatively large periodicity of SLs, they are highly constrained and the variation of BZ ratio allowed modifying strains between layers. The temperature dependence of the Raman spectra for BT{sub 0.3Λ}/BZ{sub 0.7Λ} and BT{sub 0.7Λ}/BZ{sub 0.3Λ} samples revealed giant shift of the ferroelectric phase transition. The phase transition temperature was found to be upshifted by about 300 °C with respect to BT single crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2018

Retrieval of H2O Abundance in Titan's Stratosphere: A (Re)Analysis of CIRS/Cassini and PACS/Herschel Observations

Since its first measurement 20 years ago by the Infrared Space Observatory (ISO), the water (H2O) mole fraction in Titan’s stratosphere remains uncertain due to large differences between the determinations from available measurements. More particularly, the recent measurements made from the Herschel observatory (PACS (Photodetector Array Camera and Spectrometer) and HIFI (Heterodyne Instrument for the Far Infrared)) estimated the H2O mole fraction to be 0.023 parts per billion at 12.1 millibars. A mixing ratio of 0.14 parts per billion at 10.7 millibars was, however, retrieved from nadir spatially-resolved observations of Cassini/CIRS (Composite Infrared Spectrometer). At the same pressure level (10.7 millibars), this makes a difference of a factor of 5.5 between PACS and CIRS measurements, and this has notably prevented current models from fully constraining the oxygen flux flowing into Titan’s atmosphere. In this work, we try to understand the differences between the H2O mole fractions estimated from Herschel/PACS and Cassini/CIRS observations. The strategy for this is to 1) analyse recent disc-averaged observations of CIRS to investigate if the observation geometry could explain the previous observed differences, and 2) (re)analyse the three types of observation with the same retrieval scheme to assess if previous differences in retrieval codes/methodology could be responsible for the previous discrepancies. With this analysis, we show that using the same retrieval method better reconcile the previous measurements of these instruments. However, the addition of the disc-averaged CIRS observations, instead of confirming the consistency between the different datasets, reveals discrepancies between one of the CIRS disc-averaged set of observations and PACS measurements. This raises new questions regarding the possibility of latitudinal variations of H2O, which could be triggered by seasonal changes of the meridional circulation. As it has already been shown for nitriles and hydrocarbons, this circulation could potentially impact the latitudinal distribution of H2O through the subsidence or up- welling of air rich in H2O. The possible influence of spatial/time variations of the OH/H2O input flux in Titan’s atmosphere is also discussed. The analysis of more observations will be needed in future work to address the questions arising from this work and to improve the understanding of the sources of H2O in Titan’s atmosphere.

Bauduin, S.↗

DOE OSTI2023

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

advanced gas reactors↗

DOE OSTI2023

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2009

Electron Backscatter Diffraction (EBSD) Analysis and U-Pb Geochronology of the Oldest Lunar Zircon: Constraining Early Lunar Differentiation and Dating Impact-Related Deformation

The evolution of the early moon was dominated by two processes (i) crystallization of the Lunar Magma Ocean (LMO) and differentiation of potassium-rare earth element-phosphorous-rich residual magma reservoir referred to as KREEP, and (ii) an intense meteorite bombardment referred to as lunar cataclysm . The exact timing of these processes is disputed, and resolution relies on collection and interpretation of precise age data. This study examines the microstructure and geochronology of zircon from lunar impact breccias collected during the Apollo 17 mission. A large zircon clast within lunar breccia 72215,195 shows sector zoning in optical microscopy, cathodoluminescence (CL) imaging and Raman mapping, and indicates that it was a relict fragment of a much larger magmatic grain. Sensitive high resolution ion microprobe (SHRIMP) U-Pb analysis of the zircon shows that U and Th concentration correlate with sector zoning, with darkest CL domains corresponding with high-U and Th (approx.150 and approx.100 ppm respectively), and the brightest-CL sectors containing approx.30-50 ppm U and approx.10-20 ppm Th. This indicates that variations in optical CL and Raman properties correspond to differential accumulation of alpha-radiation damage in each sector. Electron backscatter diffraction (EBSD) mapping shows that the quality of electron backscatter patterns (band contrast) varies with sector zoning, with the poorest quality patterns obtained from high-U and Th, dark-CL zones. EBSD mapping also reveals a deformation microstructure that is cryptic in optical, CL and Raman imaging. Two orthogonal sets of straight discrete and gradational low-angle boundaries accommodate approx.12 misorientation across the grain. The deformation bands are parallel to the crystallographic {a}-planes of the zircon, have misorientation axes parallel to the c-axis, and are geometrically consistent with formation by dislocation creep associated with <100>{010} slip. The deformation bands are unlike curved morphology of crystal-plastic microstructures in tectonically deformed terrestrial zircon, and geometrically similar to dislocation microstructures reported in experimentally shocked zircon. We interpret these crystal-plastic deformation microstructures to have resulted from a significant impact, either directly from impact shock, or during ductile flow directly following the impact. The deformation bands appear to continue undeflected through the non-indexed, radiation-damaged areas of the grain, which suggests that the orientation variation predates any significant mechanical weakening from radiation damage in the grain, and therefore occurred early in its history.

Timms, Nick↗