Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerator 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 937 records · Page 52

Coupling coarse mesh CFD with fine mesh CFD for modeling for Molten Salt Reactors in the Virtual Test Bed VTB

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Summary of Jefferson Lab LDRD on FFA@CEBAF Beam Dynamics Simulations

As Thomas Jefferson National Accelerator Facility (Jefferson Lab) looks toward the future, we are considering expanding our energy reach by using Fixed-Field Alternating Gradient (FFA) technology. Significant efforts have been made to design a hybrid accelerator which combines conventional recirculating electron LINAC design with permanent magnet-based FFA technology to increase the number of beam recirculations, and thus the energy. In an effort to further this progress, Jefferson Lab awarded a Laboratory Directed Research and Development (LDRD) grant to focus not on the design, but on detailed simulations of the designs created by the larger collaboration. This document will summarize the work performed during this LDRD, and direct the reader to other proceedings which describe elements of the work in greater detail.

Deitrick, K.↗

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗

Particle acceleration and plasma dynamics during magnetic reconnection in the magnetically dominated regime

Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two- and three-dimensional (3D) kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron–positron plasmas starting with a magnetically dominated, force-free current sheet (σ≡B{sup 2}/(4πn{sub e}m{sub e}c{sup 2})≫1). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra f∝(γ−1){sup −p} and approaches p = 1 for sufficiently large σ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to nonrelativistic regimes. In the 3D simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications for a variety of high-energy astrophysical problems.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-635: Quantitative Study On Beam Flatness Variation with Beam Energy Change

Purpose: Beam flatness check has been proposed for beam energy check for photon beams with flattering filters. In this work, beam flatness change with beam energy was investigated quantitatively using the Monte Carlo method and its significance was compared with depth dose curve change. Methods: Monte Carlo simulations for a linear accelerator with flattering filter were performed with different initial electron energies for photon beams of 6MV and 10MV. Dose calculations in a water phantom were then perform with the phase space files obtained from the simulations. The beam flatness was calculated based on the dose profile at 10 cm depth for all the beams with different initial electron energies. The percentage depth dose (PDD) curves were also analyzed. The dose at 10cm depth (D10) and the ratio of the dose at 10cm and 20cm depth (D10/D20) and their change with the beam energy were calculated and compared with the beam flatness variation. Results: It was found that the beam flatness variation with beam energy change was more significant than the change of D10 and the ratio between D10 and D20 for both 6MV and 10MV beams. Half MeV difference on the initial electron beam energy brought in at least 20% variation on the beam flatness but only half percent change on the ratio of D10 and D20. The change of D10 or D20 alone is even less significant. Conclusion: The beam energy impact on PDD is less significant than that on the beam flatness. If the PDD is used for checking the beam energy, uncertainties of the measurement could possibly disguise its change. Beam flatness changes more significantly with beam energy and therefore it can be used for monitoring the energy change for photon beams with flattering filters. However, other factors which may affect the beam flatness should be watched as well.

07 ISOTOPE AND RADIATION SOURCES↗

Simulation study of Type 2 counterstreaming electrons along auroral field lines

The production of counterstreaming electrons associated with parallel fields along auroral field lines is examined through the use of computer simulation. A 2 1/2-dimensional (two spatial and three velocity dimensions) electrostatic particle algorithm and auroral boundary conditions are used to set up a self-consistent V potential structure. The simulation produces signatures of counterstreaming electrons resembling those observed by the Dynamics Explorer 1 satellite. The main signatures are as follows: (1) the phase space contours of the electron distribution function are elongated along the V-parallel axis, and (2) the energy of electrons streaming in the upward direction is comparable to the energy of the accelerated electron beam. The simulation indicates that a portion of the accelerated electron beam is trapped by large amplitude electrostatic waves produced through the two-stream instability. Strong wave-particle interactions then thermalize the trapped electrons to produce suprathermal electrons streaming in the direction opposite to that of the accelerated electron beam. These results suggest a possible mechanism of producing counterstreaming electron fluxes through nonlinear processes of the two-stream instability.

Wagner, J. S.↗

Numerical Simulation of the Anomalous Transport of High-Energy Cosmic Rays in Galactic Superbubble

A continuous-time random-walk (CTRW) model to simulate the transport and acceleration of high-energy cosmic rays in galactic superbubbles has recently been put forward (Barghouty & Schnee 2102). The new model has been developed to simulate and highlight signatures of anomalous transport on particles' evolution and their spectra in a multi-shock, collective acceleration context. The superbubble is idealized as a heterogeneous region of particle sources and sinks bounded by a random surface. This work concentrates on the effects of the bubble's assumed astrophysical characteristics (cf. geometry and roughness) on the particles' spectra.

Barghouty, A. F.↗

Engineering-scale Modeling of High-Temperature Creep and Creep Crack Growth in Alloy 316H

