SEARCH · Search NASA

Results for “bin mapping”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Conclusions and recommendations

From the analyses performed thus far, it is clear that the ship data must be treated with great care in any satellite comparisons because of their high noise level. At the very least, the ship data should be screened to eliminate observations deviating unreasonably from climatology. In fact, a two pass screening is better, with the second pass screening out data deviating unreasonably from the first pass mapping attempt. Even with screening, the noise level is comparable with that of the satellite data and, particularly in spot comparisons, should not be used to evaluate the accuracy of the satellite data. Consideration also should be given to more sophisticated mapping procedures than simple binning averaging. Procedures that attempt local least squares fitting to individual observations would give better results because they would then be less sensitive to the spatial distribution of the data.

Source record↗

NASA NTRS2014

Northern Late Winter Planetary Waves: MRO/MARCI Observations and Mars Climate Model Simulations

As does Earth, Mars presents pronounced global atmospheric circulation patterns. Solar differential heating drives mean meridional overturning (Hadley) circulations which are deep and intense, are hemispherically asymmetric, and where a cross-equatorial single cell dominates. Within middle and high latitudes, thermally indirect eddy-driven (Ferrel) circulation cells have been indicated. Differently, however, large-amplitude orography on planetary and continental scales on Mars can force very non-Earth-like hemispheric circulation patterns. Recent observations from the Mars Reconnaissance Orbiter, "Mars Color Imager" (MARCI) instrument are utilized that emphasize water ice clouds in ultra-violet (UV) wavelengths, and these measurements have been binned into "daily global maps" (DGMs) of water-ice cloud optical depth. The presence of large-scale, extratropical quasi-stationary atmospheric wave disturbances in middle and late winter of the northern hemisphere have been found to be present in such DGMs. In combination with such observations, a full-physics Mars global climate model (NASA ARC marsgcm 2.1) is applied to place the observations into context. During late northern winter, it is found that strong, forced Rossby modes (i.e., planetary waves) exist, and with direct correlation to columnintegrated cloud opacity undulating spatial patterns. At this season, zonal wavenumber s = 2 dominates (in contrast to wavenumber s = 1), consistent with MGS/TES analyses at this particular season (Banfield et al., 2003). Large-scale, planetary waves dictate the "coherence" of the northern polar vortex. Fundamentally, such forced planetary waves influence the polar vortex's impermeability (wave-induced) to tracer transport (e.g., dust and water-ice aerosol) and temporal mean water vapor spatial variations. The large-scale dynamical features of such planetary waves will be highlighted and discussed.

Northern polar vortex↗

DOE OSTI2022

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2016

Effect of Satellite Formations and Imaging Modes on Global Albedo Estimation

We confirm the applicability of using small satellite formation flight for multi-angular earth observation to retrieve global, narrow band, narrow field-of-view albedo. The value of formation flight is assessed using a coupled systems engineering and science evaluation model, driven by Model Based Systems Engineering and Observing System Simulation Experiments. Albedo errors are calculated against bi-directional reflectance data obtained from NASA airborne campaigns made by the Cloud Absorption Radiometer for the seven major surface types, binned using MODIS' land cover map - water, forest, cropland, grassland, snow, desert and cities. A full tradespace of architectures with three to eight satellites, maintainable orbits and imaging modes (collective payload pointing strategies) are assessed. For an arbitrary 4-sat formation, changing the reference, nadir-pointing satellite dynamically reduces the average albedo error to 0.003, from 0.006 found in the static reference case. Tracking pre-selected waypoints with all the satellites reduces the average error further to 0.001, allows better polar imaging and continued operations even with a broken formation. An albedo error of 0.001 translates to 1.36 W/sq m or 0.4% in Earth's outgoing radiation error. Estimation errors are found to be independent of the satellites' altitude and inclination, if the nadir-looking is changed dynamically. The formation satellites are restricted to differ in only right ascension of planes and mean anomalies within slotted bounds. Three satellites in some specific formations show average albedo errors of less than 2% with respect to airborne, ground data and seven satellites in any slotted formation outperform the monolithic error of 3.6%. In fact, the maximum possible albedo error, purely based on angular sampling, of 12% for monoliths is outperformed by a five-satellite formation in any slotted arrangement and an eight satellite formation can bring that error down four fold to 3%. More than 70% ground spot overlap between the satellites is possible with 0.5deg of pointing accuracy, 2 Km of GPS accuracy and commands uplinked once a day. The formations can be maintained at less than 1 m/s of monthly (Delta)V per satellite.

