Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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

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

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Bayesian Optimization of a Free-Electron Laser

The Linac Coherent Light Source X-ray free-electron laser is a complex scienti c apparatus which changes con gurations multiple times per day, necessitating fast tuning strategies to reduce setup time for successive experiments. To this end, we employ a Bayesian approach to maximizing X- ray laser pulse energy by controlling groups of quadrupole magnets. A Gaussian process model provides probabilistic predictions for the machine response with respect to control parameters, enabling a balance of exploration and exploitation in the search for the global optimum. We show that the model parameters can be learned from archived scans, and correlations between devices can be extracted from the beam transport. Furthermore, the result is a sample-e cient optimization routine, combining both historical data and knowledge of accelerator physics to signi cantly outperform existing optimizers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Simulation of Mechanical Fractionation of Chopped Whole-Plant Corn (WPC) Using Discrete Element Method (DEM)

Fractionating whole-plant corn (WPC) in a single-pass harvesting system requires studies on the WPC-to-equipment interaction for improved property control, as well as mechanical and air-driven separation processes compared to the traditional multi-pass grain and stover harvesting system. The discrete element method (DEM) technique has the potential to simulate WPC mechanical fractionation and support simulation-based design of WPC separation processes. In this study, methods to develop DEM particle models of WPC (kernel, cob, stalk, and husk) and their material properties for simulating mass fractionation using the ASABE standard mechanical shaker were proposed. Measurement was done on the axial dimensions (major, intermediate, and minor) and mass of each WPC type (mean sample size is 56), sampled from single-pass harvesting. Applying gaussian multivariate regression and bootstrapping re-sampling techniques, a DEM particle approximate to each WPC was developed. Sensitivity analysis of the DEM Young‘s modulus, Poisson‘s ratio, and interaction parameters of coefficient of restitution, coefficient of rolling friction, and coefficient of static friction on mass fraction was performed after 156 ASABE sieve-shaking DEM simulation runs, generated using Latin Hypercube Design (LHD) design of experiment (DOE) from 19 DEM material parameters. DEM simulation using Hertz-Mindlin with flexible bond contact laws and DOE optimized material properties successfully reproduced the mass fractions retained in ASABE sieves at 9.8% mean relative error and a coefficient of determination of R2 = 0.87. Here, the DEM methodology developed for mechanical WPC mass fractionation could be deployed to perform virtual design of feedstock handling equipment and performance analysis of mechanical fraction systems.

09 BIOMASS FUELS↗
DOE OSTI2021

Waterlevel, geochemical, and borehole data from Rifle, Colorado from 2006-2016

The files included in this data package provide site wide water level data, geochemical data, and borehole information associated with the Rifle site in Colorado during the Lawrence Berkeley National Lab (LBNL) led Integrated Field Research Challenge (IFRC) and Scientific Focus Area (SFA) research projects. The data informs hydrogeochemical interactions occurring in the watershed. The readme file contains detailed metadata for the files, including naming conventions, column headers, and data validation.The borehole information is provided in the file called rifle_boreholes.csv. It provides names and location for boreholes/wells at the Rifle. The waterlevel data is provided in the folder rifle_waterlevels. The geochemical data is provided in two folders: the rifle_floodplain_geochemistry folder has data which was collected from wells across the floodplain and a variety of offsite locations that contribute solutes to the floodplain given their upgradient locations. The rifle_plotC_geochemistry folder has data which was collected primarily in plot C, which was a focused study area in which a large number of groundwater manipulation / biostimulation experiments were conducted. While not all geochemical constituents are included for each well location and sampling date, they include anions, cations, trace metals, dissolved organic and inorganic carbon, stable isotopes of water and sulfur (as sulfate), and water quality parameters including pH, groundwater electrical conductivity (EC) and dissolved oxygen (DO). All files are in csv or txt format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Geo Thermal Cloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The project is motivated by the challenges, risks, and costs associated with geothermal exploration and production. Many processes and parameters impacting geothermal conditions are poorly understood. Diverse datasets are available to help characterize subsurface geothermal conditions (public and proprietary; satellite, airborne surveys, vegetation/water sampling, geological, geophysical, etc.). Yet, it is not clear how to properly leverage these datasets for geothermal exploration due to an incomplete understanding of how physical processes impacting subsurface geothermal conditions are represented in these observations. Recent advancements in machine learning (ML) provide great promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. Our goals and work under Phases 1 and 2 (as proposed) of this project address all these needs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

