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 847 records · Page 47

Experimental study of transient dynamics of a flexible rotor

The results of an experimental program to investigate the transient response of a flexible rotor are presented. The program consisted of a series of tests conducted on a rig with a rotor designed to operate above its first bending critical speed. The purpose of the tests was to obtain experimental data on the transient behavior of a flexible rotor under conditions simulating those which might occur in a jet aircraft engine. The scope of the program included tests to measure the response of both balanced and unbalanced rotors during steady-state operation, acceleration, deceleration, and simulated blade loss.

Hibner, D. H.↗

Differential Adaptive Stress Testing of Airborne Collision Avoidance Systems

The next-generation Airborne Collision Avoidance System (ACAS X) is currently being developed and tested to replace the Traffic Alert and Collision Avoidance System (TCAS) as the next international standard for collision avoidance. To validate the safety of the system, stress testing in simulation is one of several approaches for analyzing near mid-air collisions (NMACs). Understanding how NMACs can occur is important for characterizing risk and informingdevelopment of the system. Recently, adaptive stress testing (AST) has been proposed as a way to find the most likely path to a failure event. The simulation-based approach accelerates search by formulating stress testing as a sequential decision process then optimizing it using reinforcement learning. The approach has been successfully applied to stress test a prototype of ACAS Xin various simulated aircraft encounters. In some applications, we are not as interestedin the system's absolute performance as its performance relative to another system. Such situations arise, for example, during regression testing or when deciding whether a new system should replace an existing system. In our collision avoidance application, we are interested in finding cases where ACAS X fails but TCAS succeeds in resolving a conflict. Existing approaches do not provide an efficient means to perform this type of analysis. This paper extends the AST approach to differential analysis by searching two simulators simultaneously and maximizing the difference between their outcomes. We call this approach differential adaptive stress testing (DAST). We apply DAST to compare a prototype of ACAS X against TCAS and show examples of encounters found by the algorithm.

Lee, Ritchie↗

Measuring Saturn's Electron Beam Energy Spectrum using Webb's Wedges

It is very difficult to measure the voltage of the load on the Saturn accelerator. Time-resolved measurements such as vacuum voltmeters and V-dot monitors are impractical at best and completely change the pulsed power behavior at the load at worst. We would like to know the load voltage of the machine so that we could correctly model the radiation transport and tune our x-ray unfold methodology and circuit simulations of the accelerator. Step wedges have been used for decades as a tool to measure the end - point energies of high energy particle beams. Typically, the technique is used for multi-megavolt accelerators, but we have adapted it to Saturn's modest <2 MV end-point energy and modified the standard bremsstrahlung x-ray source to extract the electron beam without changing the physics of the load region. We found clear evidence of high energy electrons >2 MV. We also attempted to unfold an electron energy spectrum using a machine learning algorithm and while these results come with large uncertainties, they qualitatively agree with PIC simulation results.

43 PARTICLE ACCELERATORS↗

Effective post-acceleration of ion bunches in foils irradiated by ultra-intense laser pulses

Two-step laser acceleration of protons with two foils and two laser pulses is modelled and optimized. It is shown that a nearly mono-energetic distribution of proton bunches can be realized by a suitable parameter choice. Two-step acceleration schemes make it possible to obtain both higher efficiency and energy as compared to the acceleration with only one laser pulse of an energy equal to the sum of the energy of the two pulses. With the aid of our analytical model, the optimal distance between the two targets, the delay between the two laser pulses, and the parameters of the laser pulses are determined. Estimates and results of the modelling are proven with 2D PIC simulations of the acceleration of proton bunches moving through the second target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A 6 degree-of-freedom Lorentz force vibration isolator with nonlinear controller

