Search NASA⌕ Search

SEARCH · Search NASA

Results for “beam 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 343 records · Page 19

Dynamics of a high-current relativistic electron beam

The dynamics of a high-current relativistic electron beam is studied experimentally and by numerical simulation. The beam is formed in a magnetically insulated diode with a transverse-blade explosive-emission cathode. It is found experimentally that the radius of a 500-keV beam with a current of 2 kA and duration of 500 ns decreases with time during the beam current pulse. The same effect was observed in numerical simulations. This effect is explained by a change in the shape of the cathode plasma during the current pulse, which, according to calculations, leads to a change in the beam parameters, such as the electron pitch angle and the spread over the longitudinal electron momentum. These parameters are hard to measure experimentally; however, the time evolution of the radial profile of the beam current density, which can be measured reliably, coincides with the simulation results. This allows one to expect that the behavior of the other beam parameters also agrees with numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomic oxygen beam source for erosion simulation

A device for the production of low energy (3 to 10 eV) neutral atomic beams for surface modification studies is described that reproduces the flux of atomic oxygen in low Earth orbit. The beam is produced by the acceleration of plasma ions onto a negatively biased plate of high-Z metal; the ions are neutralized and reflected by the surface, retaining some fraction of their incident kinetic energy, forming a beam of atoms. The plasma is generated by a coaxial RF exciter which produces a magnetically-confined (4 kG) plasma column. At the end of the column, ions fall through the sheath to the plate, whose bias relative to the plasma can be varied to adjust the beam energy. The source provides a neutral flux approximately equal to 5 x 10(exp 16)/sq cm at a distance of 9 cm and a fluence approximately equal to 10(exp 20)/sq cm in five hours. The composition and energy of inert gas beams was diagnosed using a mass spectometer/energy analyzer. The energy spectra of the beams demonstrate energies in the range 5 to 15 eV, and qualitatively show expected dependences upon incident and reflecting atom species and potential drop. Samples of carbon film, carbon-based paint, Kapton, mylar, and teflon exposed to atomic O beams show erosion quite similar to that observed in orbit on the space shuttle.

Cuthbertson, J. W.↗

Simulations of Hyperthermal Oxygen Beam Exposures

Pulsed sources of hyper-thermal O-atoms are now being extensively used to simulate low-earth orbit (LEO) surface exposure environments. The peak flux of these sources is many orders of magnitude larger than the corresponding LEO flux. Although it is desirable to accelerate the test by using higher fluxes than found in LEO, even commonly used fluxes are large enough to produce multi-collision effects by causing a build-up of gas at the sample surface. In this paper we characterize the physical consequences to the experiment using the direct simulation Monte Carlo (DSMC) method, DSMC allows us to extract the distributions of energy and impact angle for the O-atoms that reach the surface, and to record how strongly the gas build-up at the target assembly deflects flux from downstream instrumentation. By considering a range of source fluxes, we determine the onset conditions and severity of these multi-collision effects. We find that even at common experimental fluxes with a normally incident beam string a flat surface sample, the energy distribution of incident O-atoms broadens and develops a significant low-energy tail. The angular distributions also broaden significantly. The number of O-atoms that reach downstream instrumentation is decreased by approximately 50%. These simulations will aid in the calibration of ground-based O-atom measurements,a nd provide a model for the energy and angular distributions of O-atoms that actually impinge on surface samples.

Finchum, A.↗

DSMC Simulations of Hyperthermal Oxygen Beam Exposures

Pulsed sources of hyper-thermal O-atoms are now being extensively used to simulate low-earth orbit (LEO) surface exposure environments. The peak flux of these sources is many orders of magnitude larger than the corresponding LEO flux. Although it is desirable to accelerate the test by using higher fluxes than found in LEO, even commonly used fluxes are large enough to produce multi-collision effects by causing a build-up of gas at the sample surface. In this paper we characterize the physical consequences to the experiment using the direct simulation Monte Carlo (DSMC) method. DSMC allows us to extract the distributions of energy and impact angle for the O-atoms that reach the surface, and to record how strongly the gas build-up at the target assembly deflects flux from downstream instrumentation. By considering a range of source fluxes, we determine the onset conditions and severity of these multi-collision effects. We find that even at common experimental fluxes with a normally incident beam striking a flat surface sample, the energy distribution of incident O-atoms broadens and develops a significant low-energy tail. The angular distributions also broaden significantly. The number of O-atoms that reach downstream instrumentation is decreased by approximately 50%. These simulations will aid in the calibration of ground-based O-atom measurements, and provide a model for the energy and angular distributions of O-atoms that actually impinge on surface samples.