CIS Project 22359, Final Technical Report. Discretized Posterior Approximation in High Dimensions

Our primary aim in this work is to understand how to efficiently obtain reliable uncertainty quantification in automatic learning algorithms with limited training datasets. Standard approaches rely on cross-validation to tune hyper parameters. Unfortunately, when our datasets are too small, holdout datasets become unreliable—albeit unbiased—measures of prediction quality due to the lack of adequate sample size. We should not place confidence in holdout estimators under conditions wherein the sample variance is both large and unknown. More poigniantly, our training experiments on limited data (Duersch and Catanach, 2021) show that even if we could improve estimator quality under these conditions, the typical training trajectory may never even encounter generalizable models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

R-matrix Resolved Resonance Region Evaluation of 140,142 Ce

Oak Ridge National Laboratory completed the resolved resonance region (RRR) evaluation of the two most abundant cerium isotopes, 140 Ce (88.45%) and 142 Ce (11.11%), as requested by the US Nuclear Criticality Safety Program. These evaluations are based on recent high-resolution transmission and capture measurements performed on nat Ce and highly enriched 142 Ce samples at the JRC-Geel Linear Accelerator facility, as well as measured thermal constants available from the EXFOR database. Starting from the resonance parameters of the ENDF/B-VIII.0 library followed by a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information were derived by fitting these measured data using the Reich–Moore approximation of the R -matrix theory, as implemented in the SAMMY code system. The 140 Ce RRR upper energy limit was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. Updated statistical properties were obtained for the new evaluations and compared to those derived from the ENDF/B-VIII.0 nuclear data library. The new evaluation work improved some of the discrepancies found in previous work, such as the capture resonance integral and stellar Maxwellian-averaged cross sections. These integral quantities were mainly derived from the fit of the latest measured data, especially the neutron capture yield data for 142 Ce isotope.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Probing Non-Standard Interactions in NOvA with Bayesian MCMC Methods

We present a search for non-standard neutrino interactions (NSI) using NOvA's joint $\nu_e$ appearance and $\nu_\mu$ disappearance samples in both neutrino and antineutrino beam modes. A Bayesian Markov Chain Monte Carlo approach is used to map the posterior distribution over NSI parameters $\varepsilon_{e\mu}$, $\varepsilon_{e\tau}$, and $\varepsilon_{\mu\tau}$, simultaneously with standard oscillation parameters. We investigate the effect of prior choice for the complex NSI parameters and present projected sensitivity to off-diagonal NSI.

Huang, Xiaoyan [Mississippi U.]↗
DOE OSTI2022

Mercury Accumulation in a Stream Ecosystem: Linking Labile Mercury in Sediment Porewaters to Bioaccumulative Mercury in Trophic Webs

Mercury (Hg) deposition and accumulation in the abiotic and biotic environments of a stream ecosystem were studied. This study aimed to link labile Hg in porewater to bioaccumulative Hg in biota. Sediment cores, porewaters, and biota were sampled from four sites along the Fourmile Branch (SC, USA) and measured for total Hg (THg) and methyl-Hg (MHg) concentrations. Water quality parameters were also measured at the sediment–water interface (SWI) to model the Hg speciation. In general, Hg concentrations in porewaters and bulk sediment were relatively high, and most of the sediment Hg was in the solid phase as non-labile species. Surface sediment presented higher Hg concentrations than the medium and bottom layers. Mercury methylation and MHg production in the sediment was primarily influenced by sulfate levels, since positive correlations were observed between sulfate and Hg in the porewaters. The majority of Hg species at the SWI were in non-labile form, and the dominant labile Hg species was complexed with dissolved organic carbon. MHg concentrations in the aquatic food web biomagnified with trophic levels (biofilm, invertebrates, and fish), increasing by 3.31 times per trophic level. Based on the derived data, a modified MHg magnification model was established to estimate the Hg bioaccumulation at any trophic level using Hg concentrations in the abiotic environment (i.e., porewater).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

