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 1,297 records · Page 72

ADVANTG Results for a Proposed NDA System for Detecting Pin Diversion from SNF - Paper 120

This summary describes the results of applying the ADVANTG variance reduction parameter generator to accelerate continuous-energy Monte Carlo simulations of the {sup 252}Cf Interrogation with Prompt Neutron (CIPN) detection system proposed by Hu et al. for detecting the diversion of fuel elements from spent nuclear fuel (SNF) assemblies. In previous work, ADVANTG was applied to a suite of nine radiation detection, special nuclear material (SNM) movement detection, and safeguards test cases. That effort demonstrated that ADVANTG was able to provide large increases in tally figures of merit (FOMs) in nearly all of the cases without introducing bias into the Monte Carlo results. The objective of the work described here is to further evaluate ADVANTG, and also to identify gaps in capability with regard to nuclear safeguards applications, by applying the code to a challenging problem that is of a different nature from the ones previously studied. The CIPN nondestructive assay (NDA) system consists of three fission chambers embedded in a C-shaped block of polyethylene and a polyethylene-moderated {sup 252}Cf source that are placed on opposite sides of a subject SNF assembly. The idea of the detection system is to measure the difference between the background neutron count rate, due to spontaneous fission in the fuel material, and the count rate observed once the source is introduced, which is dominated by multiplication within the assembly. Based on the difference between these count rates, the fissile content of the assembly can be inferred, and this information can be used to determine if some fraction of the original fuel pins have been replaced with depleted uranium substitutes. The authors report that the system can detect the replacement of at least eight pins (3% of total mass). The ADVANTG code generates space- and energy-dependent weight-window bounds and biased source distributions based on approximate transport solutions that are efficiently generated by the Denovo 3-D discrete ordinates package. ADVANTG implements the Consistent Adjoint Driven Importance Sampling (CADIS) method for accelerating individual tallies and the Forward-Weighted CADIS (FW-CADIS) method for obtaining relatively uniform uncertainties across tallies over multiple regions and/or energy bins. The variance reduction parameters are generated in a format directly usable by standard versions of MCNP5 and MCNPX. In this work, the ADVANTG variance reduction parameter generator was applied to accelerate continuous energy MCNPX simulations of the {sup 252}Cf Interrogation with Prompt Neutron detection system for detecting the diversion of fuel elements from SNF assemblies. Variance reduction parameters generated with the existing version of ADVANTG actually slowed down tally convergence in all cases. The cause of the poor performance was found to be the lack of an induced fission source treatment in the fixed-source discrete ordinates calculations employed by ADVANTG. Because fission was treated as capture and not as a source, the importance maps generated by ADVANTG significantly underestimated the importance of thermal neutrons in and around the fuel elements. In order to obtain improved results, we developed an experimental implementation of a fission source outer iteration for the adjoint calculations. We used this prototype capability to generate new variance reduction parameters and found that the mean fission rates agreed well with those predicted by the conventional MCNPX simulations, as was expected. With these new parameters, speedup factors of 59 to 80 were obtained for the case with the background source, and factors of 17 to 24 were obtained for the cases with the interrogating source. These results demonstrate that a fission source treatment is needed to improve the effectiveness of ADVANTG in problems where the effects of induced fission are important. This includes many types of NDA problems that are relevant to safeguards applications. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Unified differentiable digital twin for the IOTA/FAST facility

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or produce insights about accelerator state. One notable example of additional information are gradients of physical observables with respect to design parameters produced by differentiable simulations. The IOTA/FAST facility has recently begun a program to implement and experimentally validate a unified start-to-end differentiable digital twin to serve as a virtual accelerator test stand, allowing for rapid prototyping of new software and experiments with minimal beam time costs. In this contribution we will discuss our plans and progress. Specifically, we will cover the selection and benchmarking of both physics and ML codes, the development of generic interfaces between device models and surrogate or physics-based sections, and the export of the parameters through either a deterministic event loop or a fully asynchronous EPICS soft input/output controller. We will also discuss challenges in model calibration and uncertainty quantification, as well as future plans to support larger proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗

Man-systems evaluation of moving base vehicle simulation motion cues

A motion cue investigation program is reported that deals with human factor aspects of high fidelity vehicle simulation. General data on non-visual motion thresholds and specific threshold values are established for use as washout parameters in vehicle simulation. A general purpose similator is used to test the contradictory cue hypothesis that acceleration sensitivity is reduced during a vehicle control task involving visual feedback. The simulator provides varying acceleration levels. The method of forced choice is based on the theory of signal detect ability.

