Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle dispersion”

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 1,135 records · Page 63

Interplanetary Magnetic Field Guiding Relativistic Particles

The origin and the propagation of relativistic solar particles (0.5 to few Ge V) in the interplanetary medium remains a debated topic. These relativistic particles, detected at the Earth by neutron monitors have been previously accelerated close to the Sun and are guided by the interplanetary magnetic field (IMF) lines, connecting the acceleration site and the Earth. Usually, the nominal Parker spiral is considered for ensuring the magnetic connection to the Earth. However, in most GLEs the IMF is highly disturbed, and the active regions associated to the GLEs are not always located close to the solar footprint of the nominal Parker spiral. A possible explanation is that relativistic particles are propagating in transient magnetic structures, such as Interplanetary Coronal Mass Ejections (ICMEs). In order to check this interpretation, we studied in detail the interplanetary medium where the particles propagate for 10 GLEs of the last solar cycle. Using the magnetic field and the plasma parameter measurements (ACE/MAG and ACE/SWEPAM), we found widely different IMF configurations. In an independent approach we develop and apply an improved method of the velocity dispersion analysis to energetic protons measured by SoHO/ERNE. We determined the effective path length and the solar release time of protons from these data and also combined them with the neutron monitor data. We found that in most of the GLEs, protons propagate in transient magnetic structures. Moreover, the comparison between the interplanetary magnetic structure and the interplanetary length suggest that the timing of particle arrival at Earth is dominantly determined by the type of IMF in which high energetic particles are propagating. Finally we find that these energetic protons are not significantly scattered during their transport to Earth.

Masson, S.↗

Electron Microscopic Examination of an Irradiated TRISO Coated Particle of AGR-2 Experiment: AGR2-222-RS019

The transmission electron microscopic (TEM) examination of an irradiated tristructural isotropic (TRISO) coated particle from advanced gas reactor (AGR) 2 fuel Compact 2-2-2 is presented in this report. Compact 2-2-2 refers to the compact in Capsule 2 at Level 2 of Stack 2. The examination on AGR2 222 RS019 is focused on the interface between the silicon carbide (SiC) to inner pyrolytic carbon (IPyC) layer, central part of SiC, and the interface between SiC and outer pyrolytic carbon layer. The SiC grain boundary distribution and associated fission product precipitates were characterized by scanning transmission electron microscopy (STEM) energy dispersive spectroscopy (EDS), and precession electron diffraction (PED) in the transmission electron microscope. Two separate locations on SiC layer were selected for on the basis of the intact of the buffer layer in the irradiated TRISO particle. The key observation in this particle is that the compositional nature of grain boundary precipitates is more complex near IPyC/SiC layer where Pd, Pd-Ag, Pd-U containing precipitates can be observed. However, the compositions of precipitate towards outer most region of SiC mainly contain Pd. Ag has been observed in the IPyC/SiC layer in one of the two locations. Comparison of the results to previously irradiated and analyzed, Variant 3, safety tested AGR-1 TRISO particles (AGR1-433-001 and AGR1 222-RS004) showed the grain boundary distributions have some noticeable differences. The differences in coincident site lattice grain boundary and misorientation angles are possibly inherent to the fabrication process. Neutron irradiation induced void distribution has been found to be nonuniform. A high concentration of smaller voids is found at stacking faults compared to those in the SiC matrix. Larger void sizes are observed in the SiC location that is near to the broken buffer layer. This void size variation with integrity of buffer layer can potentially affect the fission product retention. However, the average void size continuously decreases along the radius of SiC in both locations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigation of irradiation damage and heat deposition: a comparative analysis for HEU-to-LEU conversion in HFIR

