Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIC 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Efficient Fourier basis particle simulation

The standard particle-in-cell algorithm suffers from grid heating. There exists a gridless alternative which bypasses the deposition step and calculates each Fourier mode of the charge density directly from the particle positions. We show that a gridless method can be computed efficiently through the use of an Unequally Spaced Fast Fourier Transform (USFFT) algorithm. After a spectral field solve, the forces on the particles are calculated via the inverse USFFT (a rapid solution of an approximate linear system). We provide one and two dimensional implementations of this algorithm with an asymptotic runtime of O ( N p + N m D log ⁡ N m D ) for each iteration, identical to the standard PIC algorithm (where N p is the number of particles, N m is the number of Fourier modes, and D is the spatial dimensionality of the problem). We demonstrate superior energy conservation and reduced noise, as well as convergence of the energy conservation at small time steps.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code

As an alternative option to kinetic electrons, the gyrokinetic total-f particle-in-cell (PIC) code XGC1 has been extended to the MHD/fluid type electromagnetic regime by combining gyrokinetic PIC ions with massless drift-fluid electrons. Here, two representative long wavelength modes, shear Alfven waves and resistive tearing modes, are verified in cylindrical and toroidal magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Computational study of plasma sustainability in radio frequency micro-discharges

We apply an implicit particle-in-cell Monte-Carlo (PIC-MC) method to study a radio-frequency argon microdischarge at steady state in the glow discharge limit, in which the microdischarge is sustained by secondary electron emission from the electrodes. The plasma density, electron energy distribution function (EEDF), and electron temperature are calculated in a wide range of operating conditions, including driving voltage, microdischarge gap, and pressure. Also, the effect of gap size scaling (in the range of 50-1000 μm) on the plasma sustaining voltage and peak electron density at atmospheric pressure is examined, which has not been explored before. In our simulations, three different EEDFs, i.e., a so-called three temperature hybrid mode, a two temperature α mode, and a two temperature γ mode distribution, are identified at different gaps and voltages. The maximum sustaining voltage to avoid a transition from the glow mode to an arc is predicted, as well as the minimum sustaining voltage for a steady glow discharge. Our calculations elucidate that secondary electrons play an essential role in sustaining the discharge, and as a result the relationship between breakdown voltage and gap spacing is far away from the Paschen law at atmospheric pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effect of gas properties on the dynamics of the electrical slope asymmetry effect in capacitive plasmas: comparison of Ar, H 2 and CF 4

Tailored voltage excitation waveforms provide an efficient control of the ion energy (through the electrical asymmetry effect) in capacitive plasmas by varying the 'amplitude' asymmetry of the waveform. In this work, the effect of a 'slope' asymmetry of the waveform is investigated by using sawtooth-like waveforms, through which the sheath dynamic can be manipulated. A remarkably different discharge dynamic is found for Ar, H 2 , and CF 4 gases, which is explained by the different dominant electron heating mechanisms and plasma chemistries. In comparison to Argon we find that the electrical asymmetry can even be reversed by using an electronegative gas such as CF 4 . Phase resolved optical emission spectroscopy measurements, probing the spatiotemporal distribution of the excitation rate show excellent agreement with the results of particle-in-cell simulations, confirming the high degree of correlation between the excitation rates with the dominant heating mechanisms in the various gases. It is shown that, depending on the gas used, sawtooth-like voltage waveforms may cause a strong asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Research on effect of emission uniformity on X-band relativistic backward oscillator using conformal PIC code

Explosive emission cathodes (EECs) are adopted in relativistic backward wave oscillators (RBWOs) to generate intense relativistic electron beam. The emission uniformity of the EEC can render saturation of the power generation unstable and the output mode impure. However, the direct measurement of the plasma parameters on the cathode surface is quite difficult and there are very few related numerical study reports about this issue. In this paper, a self-developed three-dimensional conformal fully electromagnetic particle in cell code is used to study the effect of emission uniformity on the X-band RBWO; the electron explosive emission model and the field emission model are both implemented in the same cathode surface, and the local field enhancement factor is also considered in the field emission model. The RBWO with a random nonuniform EEC is thoroughly studied using this code; the simulation results reveal that when the area ratio of cathode surface for electron explosive emission is 80%, the output power is unstable and the output mode is impure. When the annular EEC does not emit electron in the angle range of 30°, the RBWO can also operate normally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling Transverse Space Charge effects in IOTA with pyORBIT

