Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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 685 records · Page 38

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗

Experimental and Modeled Assessment of Interventions to Reduce PM2.5 in a Residence during a Wildfire Event

Increasingly large and frequent wildfires affect air quality even indoors by emitting and dispersing fine/ultrafine particulate matter known to pose health risks to residents. With this health threat, we are working to help the building science community develop simplified tools that may be used to estimate impacts to large numbers of homes based on high-level housing characteristics. In addition to reviewing literature sources, we performed an experiment to evaluate interventions to mitigate degraded indoor air quality. We instrumented one residence for one week during an extreme wildfire event in the Pacific Northwest. Outdoor ambient concentrations of PM2.5 reached historic levels, sustained at over 200 μg/m3 for multiple days. Outdoor and indoor PM2.5 were monitored, and data regarding building characteristics, infiltration, and mechanical system operation were gathered to be consistent with the type of information commonly known for residential energy models. Two conditions were studied: a high-capture minimum efficiency rated value (MERV 13) filter integrated into a central forced air (CFA) system, and a CFA with MERV 13 filtration operating with a portable air cleaner (PAC). With intermittent CFA operation and no PAC, indoor corrected concentrations of PM2.5 reached 280 μg/m3, and indoor/outdoor (I/O) ratios reached a mean of 0.55. The measured I/O ratio was reduced to a mean of 0.22 when both intermittent CFA and the PAC were in operation. Data gathered from the test home were used in a modeling exercise to assess expected I/O ratios from both interventions. The mean modeled I/O ratio for the CFA with an MERV 13 filter was 0.48, and 0.28 when the PAC was added. The model overpredicted the MERV 13 performance and underpredicted the CFA with an MERV 13 filter plus a PAC, though both conditions were predicted within 0.15 standard deviation. The results illustrate the ways that models can be used to estimate indoor PM2.5 concentrations in residences during extreme wildfire smoke events.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Assessment of Technoeconomic Opportunities in Automation for Nuclear Microreactors

Achieving full decarbonization of all economic sectors remains a challenge, especially in niche markets. For example, remote communities and industrial or mining activities detached from the main electric grid heavily rely on fossil fuels, similar to urban and industrial microgrids with combined heat and power needs. A combination of renewables and energy storage is often not suitable due to cost, reliability, intermittency, and large storage requirements. Small nuclear reactors with a flexible purpose could serve these applications. Microreactors (MR) are a class of reactors that are compact, factory manufactured, transportable, and self-regulating. Typically, they generate much less power than their large reactor counterparts. The main advantages of microreactors include the versatile nature of the energy produced, the reliability of supply, and freedom from having to transport and store large quantities of fuels on-site, coupled with the absence of dependence on an electrical grid. A strong business case is needed to move from the microreactor prototype to the commercialization phase. In fact, fossil fuels are still relatively inexpensive, and in the near term, carbon credits will be available to virtually compensate for emissions. For microreactors, one of the main costs in operation and maintenance (O&M) is their staffing levels. In this study, we investigate how to optimize the number (and thus the cost) of workers, moving from a traditional, fully manned, on-site personnel approach to an unmanned, remote personnel approach. We examine four different staffing models that can be implemented as the technology matures and evolves. We estimate the staffing needs of each model and build a business case to justify the substitution of on-site personnel with adequate technologies. To do so, we propose a cost model to quantify potential cost reductions from automating O&M activities. The model accounts for both the reduction in cost derived from the reduced number of full-time-equivalent (FTE) employees and the increase in cost derived from the need to buy new control hardware as needed. Applying the cost model that we created to different scenarios, an on-site O&M cost reduction exceeding 80% can be expected. Additionally, we found that it is more impactful to focus on automating routine O&M tasks rather than attempting to automate transient management (shutdowns, restarts, monitoring condition deviations). In fact, transients typically account for less than 1% of the total FTE time spent on the reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Quantifying resolution in cosmological N- body simulations using self-similarity