This program demonstrated the technical feasibility of constructing large-stroke magnetic suspensions that can meet the active vibration isolation requirements of Space Station. These requirements include: (1) strokes over 1 cm in all directions, (2) actuator bandwidths over 100 Hz, (3) isolator roll-off frequencies below 10(exp -2) Hertz, and (4) force capability over 1 Newton in all axes. The 100 Hz actuator bandwidth allows the suspension to reject any direct force disturbances that act on the microgravity experiment, for example forces created by cable connections. The low isolator roll-off frequency and large stroke allow the magnetic suspension to isolate the microgravity experiment from Space Station vibrations above the roll-off frequency. The capability to meet these requirements was demonstrated by designing, constructing and testing a six-degree-of-freedom, prototype magnetic suspension system that featured high-performance, Lorentz-force actuators and full multi-input, multi-output control. This prototype suspension is designed to isolate large orbiter locker experiments under typical spacecraft constraints of size, weight, and power. Suspension in the full six-degrees-of-freedom was successfully demonstrated in this program while using a gravity-force unload mechanism to simulate a space environment. The prototype isolator is capable of space-based isolation service with relatively minor modification. The use of advanced, nonlinear control algorithms were investigated on a specially designed single-degree-of-freedom testbed. This low acceleration test facility simulates the Space Station vibration isolation problem in a single horizontal axis with low-friction, air-slide support. This allowed testing at the desired microgravity levels, without the gravity bias effects that are seen in a full six-degrees-of-freedom suspension. Precision components were used to reduce residual accelerations to microgravity levels so that the effects of sensor, actuator, and electronic noise could be evaluated. During the Phase 2 program, this testbed was used to demonstrate the advantages of nonlinear control algorithms to provide increased vibration isolation performance.

Fenn, Ralph↗

PARTICLE ACCELERATION AND WAVE EXCITATION IN QUASI-PARALLEL HIGH-MACH-NUMBER COLLISIONLESS SHOCKS: PARTICLE-IN-CELL SIMULATION

We herein investigate shock formation and particle acceleration processes for both protons and electrons in a quasi-parallel high-Mach-number collisionless shock through a long-term, large-scale, particle-in-cell simulation. We show that both protons and electrons are accelerated in the shock and that these accelerated particles generate large-amplitude Alfvénic waves in the upstream region of the shock. After the upstream waves have grown sufficiently, the local structure of the collisionless shock becomes substantially similar to that of a quasi-perpendicular shock due to the large transverse magnetic field of the waves. A fraction of protons are accelerated in the shock with a power-law-like energy distribution. The rate of proton injection to the acceleration process is approximately constant, and in the injection process, the phase-trapping mechanism for the protons by the upstream waves can play an important role. The dominant acceleration process is a Fermi-like process through repeated shock crossings of the protons. This process is a “fast” process in the sense that the time required for most of the accelerated protons to complete one cycle of the acceleration process is much shorter than the diffusion time. A fraction of the electrons are also accelerated by the same mechanism, and have a power-law-like energy distribution. However, the injection does not enter a steady state during the simulation, which may be related to the intermittent activity of the upstream waves. Upstream of the shock, a fraction of the electrons are pre-accelerated before reaching the shock, which may contribute to steady electron injection at a later time.

79 ASTRONOMY AND ASTROPHYSICS↗

Surface-Height Determination of Crevassed Glaciers-Mathematical Principles of an Autoadaptive Density-Dimension Algorithm and Validation Using ICESat-2 Simulator (SIMPL) Data

Glacial acceleration is a main source of uncertainty in sea-level-change assessment. Measurement of ice-surface heights with a spatial and temporal resolution that not only allows elevation-change calculation, but also captures ice-surface morphology and its changes is required to aid in investigations of the geophysical processes associated with glacial acceleration.The Advanced Topographic Laser Altimeter System aboard NASAs future ICESat-2 Mission (launch 2017) will implement multibeam micropulse photon-counting lidar altimetry aimed at measuring ice-surface heights at 0.7-m along-track spacing. The instrument is designed to resolve spatial and temporal variability of rapidly changing glaciers and ice sheets and the Arctic sea ice. The new technology requires the development of a new mathematical algorithm for the retrieval of height information.We introduce the density-dimension algorithm (DDA) that utilizes the radial basis function to calculate a weighted density as a form of data aggregation in the photon cloud and considers density an additional dimension as an aid in auto-adaptive threshold determination. The auto-adaptive capability of the algorithm is necessary to separate returns from noise and signal photons under changing environmental conditions. The algorithm is evaluated using data collected with an ICESat-2 simulator instrument, the Slope Imaging Multi-polarization Photon-counting Lidar, over the heavily crevassed Giesecke Braer in Northwestern Greenland in summer 2015. Results demonstrate that ICESat-2 may be expected to provide ice-surface height measurements over crevassed glaciers and other complex ice surfaces. The DDA is generally applicable for the analysis of airborne and spaceborne micropulse photon-counting lidar data over complex and simple surfaces.

Herzfeld, Ute C.↗

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗

Experimental concept and design of DarkLight, a search for a heavy photon