ISMIP6-based Antarctic projections to 2100: simulations with the BISICLES ice sheet model

The contribution of the Antarctic Ice Sheet is one of the most uncertain components of sea level rise to 2100. Ice sheet models are the primary tool for projecting future sea level contribution from continental ice sheets. The Ice Sheet Model Intercomparison for the Coupled Model Intercomparison Phase 6 (ISMIP6) provided projections of the ice sheet contribution to sea level over the 21st century, quantifying uncertainty due to ice sheet model, climate model, emission scenario, and uncertain parameters. We present simulations following the ISMIP6 framework with the BISICLES ice sheet model and new experiments extending the ISMIP6 protocol to more comprehensively sample uncertainties in future climate, ice shelf sensitivity to ocean melting, and their interactions. These results contributed to the land ice projections of Edwards et al. (2021), which formed the basis of sea level projections for the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (AR6). Our experiments show the important interplay between surface mass balance processes and ocean-driven melt in determining Antarctic sea level contribution. Under higher-warming scenarios, high accumulation offsets more ocean-driven mass loss when sensitivity to ocean-driven melt is low. Conversely, we show that when sensitivity to ocean warming is high, ocean melting drives increased mass loss despite high accumulation. Overall, we simulate a sea level contribution range across our experiments from 2 to 178 mm. Finally, we show that collapse of ice shelves due to surface warming increases sea level contribution by 25 mm relative to the no-collapse experiments, for both moderate and high sensitivity of ice shelf melting to ocean forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tuned oscillator atomic force microscopy methods and apparatus

Techniques for operating an atomic force microscope, the atomic force microscope comprising a cantilever and configured to image a surface of a sample using a probe tip coupled to the cantilever, the techniques comprising using a controller to perform: obtaining, based on at least one intrinsic parameter of the cantilever, a first quality factor and a first free oscillation amplitude, wherein the cantilever exhibits only one stable oscillation state when oscillating at the first free oscillation amplitude and operating at the first quality factor; and controlling the cantilever to exhibit the only one stable oscillation state by controlling the cantilever to oscillate at a fixed frequency at or near a resonance frequency of the cantilever, oscillate at the first free oscillation amplitude, and operate at the first quality factor.

Schwarz, Udo↗
DOE OSTI2026

Cross-correlation of SPT-3G D1 CMB lensing and DES Y3 galaxy lensing

Measurements of the weak lensing of galaxies and of the cosmic microwave background (CMB) provide direct probes of the cosmic matter density field, but the two observables are sensitive to different spatial scales, redshift ranges, and survey systematics. Their cross-correlation thus enables consistency checks of the theoretical model and of potential systematics in either dataset. We present measurements of the cross-correlation between CMB lensing and cosmic shear over $\sim$1,300 deg$^2$ of the sky using the SPT-3G D1 CMB lensing maps and the Dark Energy Survey Year 3 (DES Y3) shear catalogs. For the first time, we measure this cross-correlation at high significance ($\sim 14σ$) when using a polarization-only CMB lensing reconstruction that is expected to be robust against biases induced by extragalactic foregrounds. We test a variety of other CMB lensing estimators that include temperature information and exhibit different tradeoffs between foreground biases and noise, as well as a shear sample that consists of blue, star-forming galaxies and has been shown to be less impacted by galaxy intrinsic alignments. Assuming $Λ$CDM and marginalizing over uncertainties in intrinsic alignments, baryonic feedback, and various nuisance parameters, we obtain a constraint on the amplitude of matter clustering $S_8 \equiv σ_8 \sqrt{Ω_m / 0.3} = 0.833^{+0.047}_{-0.061}$, consistent with both the primary CMB results from Planck and shear-only results from DES Y3. By combining our measurement with Planck, we find mild constraints on the astrophysical processes that impact the cross-correlation. We obtain a constraint on the intrinsic alignment amplitude of the DES sample that is competitive with that from shear-only analyses, and we find a lower limit on the strength of baryonic feedback.