ABSTRACT We demonstrate that testing for self-similarity in scale-free simulations provides an excellent tool to quantify the resolution at small scales of cosmological N-body simulations. Analysing two-point correlation functions measured in simulations using abacus , we show how observed deviations from self-similarity reveal the range of time and distance scales in which convergence is obtained. While the well-converged scales show accuracy below 1 per cent, our results show that, with a small force softening length, the spatial resolution is essentially determined by the mass resolution. At later times, the lower cut-off scale on convergence evolves in comoving units as a−1/2 (a being the scale factor), consistent with a hypothesis that it is set by two-body collisionality. A corollary of our results is that N-body simulations, particularly at high red-shift, contain a significant spatial range in which clustering appears converged with respect to the time-stepping and force softening but has not actually converged to the physical continuum result. The method developed can be applied to determine the resolution of any clustering statistic and extended to infer resolution limits for non-scale-free simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

The Spatial Clustering of ROSAT All-Sky Survey Active Galactic Nuclei. V. The Evolution of Broad-line AGN Clustering Properties in the Last 6 Gyr

This is the fifth paper in a series of investigations of the clustering properties of luminous, broad-emission-line active galactic nuclei (AGNs) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this work we measure the cross-correlation function (CCF) between RASS/SDSS Data Release 14 AGNs with the SDSS CMASS galaxy sample at 0.44 < z < 0.64. We apply halo occupation distribution (HOD) modeling to the CCF along with the autocorrelation function of the CMASS galaxies. We find that X-ray-selected and optically selected AGNs at 0.44 < z < 0.64 reside in statistically identical halos with a typical dark matter halo (DMH) mass of $M$$^{typ,AGN}_{DMH}$ ~ 10 12.7 h –1 M ⊙ . The acceptable HOD parameter space for these two broad-line AGN samples have only statistically marginal differences caused by small deviations of the CCFs in the one-halo-dominated regime on small scales. In contrast to optically selected AGNs, the X-ray AGN sample may contain a larger population of satellites at MDMH ~ 1013h–1M⊙. We compare our measurements in this work with our earlier studies at lower independent redshift ranges, spanning a lookback time of 6 Gyr. The comparison over this wider redshift range of 0.07 < z < 0.64 reveals (i) no significant difference between the typical DMH masses of X-ray-selected and optically selected AGNs, (ii) weak positive clustering dependencies of $M$$^{typ,AGN}_{DMH}$ with L X and M BH , (iii) no significant dependence of $M$$^{typ,AGN}_{DMH}$ on Eddington ratio, and (iv) the same DMH masses host more-massive accreting black holes at high redshift than at low redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗

Illustrating galaxy–halo connection in the DESI era with illustrisTNG

ABSTRACT We employ the hydrodynamical simulation illustrisTNG to inform the galaxy–halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ∼ 0.8. Specifically, we model the galaxy colours of illustrisTNG and apply sliding DESI colour–magnitude cuts, matching the DESI target densities. We study the halo occupation distribution (HOD) model of the selected samples by matching them to their corresponding dark matter haloes in the illustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as ‘galaxy assembly bias’. The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modelling galaxy–halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10–15 ${{\ \rm per\ cent}}$ excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures.

79 ASTRONOMY AND ASTROPHYSICS↗

Large-scale simulation-based parametric analysis of an optimal precooling strategy for demand flexibility in a commercial office building

