Search NASA⌕ Search

SEARCH · Search NASA

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 739 records · Page 41

Simulation of Cascaded Longitudinal-Space-Charge Amplifier at the Fermilab Accelerator Science & Technology (Fast) Facility

Cascaded Longitudinal Space Charge Amplifiers (LSCA) have been proposed as a mechanism to generate density modulation over a board spectral range. The scheme has been recently demonstrated in the optical regime and has confirmed the production of broadband optical radiation. In this paper we investigate, via numerical simulations, the performance of a cascaded LSCA beamline at the Fermilab Accelerator Science & Technology (FAST) facility to produce broadband ultraviolet radiation. Our studies are carried out using elegant with included tree-based grid-less space charge algorithm.

43 PARTICLE ACCELERATORS↗

Impact of electron transport models on capillary discharge plasmas

Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on the capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing next-generation plasma accelerators. To do so, we perform two-dimensional high-resolution MHD simulations using an argon-filled capillary discharge waveguide with three different electron transport coefficients models. The models tested include (i) Davies et al., (ii) Spitzer, and (iii) Epperlein–Haines (EH). We found that the EH model overestimates the electron temperature inside the channel by over 20% while predicting a lower azimuthal magnetic field. Furthermore, the Spitzer model, often used in MHD simulations for plasma-based accelerators, predicts a significantly higher electron temperature than the other models suggest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SuperNeuro: A Fast and Scalable Simulator for Neuromorphic Computing

In many neuromorphic workflows, simulators play a vital role for important tasks such as training spiking neural networks, running neuroscience simulations, and designing, implementing, and testing neuromorphic algorithms. Currently available simulators cater to either neuroscience workflows (e.g., NEST and Brian2) or deep learning workflows (e.g., BindsNET). Problematically, the neuroscience-based simulators are slow and not very scalable, and the deep learning-based simulators do not support certain functionalities that are typical of neuromorphic workloads (e.g., synaptic delay). In this paper, we address this gap in the literature and present SuperNeuro, which is a fast and scalable simulator for neuromorphic computing capable of both homogeneous and heterogeneous simulations as well as GPU acceleration. We also present preliminary results that compare SuperNeuro to widely used neuromorphic simulators such as NEST, Brian2, and BindsNET in terms of computation times. We demonstrate that SuperNeuro can be approximately 10×--300× faster than some of the other simulators for small sparse networks. On large sparse and large dense networks, SuperNeuro can be approximately 2.2×--3.4× faster than the other simulators, respectively.

Date, Prasanna↗

Comparisons of time explicit hybrid kinetic-fluid code Architect for Plasma Wakefield Acceleration with a full PIC code

Architect, a time explicit hybrid code designed to perform quick simulations for electron driven plasma wakefield acceleration, is described. In order to obtain beam quality acceptable for applications, control of the beam-plasma-dynamics is necessary. Particle in Cell (PIC) codes represent the state-of-the-art technique to investigate the underlying physics and possible experimental scenarios; however PIC codes demand the necessity of heavy computational resources. Architect code substantially reduces the need for computational resources by using a hybrid approach: relativistic electron bunches are treated kinetically as in a PIC code and the background plasma as a fluid. Cylindrical symmetry is assumed for the solution of the electromagnetic fields and fluid equations. In this paper both the underlying algorithms as well as a comparison with a fully three dimensional particle in cell code are reported. The comparison highlights the good agreement between the two models up to the weakly non-linear regimes. In highly non-linear regimes the two models only disagree in a localized region, where the plasma electrons expelled by the bunch close up at the end of the first plasma oscillation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Prioritization of Existing Reactor Materials

The Advanced Materials and Manufacturing Technologies (AMMT) Program is aiming at the faster incorporation of new materials and manufacturing technologies into complex nuclear-related systems. An integrated approach, combining advanced characterization, high-throughput and accelerated testing, modeling and simulation, including machine learning and artificial intelligence, will be employed. Although 316H (Fe–[16–18]Cr–[10–14]Ni–[2–3]Mo–[0.04–0.1]C) has been identified as a key alloy to be integrated into the AMMT accelerated alloy qualification approach because of its relevance for many current and future nuclear energy reactors, many other alloys could be considered for the advanced fabrication of innovative, high-performance nuclear components. Argonne National Laboratory (ANL), Idaho National Laboratory (INL), Oak Ridge National Laboratory (ORNL), and Pacific Northwest National Laboratory (PNNL) are collaborating on identifying the most promising alloy candidates relevant for the AMMT Program. A selection criteria matrix was established to evaluate the alloys considering their relative importance and technological readiness levels for nuclear energy applications, with a focus on laser powder bed fusion (LPBF). Because of the broad range of potential candidate alloys, ORNL and INL focused on nickel-based alloys, and ANL and PNNL mainly evaluated iron-based alloys. PNNL previously published material scorecards reports on several key alloys, and this report provides a broader overview of iron- and nickel-based candidate alloys, expending beyond alloys well-known to the nuclear community.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Preliminary Feasibility of Printing, Microstructure Analysis and Mechanical Performance of a Down Selected Ni Alloy

