Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity 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 559 records · Page 31

Experimental, numerical and analytical evaluation of j × B -thrust for fast-liquid-metal-flow divertor systems of nuclear fusion devices

Abstract Divertor systems of fusion devices are exposed to intense heat loads from plasmas, which degrade solid plasma-facing components. Fast liquid metal (LM) flow divertors may be more advantageous for this purpose but have risk of piling due to intense magnetohydrodynamic (MHD) drag. However, severe deceleration of the flow could be countered with the injection of currents that are transverse to external magnetic fields, allowing to thrust the flow with j × B (Lorentz) forces. Given that the injection of currents as an approach to propel LM-divertor flows has remained experimentally understudied, this article focuses on the evaluation of j × B -thrust and finding its drawbacks. j × B -thrust was experimentally tested with free-surface-LM flows, a vertical magnetic field and an externally applied current. Experiments were reviewed with a theoretical model, showing agreement in the trends of theory and experiments. Full 3D-MHD-free-surface-flow simulations were also performed with FreeMHD and confirmed the sensitivity to unstable flow behavior in LM systems when applying external currents. Furthermore, excessive power requirements are expected for the implementation of j × B -thrust at the reactor scale, making these systems inefficient for commercial devices. This paper evidences that the simple operation of a LM-flow divertor with j × B -thrust, without any of the instabilities caused from reactor plasmas or parasitic currents, already presents intrinsic challenges.

Magnetohydrodynamics↗

PSA 2025 Presentation: "Modeling and Sensitivity Analysis of a Generation IV Pebble Bed Reactor Using MELCOR 2.2"

Accompanying the advancement of reactor technologies is the need for computational modeling and simulation to predict their behavior under normal operating conditions and accident scenarios. New Generation IV reactor designs which employ non-conventional fuel have a particular need for modeling the behavior and release of radionuclides and other material from the fuel. In this work, MELCOR version 2.2, a system-level safety and accident scenario code developed by Sandia National Laboratories, was used to model a 200-MWth pebble bed modular reactor and calculate the inventories of circulating and deposited graphite, metal dust, and elemental components released from the fuel elements. A base case modeling the reactor under standard operating conditions was calculated using MELCOR and the inventories were extrapolated to 30 years of operation time using a logarithmic regression fit. A sensitivity analysis was also performed in which several key parameters for the base case model were modified to explore the effect of these changes on the inventories calculated by MELCOR. A set of transient scenario simulations for a depressurized loss of forced cooling (DLOFC) accident were also performed. The results of the sensitivity analysis and transient simulations are reported and discussed in relation to the modeling techniques used for this study.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Nuclear Data Sensitivity Study for the EBR-II Fast Reactor Benchmark Using SCALE with ENDF/B-VII.1 and ENDF/B-VIII.0

The EBR-II benchmark, which was recently included in the International Handbook of Evaluated Reactor Physics Benchmark Experiments, served as a basis for assessing the performance of the SCALE code system for fast reactor analyses. A reference SCALE model was developed based on the benchmark specifications. Great agreement was observed between the eigenvalue calculated with this SCALE model and the benchmark eigenvalue. To identify potential gaps and uncertainties of nuclear data for the simulation of various quantities of interest in fast spectrum systems, sensitivity and uncertainty analyses were performed for the eigenvalue, reactivity effects, and the radial power profile of EBR-II using the two most recent ENDF/B nuclear data library releases. While the nominal results are consistent between the calculations with the different libraries, the uncertainties due to nuclear data vary significantly. The major driver of observed uncertainties is the uncertainty of the 235 U ( n,γ ) reaction. Since the uncertainty of this reaction is significantly reduced in the ENDF/B-VIII.0 library compared to ENDF/B-VII.1, the obtained output uncertainties tend to be smaller in ENDF/B-VIII.0 calculations, although the decrease is partially compensated by increased uncertainties in 235 U fission and ν ¯ .

42 ENGINEERING↗

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BEYONDPLANCK IV. On end-to-end simulations in CMB analysis — Bayesian versus frequentist statistics

End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely, how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data or as a random realization independent from the data. We refer to these as posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before BEYONDPLANCK, most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BEYONDPLANCK represents the first end-to-end CMB simulation framework that is able to generate both types of simulations and these new capabilities have brought this topic to the forefront. The BEYONDPLANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level BEYONDPLANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating the Phonon Collection Efficiency in Kinetic Inductance Phonon-Mediated Detectors

Kinetic inductance phonon-mediated (KIPM) detectors use microwave kinetic inductance detectors (MKIDs) to read out phonon signals in the substrate. They are a promising class of detectors to be used in light dark matter (DM) searches due to their potential eV-scale sensitivity and native frequency-domain multiplexability. In order to improve upon the design of these detectors, simulations are needed to understand the effects of detector design on measurable physical parameters including the phonon collection efficiency, $\eta_{ph}$, which is defined as the ratio of phonon energy detected by the sensitive target element to the incident energy deposited in the substrate. This work simulates the phonon collection efficiency for two existing detector architectures to study phonon absorption at substrate-sensor interfaces and phonon loss to non-sensitive detector elements.

