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 1,621 records · Page 90

Stable long range proton acceleration driven by intense laser pulse with underdense plasmas

Proton acceleration is investigated by 2.5-dimensional particle-in-cell simulations in an interaction of an ultra intense laser with a near-critical-density plasma. It was found that multi acceleration mechanisms contribute together to a 1.67 GeV collimated proton beam generation. The W-BOA (breakout afterburner based on electrons accelerated by a wakefield) acceleration mechanism plays an important role for the proton energy enhancement in the area far from the target. The stable and continuous acceleration maintains for a long distance and period at least several pico-seconds. Furthermore, the energy scalings are also discussed about the target density and the laser intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fingering convection induced by atomic diffusion in stars: 3D numerical computations and applications to stellar models

Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature, and detailed structure can affect the excitation process of various stellar pulsation modes and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, as well as its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection (which neglect radiative transfer and complex opacity effects) and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (within the scope of the aforementioned approximation) and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. prescription in the Toulouse-Geneva Evolution Code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al., that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases it ultimately triggers dynamical convection. The direct effects of radiative acceleration on the dynamics of fingering convection (especially in the nonlinear regime) remain to be added in the future to improve the quantitative predictions of the model.

79 ASTRONOMY AND ASTROPHYSICS↗

Coupling of chemical kinetics and thermodynamics for simulations of leaching of cement paste in ammonium nitrate solution

A novel numerical method is proposed in this paper to model the leaching of various types of cement pastes in the ammonium nitrate solution by coupling chemical kinetics with chemical thermodynamics. Arrhenius equation and Eyring equation are applied to predict the rate constant of fast dissolution of portlandite in the ammonium nitrate, and a modified rate law is developed by incorporating the saturation ratio of the reaction. The estimations of leached depths from the proposed method comply with the well-known square root-time relationship and they also agree well with the reported experimental results. Moreover, the evolutions of porosity and major hydrate compositions during leaching are obtained. It is found that a simple linear relation exists for the acceleration factor of leaching in 6 M of ammonium nitrate with respect to leaching in pure water regardless of the type of the cement clinker. - Highlights: • Chemical kinetics and thermodynamics are coupled to simulate the accelerated leaching. • The proposed method generates good estimations for leached depths and changes of porosity. • The proposed method can predict the detailed distributions of all the major cement hydrates. • A simple linear relation exists for predicting the widely scattered acceleration factors.

36 MATERIALS SCIENCE↗

RF Gap Transformation in PARMILA

The PARMILA (Phase And Radial Motion in Ion Linear Accelerators) code is used to design ion linacs and simulate the beam dynamics therein. The original code was developed in the 1960’s and later the Los Alamos Accelerator Code Group developed a PC version that has been used worldwide. Recently, in an effort to modernize the code and use it as the basis for a high-performance beam dynamics simulator with the LANSCE linac, the details of the code have come under additional scrutiny. PARMILA is a “z code” that uses various transformations to approximately represent the effect of the beam line element on the beam. In the case of the RF gap, it treats the particle dynamics inside the gap as a combination of drift + kick at the mid-plane + drift. Recently, we uncovered an undocumented keyword “off-exact” in PARMILA that enables the code to use the full modified-Bessel functions instead of low-order approximations in the RF gap transformation. This change improved the accuracy of the transformation for off-axis particles. However, a discrepancy between the analytical formula and the code was found for the transformation of the transverse coordinates. The purpose of this tech note is to document the derivation of the gap transformation and the related DTL cell length calculation used in the PARMILA code.

43 PARTICLE ACCELERATORS↗

Reconstruction of implanted marker trajectories from cone-beam CT projection images using interdimensional correlation modeling

Purpose: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. Methods: Because the superior–inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left–right and anterior–posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors’ simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. Results: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. Conclusions: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.

60 APPLIED LIFE SCIENCES↗

Nonthermally Dominated Electron Acceleration during Magnetic Reconnection in a Low-beta Plasma

This work was motivated by electron acceleration during solar flares. After some introductory remarks on proposed particle acceleration mechanisms and questions needing answers, dynamic simulations and simulation results are presented including energy spectra and the formation of the power law distribution. In summary, magnetic reconnection is highly efficient at converting the free magnetic energy stored in a magnetic shear and accelerating electrons to nonthermal energies in low-β regime. The nonthermal electrons have a dominant fraction and form power-law energy spectra with spectral index p ~ 1 in low-β regime. Electrons are preferentially accelerated along the curvature drift direction along the electric field induced by the reconnection outflow. The results can be applied to explain the observations of electron acceleration during solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. I. A Systematic Study

Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal-dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell (PIC) simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical evolution of magnetic field during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands (ultraviolet to X-ray) tend to flare earlier than lower-frequency bands (infrared to optical) in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the variation in temporal polarization depends strongly on the guide field strength. Specifically, reconnection with a significant guide field component leads to a very high polarization degree that contradicts typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly antiparallel magnetic field lines. These rotations are at least 90° and can extend to >180°, and they may be in either direction. Furthermore, these results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly antiparallel magnetic field morphology.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulation of radiation damping in rings, using stepwise ray-tracing methods

