Search NASA⌕ Search

SEARCH · Search NASA

Results for “batch simulation”

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 145 records · Page 8

Simulation of Bunch Formation for the Mu2e Experiment

The Fermilab Recycler is an 8 GeV storage ring composed of permanent magnets that was crucial to the success of the Fermilab Tevatron Collider program. It is currently being used to slip-stack protons for the high energy neutrino program and to re-bunch protons for use in the Muon g-2 and Mu2e experiments. For the latter applications, the Recycler re-bunches each 1.6 µs "batch" from the Fermilab Booster into four 2.5 MHz bunches. For the Mu2e experiment, it is crucial that beam more than 125 ns from the nominal bunch center be suppressed by at least a factor of 1E-5. While bunch formation is currently in operation for the g-2 experiment, this out of time requirement has not been met, and the reason is not understood. This work presents a simulation of bunch formation in the Recycler, in an effort to understand the reason for this excessive out of time beam and to search for a way to reduce it.

43 PARTICLE ACCELERATORS↗

An end-to-end pipeline for succinic acid production at an industrially relevant scale using Issatchenkia orientalis

Abstract Microbial production of succinic acid (SA) at an industrially relevant scale has been hindered by high downstream processing costs arising from neutral pH fermentation for over three decades. Here, we metabolically engineer the acid-tolerant yeast Issatchenkia orientalis for SA production, attaining the highest titers in sugar-based media at low pH (pH 3) in fed-batch fermentations, i.e. 109.5 g/L in minimal medium and 104.6 g/L in sugarcane juice medium. We further perform batch fermentation using sugarcane juice medium in a pilot-scale fermenter (300×) and achieve 63.1 g/L of SA, which can be directly crystallized with a yield of 64.0%. Finally, we simulate an end-to-end low-pH SA production pipeline, and techno-economic analysis and life cycle assessment indicate our process is financially viable and can reduce greenhouse gas emissions by 34–90% relative to fossil-based production processes. We expect I. orientalis can serve as a general industrial platform for production of organic acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

User guide for the Geospatial Risk Analysis Tool (V.1.0.1)

Developed by the Pacific Northwest National Laboratory (PNNL), the Geospatial Risk Analysis Tool (GRAT) links to Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT), a complete system for computing simple air parcel trajectories, as well as complex transport, dispersion, chemical transformation, and deposition simulations. The HYSPLIT model was developed by the National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory and the Australian Bureau of Meteorology Research Center in 1998. GRAT is used to batch process meteorological data into ARL format. While HYSPLIT GUI limits users to convert meteorological data up to 6 date time points per operation, GRAT enables users to convert a significantly large amount of data (For example, 10 years of meteorological data measured with an interval of 5 minutes) within one operation. Incorporated with the Potential source distribution function (PSDF), GRAT can also be used to determine the areas influenced by the emission of hazardous chemicals.

54 ENVIRONMENTAL SCIENCES↗

Formulation and solution approach for calibrating activity-based travel demand model-system via microsimulation

This study addresses the problem of calibrating utility-maximizing nested logit activity-based travel demand model-systems. After estimation, it is common practice to use aggregate measurements to calibrate the estimated model-system’s parameters prior to their application in transportation planning, policy making, and operations. However, calibration of activity-based model-systems has received much less attention. Existing calibration approaches are myopic heuristics in the sense that they do not consider the fundamental inter-dependencies among choice-models and do not have a systematic way to adjust model parameters. Also, other purely simulation-based approaches do not perform well in large-scale applications. In this study, we focus on utility-maximizing nested logit activity-based model-systems and calibrating aggregate statistics such as activity shares, mode shares, time-dependent & mode-specific OD flows, and time-dependent & mode-specific sensor counts. We formulate the calibration problem as a simulation-based optimization problem and propose a stochastic gradient-based solution procedure to solve it. The solution procedure relies on microsimulation to calculate expectations of the aggregate statistics of interest to the calibration problem. Additionally, we derive approximate analytical expressions for the gradient of the objective function —that are evaluated through microsimulation on mini-batches of the population. The proposed solution procedure is sensitive to the fundamental structure of the activity-based model-system and is non-myopic in considering the dependencies across its model components. The formulated optimization problem is non-convex, highly nonlinear, and potentially has multiple-minima. Lastly, we show —through a real-world application— that the proposed solution procedure outperforms other state-of-the-art purely simulation-based optimization approaches in terms of computational efficiency, stability, and convergence. We also compare various gradient-based solution algorithms to determine the best algorithm to update the parameters. This work has the potential to facilitate wider and easier application of activity-based model-systems.