The planned conversion of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel requires detailed evaluation of experiment-relevant parameters to ensure continued performance for materials testing and isotope production. Here, this study presents the first comprehensive assessment of displacements per atom (dpa) and heat deposition rates in target materials within the HFIR flux trap with both HEU and candidate LEU core configurations. Seven analyses were conducted to evaluate key performance metrics, including fast neutron flux distribution, cross section response functions, cross section data, and local dpa and heat deposition rates using mesh- and cell-based tallies. Simulations employed Shift, Monte Carlo N-Particle (MCNP), and the HIFR Controller (HFIRCON) tool suite for high-fidelity transport and depletion modeling. The LEU designs—using U 3 Si 2 -Al dispersion fuel and operating at 95 MW—were compared to the current 85 MW HEU configuration. Results show that while the candidate LEU cores exhibit higher dpa rates due to a harder spectrum and extended cycle lengths, they also demonstrate reduced heat deposition rates in irradiation experiments, primarily due to increased gamma self-shielding from higher 238 U content in the core. These findings confirm that LEU conversion can maintain HFIR’s materials irradiation capabilities but may require redesigning existing experimental hardware.

HEU↗

Wildfire towers drive firebrand lofting: insights from coupled fire-atmosphere model simulations

Wildfire behavior is shaped by complex fire dynamics, with firebrands playing a critical role in spot fire ignition and fire spread. While previous studies have explored firebrand generation and transport, the specific role of towers and troughs from wildland fires in the lofting of firebrands remains unquantified. This study addresses that gap by using physics-based coupled fire-atmosphere model simulations to examine how wildfire towers (updrafts) and troughs (downdrafts) influence firebrand lofting. Our results show that the majority of firebrands (78.85%) are lofted from towers, where strong updrafts drive long-range transport. In contrast, only 21.15% of firebrands are lofted within troughs, where downdrafts cause most firebrands to fall near the fireline. We also find that firebrand size significantly influences lofting behavior, with smaller particles (1 mm radius) exhibiting the strongest correlation with updraft intensity. These findings highlight the dominant role of wildfire towers in promoting long-distance firebrand dispersal—an essential factor in rapid wildfire growth and wildland-urban interface (WUI) fire risks. By quantifying the relationship between firebrand lofting and fire-induced atmospheric features, this study provides critical insights to improve spot fire modeling, support mitigation planning, and enhance firefighter and WUI community safety in spot fire-prone regions.

54 ENVIRONMENTAL SCIENCES↗

Dispersive optical model analysis of 208 Pb generating a neutron-skin prediction beyond the mean field

