Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core 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 1,765 records · Page 98

Coupled Monte Carlo neutronics and thermal hydraulics for power reactors

The availability of high performance computing resources enables more and more the use of detailed Monte Carlo models even for full core power reactors. The detailed structure of the core can be described by lattices, modeled by so-called repeated structures e.g. in Monte Carlo codes such as MCNP5 or MCNPX. For cores with mainly uniform material compositions, fuel and moderator temperatures, there is no problem in constructing core models. However, when the material composition and the temperatures vary strongly a huge number of different material cells must be described which complicate the input and in many cases exceed code or memory limits. The second problem arises with the preparation of corresponding temperature dependent cross sections and thermal scattering laws. Only if these problems can be solved, a realistic coupling of Monte Carlo neutronics with an appropriate thermal-hydraulics model is possible. In this paper a method for the treatment of detailed material and temperature distributions in MCNP5 is described based on user-specified internal functions which assign distinct elements of the core cells to material specifications (e.g. water density) and temperatures from a thermal-hydraulics code. The core grid itself can be described with a uniform material specification. The temperature dependency of cross sections and thermal neutron scattering laws is taken into account by interpolation, requiring only a limited number of data sets generated for different temperatures. Applications will be shown for the stationary part of the Purdue PWR benchmark using ATHLET for thermal- hydraulics and for a generic Modular High Temperature reactor using THERMIX for thermal- hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Visualizing a Vulnerability: Its Connections to Hardware and Software

All Hazards Analysis (AHA) is a framework developed by Idaho National Laboratory that provides capabilities to collect, store, analyze, and visualize critical infrastructure information. A core function of AHA is its ability to simulate faults or outages in networks of infrastructure originating from a plethora of causes, ranging from natural disasters to cyberattacks. AHA utilizes Hardware and Software Bills of Material (HBOM and SBOM, respectively) along with Known Exploited Vulnerabilities (KEVs) to document the potential attack vectors for each piece of infrastructure. The objective of this contribution to AHA was to create a visualization tool that could capture the small details held in each individual artifact as well as preserve the large-scale connections that link them together to aid threat modeling.

58 GEOSCIENCES↗

Experimental studies and simulations of hydrogen pellet ablation in the stellarator TJ-II

Plasma core fuelling is a key issue for the development of steady-state scenarios in large magnetically-confined fusion devices, in particular for helical-type machines. At present, cryogenic pellet injection is the most promising technique for efficient fuelling. Here, pellet ablation and fuelling efficiency experiments, using a compact pellet injector, are carried out in electron cyclotron resonance and neutral beam injection heated plasmas of the stellarator TJ-II. Ablation profiles are reconstructed from light emissions collected by silicon photodiodes and a fast-frame camera system, under the assumptions that such emissions are loosely related to the ablation rate and that pellet radial acceleration is negligible. In addition, pellet particle deposition and fuelling efficiency are determined using density profiles provided by a Thomson scattering system. Furthermore, experimental results are compared with ablation and deposition profiles provided by the HPI2 pellet code, which is adapted here for the stellarators Wendelstein 7-X (W7-X) and TJ-II. Finally, the HPI2 code is used to simulate ablation and deposition profiles for pellets of different sizes and velocities injected into relevant W7-X plasma scenarios, while estimating the plasmoid drift and the fuelling efficiency of injections made from two W7-X ports.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling the Ex-core Detector Response into a Production Option in MPACT

In the past few years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has overseen a successful activity to extend the radiation transport capability from the reactor in-core geometry to the ex-core geometry by coupling the deterministic neutronics code MPACT, with the Monte Carlo code Shift. However, transient problems such as RIAs require an efficient prediction of the ex-core detector response to determine the timeframe for a reactor trip. Therefore, in addition to the high-fidelity MPACT-Shift coupled calculation, separate efforts have been made to explore simplified methods for ex-core detector response with MPACT standalone calculation. In a previous milestone L3_RTM.MCH.P17.06, several simplified methods for ex-core detector response were developed in MPACT. In particular, the double kernel method is based on the default MPACT core geometry, so extensive modeling or transport calculation with full reflector and ex-core details are not needed in this method. Since the 1-D diffusion and neutron streaming kernels are used for extended ex-core calculations, the extra computational efforts are very small. The work implemented in the current milestone provides a fully integrated option for MPACT to approximately calculate the ex-core detector response using the double kernel method. To overcome the difficulty of obtaining the consistent diffusion coefficients (these were obtained through separate CMFD calculation and edits), the 1-D diffusion solver is replaced by a 1-D Sn transport solver to compute the neutron flux at the vessel outer surface. The new VERA inputs for ex-core detector are processed in MPACT. The detector signals are linked to the SCRAM logic as an option to trip the core during transient calculation. Two sets of 3-D quarter core problems are run with MPACT standalone and VeraShift to verify the radial and axial effects on detector response by varying the moderator density and control rod position. For the cases with moderator density changes, the detector responses from MPACT are within 2 sigma error of VeraShift. For the cases at different control rod positions, the axial offsets computed from MPACT are within 2% error of VeraShift results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mapping the sawtooth