This report is submitted as completion of a Milestone 3 deliverable under work package ORNL CT-23OR1304051 in support of the Advanced Materials and Manufacturing Technologies (AMMT) program. The AMMT program is aiming at the faster incorporation of new materials and manufacturing technologies into complex nuclear-related systems. An integrated approach, combining advanced characterization, high-throughput and accelerated testing, modeling and simulation including machine learning and artificial intelligence will be employed. While 316H (Fe-16-18Cr-10-14Ni-2-3Mo-0.04-0.1C) has been identified as a key alloy to be integrated into the AMMT accelerated alloy qualification approach due its relevance for many current and future nuclear energy reactors, many other alloys could be considered for the advanced fabrication of innovate high-performance nuclear components. ANL, INL, ORNL and PNNL are collaborating on identifying the most promising alloy candidates relevant for the AMMT program. A selection criteria matrix was established to evaluate the alloys considering their relative importance and technological readiness levels for nuclear energy applications, with a focus on laser powder bed fusion (LPBF). Due to the broad range of potential candidate alloys, ORNL and INL focused on Ni-based alloys, while ANL and PNNL mainly evaluated Fe-based alloys. PNNL previously published materials scorecards reports on several key alloys and this report is providing a broader overview of Ni-based candidate alloys, expending beyond alloys well-known to the nuclear community. Of particular interest are alloys that are currently commercially available in powder form due to the growing demand from other industries that have invested heavily in additive manufacturing. Among these alloys, Haynes 282 (Fe-20Cr-20Co-8Mo) was selected due to its superior strength at high temperature compared to the code-qualified alloy 617 (Ni-20-24Cr-10-15Co-8-10Mo). To assess the integration of this new alloy into the AMMT digital manufacturing framework, we performed the rapid optimization of alloy 282 printing parameters on a Renishaw 250 machine, the fabrication of sufficient materials for extensive characterization and mechanical testing both at ORNL and INL, and added the printing data into the digital platform via the Peregrine software. A detailed analysis of the LPBF 718 alloy was also conducted, with creep specimens being tested at 600-650°C and characterized by advanced electron microscopy. The alloy superior mechanical strength and great printability associated with the extensive database that has already been generated highlight the promising potential of LPBF 718 as a candidate alloy for the AMMT program. INL, ANL and PNNL have generated similar reports and all the information will be compiled into a final M2 milestone report to provide the AMMT leadership team with clear recommendations on the down selection of reactor materials, as well as establish a roadmap for the qualification of these selected alloys.

36 MATERIALS SCIENCE↗

Modeling MHD accretion-ejection: episodic ejections of jets triggered by a mean-field disk dynamo

We present MHD simulations exploring the launching, acceleration, and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an α{sup 2}Ω mean-field dynamo. The disk magnetization remains on a rather low level, which helps to evolve the simulations for T > 10, 000 dynamical time steps on a domain extending 1500 inner disk radii. We find the magnetic field of the inner disk to be similar to the commonly found open field structure, favoring magneto-centrifugal launching. The outer disk field is highly inclined and predominantly radial. Here, differential rotation induces a strong toroidal component, which plays a key role in outflow launching. These outflows from the outer disk are slower, denser, and less collimated. If the dynamo action is not quenched, magnetic flux is continuously generated, diffuses outward through the disk, and fills the entire disk. We have invented a toy model triggering a time-dependent mean-field dynamo. The duty cycles of this dynamo lead to episodic ejections on similar timescales. When the dynamo is suppressed as the magnetization falls below a critical value, the generation of the outflows and also accretion is inhibited. The general result is that we can steer episodic ejection and large-scale jet knots by a disk-intrinsic dynamo that is time-dependent and regenerates the jet-launching magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulator for SUPO, a Benchmark Aqueous Homogeneous Reactor (AHR)

A simulator has been developed for SUPO (Super Power) an aqueous homogeneous reactor (AHR) that operated at Los Alamos National Laboratory (LANL) from 1951 to 1974. During that period SUPO accumulated approximately 600,000 kWh of operation. It is considered the benchmark for steady-state operation of an AHR. The SUPO simulator was developed using the process that resulted in a simulator for an accelerator-driven subcritical system, which has been previously reported.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Relativistic laser driven electron accelerator using micro-channel plasma targets