Dang, Stella Q.↗

On the Effect of Historical SST Patterns on Radiative Feedback

Here, we investigate the dependence of radiative feedback on the pattern of sea-surface temperature (SST) change in fourteen Atmospheric General Circulation Models (AGCMs) forced with observed variations in SST and sea-ice over the historical record from 1871 to near-present. We find that over 1871-1980, the Earth warmed with feedbacks largely consistent and strongly correlated with long term climate sensitivity feedbacks (diagnosed from corresponding atmosphere-ocean GCM abrupt 4xCO2 simulations). Post 1980 however, the Earth warmed with unusual trends in tropical Pacific SSTs (enhanced warming in the west, cooling in the east) and cooling in the Southern Ocean that drove climate feedback to be uncorrelated with – and indicating much lower climate sensitivity than – that expected for long-term CO 2 increase. We show that these conclusions are not strongly dependent on the AMIP II SST dataset used to force the AGCMs, though the magnitude of feedback post 1980 is generally smaller in nine AGCMs forced with alternative HadISST1 SST boundary conditions. We quantify a ‘pattern effect’ (defined as the difference between historical and longterm CO 2 feedback) equal to 0.48 ± 0.47 [5-95%] W m -2 K -1 for the time-period 1871-2010 when the AGCMs are forced with HadISST1 SSTs, or 0.70 ± 0.47 [5-95%] W m -2 K -1 when forced with AMIP II SSTs. Assessed changes in the Earth’s historical energy budget agree with the AGCM feedback estimates. Furthermore satellite observations of changes in top-of-atmosphere radiative fluxes since 1985 suggest that the pattern effect was particularly strong over recent decades but may be waning post 2014.

54 ENVIRONMENTAL SCIENCES↗

Cooperative multivalent receptor binding promotes exposure of the SARS-CoV-2 fusion machinery core

The molecular events that permit the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to bind and enter cells are important to understand for both fundamental and therapeutic reasons. Spike proteins consist of S1 and S2 domains, which recognize angiotensin-converting enzyme 2 (ACE2) receptors and contain the viral fusion machinery, respectively. Ostensibly, the binding of spike trimers to ACE2 receptors promotes dissociation of the S1 domains and exposure of the fusion machinery, although the molecular details of this process have yet to be observed. We report the development of bottom-up coarse-grained (CG) models consistent with cryo-electron tomography data, and the use of CG molecular dynamics simulations to investigate viral binding and S2 core exposure. We show that spike trimers cooperatively bind to multiple ACE2 dimers at virion-cell interfaces in a manner distinct from binding between soluble proteins, which processively induces S1 dissociation. We also simulate possible variant behavior using perturbed CG models, and find that ACE2-induced S1 dissociation is primarily sensitive to conformational state populations and the extent of S1/S2 cleavage, rather than ACE2 binding affinity. These simulations reveal an important concerted interaction between spike trimers and ACE2 dimers that primes the virus for membrane fusion and entry.

59 BASIC BIOLOGICAL SCIENCES↗

Future projections of storm surge in Hurricane Katrina and sensitivity to meteorological forcing resolution

In this study, we investigated whether and how the storm surge induced by Hurricane Katrina could change if it occurs in a future warmer climate, and the sensitivity of the changes to atmospheric forcing resolution. Climate model simulations of Hurricane Katrina at 27 km, 4.5 km, and 3 km resolutions were used to drive storm surge simulations in historical and future climates using the ADvanced CIRCulation (ADCIRC) model. We found that peak surge height increased significantly in the future with all forcing resolutions. However, the future projection is 22% greater in the 3 km forcing, typical of regional climate models, compared to the 27 km forcing, typical of state-of-the-art global climate models. Additionally, the spatial extent of the future change is highly sensitive to forcing resolution, extending most broadly under the 27 km forcing. Furthermore, we found that storm surge duration decreases in the future with all forcing resolutions due to increasing TC translation speed and decreasing ocean lifetime. However, the future change in the surge duration is sensitive to the forcing resolution, decreasing by 31% in the 27 km forcing and 6% in the 3 km forcing.

54 ENVIRONMENTAL SCIENCES↗

Influence of mass ablation on ignition and burn propagation in layered fusion capsules