Cline, Jason A.↗

Uniformity simulation of multiple-beam irradiation of a spherical laser target with the inclusion of radiation absorption and refraction

We consider the physicomathematical model for the simulation of the angular laser flux distribution in a spherical target with the inclusion of radiation refraction in the corona. The irradiation uniformity is simulated using the target irradiation scheme on a megajoule multiple-beam laser facility. Calculations are made of the time-dependent angular distribution of absorbed flux and of the angular distribution of the time-integrated characteristic – the absorbed energy with allowance for the temporal shape of laser pulses – as well as of the time dependence of absorption efficiency. Different versions of target irradiation are considered. Under conventional irradiation it is assumed that the axes of all laser beams of equal energy pass through the target centre and the irradiation nonuniformity is due to its geometry only. Also examined are three unconventional irradiation versions: when there is a small spread of beam energies, when the target is displaced from the common beam pointing centre, and when the beam axes have small random deviations from the target centre. It is shown that the nonuniformity of the angular distribution of absorbed flux is greatest when the target is displaced from the beam pointing centre. (interaction of laser radiation with matter. laser plasma)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermal fluctuations from an artificial ion beam injection into the ionosphere

Simulations of argon beam experiments flown on two sounding rocket flights (ARCS 1, 2) reveal two bands of electrostatic wave emissions, one at high frequencies around the upper hybrid mode, the other at low frequencies around the lower hybrid mode. The latter is significantly enhanced by the argon beam in the simulations, and the enhancement was clearly observed on ARCS 2. This enhancement at phase velocities greater than the beam velocity is suggested to be due to an increase in the thermal fluctuation level of the plasma when the argon beam is present. Inclusion of electron dynamics and oblique angles of wave propagation with respect to B allows investigation of electron heating. Electron tail heating is observed parallel to B for k-parallel/k-perpendicular proportional to the square root of m(e)/m(0), while background ions are heated perpendicular to B.

Hudson, M. K.↗

Fully Kinetic Simulations of Proton-beam-driven Instabilities from Parker Solar Probe Observations

The expanding solar wind plasma ubiquitously exhibits anisotropic nonthermal particle velocity distributions. Typically, proton velocity distribution functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by the Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam population that drift with respect to each other. This configuration triggers a proton beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing perpendicular heating at the expense of the parallel temperature.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating the Clamped Tapered Beam Specimen Under Quasi-Static and Fatigue Loading Using Floating Node Method

As part of the NASA Advanced Composites Project (ACP), a sub-element has been designed to provide validation data for progressive damage analysis models. The clamped tapered beam is a cross-ply laminated composite specimen designed to validate the simulation of the onset of matrix cracks and their interaction with delaminations, including delamination migration. A tapered geometry was used to localize the first damage occurrence in the tapered region, without prescribing an initial crack. The boundary and loading conditions were chosen to favor delamination growth and subsequent migration after the first damage occurrence. The typical sequence of events consists of a matrix crack located at the tapered region, leading to delamination onset, followed by delamination growth and subsequent delamination migration to a different interface via a dominant matrix crack. The Clamped Tapered Beam (CTB) was tested in both quasi-static and fatigue regimes. The results obtained are used in this study to assess and validate a methodology based on the Floating Node Method (FNM) implemented as an Extended Interface Element. In this methodology, quasi-static and fatigue damage formation and development are modeled by combining FNM to represent crack networks, with Directional Cohesive Zone Elements (DCZE) and Virtual Crack Closure Technique (VCCT), respectively. Qualitatively, the methodology is capable of predicting the sequence of events and overall failure morphology. Quantitatively, the simulation results generally bound the experimental data, based on the range of the characterization data used. In this paper, the results from quasi static and fatigue simulations are compared and correlated with experimental data.

Seshadri, B. R.↗

Simulations of coherent electron cooling with two types of amplifiers

