Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coded computing”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Extracting Topological Orders of Generalized Pauli Stabilizer Codes in Two Dimensions

In this paper, we introduce an algorithm for extracting topological data from translation invariant generalized Pauli stabilizer codes in two-dimensional systems, focusing on the analysis of anyon excitations and string operators. The algorithm applies to Z d qudits, including instances where d is a nonprime number. This capability allows the identification of topological orders that differ from the Z d toric codes. It extends our understanding beyond the established theorem that Pauli stabilizer codes for Z p qudits (with p being a prime) are equivalent to finite copies of Z p toric codes and trivial stabilizers. The algorithm is designed to determine all anyons and their string operators, enabling the computation of their fusion rules, topological spins, and braiding statistics. The method converts the identification of topological orders into computational tasks, including Gaussian elimination, the Hermite normal form, and the Smith normal form of truncated Laurent polynomials. Furthermore, the algorithm provides a systematic approach for studying quantum error-correcting codes. We apply it to various codes, such as self-dual CSS quantum codes modified from the two-dimensional honeycomb color code and non-CSS quantum codes that contain the double semion topological order or the six-semion topological order. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Toward exascale whole-device modeling of fusion devices: Porting the GENE gyrokinetic microturbulence code to GPU

GENE solves the five-dimensional gyrokinetic equations to simulate the development and evolution of plasma microturbulence in magnetic fusion devices. The plasma model used is close to first principles and computationally very expensive to solve in the relevant physical regimes. In order to use the emerging computational capabilities to gain new physics insights, several new numerical and computational developments are required. Here, we focus on the fact that it is crucial to efficiently utilize GPUs (graphics processing units) that provide the vast majority of the computational power on such systems. In this paper, we describe the various porting approaches considered and given the constraints of the GENE code and its development model, justify the decisions made, and describe the path taken in porting GENE to GPUs. We introduce a novel library called gtensor that was developed along the way to support the process. Performance results are presented for the ported code, which in a single node of the Summit supercomputer achieves a speed-up of almost 15× compared to running on central processing unit (CPU) only. Typical GPU kernels are memory-bound, achieving about 90% of peak. Our analysis shows that there is still room for improvement if we can refactor/fuse kernels to achieve higher arithmetic intensity. We also performed a weak parallel scalability study, which shows that the code runs well on a massively parallel system, but communication costs start becoming a significant bottleneck.

Germaschewski, K. (ORCID:0000000284956354)↗
DOE OSTI2024

Cold Plasma Measurements

We have continued the simulation campaign in support of our ongoing magnetospheric cold plasma research project. This project aims to develop the next-generation particle instruments to measure the properties of the cold particle populations in the Earth’s magnetosphere. For this purpose, simulations have been performed with a Particle-In-Cell (PIC) code called the Curvilinear PIC (CPIC). The code is formulated in curvilinear geometry and couples the standard PIC algorithm with algorithms for the generation and adaptation of the underlaying computational mesh. It conforms to complex objects like spacecraft and it can place more grid points in regions where higher resolution is needed. The code also features a scalable solver based on the multigrid algorithm and it is fully parallelized via domain decomposition and MPI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Progress Towards Accelerating the Unified Model on Hybrid Multi-Core Systems

The cloud microphysics scheme, CASIM, and the radiation scheme, SOCRATES, are two computationally intensive parts within the Met Office's Unified Model (UM). This study enables CASIM and SOCRATES to use accelerated multi-core systems for optimal computational performance of the UM. Using profiling to guide our efforts, we refactored the code for optimal threading and kernel arrangement and implemented OpenACC directives manually or through the CLAW source-to-source translator. Initial porting results achieved 10.02x and 9.25x speedup in CASIM and SOCRATES respectively on 1 GPU compared with 1 CPU core. A granular performance analysis of the strategy and bottlenecks are discussed. These improvements will enable UM to run on heterogeneous computers and a path forward for further improvements is provided.

Zhang, Wei↗
DOE OSTI2022

An Algorithmic and Software Pipeline for Very Large Scale Scientific Data Compression with Error Guarantees

Efficient data compression is becoming increasingly critical for storing scientific data because many scientific applications produce vast amounts of data. This paper presents an end-to-end algorithmic and software pipeline for data compression that guarantees both error bounds on primary data (PD) and derived data, known as Quantities of Interest (QoI).We demonstrate the effectiveness of the pipeline by compressing fusion data generated by a large-scale fusion code, XGC, which produces tens of petabytes of data in a single day. We demonstrate that the compression is conducted by setting aside computational resources known as staging nodes, and does not impact the simulation performance. For efficient parallel I/O, the pipeline uses ADIOS2, which many codes such as XGC already use for their parallel I/O. We show that our approach can compress the data by two orders of magnitude while guaranteeing high accuracy on both the PD and the QoIs. Further, the amount of resources required by compression is a few percent of the resources required by simulation while ensuring that the compression time for each stage is less than the corresponding simulation time.This pipeline consists of three main steps. The first step decomposes the data using domain decomposition into small subdomains. Each subdomain is then compressed independently to achieve a high level of parallelism. The second step uses existing techniques that guarantee error bounds on the primary data for each subdomain. The third step uses a post-processing optimization technique based on Lagrange multipliers to reduce the QoI errors for data corresponding to each subdomain. The Lagrange multipliers generated can be further quantized or truncated to increase the compression level. All of the above characteristics of our approach make it highly practical to apply on-the-fly compression while guaranteeing errors on QoIs that are critical to the scientists.