This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.

316H↗

SU-E-T-461: Impact of the Radiation Sensitivity of Tumor and Geometric Design of the GRID Block On the Clinical Response of Grid Therapy

Purpose: The spatially fractionated radiation therapy (Grid therapy) is a novel technique that has been recently introduced for treatments of the patients with advanced bulky tumors. The purpose of this study is to investigate the effect of the radiation sensitivity of tumor and Grid block design on the clinical response of the Grid therapy. Methods: The Geant4 toolkit is used to simulate a Varian2100C linear accelerator with the block-based Grid. The dose profile has been utilized to calculate the Therapeutic Ratio (TR) and Equivalent Uniform Dose (EUD) for different types of tumors. These calculations were performed with respect to their radiation sensitivities using Linear Quadratic model (LQ). Clinical responses of few hundred patients from different publications that were treated with Grid therapy, have been used to validate correlation between tumor radiation sensitivity and clinical response of Grid therapy. Moreover, the influence of Grid design on clinical response of this radiation technique was investigated with using Monte Carlo simulation of a Grid block with different hole-diameters of 0.5 cm, 0.75 cm, 1.25 cm, and 1.5 cm. The Geometrical Sparing Factor (GSF) was considered to compare the results from different Grid designs. Results: The higher clinical response of Grid therapy for more radio-resistant tumors was confirmed in both theoretical assumptions and clinical outcomes. In addition, the results indicate that the Grid hole diameter have more effect on increasing the TR in comparison with the hole center-to-center distance effect. From this Monte Carlo study, a Grid design with the hole diameter between 1.25 cm and 1.5 cm with GSF value of 0.90 is recommended. Conclusion: Grid Therapy is more effective for more radio-resistant tumors.The optimum results of the Grid therapy depend on selection of the appropriate Grid design.This study is a first step toward clarifying some important aspects in a clinical response of spatially fractionated radiation therapy.

60 APPLIED LIFE SCIENCES↗

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. II. Connection with γ-Rays

Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LAT γ -ray bands as well as closely correlated optical and γ -ray flares. The swings in optical polarization angle are also accompanied by γ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast γ -ray flares or orphan γ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating the Performance of the NASA LaRC CMF Motion Base Safety Devices

This paper describes the initial measured performance results of the previously documented NASA Langley Research Center (LaRC) Cockpit Motion Facility (CMF) motion base hardware safety devices. These safety systems are required to prevent excessive accelerations that could injure personnel and damage simulator cockpits or the motion base structure. Excessive accelerations may be caused by erroneous commands or hardware failures driving an actuator to the end of its travel at high velocity, stepping a servo valve, or instantly reversing servo direction. Such commands may result from single order failures of electrical or hydraulic components within the control system itself, or from aggressive or improper cueing commands from the host simulation computer. The safety systems must mitigate these high acceleration events while minimizing the negative performance impacts. The system accomplishes this by controlling the rate of change of valve signals to limit excessive commanded accelerations. It also aids hydraulic cushion performance by limiting valve command authority as the actuator approaches its end of travel. The design takes advantage of inherent motion base hydraulic characteristics to implement all safety features using hardware only solutions.

Gupton, Lawrence E.↗

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗

The effects of bed rest on crew performance during simulated shuttle reentry. Volume 1: Study overview and physiological results

A centrifuge study was carried out to measure physiological stress and control task performance during simulated space shuttle orbiter reentry. Jet pilots were tested with, and without, anti-g-suit protection. The pilots were exposed to simulated space shuttle reentry acceleration profiles before, and after, ten days of complete bed rest, which produced physiological deconditioning similar to that resulting from prolonged exposure to orbital zero g. Pilot performance in selected control tasks was determined during simulated reentry, and before and after each simulation. Physiological stress during reentry was determined by monitoring heart rate, blood pressure, and respiration rate. Study results indicate: (1) heart rate increased during the simulated reentry when no g protection was given, and remained at or below pre-bed rest values when g-suits were used; (2) pilots preferred the use of g-suits to muscular contraction for control of vision tunneling and grayout during reentry; (3) prolonged bed rest did not alter blood pressure or respiration rate during reentry, but the peak reentry acceleration level did; and (4) pilot performance was not affected by prolonged bed rest or simulated reentry.

Chambers, A.↗

Neural Phase Simulation

The Neural Phase Simulation (NPS) is a package of codes for simulating microstructure evolution and accelerated molecular dynamics with deep neural-networks based surrogate models. NPS is designed to offer quantitatively accurate and computationally efficient simulation capabilities by leveraging modern machine-learning techniques. The primary intended use cases of NPS are training neural network surrogate models, though performing simulations on a single node is also supported. The NPS surrogate models can be trained from ground truth simulation methods, which are supposed to be accurate but expensive, such as molecular dynamics, phase field methods, kinetic Monte Carlo and discrete dislocation dynamics.