We present an experimental demonstration of the efficient acceleration of electrons beyond 60MeV using micro-channel plasma targets. We employed a high-contrast, 2.5 J, 32 fs short pulse laser interacting with a 5 lm inner diameter, 300 lm long microchannel plasma target. The micro-channel was aligned to be collinear with the incident laser pulse, confining the majority of the laser energy within the channel. The measured electron spectrum showed a large increase in the cut-off energy and slope temperature when compared to that from a 2 lm flat Copper target, with the cutoff energy more than doubled and the total energy in electrons >5MeV enhanced by over 10 times. Three-dimensional particle-in-cell simulations confirm efficient direct laser acceleration enabled by the novel structure as the dominant acceleration mechanism for the high energy electrons. The simulations further reveal the guiding effect of the channel that successfully explains preferential acceleration on the laser/channel axis observed in experiments. Finally, systematic simulations provide scalings for the energy and charge of the electron pulses. Our results show that the micro-channel plasma target is a promising electron source for applications such as ion acceleration, Bremsstrahlung X-ray radiation, and THZ generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A UNIFIED COMPUTATIONAL MODEL FOR SOLAR AND STELLAR FLARES

We present a unified computational framework that can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker–Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

79 ASTRONOMY AND ASTROPHYSICS↗

GPU Direct I/O with HDF5

Exascale HPC systems are being designed with accelerators, such as GPUs, to accelerate parts of applications. In machine learning workloads as well as large-scale simulations that use GPUs as accelerators, the CPU (or host) memory is currently used as a buffer for data transfers between GPU (or device) memory and the file system. If the CPU does not need to operate on the data, then this is sub-optimal because it wastes host memory by reserving space for duplicated data. Furthermore, this “bounce buffer” approach wastes CPU cycles spent on transferring data. A new technique, NVIDIA GPUDirect Storage (GDS), can eliminate the need to use the host memory as a bounce buffer. Thereby, it becomes possible to transfer data directly between the device memory and the file system. This direct data path shortens latency by omitting the extra copy and enables higher-bandwidth. To take full advantage of GDS in existing applications, it is necessary to provide support with existing I/O libraries, such as HDF5 and MPI-IO, which are heavily used in applications. In this paper, we describe our effort of integrating GDS with HDF5, the top I/O library at NERSC and at DOE leadership computing facilities. We design and implement this integration using a HDF5 Virtual File Driver (VFD). The GDS VFD provides a file system abstraction to the application that allows HDF5 applications to perform I/O without the need to move data between CPUs and GPUs explicitly. We compare performance of the HDF5 GDS VFD with explicit data movement approaches and demonstrate superior performance with the GDS method.

Ravi, J↗

NASA’s first ground-based Galactic Cosmic Ray Simulator: Enabling a new era in space radiobiology research