The ray-tracing code Zgoubi computes particle trajectories in arbitrary magnetic and/or electric field maps or analytical field models. It includes a built-in fitting procedure, spin tracking many Monte Carlo processes. The accuracy of the integration method makes it an efficient tool for multi-turn tracking in periodic machines. Energy loss by synchrotron radiation, based on Monte Carlo techniques, had been introduced in Zgoubi in the early 2000s for studies regarding the linear collider beam delivery system. However, only recently has this Monte Carlo tool been used for systematic beam dynamics and spin diffusion studies in rings, including eRHIC electron-ion collider project at the Brookhaven National Laboratory. Some beam dynamics aspects of this recent use of Zgoubi capabilities, including considerations of accuracy as well as further benchmarking in the presence of synchrotron radiation in rings, are reported here.

43 PARTICLE ACCELERATORS↗

Spinoff 2013

Topics covered include: Innovative Software Tools Measure Behavioral Alertness; Miniaturized, Portable Sensors Monitor Metabolic Health; Patient Simulators Train Emergency Caregivers; Solar Refrigerators Store Life-Saving Vaccines; Monitors Enable Medication Management in Patients' Homes; Handheld Diagnostic Device Delivers Quick Medical Readings; Experiments Result in Safer, Spin-Resistant Aircraft; Interfaces Visualize Data for Airline Safety, Efficiency; Data Mining Tools Make Flights Safer, More Efficient; NASA Standards Inform Comfortable Car Seats; Heat Shield Paves the Way for Commercial Space; Air Systems Provide Life Support to Miners; Coatings Preserve Metal, Stone, Tile, and Concrete; Robots Spur Software That Lends a Hand; Cloud-Based Data Sharing Connects Emergency Managers; Catalytic Converters Maintain Air Quality in Mines; NASA-Enhanced Water Bottles Filter Water on the Go; Brainwave Monitoring Software Improves Distracted Minds; Thermal Materials Protect Priceless, Personal Keepsakes; Home Air Purifiers Eradicate Harmful Pathogens; Thermal Materials Drive Professional Apparel Line; Radiant Barriers Save Energy in Buildings; Open Source Initiative Powers Real-Time Data Streams; Shuttle Engine Designs Revolutionize Solar Power; Procedure-Authoring Tool Improves Safety on Oil Rigs; Satellite Data Aid Monitoring of Nation's Forests; Mars Technologies Spawn Durable Wind Turbines; Programs Visualize Earth and Space for Interactive Education; Processor Units Reduce Satellite Construction Costs; Software Accelerates Computing Time for Complex Math; Simulation Tools Prevent Signal Interference on Spacecraft; Software Simplifies the Sharing of Numerical Models; Virtual Machine Language Controls Remote Devices; Micro-Accelerometers Monitor Equipment Health; Reactors Save Energy, Costs for Hydrogen Production; Cameras Monitor Spacecraft Integrity to Prevent Failures; Testing Devices Garner Data on Insulation Performance; Smart Sensors Gather Information for Machine Diagnostics; Oxygen Sensors Monitor Bioreactors and Ensure Health and Safety; Vision Algorithms Catch Defects in Screen Displays; and Deformable Mirrors Capture Exoplanet Data, Reflect Lasers.

Source record↗