The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Efficient approach to kinetic simulation in the inner magnetically insulated transmission line on Z

This project explores the idea of performing kinetic numerical simulations in the Z inner magnetically insulated transmission line (inner MITL) by reduced physics models such as a guiding center drift kinetic approximation for particles and electrostatic and magnetostatic approximation for the fields. The basic problem explored herein is the generation, formation, and evolution of vortices by electron space charge limited (SCL) emission. The results indicate that for relevant to Z values of peak current and pulse length, these approximations are excellent, while also providing tens to hundreds of times reduction in the computational load. The benefits could be enormous: Implementation of these reduced physics models in present particle-in-cell (PIC) codes could enable them to be routinely used for experimental design while still capturing essential non-thermal (kinetic) physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Modeling laser-driven electron acceleration using WARP with Fourier decomposition

WARP is used with the recent implementation of the Fourier decomposition algorithm to model laser-driven electron acceleration in plasmas. Simulations were carried out to analyze the experimental results obtained on ionization-induced injection in a gas cell. The simulated results are in good agreement with the experimental ones, confirming the ability of the code to take into account the physics of electron injection and reduce calculation time. We present a detailed analysis of the laser propagation, the plasma wave generation and the electron beam dynamics.

43 PARTICLE ACCELERATORS↗
NASA NTRS1995

Vorticity and divergence in the solar photosphere

We have studied an outstanding sequence of continuum images of the solar granulation from Pic du Midi Observatory. We have calculated the horizontal vector flow field using a correlation tracking algorithm, and from this determined three scalar field: the vertical component of the curl; the horizontal divergence; and the horizontal flow speed. The divergence field has substantially longer coherence time and more power than does the curl field. Statistically, curl is better correlated with regions of negative divergence - that is, the vertical vorticity is higher in downflow regions, suggesting excess vorticity in intergranular lanes. The average value of the divergence is largest (i.e., outflow is largest) where the horizontal speed is large; we associate these regions with exploding granules. A numerical simulation of general convection also shows similar statistical differences between curl and divergence. Some individual small bright points in the granulation pattern show large local vorticities.

Wang, YI↗
DOE OSTI2013

Role of positive ions on the surface production of negative ions in a fusion plasma reactor type negative ion source—Insights from a three dimensional particle-in-cell Monte Carlo collisions model

Results from a 3D self-consistent Particle-In-Cell Monte Carlo Collisions (PIC MCC) model of a high power fusion-type negative ion source are presented for the first time. The model is used to calculate the plasma characteristics of the ITER prototype BATMAN ion source developed in Garching. Special emphasis is put on the production of negative ions on the plasma grid surface. The question of the relative roles of the impact of neutral hydrogen atoms and positive ions on the cesiated grid surface has attracted much attention recently and the 3D PIC MCC model is used to address this question. The results show that the production of negative ions by positive ion impact on the plasma grid is small with respect to the production by atomic hydrogen or deuterium bombardment (less than 10%)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Study of the negative ion extraction mechanism from a double-ion plasma in negative ion sources

We have developed a 2D3V-PIC model of the extraction region, aiming to clarify the basic extraction mechanism of H{sup −} ions from the double-ion plasma in H{sup −} negative ion sources. The result shows the same tendency of the H{sup −} ion density n{sub H{sup −}} as that observed in the experiments, i.e.,n{sub H{sup −}} in the upstream region away from the plasma meniscus (H{sup −} emitting surface) has been reduced by applying the extraction voltage. At the same time, relatively slow temporal oscillation of the electric potential compared with the electron plasma frequency has been observed in the extraction region. Results of the systematic study using a 1D3V-PIC model with the uniform magnetic field confirm the result that the electrostatic oscillation is identified to be lower hybrid wave. The effect of this oscillation on the H{sup −} transport will be studied in the future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Extreme Scale Plasma Turbulence Simulations on Top Supercomputers Worldwide

