Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic timing analysis”

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 505 records · Page 28

Modified Eigen-Decomposition-based Interval Analysis (MEDIA) for Power System Dynamic State Estimation

The Bayesian approach has been used for the dynamic state estimation (DSE) of a power system. However, due to the complexity of noise resources, it is difficult to quantify measurement and process noise using probability density functions (PDFs). To overcome the difficulty, the authors of this paper propose a modified eigen-decomposition-based interval analysis (MEDIA) method, which employs bounds instead of PDFs to quantify the noise, and uses the eigen decomposition method to reduce the negative impact of the overestimation problem. Using the simulation data generated from IEEE 16-machine and IEEE 10-machine systems, it is shown that the proposed MEDIA method can estimate the hard boundaries of dynamic states in real time. Furthermore, comparison with the forward-backward propagation method and the extended set-membership filter also shows that the proposed MEDIA method performs better by providing narrower boundaries in the DSE.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thermodynamics of optical Bloch equations

Abstract Optical Bloch equations (OBE) describe the coherent exchange of energy between a quantum bit (qubit) and a quasi-resonant driving field in the presence of a thermal bath. Despite it being an ubiquitous process in quantum technologies, a sound thermodynamic analysis is still missing. We hereby provide such an analysis, by deriving the relevant framework from first principles. We start from a complete microscopic description of the qubit-bath system where definitions of heat, work and entropy production are unambiguous. We trace out the bath and coarse-grain the resulting expressions in time, using a methodology similar to the derivation of the dynamical master equation, to derive closed expressions for the first and second law in terms of system properties. Long coarse graining times yield the Floquet Master equation and its already known thermodynamic description. Short coarse-graining times yield instead the OBE and a novel thermodynamic framework which explicitly depends on quantum coherences in the qubit’s energy basis which produce quantum signatures in the heat and entropy production flows. This allows us to characterize a genuinely quantum non-equilibrium situation, where the coherences created by the driving field are continuously erased by the bath. Our findings can be readily extended to larger open quantum systems. They carry the seeds for future thermodynamic analyses of quantum gates and the design of quantum engines in the strong coherent driving regime.

Elouard, Cyril↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗

Evaluation of AI-enabled Digital Documented Safety Analysis

The National Reactor Innovation Center (NRIC) is leading a transformative initiative to accelerate advanced reactor deployment by fundamentally reimagining how nuclear safety basis documentation is developed, reviewed, and maintained. Traditional Documented Safety Analysis (DSA) processes for DOE-authorized facilities rely on static, document-centric workflows that consume significant time and resources, exemplified by recent major licensing efforts requiring hundreds of thousands of staff hours and millions of pages of documentation review. These conventional approaches create barriers to the rapid, cost-effective deployment of advanced reactors that America's future energy needs demand. NRIC's DOE Authorization Digital Transformation Project addresses these challenges through an innovative framework that integrates artificial intelligence (AI), digital engineering, and systems-based data management into a cohesive digital ecosystem. This white paper presents NRIC's methodology for evaluating AI-enabled document generation capabilities within this broader digital infrastructure, using the Demonstration of Microreactor Experiments (DOME) facility as a pilot case study. The evaluation will assess an AI tool's ability to generate a Preliminary Documented Safety Analysis (PDSA) through progressive integration stages—from standalone document processing to full digital thread connectivity—while maintaining rigorous verification, validation, and regulatory acceptance standards. By establishing dynamic, traceable connections between design data and safety documentation, NRIC's approach has the potential to reduce both document development time and regulatory review cycles by as much as 50%, while simultaneously improving accuracy, consistency, and traceability. This initiative represents a critical step toward establishing reusable digital infrastructure that reactor developers can leverage to accelerate their path from concept to commercial operation, directly supporting NRIC's mission to demonstrate and deploy advanced nuclear energy technologies.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Visual HPC Workflows for the Analysis of System Dynamics Models

Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