This talk gives an overview of the DarkLight experimental concept: a search for a heavy photon A′ in the 10-90 MeV/c{sup 2} mass range. After briefly describing the theoretical motivation, the talk focuses on the experimental concept and design. Topics include operation using a half-megawatt, 100 MeV electron beam at the Jefferson Lab FEL, detector design and performance, and expected backgrounds estimated from beam tests and Monte Carlo simulations.

43 PARTICLE ACCELERATORS↗

ON THE POLAR CAP CASCADE PAIR MULTIPLICITY OF YOUNG PULSARS

We study the efficiency of pair production in polar caps of young pulsars under a variety of conditions to estimate the maximum possible multiplicity of pair plasma in pulsar magnetospheres. We develop a semi-analytic model for calculation of cascade multiplicity which allows efficient exploration of the parameter space and corroborate it with direct numerical simulations. Pair creation processes are considered separately from particle acceleration in order to assess different factors affecting cascade efficiency, with acceleration of primary particles described by recent self-consistent non-stationary model of pair cascades. We argue that the most efficient cascades operate in the curvature radiation/synchrotron regime, the maximum multiplicity of pair plasma in pulsar magnetospheres is ∼few × 10{sup 5}. The multiplicity of pair plasma in magnetospheres of young energetic pulsars weakly depends on the strength of the magnetic field and the radius of curvature of magnetic field lines and has a stronger dependence on pulsar inclination angle. This result questions assumptions about very high pair plasma multiplicity in theories of pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulation of the collapse and dissipation of Langmuir wave packets

Particle-in-cell (PIC) simulations and Zakharov's partial differential equations (PDEs) are used to investigate the collapse of isolated Langmuir wave packets in two dimensions. Collapse thresholds are determined numerically, and the roles of enhanced Langmuir wave damping and nonlinearities not included in the standard Zakharov equations are discussed. The Langmuir wave and ion dynamics in PIC simulations are compared with the predictions of PDE simulations incorporating enhanced Langmuir damping. Electron heating and coherent acceleration in the PIC simulations are discussed and compared with predictions of the transit-time theory.

Newman, D. L.↗

Program For Simulation Of Trajectories And Events

Universal Simulation Executive (USE) program accelerates and eases generation of application programs for numerical simulation of continuous trajectories interrupted by or containing discrete events. Developed for simulation of multiple spacecraft trajectories with events as one spacecraft crossing the equator, two spacecraft meeting or parting, or firing rocket engine. USE also simulates operation of chemical batch processing factory. Written in Ada.

Gottlieb, Robert G.↗

Study the sensitivity of dose calculation in prism treatment planning system using Monte Carlo simulation of 6 MeV electron beam

Prism is a non-commercial Radiotherapy Treatment Planning System (RTPS) develop by Ira J. Kalet from Washington University. Inhomogeneity factor is included in Prism TPS dose calculation. The aim of this study is to investigate the sensitivity of dose calculation on Prism using Monte Carlo simulation. Phase space source from head linear accelerator (LINAC) for Monte Carlo simulation is implemented. To achieve this aim, Prism dose calculation is compared with EGSnrc Monte Carlo simulation. Percentage depth dose (PDD) and R50 from both calculations are observed. BEAMnrc is simulated electron transport in LINAC head and produced phase space file. This file is used as DOSXYZnrc input to simulated electron transport in phantom. This study is started with commissioning process in water phantom. Commissioning process is adjusted Monte Carlo simulation with Prism RTPS. Commissioning result is used for study of inhomogeneity phantom. Physical parameters of inhomogeneity phantom that varied in this study are: density, location and thickness of tissue. Commissioning result is shown that optimum energy of Monte Carlo simulation for 6 MeV electron beam is 6.8 MeV. This commissioning is used R50 and PDD with Practical length (R{sub p}) as references. From inhomogeneity study, the average deviation for all case on interest region is below 5 %. Based on ICRU recommendations, Prism has good ability to calculate the radiation dose in inhomogeneity tissue.

61 RADIATION PROTECTION AND DOSIMETRY↗

Full-f gyrokinetic simulation over a confinement time