The SPACE code is a parallel, relativistic, three-dimensional, electromagnetic, particle-in-cell code developed for the simulations of beam dynamics and interactions. The SPACE code can perform simulations with various beam distributions, different types of boundary conditions, and flexible beam line, as well as provide sufficient processing routines for data analyses and visualization. Finally, self-consistent, highly resolved and realistic simulations were performed using the SPACE code to study the coherent electron cooling (CeC) experiment with two types of amplifiers, the free electron laser (FEL) and plasma-cascade micro-bunching amplifiers.

43 PARTICLE ACCELERATORS↗

Simulation Study of Hollow Electron Beam Collimation for LHC

In this note we describe the results of numerical simulations of hollow electron beam collimation for the nominal LHC machine. The halo cleaning rates for various operating scenarios are predicted. The impact of electron lens imperfections on the collider luminosity performance is estimated.

43 PARTICLE ACCELERATORS↗

The effect of voxel size on dose distribution in Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Recently, Monte Carlo (MC) calculation method has reported as the most accurate method of predicting dose distributions in radiotherapy. The MC code system (especially DOSXYZnrc) has been used to investigate the different voxel (volume elements) sizes effect on the accuracy of dose distributions. To investigate this effect on dosimetry parameters, calculations were made with three different voxel sizes. The effects were investigated with dose distribution calculations for seven voxel sizes: 1 × 1 × 0.1 cm{sup 3}, 1 × 1 × 0.5 cm{sup 3}, and 1 × 1 × 0.8 cm{sup 3}. The 1 × 10{sup 9} histories were simulated in order to get statistical uncertainties of 2%. This simulation takes about 9-10 hours to complete. Measurements are made with field sizes 10 × 10 cm2 for the 6 MV photon beams with Gaussian intensity distribution FWHM 0.1 cm and SSD 100.1 cm. MC simulated and measured dose distributions in a water phantom. The output of this simulation i.e. the percent depth dose and dose profile in d{sub max} from the three sets of calculations are presented and comparisons are made with the experiment data from TTSH (Tan Tock Seng Hospital, Singapore) in 0-5 cm depth. Dose that scored in voxels is a volume averaged estimate of the dose at the center of a voxel. The results in this study show that the difference between Monte Carlo simulation and experiment data depend on the voxel size both for percent depth dose (PDD) and profile dose. PDD scan on Z axis (depth) of water phantom, the big difference obtain in the voxel size 1 × 1 × 0.8 cm{sup 3} about 17%. In this study, the profile dose focused on high gradient dose area. Profile dose scan on Y axis and the big difference get in the voxel size 1 × 1 × 0.1 cm{sup 3} about 12%. This study demonstrated that the arrange voxel in Monte Carlo simulation becomes important.

61 RADIATION PROTECTION AND DOSIMETRY↗

Symplectic multi-particle tracking on GPUs

A symplectic multi-particle tracking model is implemented on the Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) language. The symplectic tracking model can preserve phase space structure and reduce non-physical effects in long term simulation, which is important for beam property evaluation in particle accelerators. Though this model is computationally expensive, it is very suitable for parallelization and can be accelerated significantly by using GPUs. In this paper, we optimized the implementation of the symplectic tracking model on both single GPU and multiple GPUs. Using a single GPU processor, the code achieves a factor of 2-10 speedup for a range of problem sizes compared with the time on a single state-of-the-art Central Processing Unit (CPU) node with similar power consumption and semiconductor technology. It also shows good scalability on a multi-GPU cluster at Oak Ridge Leadership Computing Facility. In an application to beam dynamics simulation, the GPU implementation helps save more than a factor of two total computing time in comparison to the CPU implementation.

43 PARTICLE ACCELERATORS↗

Characterisation of the muon beams for the Muon Ionisation Cooling Experiment

A novel single-particle technique to measure emittance has been developed and used to characterise seventeen different muon beams for the Muon Ionisation Cooling Experiment (MICE). The muon beams, whose mean momenta vary from 171 to 281 MeV/c, have emittances of approximately 1.5--2.3 \pi mm-rad horizontally and 0.6--1.0 \pi mm-rad vertically, a horizontal dispersion of 90--190 mm and momentum spreads of about 25 MeV/c. There is reasonable agreement between the measured parameters of the beams and the results of simulations. The beams are found to meet the requirements of MICE.