The sawtooth oscillation involves a topological change of the magnetic field lines in the core of a tokamak discharge. In this work, we simulate a steady state sawtoothing discharge using the M3DC 1 -code. During the sawtooth the topological structure of the magnetic field is analyzed by tracking the fixed points of the field line map: the mapping induced by the field lines going once around the tokamak. We numerically locate the fixed points of the field line map (of which the magnetic axis is one) and use them construct a Poincaré section that resolves the core region and the (1, 1) island. Finally, by following the evolution of the magnetic field through the motion and bifurcation of the fixed points of the field line map we get an intuitive picture of the processes that take place during a Kadomtsev style sawtooth crash.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Veres, Joseph P.↗

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Jorgenson, Philip C. E.↗

TPSAS-NF1676L-10767-DND

Version 2.1ae of the NASA CFD code OVERFLOW has been used to simulate the flow past tandem cylinders. The multiblock grid system consisted of 15 overset meshes with a total grid count of 15.7 million points. The grid extended 3 cylinder diameters in the spanwise direction with periodic boundary conditions at the ends. Riemann conditions were used at the outer boundaries, and a no-slip condition on the cylinder surfaces. The 3rd order HLLC upwind spatial discretization scheme was used with a dual-time stepping algorithm. A modified version of the Delayed Detached Eddy Simulation hybrid RANS/LES turbulence model proposed by Spalart was used in a fully turbulent mode. The code was run on 53 Core 2 Duo nodes using MPI across nodes and OPENMP within each node. Simulation results will be compared with the available experimental data.

David P. Lockard↗

Sensitivity Analysis of Lbloca In Response To 10 Cfr50.46c New Rulemaking

The Nuclear Regulatory Commission (NRC) is proposing a new rulemaking on emergency core system/loss-of-coolant accident (LOCA) performance analysis. In the proposed rulemaking, designated as 10 CFR 50.46c, the US NRC puts forward an equivalent cladding oxidation criterion as a function of cladding pre-transient hydrogen content. The proposed rulemaking imposes more restrictive and burnupdependent cladding embrittlement criteria; consequently more fuel rods need to be analyzed under LOCA conditions to maintain the safety margin, in contrast to the current practice with which only one hot rod needs to be analyzed. New multi-physics analysis methods are required to provide a thorough characterization of the reactor core in order to identify the locations of the limiting rods and quantify safety margins under LOCA conditions. The U.S. DOE’s Light Water Reactor Sustainability (LWRS) Program has initiated a project to develop multi-physics analytical capabilities, called LOTUS, to support the industry in the transition to the proposed rule. An approach to uncertainty quantification and sensitivity analysis with LOTUS was developed. A typical four-loop PWR plant model was developed and simulated by RELAP5-3D with inputs generated from core design and fuel performance analyses, and uncertainty quantification and sensitivity analysis were performed with seventeen uncertain input parameters. The maximum equivalent cladding reacted ratio (ECRR) and peak clad temperature ratio (PCTR) were selected as the figures of merit (FOM). Pearson, Spearman, partial correlation coefficients, and Sobol indices were considered for all of the figures of merit in the sensitivity analysis.

10CFR50.46c↗

METHODOLOGY AND APPLICATION OF PHYSICAL SECURITY EFFECTIVENESS BASED ON DYNAMIC FORCE-ON-FORCE MODELING

This paper describes ongoing work within the Light Water Reactor Sustainability (LWRS) Program at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants (NPPs). It reviews the conservatisms in conventional physical security posture and regulations. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). The dynamic assessment methodology leverages EMRALD to process results of force-on-force (FOF) simulations and crediting safety mitigation actions from probabilistic risk assessment (PRA) models as well as diverse and flexible coping strategies (FLEX) mitigation strategies. Timing information from these simulations are compared against the available time to perform mitigations obtained from Reactor Excursion and Leak Analysis Program (RELAP5) simulations. To illustrate the methodology, a station blackout (SBO) attack scenario was modeled in commercially available FOF simulation tools. The simulation results provide valuable insights into possible attack outcomes and as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD, and were dependent on the attack outcomes by considering human operator uncertainties. RELAP5 simulations incorporating human and hardware uncertainties were performed to estimate the distribution of time-to-core damage. The results demonstrate that, even in the extreme case of a successful adversarial attack, plant mitigation strategies provide significantly high-likelihood of preventing radiological release. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the NPPs to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A satellite observational and numerical study of precipitation characteristics in western North Atlantic tropical cyclones