The goal of the extreme scale plasma turbulence studies described in this paper is to expedite the delivery of reliable predictions on confinement physics in large magnetic fusion systems by using world-class supercomputers to carry out simulations with unprecedented resolution and temporal duration. This has involved architecture-dependent optimizations of performance scaling and addressing code portability and energy issues, with the metrics for multi-platform comparisons being 'time-to-solution' and 'energy-to-solution'. Realistic results addressing how confinement losses caused by plasma turbulence scale from present-day devices to the much larger $25 billion international ITER fusion facility have been enabled by innovative advances in the GTC-P code including (i) implementation of one-sided communication from MPI 3.0 standard; (ii) creative optimization techniques on Xeon Phi processors; and (iii) development of a novel performance model for the key kernels of the PIC code. Our results show that modeling data movement is sufficient to predict performance on modern supercomputer platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

44th Annual Anomalous Absorption Conference

Conference Grant Report July 14, 2015 Submitted to the U. S. Department of Energy Attn: Dr. Sean Finnegan By the University of California, San Diego 9500 Gilman Drive La Jolla, California 92093 On behalf of the 44th Annual Anomalous Absorption Conference 8-13 June 2014, in Estes Park, Colorado Support Requested: $10,100 Amount expended: $3,216.14 Performance Period: 1 March 20 14 to 28 February 20 15 Principal Investigator Dr. Farhat Beg Center for Energy Research University of California, San Diego 9500 Gilman Drive La Jolla, California 92093-0417 858-822-1266 (telephone) 858-534-4543 (fax) fbeg@ucsd.edu Administrative Point of Contact: Brandi Pate, 858-534-0851, blpate®ucsd.edu I. Background The forty-fourth Anomalous Absorption Conference was held in Estes Park, Colorado from June 5-8, 2014 (aac2014.ucsd.edu). The first Anomalous Absorption Conference was held in 1971 to assemble experts in the poorly understood area of laser-plasma absorption. The goal of that conference was to address the anomalously large laser absorption seen in plasma experiments with respect to the laser absorption predicted by linear plasma theory. Great progress in this research area has been made in the decades since that first meeting, due in part to the scientific interactions that have occurred annually at this conference. Specifically, this includes the development of nonlinear laser-plasma theory and the simulation of laser interactions with plasmas. Each summer since that first meeting, this week-long conference has been held at unique locations in North America as a scientific forum for intense scientific exchanges relevant to the interaction of laser radiation with plasmas. Responsibility for organizing the conference has traditional rotated each year between the major Inertial Confinement Fusion (ICF) laboratories and universities including LANL, LLNL, LLE, UCLA UC Davis and NRL. As the conference has matured over the past four decades, its technical footprint has expanded beyond ICF-related laser-plasma interactions to encompass closely related technical areas including laser particle acceleration, high-intensity laser effects, short­ pulse laser interactions, PIC and Vlasov/rad-hydro modeling, inertial and magnetic fusion plasmas, advanced plasma diagnostics, alternate ignition schemes, EOS/transport/opacity, and this year, x­ ray free-electron lasers and their applications. The conference continues to be a showcase for the presentation and discussion of the latest developments in these areas. II. Meeting Report The conference was extremely successful with more than one hundred participants. There were ninety-nine (99) abstracts submitted. There were forty-four (44) presentations including eleven (11) invited talks. The following topics were covered: a) Radiation Hydrodynamics b) Implosion Plasma Kinetic Effects c) Alternate Ignition Schemes d) Astrophysical Phenomena e) Opacity/Transport/EOS f) High Power Lasers and Facilities g) High-Intensity Laser-Matter Interactions h) Hydrodynamics and Hydro-instabilities i) Hot Dense Plasma Atomic Processes j) High Energy Density Physics k) Laser Particle Acceleration Physics l) Advanced Plasma Diagnostics m) Advanced light sources and applications Despite significant advertising, there were two students who applied for the travel grants: Charlie Jarrott and Joohwan Kim. The total funds expended were $3,216.14.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

THE GEMINI NICI PLANET-FINDING CAMPAIGN: THE FREQUENCY OF GIANT PLANETS AROUND YOUNG B AND A STARS

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (≈1.5-2.5 M{sub ☉}) conducted to date and includes the planet hosts β Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58{sup +21}{sub -20} M{sub Jup} and 55{sup +20}{sub -19} M{sub Jup}, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M{sub ☉} stars can have giant planets greater than 4 M{sub Jup} between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M{sub Jup} between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M{sub Jup}, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2019

A Detect and Avoid System in the Context of Multiple-Unmanned Aircraft Systems Operations