43 PARTICLE ACCELERATORS↗

The Accelerator Design Progress for EIC Strong Hadron Cooling

The Electron-Ion Collider will achieve a luminosity of 1034 cm-2 s-1 by incorporating strong hadron cooling to counteract hadron Intra-Beam Scattering, using a coherent electron cooling scheme. An accelerator will deliver the beams with key parameters, such as 1 nC bunch charge, and 1e-4 energy spread. The paper presents the design and beam dynamics simulation results. Methods to minimize beam noise, the challenges of the accelerator design, and the R&D topics being pursued are discussed.

Wang, E.↗

SU-E-T-314: The Application of Cloud Computing in Pencil Beam Scanning Proton Therapy Monte Carlo Simulation

Purpose: Monte Carlo simulation plays an important role for proton Pencil Beam Scanning (PBS) technique. However, MC simulation demands high computing power and is limited to few large proton centers that can afford a computer cluster. We study the feasibility of utilizing cloud computing in the MC simulation of PBS beams. Methods: A GATE/GEANT4 based MC simulation software was installed on a commercial cloud computing virtual machine (Linux 64-bits, Amazon EC2). Single spot Integral Depth Dose (IDD) curves and in-air transverse profiles were used to tune the source parameters to simulate an IBA machine. With the use of StarCluster software developed at MIT, a Linux cluster with 2–100 nodes can be conveniently launched in the cloud. A proton PBS plan was then exported to the cloud where the MC simulation was run. Results: The simulated PBS plan has a field size of 10×10cm{sup 2}, 20cm range, 10cm modulation, and contains over 10,000 beam spots. EC2 instance type m1.medium was selected considering the CPU/memory requirement and 40 instances were used to form a Linux cluster. To minimize cost, master node was created with on-demand instance and worker nodes were created with spot-instance. The hourly cost for the 40-node cluster was $0.63 and the projected cost for a 100-node cluster was $1.41. Ten million events were simulated to plot PDD and profile, with each job containing 500k events. The simulation completed within 1 hour and an overall statistical uncertainty of < 2% was achieved. Good agreement between MC simulation and measurement was observed. Conclusion: Cloud computing is a cost-effective and easy to maintain platform to run proton PBS MC simulation. When proton MC packages such as GATE and TOPAS are combined with cloud computing, it will greatly facilitate the pursuing of PBS MC studies, especially for newly established proton centers or individual researchers.

07 ISOTOPE AND RADIATION SOURCES↗

WE-EF-303-02: BEST IN PHYSICS (JOINT IMAGING- THERAPY): A Comprehensive Simulation of Image Guided Beam Gating for Liver Tumor Treatments Using Scanned Proton Therapy

Purpose: To evaluate the effectiveness of image guided beam gating for PBS liver treatments under realistic breathing conditions. Methods: We have previously proposed a Beams’ Eye View (BEV) X-ray image system as an online motion monitoring device for deriving a gating signal for PBS proton therapy. Using dedicated 4D dose calculations (4DDC), in this work we have simulated gated liver treatments using three amplitude-based gating windows (10/5/3mm) based on motion extracted from BEV imaging of fiducial markers or the diaphragm. In order to improve motion mitigation, BEV guided gating has also been combined with either volumetric (VS) or layered (LS) rescanning. Nine 4DCT(MRI) liver data-sets have been used for the investigation, which not only consider realistic patient geometries but also motion variations between different breathing cycles. All 4D plans have been quantified in terms of plan homogeneity in the PTV (D5-D95), the total estimated treatment time and the beam-on duty cycle. Results: Neither gating nor rescanning can fully retrieve a comparable plan homogeneity to the static case, and considerable reductions of the duty cycle (<10%) were observed as a Result motion variations when small gating windows are used. However, once combined with rescanning, dose homogeneity within 1% of the static plan could be achieved with reasonable prolongation of the treatment time for all 9 subjects. No differences were observed between the efficacy of layered or volumetric re-scanning, or of gating signals extracted from fiducial or diaphragm motions. However, layered rescanning may be preferred over volumetric rescanning when performed in combination with gating as it is generally more time-efficient and dosimetrically robust to patient and motion variations Conclusion Combining BEV beam gating with rescanning is an efficient and effective approach to treating mobile liver tumours, and is equally effective if either the diaphragm or fiducial markers are used as motion surrogates.