BRDF↗

NASA NTRS1982

Analyzing the Broken Ridge area of the Indian Ocean using magnetic and gravity anomaly maps and geoid undulation and bathymetry data

A higher resolution anomaly map of the Broken Ridge area (2 degree dipole spacing) was produced and reduced to the pole using quiet time data for this area. The map was compared with equally scaled maps of gravity anomaly, geoid undulation, and bathymetry. The ESMAP results were compared with a NASA MAGSAT map derived by averaging data in two-degree bins. A survey simulation was developed to model the accuracy of MAGSAT anomaly maps as a function of satellite altitude, instrument noise level, external noise model, and crustal anomaly field model. A preliminary analysis of the geophysical structure of Broken Ridge is presented and unresolved questions are listed.

Lazarewicz, A. R.↗

DOE OSTI2020

Planck intermediate results

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements include: maps that retain the CMB Solar dipole, allowing for high-precision relative calibration in higher-level analyses; well-defined single-detector maps, allowing for robust CO extraction; and HFI temperature maps between 217 and 857 GHz that are binned into 0'.9 pixels (Nside = 4096), ensuring that the full angular information in the data is represented in the maps even at the highest Planck resolutions. The net effect of these improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is (3366.6 ± 2.7) μK, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of τ = 0.051 ± 0.006, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

SPT-3G D1: Maps of the millimeter-wave sky from 2019 and 2020 observations of the SPT-3G Main field

Maps of the sky in millimeter wavelengths contain rich information on cosmology through anisotropies of the cosmic microwave background (CMB). Creating multifrequency sky maps of anisotropies in the $I$, $Q$, and $U$ Stokes parameters is one of the first steps of CMB cosmology analyses. In this work, we describe the production and validation of a set of sky maps from the South Pole Telescope's third-generation camera, SPT-3G. The maps are from data taken in frequency bands centered at 95, 150, and 220 GHz and taken during the first two years, 2019 and 2020, of the SPT-3G Main survey, which covers $4\%$ of the sky. We applied high-pass filters to time series of individual detectors and binned the filtered time series samples into map pixels. After that, we calibrated and cleaned the maps to reduce known systematic errors. In addition, we searched for other systematic errors through null tests and mitigated a significant systematic error detected therein. The white noise levels of the full-depth maps of the $I$ Stokes parameter are $5.4$, $4.4$, and $16.2$$\mathrm{μK}$-$\mathrm{arcmin}$ in the 95, 150, and 220 GHz bands, respectively, and $8.4$, $6.6$, and $25.8$$\mathrm{μK}$-$\mathrm{arcmin}$ for $Q/U$. These maps are the deepest to date used for measurements of mid-to-high-$\ell$ primary temperature and $E$-mode polarization CMB anisotropies, and reconstructions of the CMB gravitational lensing potential. We make these maps and supporting data products publicly accessible.