Demonstrating the viability of Lagrangian in situ reduction on supercomputers

Performing exploratory analysis and visualization of large-scale time-varying computational science applications is challenging due to inaccuracies that arise from under-resolved data. In recent years, Lagrangian representations of the vector field computed using in situ processing are being increasingly researched and have emerged as a potential solution to enable exploration. However, prior works have offered limited estimates of the encumbrance on the simulation code as they consider “theoretical” in situ environments. Further, the effectiveness of this approach varies based on the nature of the vector field, benefitting from an in-depth investigation for each application area. With this study, an extended version of Sane et al. (2021), we contribute an evaluation of Lagrangian analysis viability and efficacy for simulation codes executing at scale on a supercomputer. We investigated previously unexplored cosmology and seismology applications as well as conducted a performance benchmarking study by using a hydrodynamics mini-application targeting exascale computing. Here, to inform encumbrance, we integrated in situ infrastructure with simulation codes, and evaluated Lagrangian in situ reduction in representative homogeneous and heterogeneous HPC environments. To inform post hoc accuracy, we conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. Additionally, our study contributes cost estimates for distributed-memory post hoc reconstruction. In all, we demonstrate viability for each application — data reduction to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 90% of our experiments.

97 MATHEMATICS AND COMPUTING↗

Fe protein docking transduces conformational changes to MoFe nitrogenase active site in a nucleotide-dependent manner

The reduction of dinitrogen to ammonia catalyzed by nitrogenase involves a complex series of events, including ATP hydrolysis, electron transfer, and activation of metal clusters for N 2 reduction. Early evidence shows that an essential part of the mechanism involves transducing information between the nitrogenase component proteins through conformational dynamics. Here, millisecond time-resolved hydrogen-deuterium exchange mass spectrometry was used to unravel peptide-level protein motion on the time scale of catalysis of Mo-dependent nitrogenase from Azotobacter vinelandii. Normal mode analysis calculations complemented this data, providing insights into the specific signal transduction pathways that relay information across protein interfaces at distances spanning 100 Å. Together, these results show that conformational changes induced by protein docking are rapidly transduced to the active site, suggesting a specific mechanism for activating the metal cofactor in the enzyme active site.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparing temperature convergence of shocked thin films of tin and iron to a bulk temperature source

An outstanding challenge in developing a complete equation of state for materials at elevated pressure and temperature is a robust method of determining the bulk temperature state under dynamic conditions. In metals, the determination of bulk temperature states by optical pyrometry is complicated by the small optical depth and thermal conduction effects. These effects lead to observed temperatures differing by 20% or more from the bulk temperature state. In this work, we show the presence of thermal conduction effects in temperature measurements of tin and iron coatings during dynamic compression experiments. We demonstrate that tin, in contrast to iron, coatings can fail to converge to a bulk temperature source over the time scale of the experiment, requiring the experimenter to modify assumptions, design, or analysis. This work bounds thermal transport at shocked conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Observation of D + → η ′ μ + ν μ and First Experimental Study of D + → η ′ ℓ + ν ℓ Decay Dynamics

Using 20.3 fb − 1 of e + e − collision data collected at the center-of-mass energy 3.773 GeV with the BESIII detector, we report the first observation of the semileptonic decay D + → η ′ μ + ν μ with significance of 8.6 σ including systematic uncertainties, and an improved measurement of D + → η ′ e + ν e . The branching fractions of D + → η ′ μ + ν μ and D + → η ′ e + ν e are determined to be ( 1.92 ± 0.2 8 stat ± 0.0 8 syst ) × 10 − 4 and ( 1.79 ± 0.1 9 stat ± 0.0 7 syst ) × 10 − 4 , respectively. The ratio of the two branching fractions is determined to be R μ / e = 1.07 ± 0.1 9 stat ± 0.0 3 syst , which agrees with the theoretical expectation of lepton flavor universality within the standard model. From an analysis of the D + → η ′ ℓ + ν ℓ decay dynamics, the product of the hadronic form factor f + η ′ ( 0 ) and the Cabibbo-Kobayashi-Maskawa matrix element | V c d | is measured for the first time, giving f + η ′ ( 0 ) | V c d | = ( 5.92 ± 0.5 6 stat ± 0.1 3 syst ) × 10 − 2 . The η − η ′ mixing angle in the quark flavor basis is determined to be ϕ P = ( 39.8 ± 0.8 stat ± 0.3 syst ) ° . Published by the American Physical Society 2025