Achieving success with grid-interactive efficient buildings (GEBs) is closely tied to the utilization of flexible loads. A valuable strategy involves the implementation of precooling techniques before high-demand events, such as peak hours, by adjusting zone air temperature setpoints. This leads to a reduction in thermal loads and peak electricity demand during these times, as the building’s thermal mass stores and subsequently releases thermal energy. However, the effectiveness of the pre-cooling optimization is highly contingent on specific conditions such as building thermal properties, weather conditions, utility rate structure, HVAC equipment sizing, etc. Therefore, investigating the impacts of these condition-specific factors is crucial, especially when considering precooling strategies that utilize thermal mass in commercial buildings. In this paper, we first devised a novel heuristic control approach that incorporates parameterized optimal precooling thermostat schedules to enhance demand flexibility in a commercial office building. Subsequently, we conducted a thorough performance evaluation of this control strategy. Here, the optimal thermostat schedule was parameterized using three optimization variables: the precooling start time, the precooling end time, and the precooling temperature setpoint. Utilizing the DOE medium-sized office building as the virtual testbed, we showed that the parameterized schedule effectively approximates model predictive control and requires drastically reduced computational overhead. In addition, we investigated the impact of different influencing factors on the optimal precooling strategy. These factors include building thermal mass, outdoor air conditions, and energy price profiles. Using high-performance computing, we simulated a total of 225 scenarios, consisting of three levels of thermal mass, five typical outdoor air temperature profiles, and fifteen time-of-use price plans. The results demonstrate that optimal thermostat scheduling could save substantial energy cost in medium-sized office buildings with heavy thermal mass but with some energy penalty. Although the potential for cost savings is lower in buildings with low and medium thermal mass, the energy penalty remains consistent in all three thermal mass scenarios. The study also highlights the need to account for zone diversity and recognize that a one-size-fits-all-zone setpoint schedule may not be suitable for all zones and can lead to unnecessary energy wastage. Furthermore, the results highlight that while outdoor air conditions play a role in cost and energy performance, the cooling load exerts a more immediate and substantial influence on cost savings in precooling strategies. Although cost savings are comparable under certain conditions with the same cooling load, observed deviations in energy penalty indicate potential disparities in the efficiency of the HVAC system during the load-shifting process. In addition, the duration of peak pricing and the ratio between peak and off-peak times exhibit clear correlations with cost savings and energy consumption, aligning with intuitive expectations. These findings offer valuable insights for optimizing precooling strategies in office buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

ZTF-SEDm Type Ia supernova sample for Twins Embedding spectrophotometric standardization

Aims. This paper has two aims: the first aim is to build a large homogeneous spectrophotometric sample of Type Ia supernovae (SNe Ia) from the second data release of the Zwicky Transient Facility (ZTF DR2). We used the spectrum sample from the low-resolution ( R ∼ 100) SEDmachine (SEDm) Integral Field Spectrograph (IFS) that gathered 3069 spectra. This is one of the largest samples of such collections that can attempt to reproduce the Twins Embedding (TE) spectrophotometric standardization method. This is our second objective. The method was developed based on high-quality spectra from 200 SNe Ia of the Nearby Supernova factory (SNfactory) and led to an exceptionally low value of 0.073 mag for the intrinsic scatter. Methods. As the SEDm is not designed as a spectrophotometric instrument, we first improved the flux-calibration accuracy of the SN Ia spectrum sample using the ZTF photometric data, which were calibrated at the percent level. We corrected the spectra for second-order polynomials, fitted by comparing the synthetic photometry in the ZTF g , r , i filters with the light-curve (LC) data. We then applied the three steps of the TE parameterization to a subset of 783 ZTF SN spectra near maximum light while comparing results from SNfactory and ZTF. We finally analyzed the standardization methods based on the TE parameters. Results. The precision of the phase-correction model, which is the first step of the TE, is estimated at 0.01 mag in g band based on ZTF data. Despite the challenge posed by the spectrum-extraction pipeline associated with the SEDm (flux calibration, leftover host signal, low signal-to-noise ratio, and low resolution), we applied a first standardization in color based on the second step of the TE, called read between the lines (RBTL), to the ZTF sample. We reached a Hubble residual scatter of 0.153 mag, all in normalized median absolute deviation, which is to be compared to the ∼0.11 mag obtained with the SNfactory data. The SALT color and stretch standardization reaches a scatter of 0.164 mag for the same ZTF SN Ia sample, and its host steps are ∼0.1 mag and zero for RBTL. When considering the scatter due to the redshift error and flux calibration error, we estimated a RBTL scatter of ∼0.129 mag for this ZTF sample as an upper limit because we identified an additional contribution from a systematic error in color. We tested the standardization based on the nonlinear TE parameters, and, as expected from the low spectrum quality, it did not improve the overall dispersion. Conclusions. We release 1897 flux calibrated spectra of 1607 SNe Ia with an estimated photometric accuracy of 0.07 mag. We further demonstrate that some amount of spectrophotometric SN Ia standardization can be done with limited-quality spectra. The RBTL standardization is more efficient than that of SALT with one parameter less, and the resulting host steps are consistent with zero. This makes it less prone to astrophysical bias. For future spectroscopic surveys, targeting the extraction pipeline for a thorough flux calibration and good signal-to-noise ratio would enable us to compute the full TE standardization, which would further reduce the scatter in the distance estimate.