Special Sensor Microwave/Imager (SSM/I) observations were used to examine the spatial and temporal changes of the precipitation characteristics of tropical cyclones. SSM/I observations were also combined with the results of a tropical cyclone numerical model to examine the role of inner-core diabatic heating in subsequent intensity changes of tropical cyclones. Included in the SSM/I observations were rainfall characteristics of 18 named western North Atlantic tropical cyclones between 1987 and 1989. The SSM/I rain-rate algorithm that employed the 85-GHz channel provided an analysis of the rain-rate distribution in greater detail. However, the SSM/I algorithm underestimated the rain rates when compared to in situ techniques but appeared to be comparable to the rain rates obtained from other satellite-borne passive microwave radiometers. The analysis of SSM/I observations found that more intense systems had higher rain rates, more latent heat release, and a greater contribution from heavier rain to the total tropical cyclone rainfall. In addition, regions with the heaviest rain rates were found near the center of the most intense tropical cyclones. Observational analysis from SSM/I also revealed that the greatest rain rates in the inner-core regions were found in the right half of fast-moving cyclones, while the heaviest rain rates in slow-moving tropical cyclones were found in the forward half. The combination of SSM/I observations and an interpretation of numerical model simulations revealed that the correlation between changes in the inner core diabetic heating and the subsequent intensity became greater as the tropical cyclones became more intense.

Rodgers, Edward B.↗

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reactivity Feedback Modeling in SAM

System Analysis Module (SAM) is being developed at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report presents the development of Point-Kinetics modeling and reactivity feedback mechanisms in SAM. It also presents a new capability for predicting the thermal expansion in various structural components by coupling SAM with an external thermomechanics module (Tensor Mechanics) from the MOOSE framework. The analysis of the transient behavior of a nuclear reactor requires the coupled simulation of reactor kinetics and thermal-hydraulics of the reactor core, especially for those unprotected transients where the reactor scram system may not function properly. The point kinetics model has been widely used for reactor safety analysis due to its simplicity to capture the transient behavior of the reactor. Various reactivity feedback models have been developed and integrated with the Point-Kinetics module, including fuel axial expansion, core radial expansion, fuel Doppler, and coolant density reactivity. The reactivity feedback models in SAM are similar to the respective models used in SAS4A/SASSYS-1. This report first presents the brief theory of the Point-Kinetics module and reactivity feedback models. A number of verification tests have been performed where the code simulations are compared to the analytical model results. The reactivity feedback due to the thermal deformation, such as the fuel axial expansion and core radial expansion, is important for SFR transient analysis. Simplified thermal expansion models for the fuel pin and reactor constrain system (e.g. grid plate) are developed and verified in SAM. Additionally, a coupling interface is developed to couple SAM with external thermomechanical analysis modules for more accurate predictions of the thermal expansion of different components during the transients. The current coupling interface has been tested with the Tensor Mechanics module from MOOSE. These point kinetics and reactivity feedback modeling capabilities have also been demonstrated by simulating the early stage of the unprotected loss-of-flow (ULOF) accident in the Advanced Burner Test Reactor (ABTR). Both the stand-alone SAM and coupled SAM and Tensor Mechanics simulations are performed. It is confirmed that the major physics phenomena in the heat transport system of the ABTR reactor are captured by SAM, and the point kinetics model, reactivity feedback models, and the coupling schemes are working well as expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Plasma Simulation Program