NASA's Unmanned Aircraft Systems Integration into the National Airspace System (UAS in the NAS) project examines the technical barriers associated with the operation of UAS in civil airspace. For UAS, the removal of the pilot from onboard the aircraft has eliminated the ability of the ground-based pilot in command (PIC) to use out-the-window visual information to make judgements about a potential threat of a loss of well clear with another aircraft. NASA's Phase 1 research supported the development of a Detect and Avoid (DAA) system that supports the ground-based pilot's ability to detect potential traffic conflicts and determine a resolution maneuver, but existing display/alerting requirements did not account for multiple UAS control (1:N). Demands for increased scalability of UAS in the NAS operations are expected to create a need for simultaneous control of UAs, and thus, a new DAA HMI design will likely be necessary. Previous research, however, has found performance degradations as the number of vehicles under operator control has increased. The purpose of the current human-in-the-loop (HITL) simulation was to examine the viability of 1:N operations with the Phase 1 DAA alerting and guidance. Sixteen UAS pilots flew three scenarios with varying number of UAs under their control (1:1, 1:3, 1:5). In addition to their supervisory and sensor mission responsibilities, pilots were to utilize the DAA system to remain DAA well clear (DWC) during scripted conflicts of mixed severity. Measured response times, separation performance, mission task data, and subjective feedback were collected to assess how the multi-UAS control configuration impacted pilots' ability to maintain DAA well clear and perform the mission tasks. Overall, the DAA system proved surprisingly adaptive to multi-UAS control for preventing losses of DAA well clear (LoDWC). The findings suggest that, while multi-UAS operators are able to maintain safe separation (DWC) from other traffic, their ability to efficiently perform missions drastically decreases with their number of controlled vehicles. Pilot feedback indicated that, for this context, the use of automation support tools for completing and managing mission tasks would be appropriate and desired, especially for ensuring efficient use of assets. Finally, human-machine interface (HMI) design considerations for multi-UAS operations are discussed.

Monk, Kevin J.↗
DOE OSTI2017

Hybrid-PIC modeling of laser-plasma interactions and hot electron generation in gold hohlraum walls

The walls of the hohlraum used in experiments at the national ignition facility are heated by laser beams with intensities ∼1015 W/cm2, a wavelength of ∼1/3 μm, and pulse lengths on the order of a ns, with collisional absorption believed to be the primary heating mechanism. X-rays generated by the hot ablated plasma at the gold walls are then used to implode a target in the hohlraum interior. In addition to the collisional absorption of laser energy at the walls, non-linear laser-plasma interactions (LPI), such as stimulated Raman scattering and two plasmon decay, are believed to generate a population of supra-thermal electrons which, if present in the hohlraum, can have a deleterious effect on target implosion. We describe results of hohlraum modeling using a hybrid particle-in-cell code. To enable this work, new particle-based algorithms for a multiple-ion magneto-hydrodynamic (MHD) treatment, and a particle-based ray-tracing model were developed. The use of such hybrid methods relaxes the requirement to resolve the laser wavelength, and allows for relatively large-scale hohlraum simulations with a reasonable number of cells. But the non-linear effects which are believed to be the cause of hot electron generation can only be captured by fully kinetic simulations with good resolution of the laser wavelength. For this reason, we employ a two-tiered approach to hohlraum modeling. Large-scale simulations of the collisional absorption process can be conducted using the fast quasi-neutral MHD algorithm with fluid particle species. From these simulations, we can observe the time evolution of the hohlraum walls and characterize the density and temperature profiles. From these results, we can transition to smaller-scale highly resolved simulations using traditional kinetic particle-in-cell methods, from which we can fully model all of the non-linear laser-plasma interactions, as well as assess the details of the electron distribution function. We find that vacuum hohlraums should be stable to both two plasmon decay and stimulated Raman scattering instabilities for intensities ≤1015 W/cm2. In gas-filled hohlraums, shocks may be induced in the blowoff gold plasma, which leads to more complex density and temperatures profiles. The resulting effect on LPI stability depends strongly on the details of the profile, and it is possible for the gas-filled hohlraum to become unstable to two plasmon decay at 1015 W/cm2 if the quarter-critical surface reaches temperatures exceeding 1 keV.

Thoma, C.↗
DOE OSTI2022

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.

97 MATHEMATICS AND COMPUTING↗