Kirkpatrick, M.↗

Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core

In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.

97 MATHEMATICS AND COMPUTING↗

Conceptual designs of two petawatt-class pulsed-power accelerators for high-energy-density-physics experiments

Here, we have developed conceptual designs of two petawatt-class pulsed-power accelerators: Z 300 and Z 800. The designs are based on an accelerator architecture that is founded on two concepts: single-stage electrical-pulse compression and impedance matching [Phys. Rev. ST Accel. Beams 10, 030401 (2007)]. The prime power source of each machine consists of 90 linear-transformer-driver (LTD) modules. Each module comprises LTD cavities connected electrically in series, each of which is powered by 5-GW LTD bricks connected electrically in parallel. (A brick comprises a single switch and two capacitors in series.) Six water-insulated radial-transmission-line impedance transformers transport the power generated by the modules to a six-level vacuum-insulator stack. The stack serves as the accelerator’s water-vacuum interface. The stack is connected to six conical outer magnetically insulated vacuum transmission lines (MITLs), which are joined in parallel at a 10-cm radius by a triple-post-hole vacuum convolute. The convolute sums the electrical currents at the outputs of the six outer MITLs, and delivers the combined current to a single short inner MITL. The inner MITL transmits the combined current to the accelerator’s physics-package load. Z 300 is 35 m in diameter and stores 48 MJ of electrical energy in its LTD capacitors. The accelerator generates 320 TW of electrical power at the output of the LTD system, and delivers 48 MA in 154 ns to a magnetized-liner inertial-fusion (MagLIF) target [Phys. Plasmas 17, 056303 (2010)]. The peak electrical power at the MagLIF target is 870 TW, which is the highest power throughout the accelerator. Power amplification is accomplished by the centrally located vacuum section, which serves as an intermediate inductive-energy-storage device. The principal goal of Z 300 is to achieve thermonuclear ignition; i.e., a fusion yield that exceeds the energy transmitted by the accelerator to the liner. 2D magnetohydrodynamic (MHD) simulations suggest Z 300 will deliver 4.3 MJ to the liner, and achieve a yield on the order of 18 MJ. Z 800 is 52 m in diameter and stores 130 MJ. This accelerator generates 890 TW at the output of its LTD system, and delivers 65 MA in 113 ns to a MagLIF target. The peak electrical power at the MagLIF liner is 2500 TW. The principal goal of Z 800 is to achieve high-yield thermonuclear fusion; i.e., a yield that exceeds the energy initially stored by the accelerator’s capacitors. 2D MHD simulations suggest Z 800 will deliver 8.0 MJ to the liner, and achieve a yield on the order of 440 MJ. Z 300 and Z 800, or variations of these accelerators, will allow the international high-energy-density-physics community to conduct advanced inertial-confinement-fusion, radiation-physics, material-physics, and laboratory-astrophysics experiments over heretofore-inaccessible parameter regimes.

42 ENGINEERING↗

International Space Station Microgravity Analytical Model Correlation And Update

The acceleration environment aboard the completed International Space Station (ISS) is a key resource for scientific and technological endeavors. Hardware verification activIties and early measurements indicate that the ISS is well on the way of meeting these "Assembly Complete" "microgravity" provisions, however, the simulation models that compute these accelerations have, to date, lacked the high degree of empirical validation typical of standard aerospace industry practices. Assembly stage, on-orbit measurements are used to address this shortcoming and to develop higher confidence in the simulation models. The Phase I correlation results show the analyses to be consistently conservative, producing higher than measured levels. The 25 to 30% greater quasi-steady computations are deemed acceptable for verification. Updates are made to localized structural dynamic and vibroacoustic parameters that reduce responses in selected one-third octave bands by almost 50%. These models are then used for the Assembly Complete verification analysis which concludes that the ISS vehicle meets the ISS microgravity requirements with minor reservations. Two of the sixteen rack are marginally non-compliant in the quasi-steady regime, and operational constraints are needed on the U. S. Lab and ESA APM vacuum resource vents, and the Russian Resistive Exercise Device in the structural dynamic regime.

DelBasso, Steve↗

gRASPA

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

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

Innovations in optimization and control of accelerators using methods of differential geometry and genetic algorithms (Final Report)

