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.
Wind tunnel, simulator and flight test results to evaluate cockpit accelerations, handling qualities and stability and control of jet transports in severe turbulence
Celeritas [1] is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications (specifically, High Lumi- nosity Large Hadron Collider (HL-LHC) simulation) on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a GPU-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.3, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling integration of Celeritas into high energy physics (HEP) experimental frameworks such as CMSSW. On the Summit supercomputer, Celeritas performs EM physics between 6 and 32 faster using the machine’s Nvidia GPUs compared to using only CPUs. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.8–2.3× on GPU and 1.2× on CPU.
As a part of the effort to expand the capabilities of CE-BAF 12 GeV (Continuous Electron Beam Accelerator Facility) at Jefferson Lab, the addition of a polarized positron source is considered. This capability would provide acceleration of high duty-cycle polarized posi-trons, with spin >60% polarization, through the same main CEBAF accelerator machine with appropriate mag-net field reversals and linac phasing to the four CEBAF experimental halls. To produce this positron beam, a high average current (3-10 mA) highly polarized electron beam with energy of 100 ? 150 MeV is required at the positron source target. The focus of this paper is the de-sign of that polarized electron beam injector. We will describe the production and delivery of a >3 mA highly polarized electron beam. We will discuss different aspects of the design, the photocathode gun, beam dynamics simulation results, spin manipulation, bunching and accelerating process and final electron beam parameters.
In negative ion electrostatic accelerators a considerable fraction of extracted ions is lost by collision processes causing efficiency loss and heat deposition over the components. Stripping is proportional to the local density of gas, which is steadily injected in the plasma source; its pumping from the extraction and acceleration stages is a key functionality for the prototype of the ITER Neutral Beam Injector, and it can be simulated with the 3D code AVOCADO. Different geometric solutions were tested aiming at the reduction of the gas density. The parameter space considered is limited by constraints given by optics, aiming, voltage holding, beam uniformity, and mechanical feasibility. The guidelines of the optimization process are presented together with the proposed solutions and the results of numerical simulations.
Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.
While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, in this work we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, as it opens up the possibility of modeling the energetic spectra without resolving microscopic diffusion regions.
The present evaluation of the basic techniques and illustrative results of charged particle-modeling numerical codes suitable for particle acceleration at oblique, fast-mode collisionless shocks emphasizes the treatment of ions as test particles, calculating particle dynamics through numerical integration along exact phase-space orbits. Attention is given to the acceleration of particles at planar, infinitessimally thin shocks, as well as to plasma simulations in which low-energy ions are injected and accelerated at quasi-perpendicular shocks with internal structure.
We present an analytical theory that reveals the importance of the longitudinal laser electric field in the course of the resonant acceleration of relativistic electrons by a tightly confined laser beam. It is shown that this laser field component always counteracts the transverse one and effectively decreases the final energy gain of electrons via the direct laser acceleration (DLA) mechanism. This effect is demonstrated by carrying out particle-in-cell simulations of the DLA of the electrons injected into the accelerating phase of the plasma wake. It is shown that the electron energy gain from the wakefield is substantially compensated by the quasiresonant energy loss to the longitudinal laser field component. The analytically obtained scalings and estimates are in good agreement with the results of the numerical simulations.
Laser photo-detachment is used as a method to measure or determine the negative ion density and temperature in electronegative plasmas. In essence, the method consists of producing an electropositive channel (negative ion free region) via pulsed laser photo-detachment within an electronegative plasma bulk. Electrostatic probes placed in this channel measure the change in the electron density. A second pulse might be used to track the negative ion recovery. From this, the negative ion density and temperature can be determined. We study the formation and relaxation of the electropositive channel via a two-dimensional Particle-In-Cell/Mote Carlo collision model. The simulation is mainly carried out in a Hydrogen plasma with an electronegativity of α = 1, with a parametric study for α up to 20. The temporal and spatial evolution of the plasma potential and the electron densities shows the formation of a double layer (DL) confining the photo-detached electrons within the electropositive channel. This DL evolves into two fronts that move in the opposite directions inside and outside of the laser spot region. As a consequence, within the laser spot region, the background and photo-detached electron energy distribution function relaxes/thermalizes via collisionless effects such as Fermi acceleration and Landau damping. Moreover, the simulations show that collisional effects and the DL electric field strength might play a non-negligible role in the negative ion recovery within the laser spot region, leading to a two-temperature negative ion distribution. The latter result might have important effects in the determination of the negative ion density and temperature from laser photo detachment diagnostic.
Numerical simulations become useful for the developing RF-ICP (Radio Frequency Inductively Coupled Plasma) negative ion sources. We are developing and parallelizing a two-dimensional three velocity electromagnetic Particle-In-Cell code. The result shows rapid increase in the electron density during the density ramp-up phase. A radial electric field due to the space charge is produced with increase in the electron density and the electron transport in the radial direction is suppressed. As a result, electrons stay for a long period in the region where the inductive electric field is strong, and this leads efficient electron acceleration and a rapid increasing of the electron density.
We present a general method to accelerate ab initio atomic structural relaxations. In this method, the conjugate gradient ab initio relaxation is preconditioned by a force field relaxation. This force field is constructed on-the-flight with the information of the ab initio forces and the current atomic configuration. At each ab initio relaxation step, the so-constructed force field is used to prerelax the system, with the relaxation direction as the preconditioned direction for the ab initio conjugated gradient method. The force field model and its parameters are rather general making this method applicable for a wide range of systems commonly used in material simulations. We report more than 80 different systems have been tested representing different cases in material simulations. Across the board, we find accelerations mostly over a factor of 2 and for some large systems with a factor of 6 ~ 9 . There is no case where the relaxation becomes worse. The code and lib for this method are provided, which can be used as a plugin in a standard ab initio atomic relaxation code.
An experimental neutron facility based on 32 MeV/18.47 kW electron linac has been studied by means of PARMELA simulation code. Beam dynamics study for a traveling wave constant gradient electron accelerator is carried out to reach the preferential operation parameters (E = 30 MeV, P = 18 kW, dE/E < 12.47% for 99% particles). The whole linac comprises mainly E-gun, pre-buncher, buncher, and 2 accelerating columns. A disk-loaded, on-axis-coupled, 2π/3-mode type accelerating rf cavity is considered for this linac. After numerous optimizations of linac parameters, 32 MeV beam energy is obtained at the end of the linac. As high electron energy is required to produce acceptable neutron flux. The final neutron flux is estimated to be 5 × 10{sup 11} n/cm{sup 2}/s/mA. Future development will be the real design of a 30 MeV electron linac based on S band traveling wave.
Spin tracking studies in the MEIC figure-8 collider ion ring are presented, based on a very preliminary design of the lattice. They provide numerical illustrations of some of the aspects of the figure-8 concept, including spin-rotator based spin control, and lay out the path towards a complete spin tracking simulation of a figure-8 ring.
It has recently been demonstrated that self-consistent particle-in-cell simulations of low-obliquity pulsar magnetospheres in flat spacetime show weak particle acceleration and no pair production near the poles. We investigate the validity of this conclusion in a more realistic spacetime geometry via general-relativistic particle-in-cell simulations of the aligned pulsar magnetosphere with pair formation. We find that the addition of the frame-dragging effect makes the local current density along the magnetic field larger than the Goldreich–Julian value, which leads to unscreened parallel electric fields and the ignition of a pair cascade. When pair production is active, we observe field oscillations in the open field bundle, which could be related to pulsar radio emission. We conclude that general-relativistic effects are essential for the existence of the pulsar mechanism in low-obliquity rotators.
Purpose: Photons from a modern high-energy therapeutic linear accelerator used in X-ray radiotherapy causes photonuclear reactions in an accelerator or patient's body. The aim of this study is to evaluate the biological effectiveness including these particles by Microdosimetric Kinetic Model (MKM) based on microdosimetry. Methods: A linear accelerator operating at 15 MV was used. CR-39 was used to obtain LET spectra of secondary ions selectively, as CR-39 is regarded insensitive to photons. CR-39 was put on the central axis of the X-ray beam at depths of 0, 5 and 10 cm in plastic phantom at a source to detector distance of 100 cm. Pits formed by the traversal of ions were etched then analyzed to obtain restricted LET distribution. Frequency-mean and dose-mean lineal energy was evaluated from the relationship between the restricted LET and the lineal energy required to evaluate the biological effectiveness by MKM. The relationship was calculated by Monte Carlo simulations with GEANT4. Results: Restricted LET distributions of secondary particles showed broad distributions that decreases exponentially with increasing LET. Frequency-mean and dose-mean lineal energy were determined uniquely within the scope of the energies of secondary particles generated from photons of 15 MeV. The frequency-mean lineal energies at the depth of 0, 5 and 10 cm were 15.1, 16.0 and 19.7 keV/μm respectively, and the dose-mean lineal energies were 18.6, 20.5 and 19.6 keV/μm, respectively. RBE of secondary particles for HSG cell evaluated by MKM was about 2.0 at all depths, and RBE of all particles including photons was evaluated 1.0. Conclusion: We investigated the biological effectiveness of secondary particles by photonuclear reactions. The method to evaluate RBE by MKM was established with measurements and simulations. However, the influence of these secondary ions on RBE was found negligible in the entire biological effectiveness of the high-energy X-ray. This study has been supported by JSPS KAKENHI Grant Number 25861144.
We quantify initial-data uncertainties on a shock accelerated heavy-gas cylinder by two-dimensional well-resolved direct numerical simulations. A high-resolution compressible multicomponent flow simulation model is coupled with a polynomial chaos expansion to propagate the initial-data uncertainties to the output quantities of interest. The initial flow configuration follows previous experimental and numerical works of the shock accelerated heavy-gas cylinder. We investigate three main initial-data uncertainties, (i) shock Mach number, (ii) contamination of SF{sub 6} with acetone, and (iii) initial deviations of the heavy-gas region from a perfect cylindrical shape. The impact of initial-data uncertainties on the mixing process is examined. The results suggest that the mixing process is highly sensitive to input variations of shock Mach number and acetone contamination. Additionally, our results indicate that the measured shock Mach number in the experiment of Tomkins et al. [“An experimental investigation of mixing mechanisms in shock-accelerated flow,” J. Fluid. Mech. 611, 131 (2008)] and the estimated contamination of the SF{sub 6} region with acetone [S. K. Shankar, S. Kawai, and S. K. Lele, “Two-dimensional viscous flow simulation of a shock accelerated heavy gas cylinder,” Phys. Fluids 23, 024102 (2011)] exhibit deviations from those that lead to best agreement between our simulations and the experiment in terms of overall flow evolution.
Machine Learning and Artificial Intelligence has been identified as an emerging priority science area within the Department of Energy. Large scale accelerator based supercomputers like Summit, while traditionally employed for modeling and simulation, provide architectures that are suitable for accelerating the ML/AI workloads at scale. With the release of Summit in 2018, there was an increase in the number of ML/AI based projects seeking time on the machine. It quickly became apparent that the allocations and job runtimes for this workload deviated from traditional large scale modeling and simulation. Accommodating this new workload requires understanding the impacts to traditional large scale modeling and simulation.
Coupling DSMC simulations with classical trajectories calculations is emerging as a powerful tool to improve predictive capabilities of computational rarefied gas dynamics. The considerable increase of computational effort outlined in the early application of the method (Koura,1997) can be compensated by running simulations on massively parallel computers. In particular, GPU acceleration has been found quite effective in reducing computing time (Ferrigni,2012; Norman et al.,2013) of DSMC-CT simulations. The aim of the present work is to study rarefied Oxygen flows by modeling binary collisions through an accurate potential energy surface, obtained by molecular beams scattering (Aquilanti, et al.,1999). The accuracy of the method is assessed by calculating molecular Oxygen shear viscosity and heat conductivity following three different DSMC-CT simulation methods. In the first one, transport properties are obtained from DSMC-CT simulations of spontaneous fluctuation of an equilibrium state (Bruno et al, Phys. Fluids, 23, 093104, 2011). In the second method, the collision trajectory calculation is incorporated in a Monte Carlo integration procedure to evaluate the Taxman’s expressions for the transport properties of polyatomic gases (Taxman,1959). In the third, non-equilibrium zero and one-dimensional rarefied gas dynamic simulations are adopted and the transport properties are computed from the non-equilibrium fluxes of momentum and energy. The three methods provide close values of the transport properties, their estimated statistical error not exceeding 3%. The experimental values are slightly underestimated, the percentage deviation being, again, few percent.