Banerjee, Tania↗
DOE OSTI2020

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Lightshow: a Python package for generating computational x-ray absorption spectroscopy input files

First-principles computational spectroscopy is a critical tool for interpreting experiment, per- forming structure refinement, and developing new physical understanding. Systematically setting up input files for different simulation codes and a diverse class of materials is a challeng- ing task with a very high barrier-to-entry, given the complexities and nuances of each individual simulation package. This task is non-trivial even for experts in the electronic structure field and nearly formidable for non-expert researchers. Lightshow solves this problem by providing a uniform abstraction for writing computational x-ray spectroscopy input files for multiple popular codes, including FEFF, VASP, OCEAN, exciting and XSpectra. Its extendable framework will also allow the community to easily add new functions and to incorporate new simulation codes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Annual IC Progress Report

The goal of this project is to computationally investigate the progenitors of astrophysical gamma-ray bursts (GRB), the most extreme explosions in the universe across length and time scales. We have developed a cutting-edge, code-bridging approach to exploring the physics behind these luminous events to test the hypothesis that GRBs with electromagnetic radio emission originate from massive star (MS) and black hole (BH) binary systems. Our multi-scale, multi-physics theoretical framework for GRBs is based on solving well-defined and testable hydrodynamics problems associated with the dynamical evolution of the system. We solve stages of the evolution in a tractable way with a suite of well-developed state-of-the-art computational tools that include the stellar evolution code MESA, the general relativistic magnetohydrodynamics (GRMHD) code Athena++, and the binary black hole (BBH) population synthesis code COSMIC, as well as a number of self-written python post-processing tools.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Understanding power and energy utilization in large scale production physics simulation codes

Power is an often-cited reason for the move to advanced architectures on the path to Exascale computing. Here, this is due to practical considerations related to delivering enough power to successfully site and operate these machines, as well as concerns about energy usage while running large simulations. Since obtaining accurate power measurements can be challenging, it may be tempting to use the processor thermal design power (TDP) as a surrogate due to its simplicity and availability. However, TDP is not indicative of typical power usage while running simulations. Using commodity and advanced technology systems at Lawrence Livermore and Sandia National Labs, we performed a series of experiments to measure power and energy usage in running simulation codes. These experiments indicate that large scale Lawrence Livermore simulation codes are significantly more efficient than a simple processor TDP model might suggest.

HPC↗
DOE OSTI2020

The Ristra Project: FY20/21 Milestone Report

The ASC Advanced Technology Development and Mitigation (ATDM) sub-program was established in 2014 to develop new simulation tools operating on exascale-class computers to serve NNSA (see Appendix B). Over the course of ATDM, LANL management have set a strategy for exascale-class application codes that follows two supportive and mutually risk-mitigating paths: evolution for established production integrated design codes (IDCs) – with a strong pedigree within the user community – based upon existing programming paradigms(MPI+X); and a new start ATDM project, Ristra, a high-risk/high-reward push for a next-generation multi-physics, multi-scale simulation toolkit based on emerging advanced programming systems(with an initial focus on data-flow task-based models exemplified by Legion). The role of Ristra as the high-risk/high-reward path for LANL’s codes was fully consistent with the goals of ATDM as described in Appendix B, in particular its emphasis on evolving ASC capabilities through novel computing programming models and computing technologies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System: Preprint

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory (NREL) has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that can be neglected for wind turbines. These include buoyant loads, added mass and inertia loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbine systems at a speed comparable to real time. While efficient for long simulations, large sets of load cases, and design studies, mid-fidelity codes cannot capture all of the potentially important physical phenomenon impacting marine turbine systems. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system designed by NREL was modeled both with OpenFAST and with the commercial CFD code, STARCCM+. The CFD model used a 3-D unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modelled with a custom coupling of the CFD solver with the open-source code, MoorDyn. This improves upon the existing quasi-static catenary solver in STARCCM+, which lacks seabed contact or line-to-line connections. Spatial and temporal convergence studies were conducted. The simulation results for a combined current and wave condition are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
DOE OSTI2022

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

An Overview of Gemma FY20 Verification Activities

Gemma verification activities for FY20 can be divided into three categories: the development of specialized quadrature rules, initial progress towards the development of manufactured solutions for code verification, and automated code-verification testing. In the method-of-moments implementation of the electric-field integral equation, the presence of a Green’s function in the four-dimensional integrals yields singularities in the integrand when two elements are nearby. To address these challenges, we have developed quadrature rules to integrate the functions through which the singularities can be characterized. Code verification is necessary to develop confidence in the implementation of the numerical methods in Gemma. Therefore, we have begun investigating the use of manufactured solutions to more thoroughly verify Gemma. Manufactured solutions provide greater flexibility for testing aspects of the code; however, the aforementioned singularities provide challenges, and existing work is limited in rigor and quantity. Finally, we have implemented automated code-verification testing using the VVTest framework to automate the mesh refinement and execution of a Gemma simulation to generate mesh convergence data. This infrastructure computes the observed order of accuracy from these data and compares it with the theoretical order of accuracy to either develop confidence in the implementation of the numerical methods or detect coding errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
DOE OSTI2022

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↗