A long time ion temperature gradient driven turbulence simulation over a confinement time is performed using the full-f gyrokinetic Eulerian code GT5D. The convergence of steady temperature and rotation profiles is examined, and it is shown that the profile relaxation can be significantly accelerated when the simulation is initialized with linearly unstable temperature profiles. In the steady state, the temperature profile and the ion heat diffusivity are self-consistently determined by the power balance condition, while the intrinsic rotation profile is sustained by complicated momentum transport processes without momentum input. The steady turbulent momentum transport is characterized by bursty non-diffusive fluxes, and the resulting turbulent residual stress is consistent with the profile shear stress theory [Y. Camenen et al., “Consequences of profile shearing on toroidal momentum transport,” Nucl. Fusion 51, 073039 (2011)] in which the residual stress depends not only on the profile shear and the radial electric field shear but also on the radial electric field itself. Based on the toroidal angular momentum conservation, it is found that in the steady null momentum transport state, the turbulent residual stress is cancelled by the neoclassical counterpart, which is greatly enhanced in the presence of turbulent fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Impact of Anomalous Resistivity in Vacuum Contaminant Plasmas on the Electrothermal Instability

This manuscript presents an assessment of the electrothermal instability (ETI) in the presence of anomalous resistivity (AR) in vacuum contaminant plasmas (VCP) when applied to a magnetized liner inertial fusion (MagLIF)-like load.Pulsed-power driven dielectrically coated metallic liners, like in MagLIF, experience the current-driven electrothermal instability which occurs when a material’s resistivity changes with temperature and is subject to ohmic heating. Large scale pulsed-power facilities that use magnetically insulated transmission lines (MITL) have been shown to generate low-density plasma which enters the target chamber and coalesces around the load. The low-density high-temperature vacuum contaminant plasmas (VCP) can parasitically divert current from the load through causing a short in the anode-cathode gap inside the target chamber. Resistive magnetohydrodynamic (MHD) simulations of these VCP experience unphysical runaway ohmic heating due to under predicting the resistivity by using a purely collisional resistivity model.AR provides a physics-based way to address this runaway heating through increasing the resistivity in a proportional way with the drift speed. In this work, 1D simulations probe the effect that AR in VCP has on the magnetic diffusion rate, and 2D simulations show how this effect manifests in the nonlinear striation form of the ETI for a MagLIF-like load. Beryllium and aluminum dielectrically coated liners are used for the 1D and 2D simulations in this work. The1D simulations show that a VCP causes a delay in the current delivery to the load by upwards of 8 ns at 60 ns into atypical current pulse for an 1 × 10 7 A scale pulsed-power accelerator. The 2D simulations show the delay observed in the 1D simulations is reduced substantially (4 ns) by 2D turbulence that disrupts the VCP layer early in time. Here, the 2Dsimulated ETI growth varies across AR models, more so for the beryllium liner than the aluminum liner because the beryllium liner shows an enhanced rate of penetration for the magnetic diffusion wave in comparison to the aluminum liner. The 2D simulations show the bulk dielectric thickness varies across AR models with the Davidson AR model beingthe largest and the Buneman AR model being the smallest, and in connection with the thickness the Rayleigh-Taylor bubble-spike distances varies correspondingly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy Production Demonstrator and Material Testing Station optimization for Megawatt proton beams

A simulation study of the Energy Production Demonstrator (EPD) concept – a solid heavy metal target irradiated by GeV-range intense proton beams and producing more energy than consuming – is carried out. Neutron production, fission, energy deposition, energy gain, testing volume and helium production are simulated using the MARS15 code for tungsten, thorium, and natural uranium targets in the proton energy range 0.5–120 GeV. This study shows that the proton energy range of 2–4 GeV is optimal for both a nat U EPD and the tungsten-based testing station for proton accelerator facilities. Simulation-based conservative estimates not including breeding and fission of plutonium suggest the proton beam current sufficient to produce 1 GW of thermal output power with the nat U EPD while supplying a relatively small fraction of that power to operate the accelerator. The thermal analysis has been performed and has shown that the EPD with a parallel proton beam feeding the core has a potential problem due to a possible core meltdown. In conclusion, a scheme has been proposed with a beam steering on the outer surface of the target; its thermal analysis indicates that slicing the target in shorter parts and maintaining a temperature of 25 °C on their surfaces would to avoid the target meltdown.

43 PARTICLE ACCELERATORS↗

PyEmittance: A General Python Package for Particle Beam Emittance Measurements With Adaptive Quadrupole Scans [Poster]

We present PyEmittance: a new Python package for general particle beam emittance measurements that offers adaptive quadrupole scans and is designed for robust and flexible integration with Python-based software, with online accelerators, and with simulation software. It is open source software, and can be installed via the command pip install pyemittance.

43 PARTICLE ACCELERATORS↗