Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder

Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the time at which L becomes resolved on the computational mesh. As a result, it is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Parallelizing Timed Petri Net simulations

The possibility of using parallel processing to accelerate the simulation of Timed Petri Nets (TPN's) was studied. It was recognized that complex system development tools often transform system descriptions into TPN's or TPN-like models, which are then simulated to obtain information about system behavior. Viewed this way, it was important that the parallelization of TPN's be as automatic as possible, to admit the possibility of the parallelization being embedded in the system design tool. Later years of the grant were devoted to examining the problem of joint performance and reliability analysis, to explore whether both types of analysis could be accomplished within a single framework. In this final report, the results of our studies are summarized. We believe that the problem of parallelizing TPN's automatically for MIMD architectures has been almost completely solved for a large and important class of problems. Our initial investigations into joint performance/reliability analysis are two-fold; it was shown that Monte Carlo simulation, with importance sampling, offers promise of joint analysis in the context of a single tool, and methods for the parallel simulation of general Continuous Time Markov Chains, a model framework within which joint performance/reliability models can be cast, were developed. However, very much more work is needed to determine the scope and generality of these approaches. The results obtained in our two studies, future directions for this type of work, and a list of publications are included.

Nicol, David M.↗
DOE OSTI2012

Wakefield Simulations for the Laser Acceleration Experiment at SLAC

Laser-driven acceleration in dielectric photonic band gap structures can provide gradients on the order of GeV/m. The small transverse dimension of the structure, on the order of the laser wavelength, presents interesting wakefield-related issues. Higher order modes can seriously degrade beam quality, and a detailed understanding is needed to mitigate such effects. On the other hand, wakefields also provide a direct way to probe the interaction of a relativistic bunch with the synchronous modes supported by the structure. Simulation studies have been carried out as part of the effort to understand the impact on beam dynamics, and to compare with data from beam experiments designed to characterize candidate structures. In this paper, we present simulation results of wakefields excited by a sub-wavelength bunch in optical photonic band gap structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluation of Portable Acceleration Solutions for LArTPC Simulation Using Wire-Cell Toolkit

The Liquid Argon Time Projection Chamber (LArTPC) technology plays an essential role in many current and future neutrino experiments. Accurate and fast simulation is critical to developing efficient analysis algorithms and precise physics model projections. The speed of simulation becomes more important as Deep Learning algorithms are getting more widely used in LArTPC analysis and their training requires a large simulated dataset. Heterogeneous computing is an efficient way to delegate computationally intensive tasks to specialized hardware. However, as the landscape of compute accelerators quickly evolves, it becomes increasingly difficult to manually adapt the code to the latest hardware or software environments. A solution which is portable to multiple hardware architectures without substantially compromising performance would thus be very beneficial, especially for long-term projects such as the LArTPC simulations. In search of a portable, scalable and maintainable software solution for LArTPC simulations, we have started to explore high-level portable programming frameworks that support several hardware backends. In this paper, we present our experience porting the LArTPC simulation code in the Wire-Cell Toolkit to NVIDIA GPUs, first with the CUDA programming model and then with a portable library called Kokkos. Preliminary performance results on NVIDIA V100 GPUs and multi-core CPUs are presented, followed by a discussion of the factors affiecting the performance and plans for future improvements.

Yu, Haiwang↗
DOE OSTI2023

TeraChem protocol buffers ( TCPB ): Accelerating QM and QM/MM simulations with a client–server model

The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. Here, we describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google’s protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs with electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ~250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI - )—a model chromophore in green fluorescent proteins—on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Simulator for an Accelerator-Driven Subcritical Fissile Solution System

LANL has developed a process to generate a progressive family of system models for a fissile solution system. This family includes a dynamic system simulation comprised of coupled nonlinear differential equations describing the time evolution of the system. Neutron kinetics, radiolytic gas generation and transport, and core thermal hydraulics are included in the DSS. Extensions to explicit operation of cooling loops and radiolytic gas handling are embedded in these systems as is a stability model. The DSS may then be converted to an implementation in Visual Studio to provide a design team the ability to rapidly estimate system performance impacts from a variety of design decisions. This provides a method to assist in optimization of the system design. Once design has been generated in some detail the C++ version of the system model may then be implemented in a LabVIEW user interface to evaluate operator controls and instrumentation and operator recognition and response to off-normal events. Taken as a set of system models the DSS, Visual Studio, and LabVIEW progression provides a comprehensive set of design support tools.

42 ENGINEERING↗
DOE OSTI2011

Dispersion And Energy Compensation in High-Gradient Linacs for Lepton Colliders