Online tuning of particle accelerators is necessary in order to achieve optimal machine performance. However, it is also a major challenge due to the large parameter space which must be searched and the fact that many of the desir- able objectives compete with one another, and so will not reach their optimal values simultaneously. In order to mitigate these issues, we have explored using dimension-reduction techniques to reduce the size of the parameter space which must be searched and multi-objective genetic algorithms to obtain the sets of tuning parameters which provide pareto-optimal values for the objectives. These methods have enabled us to obtain improved values of the vertical emittance at the Cornell Electron Storage Ring (CESR), and to do so with greater control of orbit errors. Algorithmic tuning is a multidisciplinary endeavour, requiring expertise in beam dynamics, diagnostics, control systems and computer science, and thus a key practical problem is to formulate a common language in which experts with different specialties can communicate. We developed a solution to this problem in the form of a generic accelerator software interface that allows for rapid prototyping of optimization and control algorithms. Our interface is built on the Experimental Physics and Industrial Control Systems (EPICS) and makes possible testing control code in simulation before deployment on real accelerators, as well as deployment of third-party optimization code.

43 PARTICLE ACCELERATORS↗

Particle acceleration in modified oblique non-relativistic shocks

Previous applications of the Monte Carlo technique at the quasiparallel earth bow shock has motivated the extension of this technique to oblique shock geometries typical of those found in most astrophysical shock environments. In addition, such a generalization will permit the thorough examination of theoretical predictions of rapid acceleration times at quasi-perpendicular shocks. Therefore, we have embarked on the modification of our existing Monte Carlo code and in this paper outline the major technical aspects involved in developing a simulation of cosmic-ray acceleration at modified oblique nonrelativistic shocks.

Baring, Matthew G.↗

Jupiter entry simulation with the ANAA shock tube

An annular arc accelerator (ANAA) shock tube has been built which produces shock velocities and pressures that simulate entry into the atmosphere of Jupiter. The ANAA driver deposits the energy of an arc discharge into a flowing gas, which then expands and cools without any delay for the opening of a diaphragm. A flow transducer, trigger system, and spark gap switches have been developed to coordinate the flow from a high-pressure helium driver with the discharge from a 300-kJ capacitor bank. Shock velocities up to 47 km/sec have been produced in 1.0 torr of hydrogen with the ANAA shock tube, compared with 35 km/sec velocity produced in a conical arc driver with three times the available energy. Attenuation with the ANAA shock tube is comparable to that of a conical arc driver shock tube, and initial spectroscopic measurements indicate that an impurity-free test slug is formed behind the shock wave.

Leibowitz, L. P.↗

Generalized 3D beam dynamics model for industrial traveling wave linacs design and simulations

Beam dynamics simulations in traveling wave (TW) accelerating structures with significant beam loading is a challenging problem. Some codes are capable of calculating the TW electromagnetic fields, and some can track particles through such fields, but most cannot treat both self-consistently. A few commercial codes can model the physics correctly at the expense of many processor-hour computations to obtain converged self-consistent solutions. However, simple, accurate equations of motion for intensive beam dynamics in TW accelerating structure analysis have been previously obtained by Masunov and subsequently implemented in the Hellweg code, which was developed at the Moscow Engineering-Physics Institute. Hellweg is based on equations that allow fast simulations of beam dynamics while taking into account such effects as beam loading, space charge, and external magnetic fields. In this paper, we describe in detail some recent improvements to the Hellweg physics kernel. We have generalized the Masunov results into a 3D set of equations of motion, which include all spatial components of the radiofrequency (RF) and external magnetic fields. We have also improved the Lapostolle space charge model to the general 3D ellipsoid form for any dimensions’ ratio, consideration for the particles outside the beam core, and the fields from the neighboring bunches that can exist in the real machine. These modifications allow approaching the Hellweg accuracy to the self-consistent commercial codes while keeping the simulation time short, which is essential during the linac design and optimization stage. The implementation of these new capabilities in Hellweg is carefully benchmarked against other codes and analytical calculations. The code is freely available with an open source license.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Aperture-based antihydrogen gravity experiment: Parallel plate geometry

An analytical model and a Monte Carlo simulation are presented of an experiment that could be used to determine the direction of the acceleration of antihydrogen due to gravity. The experiment would rely on methods developed by existing antihydrogen research collaborations. The configuration consists of two circular, parallel plates that have an axis of symmetry directed away from the center of the earth. The plates are separated by a small vertical distance, and include one or more pairs of circular barriers that protrude from the upper and lower plates, thereby forming an aperture between the plates. Antihydrogen annihilations that occur just beyond each barrier, within a “shadow” region, are asymmetric on the upper plate relative to the lower plate. The probability for such annihilations is determined for a point, line and spheroidal source of antihydrogen. The production of 100,000 antiatoms is predicted to be necessary for the aperture-based experiment to indicate the direction of free fall acceleration of antimatter, provided that antihydrogen is produced within a sufficiently small antiproton plasma at a temperature of 4 K.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the possibility of footprint compression with one lens in nonlinear accelerator lattice