97 MATHEMATICS AND COMPUTING↗

Accelerating science: The usage of commercial clouds in ATLAS Distributed Computing

The ATLAS experiment at CERN is one of the largest scientific machines built to date and will have ever growing computing needs as the Large Hadron Collider collects an increasingly larger volume of data over the next 20 years. ATLAS is conducting R&D projects on Amazon Web Services and Google Cloud as complementary resources for distributed computing, focusing on some of the key features of commercial clouds: lightweight operation, elasticity and availability of multiple chip architectures. The proof of concept phases have concluded with the cloud-native, vendoragnostic integration with the experiment’s data and workload management frameworks. Google Cloud has been used to evaluate elastic batch computing, ramping up ephemeral clusters of up to O(100k) cores to process tasks requiring quick turnaround. Amazon Web Services has been exploited for the successful physics validation of the Athena simulation software on ARM processors. We have also set up an interactive facility for physics analysis allowing endusers to spin up private, on-demand clusters for parallel computing with up to 4 000 cores, or run GPU enabled notebooks and jobs for machine learning applications. The success of the proof of concept phases has led to the extension of the Google Cloud project, where ATLAS will study the total cost of ownership of a production cloud site during 15 months with 10k cores on average, fully integrated with distributed grid computing resources and continue the R&D projects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

xGFabric: Coupling Sensor Networks and HPC Facilities with Private 5G Wireless Networks for Real-Time Digital Agriculture

Advanced scientific applications require coupling distributed sensor networks with centralized high-performance computing facilities. Citrus Under Protective Screening (CUPS) exemplifies this need in digital agriculture, where citrus research facilities are instrumented with numerous sensors monitoring environmental conditions and detecting protective screening damage. CUPS demands access to computational fluid dynamics codes for modeling environmental conditions and guiding real-time interventions like water application or robotic repairs. These computing domains have contrasting properties: sensor networks provide low-performance, limited-capacity, unreliable data access, while high-performance facilities offer enormous computing power through high-latency batch processing. Private 5G networks present novel capabilities addressing this challenge by providing low latency, high throughput, and reliability necessary for near-real-time coupling of edge sensor networks with HPC simulations. This work presents xGFabric, an end-to-end system coupling sensor networks with HPC facilities through Private 5G networks. The prototype connects remote sensors via 5G network slicing to HPC systems, enabling real-time digital agriculture simulation.

Digital Agriculture↗

Using Apptainer in a Pilot-based Distributed Workload

GlideinWMS is a pilot and pressure-based workload manager for distributed scientific computing. Many experiments like CMS and Fermilab’s Neutrino experiments use it to provision elastic clusters for their analysis and simulations, split into close to a million concurrent jobs. Most user jobs require containers, and the pilots use Apptainer to set up the desired platform. For the pilots that run as regular batch jobs, Apptainer is safer, lighter, and easier to use than other containerization solutions. Many images used by the pilots are expanded SIF images distributed via the CernVM-FS: this combination is very efficient. At Fermilab, for example, we store on GitHub Dockerfiles that mimic the platform in the worker nodes of local clusters. GitHub workflows build and push the images to Docker Hub, and a service periodically pulls and converts them to the expanded SIF images in the CernVM-FS, so the scientists can find a familiar environment everywhere. Apptainer has also been used to run services inside the pilot jobs, like benchmarks that characterize the worker node being used, or a Triton Inference Server that allows sharing a GPU with all the jobs that run in parallel on a node.