The shape of an rf pulse is distorted due to dispersion encountered in acceleration through traveling-wave linear accelerator structures. Simulations are made to ascertain the severity of this distortion in cavities designed to operate at various group velocities. The pulse suffers maximum degradation when propagated through accelerating cavities with a phase advance per cell in the vicinity of {pi}, where the group velocity reaches its minimum value. Several cavities are simulated to study the pulse distortion and compared with experiments performed on a high phase advance structure H60VG3, which has a phase advance of 5{pi}/6 per cell. A circuit model and a mode matching code are used to perform these simulations on accelerating structures consisting of multiple cells. Beam loading is taken into account and the implications on energy dispersion are also analyzed. Means of mitigating for this energy deviation are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

The time dependent propensity function for acceleration of spatial stochastic simulation of reaction–diffusion systems

The inhomogeneous stochastic simulation algorithm (ISSA) is a fundamental method for spatial stochastic simulation. However, when diffusion events occur more frequently than reaction events, simulating the diffusion events by ISSA is quite costly. To reduce this cost, we propose to use the time dependent propensity function in each step. In this way we can avoid simulating individual diffusion events, and use the time interval between two adjacent reaction events as the simulation stepsize. We demonstrate that the new algorithm can achieve orders of magnitude efficiency gains over widely-used exact algorithms, scales well with increasing grid resolution, and maintains a high level of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Simulation study of CO 2 laser-plasma interactions and self-modulated wakefield acceleration

3D numerical simulations of the interaction of a powerful CO 2 laser with hydrogen jets demonstrating the role of ionization on the characteristics of induced wakes are presented. Simulations using SPACE, a parallel relativistic Particle-in-Cell code, are performed in support of the plasma wakefield accelerator experiments being conducted at the Brookhaven National Laboratory (BNL) Accelerator Test Facility (ATF). A novelty of the SPACE code is its set of efficient atomic physics algorithms that compute ionization and recombination rates on the grid and transfer them to particles. The influence of ionization on the spectrum of pump laser has been studied for a range of gas densities.Simulations reproduce both Stokes and anti-Stokes shifts in the spectrum of pump laser, similar to those observed in experiments in the spectrum of probe. A good agreement has been achieved with the experiments on the effect of variation in gas density on Stokes/anti-Stokes intensity. In addition, self-injection and trapping of electrons into the self-modulated wakes has been observed and analyzed. Here, the experimentally validated code SPACE will be used for predictive simulations to guide future experiments at BNL ATF.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Hybrid Simulations of Interstellar Pickup Ion Acceleration at the Solar Wind Termination Shock

Hybrid (kinetic ion/fluid electron) simulations have been used to study self-consistently the injection and initial stages of acceleration of interstellar pickup ions at the solar wind termination shock. Results are presented from one-dimensional simulations of high Mach number oblique (40 deg-50 deg) shocks with a 10% population of interstellar pickup hydrogen. In these simulations, the pickup ions, the solar wind ions, the shock fields and the waves are all treated self-consistently. Pickup ions reflected by the shock excite large amplitude (DELTA.B / B approximately 0.3) upstream magnetosonic waves. These waves, in turn, scatter the pickup hydrogen, as expected in the diffusive shock acceleration process. The spectrum of excited waves broadens in time. We find that, for the parameters studied, the termination shock efficiently injects and accelerates the interstellar pickup hydrogen, and thus this work represents a microscopic verification of this anomalous cosmic ray injection mechanism. Comparison of injected fluxes from the simulations with Voyager fluxes extrapolated to the termination shock suggests that injection at the termination shock is probably greater than injection by interplanetary shocks. A study of accelerated pickup ion orbits shows that the energy gain comes predominantly from shock drift acceleration in the shock front, with the upstream waves aiding the acceleration by allowing multiple encounters with the shock.

Liewer, P. C.↗
DOE OSTI2014

GPUPEGAS: A NEW GPU-ACCELERATED HYDRODYNAMIC CODE FOR NUMERICAL SIMULATIONS OF INTERACTING GALAXIES

In this paper, a new scalable hydrodynamic code, GPUPEGAS (GPU-accelerated Performance Gas Astrophysical Simulation), for the simulation of interacting galaxies is proposed. The details of a parallel numerical method co-design are described. A speed-up of 55 times was obtained within a single GPU accelerator. The use of 60 GPU accelerators resulted in 96% parallel efficiency. A collisionless hydrodynamic approach has been used for modeling of stars and dark matter. The scalability of the GPUPEGAS code is shown.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Investigations on the multiple-sector hard-copper X-band accelerating structures