Ganot, C↗

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗

Ion Exchange Processing of AN-107 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AN-107 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 13.7 L of diluted and filtered supernate from Tank 241-AN-107 (hereafter referred to as AN-107) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AN-107 tank waste to meet this criterion, only 0.147% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 678. Testing with AN-107 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 40% Cs breakthrough after processing ~1700 BVs of feed; the 50% Cs breakthrough was extrapolated from the breakthrough data to occur at 1873 BVs. Testing compared to previous AP-101 and AP-107 testing at 16 °C showed ~300 BV increases in volume processed to reach the WAC limit for both lead and lag columns. The increase in capacity was determined to be due to the significantly lower K concentration in the AN-107 compared to the other tank waste matrices. A comparison in breakthrough curves for the three tests indicated slightly slower kinetic behavior in the AN-107, with variations in feed matrices (high organic complexants) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 1097 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 900 BVs. Cs breakthrough from the lag column began at 500 BVs, reaching 3.06×10 0 µCi/mL, or 2.6 % Cs breakthrough, after processing all 1700 BVs of feed. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Considerations for using Privacy Preserving Machine Learning Techniques for Safeguards

In international nuclear safeguards, the International Atomic Energy Agency (IAEA) is tasked with inspecting and verifying nuclear facilities and their activities. Data analytics and machine learning to support inspections require large amounts of data that nuclear facility operators may consider proprietary or sensitive, so the IAEA may not have full access. Allowing computation over private data without compromising its security therefore has value for safeguards inspections and analysis. Privacy-preserving machine learning (PPML) consists of security-focused techniques that allow data analytics and machine learning algorithms to run on sensitive data without revealing it. This includes ideas like homomorphic encryption (HE), secure multiparty computation (SMPC), and secure enclaves. HE allows algorithms and mathematical operations to be conducted directly on the encrypted data instead of first decrypting it. With SMPC, multiple entities collaboratively compute over distributed data such that no party is able to directly view any others’ original data. Secure enclaves allow computation to take place in a separate and heavily blocked-off section of a CPU. Techniques like these allow for several potential use cases in which the security of data is essential. With SMPC, machine learning models can be trained over the input data from multiple entities, resulting in a model that all users can benefit from without leaking the input data from any particular entity. With SMPC or a zero-knowledge proof (ZKP), an algorithm returning some single answer or truth value can be run on someone else’s data without ever needing to see that data, potentially allowing for verification or proof of some underlying question. HE can allow for outsourcing computation on data to a hostile or untrusted environment. Although most of the research in this field resides within the health and financial domains, tools from PPML may have similar applications in nuclear safeguards. Allowing the IAEA to compute over proprietary information, such as process models and raw sensor data using PPML techniques, provides the baseline for running complex analytics without needing direct unencrypted access to the underlying data, maintaining its privacy. Important limitations to consider for these techniques include the efficiency and level of security required. The security of HE and SMPC come at the cost of speed—the significant amount of overhead means that algorithms implemented in these protocols and encryption schemes are slower than when run on plaintext. Additionally, several important parameters determine what techniques or protocols are used based on the security requirements. SMPC protocols may need to be selected for resistance against a party that attempts to deviate from the protocol to distort the result or gain access to additional information, and a protocol secure against these attacks may further increase the overhead of the algorithm.

97 MATHEMATICS AND COMPUTING↗

Planck 2018 results. VI. Cosmological parameters