In this work, a nonlocal dispersive optical model analysis is carried out for neutrons and protons in 208 Pb. Elastic-scattering angular distributions, total and reaction cross sections, single-particle energies, neutron and proton numbers, the charge distribution, and the binding energy are fitted to extract the neutron and proton self-energies both above and below the Fermi energy. From the single-particle propagator derived from these self-energies, we determine the charge and matter distributions in 208 Pb. The predicted spectroscopic factors are consistent with results from the ($\textit{e,e}'\textit{p}$) reaction and inelastic-electron-scattering data to very high-spin states. Sensible results for the high-momentum content of neutrons and protons are obtained, with protons appearing more correlated, in agreement with experiment and ab initio calculations of asymmetric matter. A neutron skin of 0.25 ± 0.05 fm is deduced. An analysis of several nuclei leads to the conclusion that finite-size effects play a nonnegligible role in the formation of the neutron skin in finite nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding Platinum Ionomer Interface Properties of Polymer Electrolyte Fuel Cells

A well-designed cathode catalyst layer with optimal ionomer distribution is critical to minimizing amount of Platinum (Pt) content in polymer electrolyte fuel cells (PEFCs). The impact of Pt loading, ionomer content and carbon support types on the catalyst/ionomer interface were also investigated at dry and wet conditions. Higher Pt loadings resulted in higher double layer capacity (C dl ) and similar electrochemical surface area (ECSA) due to well dispersed ionic phase material. Higher ionomer content resulted in higher ionic conductivity but also showed similar SO 3 − group coverage. High surface area (HSA) carbon support had larger ECSA and C dl at both dry and wet conditions, as less agglomerated Pt was well dispersed in the meso-pores of the support. Lower SO 3 − group coverages were observed for HSA carbon than for Vulcan carbon due to Pt particles being buried within the porous HSA carbon support. The effect of cell conditioning and voltage recovery on the PEFC cathode catalyst layer was shown to have minimal impact on SO 3 − group coverage despite a decrease in C dl and ECSA due to the size increase of Pt particles. At dry condition, a significant increases in SO 3 − group coverage were observed for all MEAs due to higher adsorptivity of ionomer in dry conditions.

25 ENERGY STORAGE↗

Methods—Kintsugi Imaging of Battery Electrodes: Distinguishing Pores from the Carbon Binder Domain using Pt Deposition

The mesostructure of porous electrodes used in lithium-ion batteries strongly influences cell performance. Accurate imaging of the distribution of phases in these electrodes would allow this relationship to be better understood through simulation. However, imaging the nanoscale features in these components is challenging. While scanning electron microscopy is able to achieve the required resolution, it has well established difficulties imaging porous media. This is because the flat imaging planes prepared using focused ion beam milling will intersect with the pores, which makes the images hard to interpret as the inside walls of the pores are observed. It is common to infiltrate porous media with resin prior to imaging to help resolve this issue, but both the nanoscale porosity and the chemical similarity of the resins to the battery materials undermine the utility of this approach for most electrodes. In this study, a technique is demonstrated which uses in situ infiltration of platinum to fill the pores and thus enhance their contrast during imaging. Reminiscent of the Japanese art of repairing cracked ceramics with precious metals, this technique is referred to as the kintsugi method. The images resulting from applying this technique to a conventional porous cathode are presented and then segmented using a multi-channel convolutional method. We show that while some cracks in active material particles were empty, others appear to be filled (perhaps with the carbon binder phase), which will have implications for the rate performance of the cell. Energy dispersive X-ray spectroscopy was used to validate the distribution of phases resulting from image analysis, which also suggested a graded distribution of the binder relative to the carbon additive. The equipment required to use the kintsugi method is commonly available in major research facilities and so we hope that this method will be rapidly adopted to improve the imaging of electrode materials and porous media in general.

25 ENERGY STORAGE↗

The oblique behavior of low-frequency electromagnetic waves excited by newborn cometary ions

The free energy in oxygen or hydrogen ions freshly created in the solar wind stimulates low-frequency electromagnetic waves whose growth does not always maximize at parallel propagation. Exploration of the wave vector plane discloses the frequent occurrence of islets of oblique growth unconnected to the unstable parallel modes. Contour plots of the growth rate, real frequency, polarization, and magnetic compression characterize the oblique wave behavior for large values of the initial pitch angle of the cometary particles. Although wave-particle (Landau and cyclotron) resonances feed most of the surveyed oblique instabilities, some are seemingly fluidlike. The results, obtained from the numerical solution of the kinetic dispersion and wave equations, imply that newborn ions can easily excite significant oblique hydromagnetic wave activity. Cometary environments provide the adopted plasma model, but the study is helpful in the interpretation of other low-frequency wave observations in space.

Brinca, Armando L.↗

Compositional analysis of projectile residues on LDEF instrument AO187-1

Impact craters greater than 30 microns and associated projectile residues were analyzed by Scanning Electron Microscopy (SEM) and Energy Dispersive X-ray Analysis (EDXA). Objectives were to analyze a statistically significant number of projectiles to evaluate their chemical variability and possible clustering into discrete particle types. Bay A11 exposed six collector surfaces of anodized 1100-T4 (greater than 99 percent pure) aluminum sheets, 0.32 cm thick, yielding an exposed surface area of 1.1 sq. m. Four of the six panels have been retained at JSC, and were optically scanned, one (A11E00E) was prepared for SEM/EDX analysis. Bay A03 was occupied by high purity (99.99 percent) gold sheets, 0.5 mm thick, yielding an exposed surface area of 0.85 sq. m. Sample processing included the optical scanning (6X), labeling, and dislodging (by a punch-die device) of each individual impact greater than 75 microns for the aluminum and 30 microns for the gold. The 209 craters were dislodged form A11E00E, having crater diameters up to 3500 microns. Optical examination of the gold surfaces detected 238 craters, 198 of which were retained at JSC and analyzed via SEM/EDX. The analytical procedures included maximizing the geometric efficiency (take-off angles), using relatively long count times (500-1000 sec) and sufficiently high accelerating currents (25-30Kev). Despite diligent examination, a large number of craters did not exhibit measurable signals above background. Detectable resides were classified as either micrometeoritic or as man-made debris.

Bernhard, Ronald P.↗

Particulate multi-phase flowfield analysis for advanced solid rocket motor

Particulate multi-phase flowfield with chemical reaction for a 2D advanced solid rocket motor (ASRM) is analyzed using the finite difference Navier-Stokes (FDNS) code. The flowfield in the aft dome cavity of the ASRM is examined and its significant impact on the motor operation and performance is demonstrated. Chemical reaction analysis is performed for H2O, O2, H2, O, H, OH, CO, CO2, Cl, Cl2, HCl, and N2. The turbulent dispersion effect is calculated with the Monte Carlo method. Result show that a recirculation zone exists at the entry of the aft-dome cavity. The particle impingement could cause the erosion and damage nozzle wall. Accumulating in the impingement area the particles change the wall shape and affect the motor performance.

Liaw, Paul↗

Plasma Sheet Source and Loss Processes

Data from the TIMAS ion mass spectrometer on the Polar satellite, covering 15 ev/e to 33 keV/e in energy and essentially 4(pi) in view angles, are used to investigate the properties of earthward (sunward) field-aligned flows of ions, especially protons, in the plasma sheet-lobe transition region near local midnight. A total of 142 crossings of this region are analyzed at 12-sec time resolution, all in the northern hemisphere, at R(SM) approx. 4 - 7 R(sub E), and most (106) in the poleward (sunward) direction. Earthward proton flows are prominent in this transition region (greater than 50% of the time), typically appearing as sudden "blasts" with the most energetic protons (approx. 33 keV) arriving first with weak flux, followed by protons of decreasing energy and increasing flux until either: (1) a new "blast" appears, (2) the flux ends at a sharp boundary, or (3) the flux fades away within a few minutes as the mean energy drops to a few keV. Frequent step-like changes (less than 12 sec) of the flux suggest that perpendicular gradients on the scale of proton gyroradii are common. Peak flux is similar to central plasma sheet proton flux (10(exp 5) - 10(exp 6)/[cq cm sr sec keV/e] and usually occurs at E approx. 4 - 12 keV. Only the initial phase of each "blast" (approx. 1 min) displays pronounced field-alignment of the proton velocity distribution, consistent with the time-of-flight separation of a more or less isotropic source distribution with df/d(nu) less than 0. The dispersive signatures are often consistent with a source at R(SM) less than or equal to 30 R(sub E). No systematic latitudinal velocity dispersion is found, implying that the equatorial plasma source is itself convecting. In short, the proton "blasts" appear as sudden local expansions of central plasma sheet particles along reconfigured ("dipolarized") magnetic field lines.

Lennartsson, O. W.↗

(abstract) Cometary Particles as a Tracer of Jupiter's Stratospheric Circulation

The impact of fragments of comet Shoemaker-Levy 9 on Jupiter's atmosphere in July 1994 may provide an unprecedented opportunity to study Jupiter's stratospheric circulation. Recent calculations by Z. Sekanina predict that much of the comet material will be deposited in Jupiter's stratosphere. If so, and if the material is deposited in a confined region (10 000 km or less, horizontally) we can expect a situation analogous to an El Chichon or Pinatubo event for the terrestrial stratosphere. Initially the volatile material will be vaporized and will rapidly recondense. The large ice crystals and dust particles will rain out and be lost to the troposphere. The cloud of small particles which remain may have settling times of more than a year. These submicron to micron particles would probably be easily seen in methane filter images in the near-IR, and possibly in the ultraviolet. An observational program to monitor the dispersal of this cloud or clouds would reveal much about the nature of the circulation. Some predictions about the meridional evolution of the clouds can be made already, based on the meridional circulation model of West et al. unless the impact itself significantly disrupts the annual average circulation well after the initial transients die away.

comet Shoemaker-Levy 9 impact Jupiter stratosphere↗

DED Additively Manufactured HEAs Optimized via Parametric Study of Functionally Graded Materials

Complex concentrated alloys (CCAs) are a system of alloys composition of nearly equiatomic elements, creating a state of high configurational entropy in the material. These alloys are of interest for various industrial applications due to their high strength and hardness, wide operational temperature range, creep and diffusion resistance, and radiation resistance. In the current work CoCrFeNi-base CCAs are synthesize by combining three commonly industrial alloys: IN718, SS316L, and 70Co30Cr (commercially known as Stellite 21), in-situ utilizing directed energy deposition (DED) additive manufacturing (AM). One bulk specimen is fabricated with a 1:1:1 ratio of IN718, SS316L, and 70Co30Cr where IN718 is a source of nickel and chrome, SS316 is a source of iron and chrome, and 70Co30Cr is a source of cobalt and supplemental chrome. Additionally, a functionally graded material (FGM) of the identified CCA system from SS316 was fabricated with 10 layers, where the initial five layers are the base material of SS316L and graded by steps of 20% to transition to the CCA system. Overall build quality, microstructure, compositional distribution, phases analysis, and microhardness were studied by laser optical microscopy, scanning election microscopy (SEM), energy dispersive spectroscopy (EDS), X-ray diffraction analysis (XRD), and Vickers microhardness techniques. Overall build quality showed a successful bulk build with acceptable density, however unmelted particles are observed in the as-built alloy. Compositional results show that the system meets the Boltzmann’s hypothesis definition of an HEA, and that an FCC CCA was produced with increased hardness, small grain size, and superior hardness properties to the SS316L substrate. While initial XRD results show predominately FCC structure, EDS analysis of interdendritic regions show the existence of Laves phases in the material. The fabrication of an FGM CCA allowed for potential screening of additional alloy mixes which may be of interest to explore in future research, and the resulting composition is compared to predicted values. This novel method of CCA fabrication results in substantial cost savings of these CCA systems over traditional methods.

36 MATERIALS SCIENCE↗

Unusually Long Path Length for a Nearly Scatter-Free Solar Particle Event Observed By Solar Orbiter at 0.43 au

Context. After their acceleration and release at the Sun, solar energetic particles (SEPs) are injected into the interplanetary medium and are bound to the interplanetary magnetic field (IMF) by the Lorentz force. The expansion of the IMF close to the Sun focuses the particle pitch-angle distribution, and scattering counteracts this focusing. Solar Orbiter observed an unusual solar particle event on 9 April 2022 when it was at 0.43 astronomical units (au) from the Sun. Aims. We show that the inferred IMF along which the SEPs traveled was about three times longer than the nominal length of the Parker spiral and provide an explanation for this apparently long path. Methods. We used velocity dispersion analysis (VDA) information to infer the spiral length along which the electrons and ions traveled and infer their solar release times and arrival direction. Results. The path length inferred from VDA is approximately three times longer than the nominal Parker spiral. Nevertheless, the pitch-angle distribution of the particles of this event is highly anisotropic, and the electrons and ions appear to be streaming along the same IMF structures. The angular width of the streaming population is estimated to be approximately 30 degrees. The highly anisotropic ion beam was observed for more than 12 h. This may be due to the low level of fluctuations in the IMF, which in turn is very probably due to this event being inside an interplanetary coronal mass ejection The slow and small rotation in the IMF suggests a flux-rope structure. Small flux dropouts are associated with very small changes in pitch angle, which may be explained by different flux tubes connecting to different locations in the flare region. Conclusions. The unusually long path length along which the electrons and ions have propagated virtually scatter-free together with the short-term flux dropouts offer excellent opportunities to study the transport of SEPs within interplanetary structures. The 9 April 2022 solar particle event offers an especially rich number of unique observations that can be used to limit SEP transport models.

Robert F. Wimmer-Schweingruber↗

The measurement of trace elements in interplanetary dust and cometary particles by ultra-high sensitivity INAA

Today the major elemental composition of interplanetary dust particles (IDPs) is routinely determined in many laboratories. These and mineralogical studies have revealed the presence of at least two major types of IDPs, chondritic and refractory. Preliminary results of a successful attempt to determine abundances of a large suite of trace elements from both chondritic and refractory IDPs are reported. The analytical procedure can be used in the grain-by-grain analysis of returned cometary samples. Chondritic and refractory IDPs are characterized by standard scanning electron microscopy and energy dispersive x ray spectroscopy (SEM-EDX) techniques. With this system, detection limits for many elements are well below picogram levels, and some approach femtogram levels. This technique is non-destructive, although some sample handling is required, so particles can be analyzed by other techniques after instrument neutron activation analysis (INAA) is completed. Data is presently being reduced from the analyses of 7 IDPs. These are U2015E10, U2015F1, W7029-A2, W7029-A3, W7013A8, LACl (all chondritic) and 705 (refractory). So far, 17 different major and trace elements were detected and measured in these particles, including rare earths and some very volatile elements (Br and Zn).

Zolensky, M. E.↗

Large‐Area, High‐Numerical‐Aperture, Freeform Metasurfaces

Abstract Nanophotonic devices are optical platforms capable of unprecedented wavefront control. To push the limits of experimental device performance, scalable design methodologies that combine the simplicity and fabricability of conventional design paradigms with the extended capabilities of freeform optimization are required. A novel gradient‐based design framework for large‐area freeform metasurfaces is introduced in which nonlocal interactions between simply shaped nanostructures, placed on an irregular lattice, are tailored to produce high‐order hybridized modes that support customizable large‐angle scattering profiles. Utilizing this approach, multifunctional super‐dispersive metalenses are designed and experimentally demonstrated. The approach to high‐numerical‐aperture radial metalenses capable of diffraction limited focusing and the generation of donut‐shaped point spread functions is also extended. It is anticipated that these concepts will have utility in super‐resolution microscopy, particle trapping, additive manufacturing, and metrology applications that require ultra‐high numerical apertures.

Optics↗

Microstructural and Micro-Chemical Evolutions in the Irradiated UCO Fuel Kernels of AGR-1 and AGR-2 TRISO Fuel Particles

AGR-1 and AGR-2 TRISO fuel particles were fabricated with slightly different fuel kernel chemical compositions, modified fabrication processes, different fuel kernel diameters, and changed 235U enrichments. To correlate those differences with the fuel kernel responses to neutron irradiations in terms of irradiated fuel microstructure, fission products chemical and physical states, and fission gas bubble evolutions, extensive microstructural and analytical characterizations were conducted. The studies used a state of art transmission electron microscopy (TEM) equipped with Energy-dispersive X-ray spectroscopy (EDS) of four silicon solid-state detectors which have super sensitivity and fast speed. The TEM specimens were prepared from selected AGR-1 and AGR-2 irradiated fuel kernels exposed to safety testing after irradiation. The particles were chosen to represent a representative irradiation conditions with a fuel burnup within the range from 10.8 to 18.6% FIMA, and the time-average volume-average temperatures vary from 1070 to 1287°C. The 235U enrichment was 19.74 wt.% for the AGR1 fuel kernels and 14.03 wt.% for the AGR-2 fuel kernels. The TEM results show that there were significant microstructural reconstructions in the irradiated fuel kernels for both the AGR-1 and AGR-2 fuels. There are four major phases including fuel matrix of UO2 and UC, U2RuC2, and UMoC2 in the irradiated AGR2 fuel kernel. Zr and Nb form solid solution in the UC phase. UMoC2 phase often shows a detectable concentration of Tc. Pd was found to mainly locate in the buffer layer or to be associated with fission gas bubble within the UMoC2 phase. The EDS maps qualitatively show that the rare-earth fission products (Nb, et al.) preferentially reside in the UO2 phase. In contrast, in the irradiated AGR1 fuel kernel, no U2RuC2 or UMoC2 precipitates were positively identified. Instead, there is a high number of rod-shape precipitates enriched with Ru, Tc, Rh, and Pd observed in the fuel kernel center and edge zone. The difference of microstructural and micro-chemical evolutions in irradiated fuel kernels between the AGR-1 and AGR-2 TRISO fuel particle may result from a combined factor of irradiation temperature, fuel geometry and chemical composition. However, the irradiation temperature probably play a more deterministic role. Limited electron energy loss spectroscopy (EELS) characterizations on the AGR2 fuel kernel show that there is nearly no carbon in the UO2 phase while a small fraction of oxygen was detected in the UC/UMoC2 phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Verification of a fully implicit particle-in-cell method for the <!--${MathJax: TeX-AMS-MML_HTMLorMML}--> v &#x2225; -formalism of electromagnetic gyrokinetics in the XGC code

A fully implicit particle-in-cell method for handling the v ∥ -formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the v ∥ -formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the p ∥ -formalism. The v ∥ -formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