Ablikim, M.↗

Progressing Analysis of Variable Electric Rates (PAVER) Study

The Progressing Analysis of Variable Electric Rates (PAVER) study analyzed the impact of a range of time-varying electric rates on the performance of a regional electric grid and the resulting costs for participating and non-participating customers. This analysis leveraged and extended the work of PNNL’s Distribution System Operator with Transactive (DSO+T) study. Five different rate designs were included: a flat volumetric energy charge, a typical Time of Use (TOU) rate, a dynamic energy (DE) rate (based on wholesale locational marginal prices), a dynamic energy and capacity (DE+C) rate, and, finally, a Block and Swing (B&S) rate that billed customers based on their average load profile at constant pricing, but used the DE+C dynamic price for load deviations from their average profile. These rates were analyzed in a large-scale co-simulation of an entire regional grid with a customer population representative of the current state. A large fraction (80%) of residential and commercial customers were assumed to participate in these time-varying rates with automatically controlled HVAC, water heaters, electric vehicles, and batteries. This study assumed no industrial sector participation. The DE and DE+C rates saw system peak loads reduced by 6-7%, while the large participation in the TOU rate case saw a significant rebound effect and a resulting peak load increase of >5%. The impacts to the annual and peak system demand impacted system wholesale prices and the overall grid operating costs. This cost structure determined the revenue needed to be collected from customers by each rate design. Participating customers on the DE and DE+C rates (located in one of the modeled DSOs) saw reductions in average annual electricity bills of 11-17% with average increases in monthly bill variation of no more than 13%. At such high participation levels, TOU customers saw 10% higher average annual bills (due to system-wide rebound effects) and average increased monthly bill variation of 16%. Residential owners of large flexible loads (such as electric vehicles) saw larger bill savings (17-20%) when on a fully dynamic rate. The presence of on-site generation (such as rooftop solar) did not appear to appreciably change customer outcomes. Customers on the Block and Swing rate did see 6% lower monthly bill variation (as intended) than the flat rate case, but at the expense of appreciable bill savings, which were only 3%, comparable to the savings seen by non-participants. Given this finding we recommend that additional research be conducted into how best various bill protection mechanisms can balance minimizing customer bill variation with providing financial incentives commensurate with the flexibility customers provide. We also recommend that customer outcomes be explored across a range of regions using current actual customer and system cost data.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

From Ensemble Climate to Ensemble Impacts

Many climate-risk tools rely on ensemble mean projections or endpoint climate snapshots to characterize future hazards. Although convenient for communication, these representations remove the statistical, temporal, and physical information that real infrastructure systems respond to. Infrastructure degradation and failure arise from extremes, sequences, cumulative stress, compound hazards, and nonlinear fragility relationships, none of which survive ensemble averaging or temporal compression. Power-system failure statistics and cascading failure models further show that infrastructure risk is dominated by tail events and path-dependent dynamics rather than by mean conditions. This paper demonstrates why ensemble mean or endpoint-only climate representations are mathematically and physically inconsistent with engineering-grade risk analysis. We outline a model-resolved, time-series-based workflow that preserves extremes, variability, and sequencing by propagating each climate-model realization independently through hazard formation, exposure, fragility, and cascading failure mechanisms. Taking the ensemble of impacts—rather than the ensemble of climate—provides a defensible, physically coherent foundation for infrastructure resilience planning, regulatory compliance, and long-term investment decisions.