Ouellette, A. [Illinois U., Urbana] (ORCID:0000000↗
DOE OSTI2021

Scheduling strategies for the ESPRESSO follow-up of TESS targets

ABSTRACT Radial-velocity follow-up of stars harbouring transiting planets detected by TESS is expected to require very large amounts of expensive telescope time in the next few years. Therefore, scheduling strategies should be implemented to maximize the amount of information gathered about the target planetary systems. We consider myopic and non-myopic versions of a novel uniform-in-phase scheduler, as well as a random scheduler, and compare these scheduling strategies with respect to the bias, accuracy and precision achieved in recovering the mass and orbital parameters of transiting and non-transiting planets. This comparison is carried out based on realistic simulations of radial-velocity follow-up with ESPRESSO of a sample of 50 TESS target stars, with simulated planetary systems containing at least one transiting planet with a radius below 4R⊕. Radial-velocity data sets were generated under reasonable assumptions about their noise component, including that resulting from stellar activity, and analysed using a fully Bayesian methodology. We find the random scheduler leads to a more biased, less accurate, and less precise, estimation of the mass of the transiting exoplanets. No significant differences are found between the results of the myopic and non-myopic implementations of the uniform-in-phase scheduler. With only about 22 radial velocity measurements per data set, our novel uniform-in-phase scheduler enables an unbiased (at the level of 1 per cent) measurement of the masses of the transiting planets, while keeping the average relative accuracy and precision around 16 per cent and 23 per cent, respectively. The number of non-transiting planets detected is similar for all the scheduling strategies considered, as well as the bias, accuracy and precision with which their masses and orbital parameters are recovered.

Cabona, L.↗
DOE OSTI2022

Low Cost, High Volume, Carbon Fiber Precursor for Plasma Oxidation

Light weighting with carbon fiber is critical to improving the energy efficiency of wind turbines, airplanes, and automobiles, but the high cost of carbon fiber continues to limit market growth and applications. Half of this cost is in the raw materials. For the carbon fiber market to reach the next level of adoption, the cost of both the raw materials and conversion must be reduced. This project considered the combination of low-cost, textile-grade precursors combined with the increased efficiency of plasma oxidation to greatly reduce the cost of carbon fiber. The primary objective of this project was to make industrial-grade carbon fiber (550 ksi tensile strength, 35 Msi tensile modulus) from textile-grade PAN precursors, something that cannot be achieved with conventional oxidation and to date has not been achieved with any technology. To accomplish this, the project team iteratively optimized the process conditions in the plasma oxidation and carbonization stages to maximize the tensile properties of the resultant carbon fibers. This involved implementing experimental designs for both stages, whereby process conditions (e.g. temperature) were varied and correlated to fiber properties. Samples of oxidized PAN fibers (OPFs) were prepared using a small pilot-scale plasma oxidation oven at 4XTechnologies (4XT). The thermal, chemical, and physical properties of OPF samples were analyzed by the University of Tennessee (UTK) and Oak Ridge National Laboratory (ORNL). The OPF samples were then converted to carbon fiber by subsequent low- and high-temperature carbonization using a lab-scale furnace at ORNL. The resultant carbon fibers were tested for tensile strength also at ORNL, and the results were fed back into the experimental design. Of the hundreds of carbon fiber samples produced over the course of the project, the best sample had an average break strength of 479 ksi and modulus of 33 Msi. While the best tensile properties achieved fell short of the final project targets of 550 ksi for break strength and 35 Msi for modulus, the team demonstrated the potential to produce industrial-grade carbon fiber using textile-grade precursors. The best performing fiber samples had tensile properties that are 87% of the goal for tensile strength and 94% of the goal for tensile modulus. With additional optimization work, it is likely that the targets could be achieved, since only a fraction of the process parameter space was explored during this project. 4XT and its affiliate (4M) have already secured funding through private equity, and there is a CRADA (Cooperative Research and Development Agreement) project currently underway with the Carbon Fiber Technology Center (CFTF) at ORNL to continue this research. Additionally, while not an original goal of the project, the carbon fibers for most samples had diameters >8.5μm. This is considerably larger than commercially available carbon fiber products, which can impart unique properties to the composites (e.g. increased compressive strength) and decrease manufacturing costs even further. Moreover, these fibers were oxidized with plasma oxidation using residence times ranging <60mins, which would be considered fast for normal diameter fibers, but it is as much as 5x faster when factoring in the diameter of the fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultrafast formation of domain walls of a charge density wave in SmTe 3

We study ultrafast x-ray diffraction on the charge density wave (CDW) of SmTe 3 using an x-ray free-electron laser. The high momentum and time resolution afforded by the x-ray laser enabled capturing fine wave-vector and time-dependent features of the CDW that originate from fast (in time) and sharp (in real space) variations of the CDW lattice distortion, which we attribute to an inversion of the order parameter. These domain inversions occur near the surface and are caused by the short penetration depth of the near-infrared pump with the wavelength centered at 800 nm, resulting in CDW domain walls perpendicular to the sample surface. These domain walls break the CDW long-range order on the scale of the x-ray probe depth, controlled experimentally by the x-ray incidence angle and suppress the diffraction intensity of the CDW for times much longer than the ~1–ps recovery of the electronic gap observed in time and angle-resolved photoemission spectroscopy. We model the spatial and temporal dependences of the order parameter using a simple Ginzburg-Landau model with all the parameters obtained from the published literature. We find reasonable agreement between the calculated and the measured diffraction across the momentum, time, fluence, and incidence angle dependence without adjusting any parameters. In conclusion, we reconstruct the spatial and temporal dependences of the lattice order parameter and find that at long times, depending on the pump fluence, multiple domain walls remain at distances of a few nanometers from the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗
DOE OSTI2024

Learning of networked spreading models from noisy and incomplete data

Recent years have seen a lot of progress in algorithms for learning parameters of spreading dynamics from both full and partial data. Some of the remaining challenges include model selection under the scenarios of unknown network structure, noisy data, missing observations in time, as well as an efficient incorporation of prior information to minimize the number of samples required for an accurate learning. Here, in this work, we introduce a universal learning method based on a scalable dynamic message-passing technique that addresses these challenges often encountered in real data. The algorithm leverages available prior knowledge on the model and on the data, and reconstructs both network structure and parameters of a spreading model. We show that a linear computational complexity of the method with the key model parameters makes the algorithm scalable to large network instances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Designing an Optimal LSST Deep Drilling Program for Cosmology with Type Ia Supernovae

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) is forecast to collect a large sample of Type Ia supernovae (SNe Ia) expected to be instrumental in unveiling the nature of dark energy. The feat, however, requires accurately measuring the two components of the Hubble diagram, distance modulus and redshift. Distance is estimated from SN Ia parameters extracted from light-curve fits, where the average quality of light curves is primarily driven by survey parameters. An optimal observing strategy is thus critical for measuring cosmological parameters with high accuracy. We present in this paper a three-stage analysis to assess the impact of the deep drilling (DD) strategy parameters on three critical aspects of the survey: redshift completeness, the number of well-measured SNe Ia, and cosmological measurements. We demonstrate that the current DD survey plans (internal LSST simulations) are characterized by a low completeness (z ~ 0.55–0.65), and irregular and low cadences (several days), which dramatically decrease the size of the well-measured SN Ia sample. We propose a method providing the number of visits required to reach higher redshifts. We use the results to design a set of optimized DD surveys for SN Ia cosmology taking full advantage of spectroscopic resources for host galaxy redshift measurements. The most accurate cosmological measurements are achieved with deep rolling surveys characterized by a high cadence (1 day), a rolling strategy (at least two seasons of observation per field), and ultradeep (z ≳ 0.8) and deep (z ≳ 0.6) fields. A deterministic scheduler including a gap recovery mechanism is critical to achieving a high-quality DD survey.

79 ASTRONOMY AND ASTROPHYSICS↗