Quan, W. [Argonne (main); Chicago U., EFI; Chicago↗

NASA NTRS1995

Synoptic maps of solar wind parameters from in situ spacecraft observations

Solar wind observations from the Interplanetary Monitoring Platform-8 (IMP-8) and Pioneer Venus Orbiter (PVO) spacecraft from 1982 until 1988 are combined to construct synoptic maps of solar wind parameters near 1 AU. Each map consists of 6 months of hourly averaged solar wind data, binned by heliographic latitude and Carrington longitude and projected back to the Sun. These maps show the structure and time evolution of solar wind streams near 1 AU in the heliographic latitudes of +/- 7.25 deg and provide and explicit picture of several phenomena, such as gradients, changes in the inclination of the heliospheric current sheet, and the relative positions of various structures in the inner heliosphere, that is difficult to obtain from single-spacecraft observations. The stream structure varied significantly during the last solar cycle. Between 1982 and early 1985, solar wind parameters did not depend strongly on heliographic latitude. During the last solar minimum, the solar wind developed significant latitudinal structure, and high-speed streams were excluded from the vicinity of the solar equator. The interplanetary magnetic field was strongly correlated with the coronal field, and the current sheet tended to coincide with the coronal neutral line. The solar wind speed showed the expected correlations with temperature, interplanetary magnetic field, and distance from the current sheet. The solar wind speed was anticorrelated with density, but the regions of highest density occurred east of the heliospheric current sheet and the regions of lowest solar wind speed. This is consistent with compression at the leading edge of high-speed streams.

Gazis, P. R.↗

DOE OSTI2022

Constraining low redshift [C II ] emission by cross-correlating FIRAS and BOSS data

ABSTRACT We perform a tomographic cross-correlation analysis of archival FIRAS data and the BOSS galaxy redshift survey to constrain the amplitude of [C II] 2P3/2 → 2P1/2 fine structure emission. Our analysis employs spherical harmonic tomography (SHT), which is based on the angular cross-power spectrum between FIRAS maps and BOSS galaxy over-densities at each pair of redshift bins, over a redshift range of 0.24 < z < 0.69. We develop the SHT approach for intensity mapping, where it has several advantages over existing power spectral estimators. Our analysis constrains the product of the [C II] bias and [C II] specific intensity, $b_{\rm [C \small{\rm II}]}I_{\rm [C \small{\rm II}]}$, to be <0.31 MJy/sr at z ≈ 0.35 and <0.28 MJy/sr at z ≈ 0.57 at $95{{\ \rm per\ cent}}$ confidence. These limits are consistent with most current models of the [C II] signal, as well as with higher-redshift [C II] cross-power spectrum measurements from the Planck satellite and BOSS quasars. We also show that our analysis, if applied to data from a more sensitive instrument such as the proposed PIXIE satellite, can detect pessimistic [C II] models at high significance.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Cosmology from Planck CMB lensing and DESI DR1 quasar tomography

We present a measurement of the amplitude of matter fluctuations over the redshift range 0.8 ≤ z ≤ 3.5 from the cross correlation of over 1.2 million spectroscopic quasars selected by the Dark Energy Spectroscopic Instrument (DESI) across 7,200 deg 2 (∼ 170 quasars/deg 2 ) and Planck PR4 (NPIPE) cosmic microwave background (CMB) lensing maps. We perform a tomographic measurement in three bins centered at effective redshiftsz=1.44, 2.27 and 2.75, which have ample overlap with the CMB lensing kernel. From a joint fit using the angular clustering of all three redshift bins (auto and cross-spectra), and including an Q m prior from DESI DR1 baryon acoustic oscillations to break the $Ω_{m}-σ_{8}$ degeneracy, we constrain the amplitude of matter fluctuations in the matter-dominated regime to be $σ_{8}=0.929^{+0.059}_{-0.074}$ and $S_{8}≡σ_{8}(Ω_m/0.3)^{0.5} = 0.922^{+0.059}_{-0.073}$. We provide a growth of structure measurement with the largest spectroscopic quasar sample to date at high redshift, which is ∼ 1.5σ higher than predictions from ΛCDM fits to measurements of the primary CMB from Planck PR4. The cross-correlation between PR4 lensing maps and DESI DR1 quasars is detected with a signal-to-noise ratio of 21.7 and the quasar auto-correlation at 27.2 for the joint analysis of all redshift bins. We combine our measurement with the CMB lensing auto-spectrum from the ground-based Atacama Cosmology Telescope (ACT DR6) and Planck PR4 to perform a sound-horizon-free measurement of the Hubble constant, yielding $H_{0}=69.1^{+2.2}_{-2.6} {\rm \ km\ s}^{-1}{\rm Mpc}^{-1}$ through its sensitivity to the matter-radiation equality scale.

baryon acoustic oscillations↗

DOE OSTI2025

Dark Energy Survey Year 3 results: $w$CDM cosmology from simulation-based inference with persistent homology on the sphere

We present cosmological constraints from Dark Energy Survey Year 3 (DES Y3) weak lensing data using persistent homology, a topological data analysis technique that tracks how features like clusters and voids evolve across density thresholds. For the first time, we apply spherical persistent homology to galaxy survey data through the algorithm TopoS2, which is optimized for curved-sky analyses and HEALPix compatibility. Employing a simulation-based inference framework with the Gower Street simulation suite, specifically designed to mimic DES Y3 data properties, we extract topological summary statistics from convergence maps across multiple smoothing scales and redshift bins. After neural network compression of these statistics, we estimate the likelihood function and validate our analysis against baryonic feedback effects, finding minimal biases (under $0.3σ$) in the $Ω_\mathrm{m}-S_8$ plane. Assuming the $w$CDM model, our combined Betti numbers and second moments analysis yields $S_8 = 0.821 \pm 0.018$ and $Ω_\mathrm{m} = 0.304\pm0.037$-constraints 70% tighter than those from cosmic shear two-point statistics in the same parameter plane. Our results demonstrate that topological methods provide a powerful and robust framework for extracting cosmological information, with our spherical methodology readily applicable to upcoming Stage IV wide-field galaxy surveys.

Prat, J. [Nordita; Royal Inst. Tech., Sodertalje; ↗

DOE OSTI2024

The star formation, dust, and abundance of galaxies with unWISE -CIB cross-correlations

Here, the cosmic infrared background (CIB) is the accumulated infrared (IR) radiation mainly from interstellar dust heated up by early stars. In this work, we measure the cross-correlation between galaxies from the unWISE catalog and the CIB maps from the Planck satellite to simultaneously constrain the cosmic star formation rate (SFR), dust spectral energy distribution (SED), and the halo occupation distribution (HOD). The unWISE galaxy catalog is divided into three tomographic bins centered at z ~ 0.6, 1.1, 1.5, and the CIB maps are at 353, 545, and 857 GHz. We measure the cross-correlations between these galaxy samples and CIB maps and get a 194σ signal within an angular scale 100 < ℓ < 2000, from which we constrain two CIB halo models from previous literature and one new model. The SFR, SED, and HOD model parameters are constrained consistently among the three models. Specifically, the dust temperature at z = 0 is constrained T 0 = $21.14^{+1.02}_{-1.34}$ K, which is slightly lower than T 0 = 24.4±1.9 K measured by the Planck collaboration. The halo mass that gives the most efficient star formation is around 10 11.79 $^{+0.73}_{-0.86}$ M ⊙ . From the model parameters, combined with the SFR density at z = 0 synthesized from multi-wavelength observations, we break the degeneracy between SED and SFR and recover the cosmic star formation history that is consistent with multi-wavelength surveys. We also constrain the graybody SED model in agreement with previous measurements from infrared flux stacking. From the HOD constraints, we derive an increasing trend of galaxy linear bias along redshifts that agrees with the results from cross- and auto-correlation with unWISE galaxies. This study indicates the power of using CIB-galaxy cross-correlation to study star formation, dust, and abundance of galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Introducing a New Spitzer Master BLISS Map to Remove the Instrument Systematic Phase-curve-parameter Degeneracy, as Demonstrated by a Reanalysis of the 4.5 μm WASP-43b Phase Curve

While Spitzer Infrared Array Camera (IRAC) systematics are generally well understood, each data set can provide its own challenges that continue to teach us about the underlying functional form of these systematics. Multiple groups have analyzed the phase curves of WASP-43b with varying detrending techniques, each obtaining different results. In this work, we take another look at WASP-43b while exploring the degenerate relation between Bilinearly Interpolated Subpixel Sensitivity (BLISS) mapping, point response function (PRF)–FWHM detrending, and phase curve parameters. We find that there is a strong correlation between the detrending parameters in the two models, and best-fit phase curve amplitudes vary strongly when the data are temporally binned. To remove this degeneracy, we present a new Gaussian centroided intrapixel sensitivity map (hereafter fixed sensitivity map), generated using 3,712,830 exposures spanning 5 yr, for a variety of aperture sizes at 4.5 μm. We find evidence for time variability in the sensitivity at 3.6 μm and do not generate a visit-independent map for this channel. With the fixed 4.5 μm intrapixel sensitivity map, the best fits for WASP-43b no longer vary strongly with bin size and PRF–FWHM detrending is no longer required to remove correlated noise. For data sets that do not fall completely within the sweet spot, temporal binning should not be used in the analysis of Spitzer phase curves. We confirm nightside emission for WASP-43b with a disk-integrated nightside temperature of 806 ± 48 K at 4.5 μm. The 4.5 μm maps are available at github.com/kevin218/POET.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

Alternate Inversion Paradigm for Spectroheliogram Data

Over the past five years, new methods to reconstruct spectrally pure maps of the Sun from spectroheliogram data have emerged, essentially unlocking this long-abandoned method of obtaining both spatial and spectral information over a large field of view simultaneously. The original inversion method determined the plasma’s emission measure distribution as a function of temperature at every spatial location in the field of view. To complete this inversion, a response matrix had to be created mapping the emission measure in different (temperature, space) bins to detector, requiring assumptions on the thermal and ionization equilibrium and abundance state of the plasma. We have since derived a new method of the inversion that does not require these atomic assumptions to be made. Instead, we use only the different locations of spectral lines from the same ion species and the possible ratios of those single-species spectral lines, removing the need for a priori knowledge on the state of the emitting plasma. In this presentation, we demonstrate this method using observed data from the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) and simulated data from the EUV CME and Coronal Connectivity Observatory (ECCCO) investigation.

Amy Winebarger↗

DOE OSTI2021

Towards a reliable assessment of charging effects during surface analysis: Accurate spectral shapes of ZrO 2 and Pd/ZrO 2 via X-ray Photoelectron Spectroscopy

X-ray Photoelectron Spectroscopy of large bandgap or insulating material surfaces relies on an effective mechanism that compensates for the emission (loss) of electrons by maintaining the material surface at a steady-state uniform potential. While a steady-state may be attained by utilizing an active compensation, such as low power electron emitting filament, there is the possibility that the surface potential is not uniform over the area analyzed, leading to peak shifts and incorrect spectral interpretation. Here, in this work, a spectral data processing method based on mapping the ZrO 2 and Pd/ZrO 2 surfaces utilizing photoemission peak binding energy is proposed, which provides information about the response of specific material surfaces to charge compensation. Spectromicroscopy of ZrO 2 and Pd/ZrO 2 surfaces without spatial information is used to monitor the efficacy of charge compensation. Exploiting counts distributed over many bins require the use of procedures and algorithms essential to practical mapping peak positions. Iterative singular value decomposition is therefore introduced and utilized as a means of efficiently delivering spatially resolved spectra from which binding energy for peaks is computed. The concepts developed in this work result in robust and accurate peak models of ZrO 2 and Pd/ZrO 2 that can be applied in XPS analysis of not only ZrO 2 but other large bandgap or insulating material surfaces. Supporting arguments for a peak model representing signal from Zr 3p and Pd 3d are developed within this work are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2024

The anomaly of the CMB power with the latest Planck data

Abstract The lack of power anomaly is an unexpected feature observed at large angular scales in the maps of Cosmic Microwave Background (CMB) produced by the COBE, WMAP andPlancksatellites. This signature, which consists in a missing of power with respect to that predicted by the ΛCDM model, might hint at a new cosmological phase before the standard inflationary era.The main point of this paper is taking into account the latestPlanckpolarisation data to investigate how the CMB polarisation improves the understanding of this feature. With this aim, we apply to the latestPlanckdata, both PR3 (2018) and PR4 (2020) releases, a new class of estimators capable of evaluating this anomaly by considering temperature and polarisation data both separately and in a jointly way. This is the first time that the PR4 dataset has been used to study this anomaly. To critically evaluate this feature, taking into account the residuals of known systematic effects present in thePlanckdatasets, we analyse the cleaned CMB maps using different combinations of sky masks, harmonic range and binning on the CMB multipoles.Our analysis shows that the estimator based only on temperature data confirms the presence of a lack of power with a lower-tail-probability (LTP), depending on the component separation method, ≤ 0.33% and ≤ 1.76% for PR3 and PR4, respectively. To our knowledge, the LTP≤ 0.33% for the PR3 dataset is the lowest one present in the literature obtained fromPlanck2018 data, considering thePlanckconfidence mask. We find significant differences between these two datasets when polarisation is taken into account most likely due to a different level of systematics. Especially, the analysis with PR3 data, unlike that with PR4, seems to point towards a lack of power at large scales also for polarisation.Moreover, we also show that for the PR3 dataset the inclusion of the subdominant polarisation information provides estimates that are less likely accepted in a ΛCDM cosmological model than the only-temperature analysis over the entire harmonic-range considered. In particular, at ℓ max = 26, we found that no simulation has a value as low as the data for all the pipelines.

Astronomy & Astrophysics↗

DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