SEARCH · Engineering Papers
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.
Proton accelerator, electron accelerators, and beam transport system of NASA space radiation effects laboratory - transient, permanent, surface, and chemical radiation effects
We present the first microscopic electron cooling simulations from first principles with accurate prediction of cooling time. These simulations were performed using our previously developed numerical method, PHAD, which is the first efficient large-scale collisional numerical method in beam physics. The simulation results are benchmarked with the experimental data of the low energy bunched electron cooling of ion beams at the storage ring CSRm at the IMP facility in China. Here, we have accurately considered the nonlinear dynamics in the whole accelerator system in addition to the electron cooling section. As a result, our simulations correctly reproduced cooling times of the experiments from first principles and without any tuning or fitting parameters in the code.
A realistic electrostatic simulation model is used to investigate the plasma environment in the near vicinity of a spacecraft during the injection of electon beams from the spacecraft. The model is described, and the cases of injection into vacuum and into a low-density plasma are addressed. The relationship of the two-dimensional results to the purely field-aligned one-dimensional simulations is discussed. The dependence of the results on the ambient plasma density is investigated. Then, the properties of a charge-neutral beam and the plasma response are examined for the case where the beam has nonzero velocity components parallel and perpendicular to the magnetic field. The parameters used in the simulation are described, and the properties of the beam injection into the vacuum and in the case where a plasma is present are given.
Using models for how energetic electrons heat the atmosphere during flares, we simulate the radiative-hydrodynamic response of the lower solar atmosphere to flare heating. The simulations account for much of the non-LTE, optically thick radiative transfer that occurs in the chromosphere. Our models predict an increase in white light continuum during the flare on the order of 20%, but this is highly sensitive to the electron beam flux used in the simulation. We find that a majority of the white light continuum originates in the chromosphere as a result of Balmer and Paschen recombinations, but a significant portion also forms in the photosphere which has been heated by radiative backwarming.
This work examines the multi-pass steering of six electron beams in an FFA arc ranging from approximately 10.5 GeV to 22 GeV. Shown here is an algorithm based on singular value decomposition (SVD) to successfully steer all six beams through the arc given precise knowledge of all beam positions at each of one hundred and one diagnostic locations with one hundred individual corrector magnets: that is successive application of SVD to different 100 × 101 response matrices—one for each beam energy. Further, a machine learning scheme is developed which only requires knowledge of the energy-averaged beam position at each location to provide equivalent steering. Extension of this scheme to other beam optics quantities as well as transverse and longitudinal coupling is explored.
GaAs-based DC high voltage photo-guns used at accelerators with extensive user programs must exhibit long photocathode operating lifetime. Achieving this goal represents a significant challenge for proposed high average current facilities that must operate at tens of milliamperes or more. Specifically, the operating lifetime is dominated by ion back-bombardment of the photocathode from ionized residual gas. While numerous experiments have been performed to characterize the operating lifetime under various conditions, detailed simulations of the ion back-bombardment mechanism that explains these experiments are lacking. Recently, a new user routine was implemented using the code General Particle Tracer (GPT) to simulate electron impact ionization of residual beam line gas and simultaneously track the incident electron, the ejected electron, and the newly formed ion. This new routine was benchmarked against analytical calculations and then applied to experiments performed at the CEBAF injector at the Thomas Jefferson National Accelerator Facility to study the effectiveness of limiting ions from entering the cathode-anode gap using a positively biased anode. These simulations were performed using detailed 3D field maps produced with CST Microwave Studio describing the photo-gun electrostatics. Discussion of the experiment and the application of this new GPT routine to model the experiments will be presented at the workshop.
We discuss a method for recovering complex structure factors from many simultaneously excited Bragg beam intensities is described. The method is applied to simulated transmission electron diffraction data over a wide range of crystal thickness and beam energies. The method is based on iterated projections between structure and scattering matrices, which are related by a matrix unitary transformation, exponential, which we invert. The algorithm removes multiple-scattering perturbations from diffraction data and might be extended to other fields, including X-ray and neutron diffraction and cryo-electron microscopy. Because coherent multiple scattering involves interference between Bragg beams, the method also solves the phase problem. Unlike dynamical inversion from electron microscope images or ptychography data, the method, which starts with Bragg beam intensities, provides complex structure factors unaffected by focusing errors or resolution limitations imposed by lenses. We provide inversions from simulated data with 441 simultaneously excited Bragg beams over a range of thickness and beam energy. We discuss the retrieval of chirality information from enantiomorphs, the efficient incorporation of symmetry information using the irreducible representation of the group of structure matrices, and the effect of HOLZ lines to provide three-dimensional information.
Abstract A new method to obtain the mean energy of fast electron losses in fusion plasmas using a versatile multi-energy hard x-ray (HXR) detector is presented. The method is based on measuring the thick-target emission of tungsten in the divertor region produced by fast electron losses interacting with the target and modeling the tungsten spectra by a Monte Carlo code which simulates the interaction between a beam of electrons and a solid target. The mean energy of the fast electron losses is determined through the comparison between the experimental and synthetic emission. The results show that fast electron losses during lower hybrid current drive discharges at WEST have a mean energy of 90–140 keV and represent only 2% of the total heat flux at the target. Additionally, anisotropic HXR emission has been detected for the first time at the WEST core and edge plasma, with opposite directions. It is due to the forward-peak emission of two distinctive populations of fast electrons: co-current fast electrons in the core and counter-current fast electron losses at the inner strike point. In view of future experiments like ITER where electron cyclotron current drive will generate a fast electron population, this technique could serve as a real-time monitor of fast electron losses and eventually feed an actuator on the current drive generation.
Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.
In this report, we use tracking simulations to investigate synchrotron phase space injection for electron accumulation in the electron storage ring of the Electron Ion Collider. Our simulation model accounts for both beam-beam interactions and lattice nonlinearities. Specifically, we examine how particle loss is influenced by various parameters. Additionally, we conduct a theoretical analysis and derive an analytical formula for the rapid evaluation of particle loss. Our results demon strate the feasibility of synchrotron phase-space injection for the electron storage ring and provide insights to guide parameter selection for the design of the injection line.
Here, an electro-optic sampling beam position monitor for installation at the SLAC National Accelerator Laboratory FACET-II facility is described. Simulations of the detector signal were performed using realistic electron beam current profiles, demonstrating that this single-shot, non-destructive, compact detector is capable of resolving the relative transverse offset and longitudinal separation of an ultra-relativistic, two-bunch electron beam to order 1 μm and 10 fs, respectively. A design study was performed to optimize the detector’s response to the ultra-relativistic two-bunch electron beam utilized in plasma wakefield accelerator experiments. In addition, the expected sensitivity to beam tilt is studied for typical electron beams used in free electron lasers.
Electron diffraction through a thin patterned silicon membrane can be used to create complex spatial modulations in electron distributions. By precisely varying parameters such as crystallographic orientation and wafer thickness, the intensity of reflections in the diffraction plane can be controlled and by placing an aperture to block all but one spot, we can form an image with different parts of the patterned membrane, as is done for bright-field imaging in microscopy. The patterned electron beams can then be used to control phase and amplitude of subsequent x-ray emission, enabling novel coherent x-ray methods. Furthermore, the electrons themselves can also be used for femtosecond time resolved diffraction and microscopy. As a first step toward patterned beams, we demonstrate experimentally and through simulation the ability to accurately predict and control diffraction spot intensities. We simulate MeV transmission electron diffraction patterns using the multislice method for various crystallographic orientations of a single crystal Si(001) membrane near beam normal. The resulting intensity maps of the Bragg reflections are compared to experimental results obtained at the Accelerator Structure Test Area Ultrafast Electron Diffraction (ASTA UED) facility at SLAC. Furthermore, the fraction of inelastic and elastic scattering of the initial charge is estimated along with the absorption of the membrane to determine the contrast that would be seen in a patterned version of the Si(001) membrane.
The design luminosity goal for the Electron-Ion Collider (EIC) is 1e34 cm⁻²s⁻¹. To achieve such a high luminosity, the EIC design adopts high bunch intensities, flat beams at the interaction point (IP) with a small vertical β^{*}-function, and a high collision frequency, together with crab cavities to compensate the geometrical luminosity loss due to the large crossing angle of 25mrad. In this article, we present our strategies and approaches to obtain the design luminosity by optimizing some key beam-beam related design parameters. Through our extensive strong-strong and weak-strong beam-beam simulations, we found that beam flatness, electron and proton beam size matching at the IP, electron and proton working points, and synchro-betatron resonances arising from the crossing angle collision play a crucial role in proton beam size growth and luminosity degradation. After optimizing those parameters, we found a set of beam-beam related design parameters to reach the design luminosity with an acceptable beam-beam performance.
Introduction. The Moon’s south pole is a high-priority target for human exploration and scientific study. This interest is, in part, due to the presence of Permanently Shadowed Regions (PSRs), which could contain high concentrations of unique volatiles at cryogenic temperatures [1]. Returned samples from PSRs may include a unique combination of rocks, regolith, and volatile species, providing unprecedented insights into the history of the Solar System and the potential for resource utilization on the Moon. However, because PSR samples are cryogenic up-on collection, lunar polar sample return will eventually require cold stowage for the journey from the Moon to Earth. Without cold stowage, PSR sample return will likely result in phase changes and chemical reactions within the volatile component of the sample, which would negatively impact the resulting scientific studies of those samples. This abstract summarizes the initial results from an ongoing characterization of analog PSR samples at a range of temperatures, with the goal of defining the temperatures needed for a flight cold stowage freezer. Background. Based on remote sensing observations of the Moon [2], south polar PSRs range in temperature from ~120K for small and/or shallow PSRs to ~20K at the most extreme locations in large, deep PSRs. A range of volatiles have been hypothesized to exist at the surface or subsurface of the lunar poles [3-5 and others]. This hypothesis was verified when the LCROSS mission impacted the <50-K PSR in the crater Cabeus, detecting a range of volatiles from water to low condensation temperature species such as H2S and methane [6]. Species such as H2S and ammonia (also detected by LCROSS) are also highly reactive, and increase the likelihood of chemical reactions at elevated (non-cryogenic) temperatures. At the Johnson Space Center’s Planetary Exploration and Astromaterials Research Laboratory (JSC-PEARL), we have developed a volatile-bearing lunar simulant that incorporates several of the species detected by LCROSS [Table 1] mixed cryogenically with the USGS Lunar Highlands Type (LHT) regolith simulant. The new volatile-regolith simulant will be used to assess the degree of sample alteration at room temperature, -20°C, -80°C, and -196°C (liquid nitro-gen), over a two-week period. Room temperature samples represent those likely to be returned during initial missions without cold stowage, -20°C provides an analog to Apollo cold curated samples, -80°C is the temperature of multiple flight payload freezers (e.g., MELFI), and -196°C is analogous to lunar PSRs. Two weeks is an approximation of the time between sample collection and Earth return for initial Artemis missions. Over this period of time, sample head-space gases will be analyzed using a Universal Gas Analyzer (UGA, a type of mass spectrometer) coupled with a Baratron pressure sensor. After testing, the regolith component of the simulant will be purged of volatiles and preserved for future electron beam and/or FTIR analysis. Experimental Procedure. Volatile-regolith simulants will be produced as an initial homogenous batch; this batch will then be distributed into aliquots (gas chromatography/GC vials or cryo vials), ensuring that each sample has the same starting composition and conditions [Fig. 1]. In addition to the “full” simulant shown in Table 1, less complex simulant compositions will be used as baseline and control samples [Table 2]. Aliquots will be produced in triplicate for each simulant composition, storage temperature, and date of sampling. Headspace gases in all Day 0 samples will be analyzed by the UGA immediately. Cold storage samples for future analytical days will be placed in freezers appropriate to their target temperatures (-20°C, -80°C, -196°C). For ambient-temperature samples, regolith and regolith-water samples will be stored in a fume hood, while the full simulant will be stored in a sealed Parr vessel for safety; no other simulants (RWCM/ RWCM+) will be stored at ambient temperature for this test. Samples will be analyzed by UGA in this manner on each Analysis Day outlined in Table 2. Analytical Data. The UGA measures the partial pressures in a single sample aliquot over a set mass range of 0-105 atomic mass units (AMU) [Figure 2]; this set range was selected to slightly exceed the mass of the highest-mass expected reaction product (H2SO4). The Baratron complements the UGA by measuring the total pressure in the headspace of a sample vial. Coupled together, the quantitative abundances of gases will be monitored throughout the test. UGA analyses of the triplicate samples for each storage temperature, day, and simulant composition will be averaged, and standard deviations for each will be calculated. The compositions of starting species (shown in Table 1) will be characterized as a function of time, and the presence of any new compounds (reaction products) will be monitored as well. Total pressures will be recorded for each sample analysis, and any samples that show signs of leakage (e.g., a significant reduction in pressure or simulant volatiles) will be discarded. Anticipated Results. Testing is planned to begin in January 2022. The resulting data will allow compositional and phase changes in the volatile component of the simulants to be determined. Both the reduction in initial compounds and the addition of reaction products are expected to be observed. In addition, the relative efficacy of the different temperatures at pre-serving the initial composition of the simulants will be quantified. Finally, the regolith component of each sample will be argon-purged and stored in a controlled environment for future laboratory analysis. Compositional and morphological changes in the regolith are expected for samples above 0°C. This test will be repeated three times over the course of 2022. Understanding the effect of temperature on both the volatile and regolith components of analog lunar materials will allow requirements for a cold stowage freezer to be developed. The implementation of cold stowage for lunar polar missions will maximize the preservation of returned samples, enabling ground-breaking lunar and Solar System volatiles science for decades to come.
The proton is the major component of the visible matter of our universe, so being able to determine its fundamental properties, such as its charge radius, is of vital importance. In recent years, several experiments have been performed for the accurate measurement of the proton radius with different experimental techniques. The main and most common are based on two different techniques: ep scattering and Lamb shift spectroscopy of regular hydrogen atoms. The latter can also be performed on muonic hydrogen atoms with a significant improvement in accuracy. Recent measurements that use this method have highlighted a discrepancy between the results obtained with the diffusion of electrons and those obtained with the spectroscopic technique leading to the so-called ?proton radius puzzle?. The explanation of why this difference between the measured values exists is still a matter of discussion. The last attempt to address the puzzle was the PRad Experiment. It was performed with a unique experimental setup that gave the scientists an excellent control on the systematic errors that are usually present in electron-proton scattering experiments. This has significantly increased the precision of the measurement, bringing it in the same area as the spectroscopic ones done on muonic hydrogen atoms, which are considered the most precise. However, to solve the puzzle once and for all, even more precise electron scattering measurements are needed. Alternatively, new methods can be used, such as eg. muon-proton scattering at low transferred momenta in order to lower the systematic errors that comes with the measurement. This thesis tries to address this problem by studying the set up of an intense muon beam at the Jefferson Laboratory (JLab) that can be used to carry out a high-precision m ? p scattering experiment. Using Montecarlo simulations I studied the muon production resulting from the interaction of the primary electron beam with the beam-dump at JLab Hall-A. The simulations were performed using the FLUKA package via the FLAIR graphical interface. In this study, I modified an existing geometry (built for the ?Beam Dump eXperiment?) and performed several studies in order to optimize the muon beam parameters. In order to estimate the beam intensity, I studied the muon generation and transportation through the concrete bunker surrounding the Hall-A beam-dump. Detailed studies on muon attenuation and dispersion at different depth allowed me to track the beam profile in space and momentum as a function of the distance from the beam-dump. Last, but not least, I studied a possible tagging system for the muon beam. The tagging system will allow us to know with good precision (O(1=1000)) the energy of each muon of the beam. Based on the results obtained, I created an event generator that incorporates the detailed description of the muon beam (energy, position, emittance) useful to define and optimize the experimental setup for a future measurement of the proton radius at the JLab.
Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.
Spectroscopic observations of optical emissions from the beam-plasma discharge (BPD) phenomenon were made with NASA's vacuum chamber facility, at the Johnson Space Center, configured to simulate the physical conditions of magnetospheric electron beam injection into the ionospheric/upper-atmospheric environment. Nonlinear N2 and N2(+) optical emission growth rates (with respect to incremental electron beam current values) were observed from the chamber gas during transition to the BPD state. For electron-beam currents (I) near the BPD transition value (I(sub c)), the band emissions from the chamber gas produced by relatively low energy (less than or equal to 50 eV) electrons interacting with N2 were anomalously more intense than those requiring higher energy (greater than 100 eV) electrons to excite them. For I more greater than I(sub c), the optical emissions increased linearly with I (as was the case for I less than I(sub c)) and their ratios decreased significantly from the peak values attained when I approximately equals I(sub c). These observations suggest that during BPD some of the energy of the primary electron beam is efficiently transferred, via wave-particle interactions, to local electrons produced through ionization of the chamber gas; the resulting suprathermal electrons provide an additional source of excitation for the relatively low energy states (A, B and C) of N2. Such nonlinear excitation of upper atmospheric gas may occur in certain auroral events wherein the current due to the precipitating electrons approaches a value close to I(sub c).
Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.