60 APPLIED LIFE SCIENCES↗

Commissioning of a Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Monte Carlo modelling of a linear accelerator is the first and most important step in Monte Carlo dose calculations in radiotherapy. Monte Carlo is considered today to be the most accurate and detailed calculation method in different fields of medical physics. In this research, we developed a photon beam model for Varian Clinac iX 6 MV equipped with MilleniumMLC120 for dose calculation purposes using BEAMnrc/DOSXYZnrc Monte Carlo system based on the underlying EGSnrc particle transport code. Monte Carlo simulation for this commissioning head LINAC divided in two stages are design head Linac model using BEAMnrc, characterize this model using BEAMDP and analyze the difference between simulation and measurement data using DOSXYZnrc. In the first step, to reduce simulation time, a virtual treatment head LINAC was built in two parts (patient-dependent component and patient-independent component). The incident electron energy varied 6.1 MeV, 6.2 MeV and 6.3 MeV, 6.4 MeV, and 6.6 MeV and the FWHM (full width at half maximum) of source is 1 mm. Phase-space file from the virtual model characterized using BEAMDP. The results of MC calculations using DOSXYZnrc in water phantom are percent depth doses (PDDs) and beam profiles at depths 10 cm were compared with measurements. This process has been completed if the dose difference of measured and calculated relative depth-dose data along the central-axis and dose profile at depths 10 cm is ≤ 5%. The effect of beam width on percentage depth doses and beam profiles was studied. Results of the virtual model were in close agreement with measurements in incident energy electron 6.4 MeV. Our results showed that photon beam width could be tuned using large field beam profile at the depth of maximum dose. The Monte Carlo model developed in this study accurately represents the Varian Clinac iX with millennium MLC 120 leaf and can be used for reliable patient dose calculations. In this commissioning process, the good criteria of dose difference in PDD and dose profiles were achieve using incident electron energy 6.4 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗

CEBAF Injector for K Long Beam Conditions

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab concurrently operates four experimental Halls with distinct bunch charge specifications and repetition rates. Numerous critical beam parameters within CEBAF are configured in the injector, some remaining unchanged throughout the accelerator. Consequently, the injector plays a crucial role in determining final beam characteristics, including bunch structure, beam sizes, bunch lengths, energy spread, and beam transmission. The Jefferson Lab KL experiment is scheduled to take place at CEBAF in Hall D, featuring a much lower bunch repetition rate of 7.80 MHz or 15.59 MHz, below the nominal values of 249.5 MHz or 499 MHz. Although the proposed average current of 5 ?A or 10 ?A is low compared to the maximum CEBAF cur- rent of approximately 180 ?A, the corresponding bunch charge is unusually high for CEBAF injector operation. This study focuses on the behavior of low-repetition-rate, high-bunch- charge (0.32 to 0.64 pC) beams in the CEBAF injector. We investigated the evolution and transmission of low-charge beams to space-charge dominated high-charge beams in the front end of the CEBAF injector for two configurations: the pre-existing CEBAF Phase 1 injector upgrade, operated at 130 kV, and the existing CEBAF Phase 2 injector upgrade, operated at 140 kV, 180 kV, and to be operated at 200 kV. The electron beam through the CEBAF injector is characterized using beam dynamics simulations and comparisons with the available measurements performed at 130 kV. Multi-objective genetic optimizations of the CEBAF injector were performed to determine the operating magnetic elements and RF settings for the evolution and transmission of low, moderate, and high charge beams in the CEBAF injector at 180 kV and 200 kV DC gun voltages. Subsequently, simulations at the same voltages were conducted to obtain the beam characteristics at the front end of the CEBAF injector. The laser spot size and laser pulse length at the cathode were varied to observe their effects on beam transmission in the injector at different voltages (130 kV, 180 kV, and 200 kV). Experimental studies at 130 kV, 140 kV, and 180 kV validate the simulations. Beam study measurements are carried out using EPICS tools, while optimizations and simulations are facilitated by General Particle Tracer. Based on the findings, optimal parameters for the upcoming Jefferson Lab KL experiment are proposed, utilizing a lower repetition rate and higher bunch charge

Pokharel, Sunil↗