Zhou, Fei↗

THE MECHANISMS OF ELECTRON ACCELERATION DURING MULTIPLE X LINE MAGNETIC RECONNECTION WITH A GUIDE FIELD

The interactions between magnetic islands are considered to play an important role in electron acceleration during magnetic reconnection. In this paper, two-dimensional particle-in-cell simulations are performed to study electron acceleration during multiple X line reconnection with a guide field. Because the electrons remain almost magnetized, we can analyze the contributions of the parallel electric field, Fermi, and betatron mechanisms to electron acceleration during the evolution of magnetic reconnection through comparison with a guide-center theory. The results show that with the magnetic reconnection proceeding, two magnetic islands are formed in the simulation domain. Next, the electrons are accelerated by both the parallel electric field in the vicinity of the X lines and the Fermi mechanism due to the contraction of the two magnetic islands. Then, the two magnetic islands begin to merge into one, and, in such a process, the electrons can be accelerated by both the parallel electric field and betatron mechanisms. During the betatron acceleration, the electrons are locally accelerated in the regions where the magnetic field is piled up by the high-speed flow from the X line. At last, when the coalescence of the two islands into one big island finishes, the electrons can be further accelerated by the Fermi mechanism because of the contraction of the big island. With the increase of the guide field, the contributions of the Fermi and betatron mechanisms to electron acceleration become less and less important. When the guide field is sufficiently large, the contributions of the Fermi and betatron mechanisms are almost negligible.

79 ASTRONOMY AND ASTROPHYSICS↗

Longitudinal Beam-Shaping Simulation for Enhanced Transformer Ratio in Beam-Driven Accelerators

A beam-driven structure wakefield accelerator (SWFA) is currently under consideration for future multi-user Xray free-electron-laser (FEL) light sources. One of the challenges associated with SWFA is the efficient energy transfer from the drive to the main bunch that can be enhanced via a precise control over the drive-bunch current profile. This paper discusses the development and initial beam-dynamics simulation of the drivebunch longitudinal dynamics using the fast-tracking program TWICE. This code, written in PYTHON, allows for rapid prototyping of the beam shaping process while taking into account the collective effects using simplified models. This paper describes the algorithm implementation and its capabilities of including collective effects into both forward and backward tracking. Initial results of optimizing beam shape with a conventional linac to produce a drive beam with a shape optimized for SFWA are presented.

43 PARTICLE ACCELERATORS↗

Particle acceleration by propagating interplanetary shocks

The process of charged particle acceleration in interplanetary shocks has been simulated for typical parameters. Since space probe observations of charged particle fluxes are the principle means of inferring the action of an acceleration mechanism, the simulation was designed to follow particles backwards in time from a given observing point, usually 1 AU. This allowed a full diagnosis of which observed particles had interacted with an oncoming shock, where the interacting occurred, and by how much the energy had been changed by the shock interaction. Simple assumptions about the preshock energy spectrum allow the construction of a full temporal profile of the expected intensity, anisotropy, and energy spectrum. The simulations are apparently capable of reproducing the main features of the less than 10 MeV/nuc ion flux enhancements observed at interplanetary shocks using only a laminar interplanetary magnetic field, Archimidean spiral, and a spherical oblique shock with typical speed, strength, and shock normal-magnetic field angle.

Armstrong, T. P.↗

Experimental Setup for the Laboratory Investigation of Micrometeoroid Ablation Using a Dust Accelerator

A facility has been developed to simulate the ablation of micrometeoroids in laboratory conditions. An electrostatic dust accelerator is used to generate iron particles with velocities of 10-70 kilometers. The particles are then introduced into a chamber pressurized with a target gas, where the pressure is adjustable between 0.01 and 0.5 Torr, and the particle partially or completely ablates over a short distance. An array of biased electrodes above and below the ablation path is used to collect the generated ions/electrons with a spatial resolution of 2.6 centimeters along the ablating particles path, thus allowing the study of the spatiotemporal evolution of the process. For completely ablated particles, the total collected charge directly yields the ionization coefficient of a given dust material-target gas combination. The first results of this facility measured the ionization coefficient of iron atoms with N2, air, CO2, and He target gases for impact velocities greater than 20 kilometers per second, and are reported by Thomas et al. The ablation chamber is also equipped with four optical ports that allow for the detection of the light emitted by the ablating particle. A multichannel photomultiplier tube system is used to observe the ablation process with a spatial and temporal resolution of 0.64 centimeters and 90 nanoseconds. The preliminary results indicate that it is possible to calculate the velocity of the ablating particle from the optical observations, and in conjunction with the spatially resolved charge measurements allow for experimental validation of ablation models in future studies.

Thomas, Evan↗