54 - ENVIRONMENTAL SCIENCES/GLOBAL CLIMATE CHANGE ↗

Direct Study of Changes in Catalyst Structure-Kinetic Properties During Redox Transitions

The Temporal Analysis of Products (TAP) pulse response methodology is a transient technique that provides the time resolution needed to deconvolve reaction steps from the complex networks typical in industrial catalytic processes. Traditionally, TAP measurements observe gas phase dynamics at the reactor exit but lack direct measurements of changes in the catalyst itself. Recently, a new operando technique was developed that couples gas phase transients to dynamic changes in metal centers with the precise TAP methodology for nanomole titration. Using an industrial CrOx/Al2O3 catalyst used for propane dehydrogenation, we demonstrate the capabilities of this unique device to reveal key catalytic processes: 1) total propane oxidation not associated with chromia centers, 2) reduction of Cr6+ to Cr3+ correlated with selective product formation, and 3) subsequent carbon accumulation. By utilizing incremental pulsing in a diffusion-only transport regime, the spectrokinetic device allows us to resolve detailed changes in catalyst structure, composition, and kinetic function that are otherwise indistinguishable in conventional operando devices. The unification of the TAP methodology with time-resolved spectroscopic measurements offers new and unique insights into the complex kinetic phenomena regulated by solid catalyst surfaces.

03 - NATURAL GAS↗

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗

Statistical analysis of correlations in the x-ray induced Coulomb explosion of iodopyridine

Coulomb explosion imaging is a promising experimental tool to study individual molecules. Here, in this work simulation data for the x-ray induced Coulomb explosion of 2-iodopyridine is analyzed and the involved fragmentation dynamics are described. It is found that particular final ion momenta show correlations that reflect a collision of two atoms during the explosion. Variations of the proton kinetic energies can be associated with variations in the charge build-up time. Covariances of the forces during the explosion can be utilized to simplify the description of the dynamics in reduced dimensionality using only four collective coordinates.

74 ATOMIC AND MOLECULAR PHYSICS↗

Approximate entropy analysis for nonlinear beam dynamics

In this paper, we applied an approximate entropy (ApEn) analysis to nonlinear beam dynamics. Due to the presence of strong nonlinear magnets, the degree of chaos in beam motion gradually increases proportionally with its amplitude. Such chaos can be quantitatively characterized with ApEn of beam turn-by-turn readings. ApEn analysis is a technique used to quantify the amount of regularity and the unpredictability of fluctuations with respect to time-series data. The ApEn, when applied as a chaos indicator, can then be used for nonlinear lattice optimization and analysis. The National Synchrotron Light Source II (NSLS-II) electron storage ring lattice was used as an example of a real-world application of our technique. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrating the PanDA Workload Management System with the Vera C. Rubin Observatory

The Vera C. Rubin Observatory will produce an unprecedented astronomical data set for studies of the deep and dynamic universe. Its Legacy Survey of Space and Time (LSST) will image the entire southern sky every three to four days and produce tens of petabytes of raw image data and associated calibration data over the course of the experiment’s run. More than 20 terabytes of data must be stored every night, and annual campaigns to reprocess the entire dataset since the beginning of the survey will be conducted over ten years. The Production and Distributed Analysis (PanDA) system was evaluated by the Rubin Observatory Data Management team and selected to serve the Observatory’s needs due to its demonstrated scalability and flexibility over the years, for its Directed Acyclic Graph (DAG) support, its support for multi-site processing, and its highly scalable complex workflows via the intelligent Data Delivery Service (iDDS). PanDA is also being evaluated for prompt processing where data must be processed within 60 seconds after image capture. This paper will briefly describe the Rubin Data Management system and its Data Facilities (DFs). Finally, it will describe in depth the work performed in order to integrate the PanDA system with the Rubin Observatory to be able to run the Rubin Science Pipelines using PanDA.

79 ASTRONOMY AND ASTROPHYSICS↗