Mambelli, Marco [Fermilab] (ORCID:0000000294892681↗

Recursive Gaussian Process over graphs for Integrating Multi-timescale Measurements in Low-Observable Distribution Systems

The transition to a smarter grid is empowered by enhanced sensor deployments and smart metering infrastructure in the distribution system. Measurements from these sensors and meters can be used for many applications, including distribution system state estimation (DSSE). However, these measurements are typically sampled at different rates and could be intermittent due to losses during the aggregation process. These multi timescale measurements should be reconciled in real-time to perform accurate grid monitoring. This paper tackles this problem by formulating a recursive multi-task Gaussian process (RGP-G) approach that sequentially aggregates sensor measurements. Specifically, we formulate a recursive multi-task GP with and without network connectivity information to reconcile the multi time-scale measurements in distribution systems. Here, the proposed framework is capable of aggregating the multi-time scale measurements batch-wise or in real-time. Following the aggregation of the multi time-scale measurements, the spatial states of the consistent time-series are estimated using matrix completion based DSSE approach. Simulation results on IEEE 37 and IEEE 123 bus test systems illustrate the efficiency of the proposed methods from the standpoint of both multi time-scale data aggregation and DSSE.

42 ENGINEERING↗

aphBO-2GP-3B: a budgeted asynchronous parallel multi-acquisition functions for constrained Bayesian optimization on high-performing computing architecture

High-fidelity complex engineering simulations are often predictive, but also computationally expensive and often require substantial computational efforts. The mitigation of computational burden is usually enabled through parallelism in high-performance cluster (HPC) architecture. Optimization problems associated with these applications is a challenging problem due to the high computational cost of the high-fidelity simulations. In this paper, an asynchronous parallel constrained Bayesian optimization method is proposed to efficiently solve the computationally expensive simulation-based optimization problems on the HPC platform, with a budgeted computational resource, where the maximum number of simulations is a constant. The advantage of this method are three-fold. Firstly, the efficiency of the Bayesian optimization is improved, where multiple input locations are evaluated parallel in an asynchronous manner to accelerate the optimization convergence with respect to physical runtime. This efficiency feature is further improved so that when each of the inputs is finished, another input is queried without waiting for the whole batch to complete. Second, the proposed method can handle both known and unknown constraints. Third, the proposed method samples several acquisition functions based on their rewards using a modified GP-Hedge scheme. The proposed framework is termed aphBO-2GP-3B, which means asynchronous parallel hedge Bayesian optimization with two Gaussian processes and three batches. The numerical performance of the proposed framework aphBO-2GP-3B is comprehensively benchmarked using 16 numerical examples, compared against other 6 parallel Bayesian optimization variants and 1 parallel Monte Carlo as a baseline, and demonstrated using two real-world high-fidelity expensive industrial applications. The first engineering application is based on finite element analysis (FEA) and the second one is based on computational fluid dynamics (CFD) simulations.

97 MATHEMATICS AND COMPUTING↗

MODFLOW6 models used to evaluate potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, Southern Nevada

The hydrograph for well ER-5-3-2 in Frenchman Flat, southern Nevada, has previously unexplained water-level fluctuations. Four, three-dimensional, groundwater models (MODFLOW 6) were developed to evaluate potential stresses and hydrologic conditions affecting the well ER-5-3-2 hydrograph. Four model scenarios were developed that simulated: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) shallow (low transmissivity) and deep (high transmissivity) carbonate rocks, and (4) lateral heterogeneity of carbonate rocks. Input and output files for the four model scenarios are in the model and output directories, respectively. Hydraulic conductivity, specific storage, and wellbore-leakage rates (when simulated) were estimated with parameter estimation (PEST) by minimizing a weighted composite, sum-of-squares objective function. The objective function was informed by measurement and Tikhonov regularization observations. Measurement observations included drawdowns from the constant-rate aquifer test and water-level altitudes measured in well ER-5-3-2 from 2001-2021. Tikhonov regularization informed hydraulic conductivity and specific storage parameters that were insensitive to measurement observations, where homogeneity was the preferred relation. Batch files, executables, and MODFLOW 6, PEST, and post-processing utilities are in the ancillary directory. Supplementary data also are included in the ancillary directory, including site information, high-frequency water-level and aquifer-test data, transmissivity estimates, water-chemistry data, and water-temperature analyses. This USGS data release contains data, analyses, and model files for the simulations and analysis results described in U.S. Geological Survey Scientific Investigations Report (https://doi.org/10.3133/sir20225132).

54 ENVIRONMENTAL SCIENCES↗

Dynamic Parameter Estimation with Physics-based Neural Ordinary Differential Equations

Accurate estimation of dynamic parameters of gen-erators is crucial to building a reliable model for dynamical studies and reliable operation of the power system. This paper develops a physics-based neural ordinary differential equations (ODE) approach to learn the parameters of generator dynamic model using phasor measurement units (PMU) data. We design a physics-based neural network to represent the swing equations of the power system dynamics. A loss function is defined as the difference between dynamic simulation results from the physics-based neural networks and pseudo PMU measurements. The parameters of generator dynamic model are iteratively updated using the neural ODEs and the adjoint method. By exploiting the mini-batch scheme in neural ODE training, the parameter estimation performance is significantly improved. Numerical study results on a 3-machine 9-bus system show that the proposed algorithm outperforms state-of-the-art baseline method in both computation time and dynamic parameter estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sequential Recovery of Critical Metals from Leached Liquor of Processed Spent Lithium-Ion Batteries

The processing and extraction of critical metals from black mass is important to battery recycling. Separation and recovery of critical metals (Co, Ni, Li, and Mn) from other metal impurities must yield purified metal salts, while avoiding substantial losses of critical metals. Solvent extraction in batch experiments were conducted using mixed metal sulphates obtained from the leach liquor obtained from spent and shredded lithium-ion batteries. Selective extraction of Mn 2+ , Fe 3+ , Al 3+ and Cu 2+ from simulated and real leached mixed metals solution was carried out using di-2-ethylhexylphophoric acid (D2EPHA) and Cyanex-272 at varying pH. Further experiments with the preferred extractant (D2EPHA) were performed under different conditions: changing the concentration of extractant, organic to aqueous ratio, and varying the diluents. At optimum conditions (40% v/v D2EPHA in kerosene, pH 2.5, O:A = 1:1, 25 °C, and 20 min), 85% Mn 2+ , 98% Al 3+ , 100% Fe 3+ , and 43% Cu 2+ were extracted with losses of only trace amounts (<5.0%) of Co 2+ , Ni 2+ , and Li + . The order of extraction efficiency for the diluents was found to be kerosene > Exxal-10 >>> dichloromethane (CH 2 Cl 2 ) > toluene. Four stages of stripping of metals loaded on D2EPHA were performed as co-extracted metal impurities were selectively stripped, and a purified MnSO 4 solution was produced. Spent extractant was regenerated after Fe 3+ and Al 3+ were completely stripped using 1.0 M oxalic acid (C 2 H 2 O 4 ).

25 ENERGY STORAGE↗

EdgeAI: Machine learning via direct attached accelerator for streaming data processing at high shot rate x-ray free-electron lasers

We present a case for low batch-size inference with the potential for adaptive training of a lean encoder model. We do so in the context of a paradigmatic example of machine learning as applied in data acquisition at high data velocity scientific user facilities such as the Linac Coherent Light Source-II x-ray Free-Electron Laser. We discuss how a low-latency inference model operating at the data acquisition edge can capitalize on the naturally stochastic nature of such sources. We simulate the method of attosecond angular streaking to produce representative results whereby simulated input data reproduce high-resolution ground truth probability distributions. By minimizing the mean-squared error between the decoded output of the latent representation and the ground truth distributions, we ensure that the encoding layers and resulting latent representation maintains full fidelity for any downstream task, be it classification or regression. We present throughput results for data-parallel inference of various batch sizes, some with throughput exceeding 100 k images per second. We also show in situ training below 10 s per epoch for the full encoder–decoder model as would be relevant for streaming and adaptive real-time data production at our nation’s scientific light sources.

97 MATHEMATICS AND COMPUTING↗

M-Star ® Software Test and Verification

Savannah River Mission Completion (SRMC) currently manages the risk for retained hydrogen in the Defense Waste Processing Facility (DWPF) vessels by implementing a Retained Hydrogen Program. The current program relies on Sludge Batch (SB) 8 Gas Chromatograph data and on conservative assumptions concerning gas release and retention to determine allowable vessel Quiescent time (Q-time). The authors identified M-Star ® CFD as a software that could simulate processes such as fluid flow, heat transfer, species transport, chemical reactions, particle transport, and retained hydrogen gas release. Preliminary simulation results suggest that more realistic assumptions on gas retention and release may be feasible for the DWPF retained hydrogen program. Because of the desire to use the M-Star ® software to perform analyses that support nuclear safety, SRMC has requested Savannah River National Laboratory (SRNL) to upgrade the software classification level of M-Star ® CFD from level D to level A to perform analyses that support nuclear safety.

08 HYDROGEN↗

FY25 MOOSE Usability Improvements: 3D Meshing Capabilities, Initiation of Geometry Support for Monte Carlo Tools, and Enhancement of MOOSE/Workbench User Input Interactions

Usability improvements have been made to MOOSE and Workbench in FY25 to enhance usability and user workflows. Assorted enhancement have been made to MOOSE’s intrinsic meshing capabilities in order to enable more flexible and complex meshing of nuclear reactor systems, in particular for 3D applications. Mesh generators have been added to perform operations such as batch mesh generation, surface mesh generation, and creation of 3D transition layers. These mesh generation capabilities make it much easier to generate high quality non-extruded 3D meshes. Additionally, work to integrate Monte Carlo reactor physics simulations into MOOSE-based multi-physics workflows has reached another milestone with the implementation of the Constructive Solid Geometry (CSG) base framework. This framework lays the foundation for mesh generators to offer the user a generic CSG output option (as opposed to a finite element mesh). To support users, workshop on the MOOSE Reactor Module was delivered which featured hands-on examples using the NEAMS Workbench on INL’s High Performance Computing system. Recent updates to the NEAMS Workbench, WASP, and the MOOSE language server have introduced several improvements aimed at making MOOSE-based simulation setup and input management faster, more accurate, and easier to use. Key capabilities that have been added include multi-tab-stop autocompletion, visual input diagnostics, developer-directed data visualizations, upgraded ParaView integration, and Workspace-level file tracking. Together, these changes make it easier for users to build, validate, and manage complex MOOSE-based simulation models — especially those involving reusable components, included files, and datasets. The improvements are designed to save time, reduce input errors, and help users get to a successful simulation run faster, with more confidence in the results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accumulation mechanisms for contaminants on weak-base hybrid ion exchange resins

Mechanism of hexavalent chromium removal (Cr(VI) as CrO 4 2- ) by the weak-base ion exchange (IX) resin ResinTech® SIR-700-HP (SIR-700) from simulated groundwater is assessed in the presence of radioactive contaminants iodine-129 (as IO 3 - ), uranium (U as uranyl UO 2 2+ ), and technetium-99 (as TcO 4 - ), and common environmental anions sulfate (SO 4 2- ) and chloride (Cl-). Batch tests using the acid sulfate form of SIR-700 demonstrated Cr(VI) and U(VI) removal exceeded 97%, except in the presence of high SO 4 2- concentrations (536 mg/L) where Cr(VI) and U(VI) removal decreased to ≥ 80%. However, Cr(VI) removal notably improved with co-mingled U(VI) that complexes with SO 4 2- at the protonated amine sites. These U–SO 4 2- complexes are integral to U(VI) removal, as confirmed by the decrease in U(VI) removal (<40%) when the acid chloride form of SIR-700 was used instead. Solid phase characterization revealed that CrO 4 2- is removed by IX with SO 4 2- complexes and/or reduced to amorphous Cr(III)(OH) 3 at secondary alcohol sites. Tc(VII)O 4 - and I(V)O 3 - also undergo chemical reduction, following a similar removal mechanism. Oxyanion removal preference is determined by the anion reduction potential (CrO 4 2 - >TcO 4 - >IO 3 - ), geometry, and charge density. For these reasons, 39% and 69% of TcO 4 - and 17% and 39% of IO 3 - are removed in the presence and absence of Cr(VI), respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid Sorption-Enhanced Haber–Bosch Process

The use of a liquid sorbent in a traditional Haber-Bosch process enables significant improvements in energy efficiency and potential cost savings for arguably the most important chemical process on the planet. The approach presented in this report employs an incompressible liquid sorbent that absorbs and releases ammonia (NH 3 ) under specific conditions. To achieve this, we investigate reactions of ammonia and pure phosphoric acid (H 3 PO 4 , PA), which rapidly neutralize to form an equilibrated solution of monoammonium phosphate (MAP) and diammonium phosphate (DAP) that functions as a reversible and regenerable sorbent. Through intimate contact of the gas-phase Haber-Bosch reaction mixture with this liquid absorbent, complete equilibrium uptake may be achieved in an appropriately sized separator, and facile separation occurs through the use of independent liquid and gas phases. Following depressurization and release of the ammonia product, only the incompressible fluid needs to be repressurized and returned to the reactor. This study documents proof-of-concept absorption and desorption experiments carried out in 75 mL batch reactors, predominantly charged with precise MAP and DAP mixtures that equilibrate at process-relevant temperatures and pressures. We then assemble the first thermodynamic relationships that underlie this advantaged separation strategy, validated by reactive force field (ReaxFF) interatomic potential simulations, and benchmarked with traditional separation routes via process modeling and technoeconomic analysis. The scale of energy consumption in the century-old Haber-Bosch process is massive, and the elegant liquid sorption approach reported here offers opportunities to enhance its energy efficiency for the next frontier of ammonia synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geant4 Visualization

Geant4 offers versatile visualization tools to study detector geometry, particle trajectories, and hits. These tools support both interactive and batch modes (e.g., vis.mac). With a unified interface, Geant4 integrates various graphics systems such as G4VisExecutive, OpenGL, Qt, HepRApp, DAWN, VTK, and ToolsSG. Advanced examples, including the ICRP145 Human Phantoms and the gMocren tool, highlight its applications in medical research and simulation and as well as radiation protection and shielding studies. These features make Geant4 a robust framework for exploring particle interactions and supporting domain-specific use cases.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