In this paper, we present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5 σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density Ω c h 2 = 0.120 ± 0.001, baryon density Ω b h 2 = 0.0224 ± 0.0001, scalar spectral index n s = 0.965 ± 0.004, and optical depth τ = 0.054 ± 0.007 (in this abstract we quote 68% confidence regions on measured parameters and 95% on upper limits). The angular acoustic scale is measured to 0.03% precision, with 100 θ * = 1.0411 ± 0.0003. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: Hubble constant H 0 = (67.4 ± 0.5) km s –1 Mpc –1 ; matter density parameter Ω m = 0.315 ± 0.007; and matter fluctuation amplitude σ 8 = 0.811 ± 0.006. We find no compelling evidence for extensions to the base-ΛCDM model. Combining with baryon acoustic oscillation (BAO) measurements (and considering single-parameter extensions) we constrain the effective extra relativistic degrees of freedom to be N eff = 2.99 ± 0.17, in agreement with the Standard Model prediction N eff = 3.046, and find that the neutrino mass is tightly constrained to Σ m ν < 0.12 eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base ΛCDM at over 2 σ , which pulls some parameters that affect the lensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. The joint constraint with BAO measurements on spatial curvature is consistent with a flat universe, Ω K = 0.001 ± 0.002. Also combining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w 0 = –1.03 ± 0.03, consistent with a cosmological constant. We find no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r 0.002 < 0.06. Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations. The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey’s combined-probe results including galaxy clustering (which prefers lower fluctuation amplitudes or matter density parameters), and in significant, 3.6 σ , tension with local measurements of the Hubble constant (which prefer a higher value). Simple model extensions that can partially resolve these tensions are not favoured by the Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion Exchange Processing of AW-105 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AW-105 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 9.2 L of diluted and filtered supernate from Tank 241-AW-105 (hereafter referred to as AW-105) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AW-105 tank waste to meet this criterion, only 0.225% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 445. Testing with AW-105 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 83% Cs breakthrough after processing ~1500 BVs of feed; the 50% Cs breakthrough was interpolated from the breakthrough data and occurred at 1041 BVs. Despite the AW-105 having a significantly higher K concentration (0.55 M compared to 0.10 M), testing compared to previous AP-107 ion exchange column testing at 16 °C showed no difference in BVs processed to reach the WAC on the lead column and only an approximate ~20 BV decrease in volume processed to reach the WAC limit on the lag column. The negligible differences in capacity despite the 5x concentration differences in K was determined to be due to the significantly lower NO3 concentration in the AW-105 supernate compared to the AP-107 tank waste matrix. A comparison in breakthrough curves for the two tests also indicated slightly faster kinetic behavior in the AW-105, with the variations in feed matrices (lower NO3 concentration) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 772 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 675 BVs. Cs breakthrough from the lag column began at 300 BVs, reaching 1.10×10 1 µCi/mL, or 14.13 % Cs breakthrough, after processing all 1500 BVs of feed. The polish column processed nominally 830 BVs and reached 2.10×10 -1 µCi/mL, or 0.27 % Cs breakthrough at the conclusion of the test. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Shallow Cumulus Properties as Captured by Adiabatic Fraction in High-Resolution LES Simulations

Abstract Shallow convective clouds are important players in Earth’s energy budget and hydrological cycle, and are abundant in the tropical and subtropical belts. They greatly contribute to the uncertainty in climate predictions due to their unresolved, complex processes that include coupling between the dynamics and microphysics. Analysis of cloud structure can be simplified by considering cloud motions as a combination of moist adiabatic motions like adiabatic updrafts and turbulent motions leading to deviation from adiabaticity. In this work, we study the sizes and occurrence of adiabatic regions in shallow cumulus clouds during their growth and mature stages, and use the adiabatic fraction (AF) as a continuous metric to describe cloud processes and properties from the core to the edge. To do so, we simulate isolated trade wind cumulus clouds of different sizes using the System of Atmospheric Modeling (SAM) model in high resolution (10 m) with the Hebrew University spectral bin microphysics (SBM). The fine features in the clouds’ dynamics and microphysics, including small near-adiabatic volumes and a thin transition zone at the edge of the cloud (∼20–40 m in width), are captured. The AF is shown to be an efficient measure for analyzing cloud properties and key processes determining the droplet-size distribution formation and shape during the cloud evolution. Physical processes governing the properties of droplet size distributions at different cloud regions (e.g., core, edge) are analyzed in relation to AF. Significance Statement 1) This study investigates the evolution of cumulus clouds (Cu) using a 10-m-resolution LES model with spectral bin microphysics. 2) The study improves the understanding of the mutual effects of adiabatic updrafts and lateral entrainment and mixing. 3) The study demonstrates the existence of an adiabatic core in nonprecipitating Cu. 4) Shapes of the droplet size distributions are closely related to the adiabatic fraction values. 5) Utilization of high resolution reveals the existence of physically significant small features in the cloud structure, such as a narrow cloud interface zone and small adiabatic volumes.