With exciting new NASA plans for a sustainable return to the moon, astronauts will once again leave Earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic radiation (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPEs with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include carcinogenesis, central nervous system (CNS) effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including circulatory and heart disease, as well as potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy-ion exposures and the validation of countermeasures in a relevant space environment. Historically, most research on understanding space radiation-induced health risks has been performed using acute exposures of monoenergetic single-ion beams. However, the space radiation environment consists of a wide variety of ion species over a broad energy range. Using the fast beam switching and controls systems technology recently developed at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory, a new era in radiobiological research is possible. NASA has developed the “GCR Simulator” to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (approximately 65%–75%) and helium ions (approximately 10%–20%) with heavier ions (Z ≥ 3) contributing the remainder. The GCR simulator exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the 5 heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader system is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low-energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires approximately 75 minutes. To more closely simulate the low-dose rates found in space, sequential field exposures can be divided into daily fractions over 2 to 6 weeks, with individual beam fractions as low as 0.1 to 0.2 mGy. In the large beam configuration (60 × 60 cm 2 ), 54 special housing cages can accommodate 2 to 3 mice each for an approximately 75 min duration or 15 individually housed rats. On June 15, 2018, the NSRL made a significant achievement by completing the first operational run using the new GCR simulator. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NSRL to accelerate our understanding and mitigation of health risks faced by astronauts. Ultimately, the GCR simulator will require validation across multiple radiogenic risks, endpoints, doses, and dose rates.

59 BASIC BIOLOGICAL SCIENCES↗

Beam-contamination-induced compositional alteration and its neutron-atypical consequences in ion simulation of neutron-induced void swelling

Although accelerator-based ion irradiation has been widely accepted to simulate neutron damage, neutron-atypical features need to be carefully investigated. In this study, we have shown that Coulomb force drag by ion beams can introduce significant amounts of carbon, nitrogen, and oxygen into target materials even under ultra-high vacuum conditions. The resulting compositional and microstructural changes dramatically suppress void swelling. By applying a beam-filtering technique, introduction of vacuum contaminants is greatly minimized and the true swelling resistance of the alloys is revealed and matches neutron behavior closely. These findings are a significant step toward developing standardized procedures for emulating neutron damage.

36 MATERIALS SCIENCE↗

Temperature effects in plasma-based positron acceleration schemes using electron filaments

Preserving the quality of a positron beam in a plasma-based accelerator, where a wakefield suitable for positron transport and acceleration is generated by means of an electron filament, is challenging. This is due to the nature of the wakefields, characterized by focusing fields that vary nonlinearly in the transverse direction, and by accelerating fields that are non-uniform. These fields also change slice-by-slice along the beam. Maintaining a high beam quality is pivotal for application of positron beams in a plasma-based collider. In this paper, we show that an initial background plasma temperature can help mitigate the positron beam quality degradation in plasma-based accelerators that rely on electron filaments. We show that temperature effects broaden the electron filament and smooth radially both the non-linear transverse and the non-uniform longitudinal wakefields. Using warm plasmas opens up new possibilities to improve beam quality in several plasma-based positron acceleration concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗

Do Arctic mixed-phase clouds sometimes dissipate due to insufficient aerosol? Evidence from comparisons between observations and idealized simulations

Abstract. Mixed-phase clouds are ubiquitous in the Arctic. These clouds can persist for days and dissipate in a matter of hours. It is sometimes unknown what causes this sudden dissipation, but aerosol–cloud interactions may be involved. Arctic aerosol concentrations can be low enough to affect cloud formation and structure, and it has been hypothesized that, in some instances, concentrations can drop below some critical value needed to maintain a cloud. We use observations from a Department of Energy ARM site on the northern slope of Alaska at Oliktok Point (OLI), the Arctic Summer Cloud Ocean Study (ASCOS) field campaign in the high Arctic Ocean, and the Integrated Characterisation of Energy, Clouds, Atmospheric state, and Precipitation at Summit – Aerosol Cloud Experiment (ICECAPS-ACE) project at the NSF (National Science Foundation) Summit Station in Greenland (SMT) to identify one case per site where Arctic boundary layer clouds dissipated coincidentally with a decrease in surface aerosol concentrations. These cases are used to initialize idealized large eddy simulations (LESs) in which aerosol concentrations are held constant until, at a specified time, all aerosols are removed instantaneously – effectively creating an extreme case of aerosol-limited dissipation which represents the fastest a cloud could possibly dissipate via this process. These LESs are compared against the observed data to determine whether cases could, potentially, be dissipating due to insufficient aerosol. The OLI case's observed liquid water path (LWP) dissipated faster than its simulation, indicating that other processes are likely the primary drivers of the dissipation. The ASCOS and SMT observed LWP dissipated at similar rates to their respective simulations, suggesting that aerosol-limited dissipation may be occurring in these instances. We also find that the microphysical response to this extreme aerosol forcing depends greatly on the specific case being simulated. Cases with drizzling liquid layers are simulated to dissipate by accelerating precipitation when aerosol is removed while the case with a non-drizzling liquid layer dissipates quickly, possibly glaciating via the Wegener–Bergeron–Findeisen (WBF) process. The non-drizzling case is also more sensitive to ice-nucleating particle (INP) concentrations than the drizzling cases. Overall, the simulations suggest that aerosol-limited cloud dissipation in the Arctic is plausible and that there are at least two microphysical pathways by which aerosol-limited dissipation can occur.

54 ENVIRONMENTAL SCIENCES↗

Dependence of Lamellar Eutectic Growth with Convection on Boundary Conditions and Geometric Confinement: A Phase-Field Lattice-Boltzmann Study

Simulating eutectic growth with convection is challenging because of the enormous computing demand resulting from the required domain size compared with the eutectic scale. In this study, a phase-field lattice-Boltzmann approach was employed to reproduce the eutectic growth with convection (both natural and forced). To enable and accelerate large-scale simulations, a parallel-adaptive mesh refinement algorithm was developed to combine with this numerical approach. The effects of the boundary conditions and the geometric confinement on the eutectic growth with convection were systematically investigated in the Al-Cu eutectic alloy. Results showed that both boundary conditions and the geometric confinement altered the eutectic evolution by affecting the distributions of the solute and velocity ahead of the eutectic solid/liquid interface. As the first attempt to explore the two factors on the eutectic growth with convection, our thorough investigation lays a foundation for modeling the coupled convection–solute eutectic evolution.

36 MATERIALS SCIENCE↗