The development of advanced, high gradient accelerating structures is one of the leading activity of the particle accelerator community. In the technological research of new construction methods for these devices, high-power testing is a critical step for the verification of their viability. Recent experiments showed that accelerating cavities made out of hard copper, fabricated without high-temperature processes, can achieve better performance as compared with soft copper ones. Recently, we have built cavities using Tungsten Inert Gas welding and the high-power experiments confirmed that this joining process is a robust and low-cost alternative to brazing or diffusion bonding. This is a good solution for high-gradient operation, with a gradient of about 150 MV/m in X-band, at a breakdown rate of 10 -3 /pulse/meter using a shaped RF pulse with a 150 ns flat part. We continue the design, construction and high power tests of three-cells standing-wave X-band accelerating cavities fabricated with a split-open geometry, made of hard copper and vacuum-sealed by welding. Our aim is the fabrication of accelerating structures made out of hard copper alloys by using innovative cost-effective technologies. Moreover, our method of multiple-sector structures opens the way to new technological approaches and design methods, for example providing the possibility of placing dampers of the parasitic higher-order modes. Such structures with dampers allow for acceleration of multi-bunch beams, both in standing and in travelling-wave accelerating modality. Here, this paper describes the design of two cavities made of multiple sectors that were fabricated at the Italian Institute for Nuclear Physics in Frascati, Italy (INFN-LNF), assembled by means of clamping and then joined with the Tungsten Inert Gas welding in order to preserve the hardness of the metal. These are three-cell standing-wave X-band accelerating cavities, to be operated in the π-mode. We report the low-power RF tests of the first one, made of two halves, and the RF design of the second one, made of four quadrants. Both structures were built for high-gradient tests at SLAC National Accelerator Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Calculations of detuning with amplitude for the McMillan electron lens in the Fermilab Integrable Optics Test Accelerator (IOTA)

The McMillan electron lens will be one of the experimental implementations of nonlinear integrable lattices in the Fermilab Integrable Optics Test Accelerator (IOTA). We describe the physics of the nonlinear McMillan lens and calculate the tune dependence with amplitude from the experimental parameters. Some of the implications for the design of experiments in IOTA are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Beam physics research with the IOTA electron lens

The electron lens in the Fermilab Integrable Optics Test Accelerator (IOTA) will enable new research in nonlinear integrable optics, space-charge compensation, electron cooling, and the stability of intense beams. This research addresses scientific questions on high-brightness beams and operational challenges of high-power accelerators for nuclear and particle physics. We review the roles that electron lenses play in this field and the physical principles behind their applications. The design criteria and specifications for the IOTA storage ring and electron lens are then discussed. We conclude with a description of the components of the apparatus.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Numerical modeling of a proton spin-flipping system in the spin transparency mode at an integer spin resonance in JINR's Nuclotron

In this paper we propose a lattice insertion for the Nuclotron ring called a “spin navigator” that can adjust any direction of the proton polarization in the orbital plane using weak solenoids. The polarization control is realized in the spin transparency mode at the energy of 108 MeV, which corresponds to the integer spin resonance γ G = 2. The requirements on the navigator solenoid fields are specified considering the criteria for stability of the spin motion during any manipulation of the polarization direction in an experiment. Additionally, this paper presents the results of numerical modeling of the proton spin dynamics in the Nuclotron ring operated in the spin transparency mode. The verified spin navigator is aimed at an experimental study of a spin-flipping system using the Nuclotron ring. The results are relevant to the NICA (JINR), EIC (BNL) and COSY (FZJ) facilities where the spin transparency mode can be applied for polarization control.

47 OTHER INSTRUMENTATION↗
NASA NTRS1975

New two-stage accelerator for hypervelocity impact simulation

The present knowledge of the micrometeoroid environment and the simulation requirements for laboratory devices show that the simulated mass-velocity range must be extended to higher values in both mass and velocity. This is accomplished with the new two-stage accelerator, which consists of a light gas gun as the first stage and a coaxial accelerator with a compressor coil as the second stage. The operation of the new device is described together with a theoretical investigation of the timing of the first and the second stage, which leads to an appropriate adjustment of the two stages. The theoretical and experimental results agree well and show that 600-micron diameter glass beads have been accelerated to velocities near 20 km/sec.

Igenbergs, E. B.↗
DOE OSTI2022

Fast Quasi-Static Time-Series Simulation for Accurate PV Inverter Semiconductor Fatigue Analysis with a Long-Term Solar Profile

Power system simulations with long-term data typically have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stress cycle occurs over the fundamental grid frequency (50 or 60 Hz). This requires the time step of the fatigue simulation to be around 100 µs. This small time step requires long computation times to process yearly power production profiles. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The proposed simulation calculates the steady state of the semiconductor junction temperature by using a Fast Fourier Transform (FFT). The small thermal cycling during a switching period and even over the fundamental waveform is disregarded to further accelerate the simulation speed. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset. The error of the proposed simulation is 0.16% compared to the fatigue simulation results using the complete thermal stress profile. A PV inverter that responds to a Transactive Energy System (TES) is simulated to demonstrate the use of the proposed fatigue simulation. The proposed simulation has the potential to co-simulate with system level simulation tools that also adopt the QSTS concept.

Liu, Yunting↗