54 ENVIRONMENTAL SCIENCES↗

The SRG/eROSITA All-Sky Survey: Dark Energy Survey year 3 weak gravitational lensing by eRASS1 selected galaxy clusters

Context. Number counts of galaxy clusters across redshift are a powerful cosmological probe if a precise and accurate reconstruction of the underlying mass distribution is performed – a challenge called mass calibration. With the advent of wide and deep photometric surveys, weak gravitational lensing (WL) by clusters has become the method of choice for this measurement. Aims. We measured and validated the WL signature in the shape of galaxies observed in the first three years of the Dark Energy Survey (DES Y3) caused by galaxy clusters and groups selected in the first all-sky survey performed by SRG (Spectrum Roentgen Gamma)/eROSITA (eRASS1). These data were then used to determine the scaling between the X-ray photon count rate of the clusters and their halo mass and redshift. Methods. We empirically determined the degree of cluster member contamination in our background source sample. The individual cluster shear profiles were then analyzed with a Bayesian population model that self-consistently accounts for the lens sample selection and contamination and includes marginalization over a host of instrumental and astrophysical systematics. To quantify the accuracy of the mass extraction of that model, we performed mass measurements on mock cluster catalogs with realistic synthetic shear profiles. This allowed us to establish that hydrodynamical modeling uncertainties at low lens redshifts (z < 0.6) are the dominant systematic limitation. At high lens redshift, the uncertainties of the sources’ photometric redshift calibration dominate. Results. With regard to the X-ray count rate to halo mass relation, we determined its amplitude, its mass trend, the redshift evolution of the mass trend, the deviation from self-similar redshift evolution, and the intrinsic scatter around this relation. Conclusions. The mass calibration analysis performed here sets the stage for a joint analysis with the number counts of eRASS1 clusters to constrain a host of cosmological parameters. We demonstrate that WL mass calibration of galaxy clusters can be performed successfully with source galaxies whose calibration was performed primarily for cosmic shear experiments, opening the way for the cluster cosmological exploitation of future optical and NIR surveys like Euclid and LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Energy Survey year 3 results: covariance modelling and its impact on parameter estimation and quality of fit

ABSTRACT We describe and test the fiducial covariance matrix model for the combined two-point function analysis of the Dark Energy Survey Year 3 (DES-Y3) data set. Using a variety of new ansatzes for covariance modelling and testing, we validate the assumptions and approximations of this model. These include the assumption of Gaussian likelihood, the trispectrum contribution to the covariance, the impact of evaluating the model at a wrong set of parameters, the impact of masking and survey geometry, deviations from Poissonian shot noise, galaxy weighting schemes, and other sub-dominant effects. We find that our covariance model is robust and that its approximations have little impact on goodness of fit and parameter estimation. The largest impact on best-fitting figure-of-merit arises from the so-called fsky approximation for dealing with finite survey area, which on average increases the χ2 between maximum posterior model and measurement by $3.7{{\ \rm per\ cent}}$ (Δχ2 ≈ 18.9). Standard methods to go beyond this approximation fail for DES-Y3, but we derive an approximate scheme to deal with these features. For parameter estimation, our ignorance of the exact parameters at which to evaluate our covariance model causes the dominant effect. We find that it increases the scatter of maximum posterior values for Ωm and σ8 by about $3{{\ \rm per\ cent}}$ and for the dark energy equation-of-state parameter by about $5{{\ \rm per\ cent}}$.

79 ASTRONOMY AND ASTROPHYSICS↗