Many others in the fusion energy and advanced scientific computing communities participated in the development of this plan. The core planning team is grateful for their important contributions. This summary is meant as a quick overview the Fusion Simulation Program's (FSP's) purpose and intentions. There are several additional documents referenced within this one and all are supplemental or flow down from this Program Plan. The overall science goal of the DOE Office of Fusion Energy Sciences (FES) Fusion Simulation Program (FSP) is to develop predictive simulation capability for magnetically confined fusion plasmas at an unprecedented level of integration and fidelity. This will directly support and enable effective U.S. participation in International Thermonuclear Experimental Reactor (ITER) research and the overall mission of delivering practical fusion energy. The FSP will address a rich set of scientific issues together with experimental programs, producing validated integrated physics results. This is very well aligned with the mission of the ITER Organization to coordinate with its members the integrated modeling and control of fusion plasmas, including benchmarking and validation activities. [1]. Initial FSP research will focus on two critical Integrated Science Application (ISA) areas: ISA1, the plasma edge; and ISA2, whole device modeling (WDM) including disruption avoidance. The first of these problems involves the narrow plasma boundary layer and its complex interactions with the plasma core and the surrounding material wall. The second requires development of a computationally tractable, but comprehensive model that describes all equilibrium and dynamic processes at a sufficient level of detail to provide useful prediction of the temporal evolution of fusion plasma experiments. The initial driver for the whole device model will be prediction and avoidance of discharge-terminating disruptions, especially at high performance, which are a critical impediment to successful operation of machines like ITER. If disruptions prove unable to be avoided, their associated dynamics and effects will be addressed in the next phase of the FSP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Aerosol Model Benchmarking Repository: A toolkit for model intercomparison

The Aerosol Model Benchmarking Repository and Standards (AMBRS) project was initiated to provide tools and to establish community standards for benchmarking aerosol models. This report describes a set of open-source tools for building, running, and analyzing aerosol box model simulations in a standardized framework. The framework consists of three core components: AMBuilder, a CMake-based build system that compiles supported models consistently; AMBRS, a Python module that defines unified numerical experiments and executes them with aligned inputs; and PyParticle, an aerosol analysis package that standardizes output, computes diagnostics, and visualizes simulation results. Together, these tools enable reproducible intercomparison of aerosol schemes and support process-level evaluation of how model simplifications affect predictions of size distributions, cloud condensation nuclei activity, and other relevant properties relevant for the Earth-Energy system. Beyond its role in benchmarking, AMBRS provides a platform for studying aerosol processes across scales and can be used to generate training data for AI/ML applications in support of a broader hierarchical aerosol modeling strategy.

54 ENVIRONMENTAL SCIENCES↗

Assignment and calibration of relative permeability by hydrostratigraphic units for multiphase flow analysis, case study: CO 2 -EOR operations at the Farnsworth Unit, Texas

Included in the most critical factors for geological CO 2 storage site screening, selection, and operation is effective simulations of multiphase flow and transport. Relative permeability is probably the greatest source of potential uncertainty in multiphase flow simulation, second only to intrinsic permeability heterogeneity. The specific relative permeability relationship assigned greatly impacts forecasts of CO 2 trapping mechanisms, phase behavior, and long-term plume movement. A primary goal of this study is to evaluate the impacts and implications of different methods of assigning relative permeability relationships for CO 2 -EOR model forecasts. Most simulation studies published in the literature base selection of relative permeability functions on the geologic formation or rock type alone. In this study, we initially implemented reservoir model grids with previously-identified hydrostratigraphic units based on porosity and permeability relationship of the Morrow ‘B’ Sandstone, then assigned relative permeability functions for those hydrostratigraphic units. Specific, constrained relative permeability relationships were created and assigned to each hydrostratigraphic unit using petrophysical data and Mercury Intrusion Capillary Pressure (MICP) measurements, from core samples of each hydrostratigraphic unit. Results of forward simulations with the newly-calibrated models will be compared to those of previous models as well as to simulation results for a range of different relative permeability relationships. The study site is the Farnsworth Unit (FWU) in the northeast Texas Panhandle, an active CO 2 -EOR operation. The target formation is the Morrow ‘B’ Sandstone, a clastic formation composed of medium to course sands.

58 GEOSCIENCES↗

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Predict-first experiments and modeling of perturbative cold pulses in the DIII-D tokamak

Cold pulses are presented in Ohmic DIII-D tokamak plasmas via injection of impurities with a laser blow-off system, revealing for the first time in this machine a quick increase in core electron temperature shortly after the edge cold-pulse injection at low collisionality. The experimental results are consistent with predict-first simulations of heat transport enabled by the Trapped Gyro-Landau-Fluid transport model. Measurements of electron density evolution during the cold-pulse propagation are enabled by a high time resolution density profile reflectometer. The density evolution reveals the quick propagation of a pulse from edge to core, which is a mechanism to transiently increase core temperature in low-collisionality plasmas. Local transport simulations with measured density evolution demonstrate that the core temperature response can indeed be explained by the stabilization of Trapped Electron Mode turbulence at low collisionality, hence providing confidence that local transport modeling is enough to explain cold-pulse propagation and associated phenomenology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