Constraints on Cosmic Rays, Magnetic Fields, and Dark Matter from Gamma-ray Observations of the Coma Cluster of Galaxies with VERITAS and FERMI

Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated and, on account of weak energy losses, can accumulate, suggesting that clusters may also be sources of very high energy (VHE; E greater than100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cerenkov telescopes, with complementing Fermi Large Area Telescope observations at GeV energies. No significant gamma-ray emission from the Coma Cluster was detected. Integral flux upper limits at the 99 confidence level were measured to be on the order of (2-5) x 10(sup -8) photons m(sup -2) s(sup -1) (VERITAS,greater than 220 GeV) and approximately 2 x 10(sup -6) photons m(sup -2) s(sup -1) (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain cosmic rays (CRs) and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be less than 16% from VERITAS data and less than 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be 50. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of approximately (2-5.5)microG, depending on the radial magnetic field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this renders the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, (sigma upsilon)

cosmic↗

TU-H-CAMPUS-TeP1-02: Seated Treatment: Setup Uncertainty Comparable to Supine

Purpose: For some head and neck patients, positioning in the supine position is not well tolerated. For these patients, treatment in a seated position would be preferred. We have evaluated inter- and intra- fraction uncertainty of patient set-up in a novel treatment chair which is compatible with modern linac designs. Methods: Five head-and-neck cancer patients were positioned in the chair, fitted with immobilization devices, and imaged with orthogonal X-rays. The couch (with chair attached) was rotated to simulate delivery (without actual treatment), another set of images were acquired, providing a measure of intra-fraction displacement. The patient then got off of and back onto the chair and the process was repeated, thus providing a measure of inter-fraction set-up uncertainty. Six sub-regions in the head-and-neck were rigidly registered to evaluate local intra- and interfraction displacement. Image guidance was simulated by first registering one sub-region; the residual displacement of other sub-regions was then measured. Additionally, a patient questionnaire was administered to evaluate tolerance of the seated position. Results: The chair design is such that all advantages of couch motions may be utilized. Average inter- and intrafraction displacements of all sub-regions in the seated position were less than 2 and 3 mm, respectively. When image guidance was simulated, interfraction displacements were reduced by an average of 4 mm, providing comparable setup to the supine position. The enrolled patients, who had no indication for a seated treatment position, reported no preference for the seated or the supine position. Conclusion: The novel chair design provides acceptable inter- and intra-fraction displacement, with reproducibility similar to that observed for patients in the supine position. Such a chair will be utilized for patients who cannot tolerate the supine position and use with CBCT images for planning, in a fixed-beam linac system, and for other treatment sites is under investigation. Funding: Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

Discrete Element Method Simulation of a Boulder Extraction From an Asteroid

The force required to pull 7t and 40t polyhedral boulders from the surface of an asteroid is simulated using the discrete element method considering the effects of microgravity, regolith cohesion and boulder acceleration. The connection between particle surface energy and regolith cohesion is estimated by simulating a cohesion sample tearing test. An optimal constant acceleration is found where the peak net force from inertia and cohesion is a minimum. Peak pulling forces can be further reduced by using linear and quadratic acceleration functions with up to a 40% reduction in force for quadratic acceleration.

Kulchitsky, Anton K.↗

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗

The effect of the near earth micrometeoroid environment on a highly reflective mirror surface

A resurgence of interest in placing large solar concentrator solar dynamic systems in space for power generation has brought up again a concern for maintaining the integrity of the optical properties of highly specular reflecting surfaces in the near earth space environment. One of the environmental hazards needing evaluation is the micrometeoroid environment. It has been shown that highly reflective polished metals and thin film coatings degrade when exposed to simulated micrometeoroids in the lab. At NASA-Lewis, a shock tube was used to simulate the phenomenon of micrometeoroid impact by accelerating micron sized particles to hypervelocities. Any changes in the optical properties of surfaces exposed to this impact were then evaluated. The degradation of optical properties of polished metals and thin metallic films after exposure to simulated micrometeoroids was determined as a function of impacting kinetic energy area of the particles. A calibrated sensor was developed to not only detect the micrometeoroid environment, but also to evaluate the degradation of the optical properties of thin aluminum films in space. Results of the simulation are presented and discussed.

Mirtich, Michael J.↗

The effect of the near earth micrometeoroid environment on a highly reflective mirror surface

A resurgence of interest in placing large solar concentrator solar dynamic systems in space for power generation has brought up again a concern for maintaining the integrity of the optical properties of highly specular reflecting surfaces in the near earth space environment. One of the environmental hazards needing evaluation is the micrometeoroid environment. It has been shown that highly reflective polished metals and thin film coatings degrade when exposed to simulated micrometeoroids in the lab. At NASA-Lewis, a shock tube was used to simulate the phenomenon of micrometeoroid impact by accelerating micron sized particles to hypervelocities. Any changes in the optical properties of surfaces exposed to this impact were then evaluated. The degradation of optical properties of polished metals and thin metallic films after exposure to simulated micrometeoroids was determined as a function of impacting kinetic energy area of the particles. A calibrated sensor was developed to not only detect the micrometeoroid environment, but also to evaluate the degradation of the optical properties of thin aluminum films in space. Results of the simulation are presented and discussed.

Mirtich, Michael J.↗

Simulation of automatic rotorcraft nap-of-the-earth flight in graphics workstation environment

This paper describes a three-dimensional (3D) helicopter flight simulation system. The simulation is designed to be a readily available tool for concept verification and tuning of automatic obstacle-avoidance guidance algorithms. The system has been implemented on networked workstations capable of interactive 3D graphics simulation. The simulation uses realistic terrain and obstacle models. The dynamics of the rotorcraft and the functional capabilities of the range sensors are simulated to provide all the components required to evaluate the guidance function. Standard graphics hardware available on the workstation is utilized to accelerate the range-data calculations for sensor simulation at the guidance rate. An example is given to demonstrate the performance of the obstacle-avoidance capability.

Lam, T.↗