Electromagnetic interaction of colliding beams along with other nonlinear fields often limits the beams' lifetimes and luminosities. Nonlinearities result in the spread of betatron frequencies (footprint) and, thus, may enhance dynamic diffusion of particles due to high order resonances. One of the possible ways to eliminate nonlinearities and overcome the corresponding difficulties is compensation of nonlinear forces, but, in practice, it is hardly possible to obtain exact linearity of the system. The compensation with a single nonlinear lens cannot cope with distributed nonlinearities, nonlinearities due to parasitic crossings, etc. Here, we present a method to compute parameters of nonlinear element (lens) that eliminates both the footprint and resonance strength without achieving full compensation.

43 PARTICLE ACCELERATORS↗

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗

Accelerator Test of an Imaging Calorimeter

The Imaging Calorimeter for ACCESS (ICA) utilizes a thin sampling calorimeter concept for direct measurements of high-energy cosmic rays. The ICA design uses arrays of small scintillating fibers to measure the energy and trajectory of the produced cascades. A test instrument has been developed to study the performance of this concept at accelerator energies and for comparison with simulations. Two test exposures have been completed using a CERN test beam. Some results from the accelerator tests are presented.

Christl, Mark J.↗

Study of changes in properties of solar sail materials from radiation exposure

Techniques for monitoring changes in preparation of solar sail materials resulting from space radiation simulation, stressing (e.g., thermal, mechanical) and exposure to terrestrial environments are developed. The properties of interest are: metallic coating deterioration, polymeric film deterioration, interfacial debonding and possible metallic coating diffusion into the polymeric film. Four accelerated tests were devised to simulate the possible degradation processes mentioned above. These four tests are: a thermal shock test to simulate the wide variation of temperature expected in space (260 C to -100 C), a cyclic temperature test to stimulate the 6 minute temperature cycle anticipated in space, a mechanical vibration test to simulate mechanical bonding, folding and handling, and a humidity test to simulate terrestrial environment effects. The techniques for monitoring property changes are: visual and microscopic examination, ellipsometry, surface potential difference (SPD), photoelectron emission (PEE), and water contact angles.

Smith, T.↗

Pronghorn: A Porous Media Thermal-Hydraulics Core Simulator and its Validation with the SANA Experiments

Pebble bed High Temperature Reactors (HTRs) are characterized by many advantageous design features, such as excellent passive heat removal in accidents, large margins to fuel failure, and online refueling potential. However, a significant challenge in the core modeling of pebble bed reactors is the complex fuel-coolant structure. This paper presents a new porous media simulation code, Pronghorn, that aims to alleviate modeling challenges for pebble bed reactors by providing a fast-running, medium- fidelity core simulator. Pronghorn is intended to accelerate the design and analysis cycle for pebble bed and prismatic HTRs by permitting fast scoping studies and providing boundary conditions for systems-level analysis. Pronghorn is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) using modern software practices and a thorough testing framework. This paper describes the physical models used in Pronghorn and demonstrates Pronghorn’s capability for modeling gas-cooled pebble bed HTRs by presenting simulation results obtained for the German SANA pebble bed decay heat experiments. Within the limitations of the porous media approximation and existing available closure relationships, Pronghorn predicts the SANA experimental pebble temperatures well, expanding the code’s validation base. A brief code-to-code comparison shows a similar level of accuracy with other porous media simulation tools; this enhances Pronghorn’s advantages of an arbitrary equation of state, unstructured mesh capabilities, compressible flow models, high coupling potential to MOOSE fuels performance and systems-level thermal-hydraulics codes, and modern software design.

42 ENGINEERING↗

GEANT4 simulations for beam emittance in a linear collider based on plasma wakefield acceleration

Alternative acceleration technologies are currently under development for cost-effective, robust, compact, and efficient solutions. One such technology is plasma wakefield acceleration, driven by either a charged particle or laser beam. However, the potential issues must be studied in detail. In this paper, the emittance evolution of a witness beam through elastic scattering from gaseous media and under transverse focusing wakefields is studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