After decades of research, recent laser-driven inertial fusion experiments have demonstrated rapid progress toward achieving thermonuclear ignition using capsule designs with cryogenic fuel layers. The ignition physics for these layered capsules involves a complex interplay between the dynamically forming hot spot and the dense surrounding fuel. Using analytic theory and numerical simulations, we demonstrate that the mass ablation rate into the hot spot depends sensitively upon the temperature of the dense fuel, resulting in ablative inflows up to [Formula: see text] faster than previous estimates. This produces an enthalpy flux into the hot spot that plays a critical role in controlling the hot spot temperature, the ignition threshold, and the subsequent burn propagation. The net influence of mass ablation on the ignition threshold is regulated by a dimensionless parameter that depends upon the temperature of the dense fuel. As a consequence, the ignition threshold is sensitive to any mechanism that heats the dense fuel, such as neutrons or radiation emitted from the hot spot. These predictions are confirmed using radiation-hydrodynamic simulations for a series of capsules near ignition conditions. This analysis may have relevance for understanding the variable performance of recent experiments and for guiding new capsule designs toward higher fusion yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A full-ring variable-aperture cadmium zinc telluride system for whole-body single photon emission computed tomography: realistic simulations with phantoms

Single photon emission computed tomography (SPECT) is an imaging modality that has demonstrated its utility in a number of clinical indications. Despite this progress, a high sensitivity, high spatial resolution, multi-tracer SPECT with a large field of view suitable for whole-body imaging of a broad range of radiotracers for theranostics is not available. Purpose We have designed a cadmium zinc telluride (CZT) variable-aperture full-ring SPECT scanner instrumented with a broad-energy tungsten collimator intended to fill this technological gap. The final purpose is to provide a multi-tracer solution for brain and whole-body imaging. Our static SPECT scanner breaks the paradigm of the standard dual- and triple-head rotational SPECT systems, utilizing a larger detector area in each scan increasing the sensitivity. We provide a demonstration of the performance of our design using a realistic model of our detector with simulated body-sized 99mTc phantoms. Methods We developed a realistic model of our detector by using a combination of a Geant4 Monte Carlo simulation and a CZT detector response model based on a finite element model. Our approach models the characteristic low-energy tail effect in CZT that noticeably affects the sensitivity and the quality of the scatter correction in CZT detectors. We implement a modified dual energy window scatter correction adapted to include the CZT low-energy tail effect. A dedicated correction is also developed to eliminate the undesirable truncation observed in images given the presence of detector edges and gaps between detectors, due to the non-rotational nature of our device. Corrections for the attenuation, detector response and the presence of collimators are also included. The images are reconstructed using the maximum-likelihood expectation-maximization algorithm implemented in the reconstruction open software STIR. Detector and reconstruction performance are characterized with a Derenzo phantom and a body-sized National Electrical Manufacturers Association (NEMA) Image Quality (IQ) phantom containing 99mTc. Results Our SPECT design can resolve 6.4mm rods in a Derenzo phantom and obtain a good image contrast with the IQ phantom. Explicit testing of the gap and edge correction is provided, showing an excellent performance in eliminating the image truncation artifacts. Our modified scatter correction shows no overestimation of the contrast-recovery ratio for our realistic CZT detector model, as opposed to the cases without correction and with a standard dual-energy window scatter correction. Conclusions In this paper, we further demonstrate the performance of our design for whole-body imaging purposes. This adds to our previous demonstration of improved qualitative and quantitative 99mTc imaging for brain perfusion and 123I imaging for dopamine transport with respect to state-of-the-art NaI dual-head cameras. We show that our design performs similarly to the VERITON SPECT from Spectrum Dynamics, a commercial full-ring CZT SPECT camera, with the potential advantage of the broader energy range of application given by our custom-design tungsten collimators. Here, our device combines high sensitivity and image resolution with a broad-energy imaging application for the purpose of clinical imaging and theranostics of emerging radionuclides.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Simulating the Phonon Collection Efficiency in KIPMDs

Kinetic inductance phonon-mediated (KIPM) detectors use microwave kinetic inductance detectors (MKIDs) to read out phonon signals in the substrate. They are a promising class of detectors to be used in light dark matter (DM) searches due to their potential eV-scale sensitivity and native frequency-domain multiplexability. In order to improve upon the design of these detectors, simulations are needed to understand the effects of detector design on measurable physical parameters including the phonon collection efficiency, $\eta_{ph}$, which is defined as the ratio of phonon energy detected by the sensitive target element to the incident energy deposited in the substrate. This work simulates the phonon collection efficiency for a KIPMD currently operated at the Northwestern EXperimental Underground Site (NEXUS).

Dang, Stella Q.↗

Neutron-antineutron oscillation search with MicroBooNE and DUNE

The Deep Underground Neutrino Experiment (DUNE) is an international project aiming at neutrino physics and astrophysics and a search for phenomena predicted by theories beyond the standard model (BSM). The excellent imaging capability of Liquid Argon Time Projection Chamber (LArTPC) technology, particle tracking and identification utilized in the Far Detector, as well as the Far Detector size and underground placement, allow the experiment to achieve high sensitivity to various rare processes. BSM theories predict the existence of baryon number non-conservation effects, in particular when the baryon number changes by 2. Here we discuss the sensitivity of DUNE to neutron-antineutron oscillation. With full event simulation and reconstruction using the LArSoft package, we have investigated the background to potential signal events from atmospheric neutrino interactions and particle misidentification, and utilized machine learning techniques to enhance the discrimination between signal and background. The methodologies being developed for a high-sensitivity search for neutron-antineutron oscillation in DUNE can also be demonstrated with the currently running MicroBooNE LArTPC. We discuss progress on demonstrating the developed techniques with the first-ever search for neutron-antineutron oscillation in a LArTPC using MicroBooNE data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Design of Premoderated Pulsed-Neutron Die-Away Experiments for Validating Beryllium Thermal Neutron Scattering Laws

Accurate thermal scattering laws (TSLs) are essential for reliable neutron transport simulations, particularly for systems with low-energy neutrons. Pulsed-neutron die-away (PNDA) experiments provide a highly sensitive platform for TSL validation but face challenges from fast neutron leakage and material constraints in nonhydrogenous targets such as beryllium. This work investigates the use of high-density polyethylene (HDPE) as a premoderator to reduce leakage, improve thermalization, and shorten experimental run times while preserving sensitivity to the Be TSL. Here, three premoderated configurations, encapsulation, slab, and interstitial, were evaluated using MCNP6.3® simulations with ENDF/B-VIII.0 nuclear data. The effects of the Be TSL was quantified by comparing decay constants from simulations with and without 𝑆⁡(𝛼,𝛽) treatments for both beryllium and HDPE. The results show that premoderation enables the use of as little as 5.20% of the beryllium volume required for an unreflected geometry, achieving uncertainties below 0.5% and reducing run time by up to 72.2% in the studied configurations. The interstitial configuration achieved the best balance between low statistical uncertainty and high TSL sensitivity, outperforming the slab and encapsulation designs. These findings demonstrate that carefully optimized premoderation can significantly enhance the efficiency, feasibility, and precision of PNDA experiments for Be TSL validation.

Integral Experiments↗

Sensitivity and uncertainty of the IFR-1 BISON benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered for adoption by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This paper summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment.For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. Here, the resulting variations in benchmark predictions were normalized to calculate sensitivities. These sensitivities were then used to inform input selections for the UQ.The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of the benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because the confidence intervals associated with benchmark predictions were wide. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

MRT 7365: Power flow physics and key physics phenomena

The Z accelerator at Sandia National Laboratories conducts z-pinch experiments at 26 MA in support of DOE missions in stockpile stewardship, dynamic materials, fusion, and other basic sciences. Increasing the current delivered to the z-pinch would extend our reach in each of these disciplines. To achieve increases in current and accelerator efficiency, a fraction of Z’s shots are set aside for research into transmission-line power flow. These shots, with supporting simulations and theory, are incorporated into this Advanced Diagnostics milestone report. The efficiency of Z is reduced as some portion of the total current is shunted across the transmission-line gaps prior to the load. This is referred to as “current loss”. Electrode plasmas have long been implicated in this process, so the bulk of dedicated power-flow experiments are designed to measure the plasma environment. The experimental analyses are enhanced by simulations conducted using realistic hardware and Z voltage pulses. In the same way that diagnostics are continually being improved for sensitivity and resolution, the modeling capability is continually being improved to provide faster and more realistic simulations. The specifics of the experimental hardware, diagnostics, simulations, and algorithm developments are provided in this report. The combined analysis of simulation and data confirms that electrode plasmas have the most detrimental impact on current delivery. Experiments over the last three years have tested the theoretical current-loss mechanisms of enhanced ion current, plasma gap closure, and Hall-related current. These mechanisms are not mutually exclusive and may be coincident in the final feed as well as in upstream transmission lines. The final-feed geometries tested here, however, observe lower-density plasmas without dominant ion currents which is consistent with a Hall-related current. The picture of plasma formation and transport formed from experiment and simulation is informing hardware designs being fielded on Z now and being proposed for the Next-Generation Pulsed Power (NGPP) facility. In this picture, the strong magnetic fields that heat the electrodes above particle emission thresholds also confine the charged particles near the surface. Some portion of the plasmas thus formed is transported into the transmission-line gap under the force of the electric field, with aid from plasma instabilities. The gap plasmas are then transported towards the load by a cross-field drift, where they accumulate and contribute to a likely Hall-related cross-gap current. The achievements in experimental execution, model validation, and physical analysis presented in this report set the stage for continued progress in power flow and load diagnostics on Z. The planned shot schedule for Z and Mykonos will provide data for extrapolation to higher current to ensure the predicted performance and efficiency of a NGPP facility.

43 PARTICLE ACCELERATORS↗

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