Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fusion Simulation Project”

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 451 records · Page 25

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and numerical investigation of non-neutral complex plasmas

A plasma of particles with the same sign of charge, can be easily confined under ultra-high vacuum conditions in Penning-Malmberg traps, where the time evolution of the system is monitored for very long times by means of electrostatic and optical diagnostic systems. Complex (dusty) plasmas are ionized gases that contain a distribution of micrometer-sized particles with a surface charge of the order of a few thousand electron charges. The interplay between a wide range of scales in time and space gives rise to new characteristic physical phenomena. Laboratory complex plasmas generally satisfy a global (quasi-)neutrality condition. A different concept is represented by a non-neutral complex plasma. To investigate the dynamics of this system, we are currently developing the DuEl (Dust-Electron) device, where negatively charged dust particles will be present together with a population of electrons. The experimental set-up will include a dust injection system and a Penning-Malmberg trap for the confinement of the dust-contaminated electron plasma. We describe here the main physical aims of the project and the present design of the apparatus. To support the experimental project, we have been developing a specifically tailored two-dimensional 'hybrid' Particle-In-Cell code. Using polar cylindrical coordinates, the code aims to investigate the transverse dynamics of a magnetized electron plasma contaminated by a massive, charged species. A mass-less fluid approximation for the electron population is exploited, while the dust component is treated with a kinetic description, also including the gravitational force. The preliminary results of systematic studies on the effects of heavy (magnetized or non-magnetized) dust grains on the equilibrium and stability properties of the electron fluid are presented. The implementation of other characteristic phenomena of interest, e.g. residual gas friction and dust charge fluctuations, is also under development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Report for "SciDAC Center for Simulation of Wave-Plasma Interactions (CSWPI)"

This project has provided for the participation of the PI (Dr. David Smithe) and other Tech-X Corporation personnel in the larger CSWPI (Center for the Study of Wave Plasma Interactions) SciDAC project. In addition to collaborative efforts in several areas, the PI's thrust research area has focused on the tokamak edge, in the vicinity of the antenna structure, and on the specific physics of interest to the edge, including RF sheath formation. This interest is being driven by a desire to understand measured anomalous power losses that are associated with the edge region. Realistic true 3D geometry simulation capability of the operating antenna, immersed in the edge plasma, has been developed. Geometry accuracy is at the same degree as CAD drawings themselves. In addition, the plasma model includes an RF sheath potential calculation for cells where the plasma is in contact with the metallic structure of the antenna. This capability, in a time-domain simulation, is presently unique. This model can be used to immediately estimate local sheaths (and possibly sputtering damage caused by them), and will be used in the future to model long-range sheaths. Modeling on this scale, in 3D, clearly requires very large scale parallel processing and advanced computing. Ultimately, the desire is to provide quantitative prediction of parasitic edge loss, and edge plasma behavior, analogous to the quantitative core models that have long been a staple of the RF modeling regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance Evaluation of K-DEMO Cable-in-conduit Conductors Using the Florida Electro-Mechanical Cable Model

The United States ITER Project Office (USIPO) is responsible for design of the Toroidal Field (TF) insert coil, which will allow validation of the performance of significant lengths of the conductors to be used in the full scale TF coils in relevant conditions of field, current density and mechanical strain. The Japan Atomic Energy Agency (JAEA) will build the TF insert which will be tested at the Central Solenoid Model Coil (CSMC) Test facility at JAEA, Naka, Japan. Three dimensional mathematical model of TF Insert was created based on the initial design geometry data, and included the following features: orthotropic material properties of superconductor material and insulation; external magnetic field from CSMC, temperature dependent properties of the materials; pre-compression and plastic deformation in lap joint. Major geometrical characteristics of the design were preserved including cable jacket and insulation shape, mandrel outline, and support clamps and spacers. The model is capable of performing coupled structural, thermal, and electromagnetic analysis using ANSYS. Numerical simulations were performed for room temperature conditions; cool down to 4K, and the operating regime with 68kA current at 11.8 Tesla background field. Numerical simulations led to the final design of the coil producing the required strain levels on the cable, while simultaneously satisfying the ITER magnet structural design criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

The study of local overheating and plasma formation on stainless steel z-pinch targets

Plasma formation from intensely ohmically heated conductors is known to be highly non-uniform, as local overheating can be driven by micron-scale imperfections. Detailed understanding of plasma formation is required to predict the performance of magnetically driven physics targets and magnetically-insulated transmission lines (MITLs). Previous LDRD-supported work (projects 178661 and 200269) developed the electrothermal instability (ETI) platform, on the Mykonos facility, to gather high-resolution images of the self-emission from the non-uniform ohmic heating of z-pinch rods. Experiments studying highly inhomogeneous alloyed aluminum captured complex heating topography. To enable detailed comparison with magnetohydrodynamic (MHD) simulation, 99.999% pure aluminum rods in a z-pinch configuration were diamond-turned to ~10nm surface roughness and then further machined to include well-characterized micron-scale "engineered" defects (ED) on the rod's surface (T.J. Awe, et al., Phys. Plasmas 28, 072104 (2021)). In this project, the engineered defect hardware and diagnostic platform were used to study ETI evolution and non-uniform plasma formation from stainless steel targets. The experimental objective was to clearly determine what, if any, role manufacturing, preparation, or alloy differences have in encouraging nonuniform heating and plasma formation from high-current density stainless steel. Data may identify improvements that may be implemented in the fabrication/preparation of electrodes used on the Z machine. Preliminary data shows that difference in manufacturer has no observed effect on ETI evolution, stainless alloy 304L heated more uniformly than alloy 310 at similar current densities, and that stainless steel undergoes the same evolutionary ETI stages as ultra-pure aluminum, with increased emission tied to areas of elevated surface roughness.

36 MATERIALS SCIENCE↗

Investigation of dynamics of ELM crashes and their mitigation techniques

The accurate prediction of H-mode pedestal dynamics is critical for planning experiments in existing tokamaks and in the design of future tokamaks such as ITER and DEMO. The main objective of the proposed research is to advance the understanding of the physics of H-mode pedestal. Through advances in coupled kinetic-MHD simulations, a new model for H-mode pedestal and ELM crashes as well as an improved model for the bootstrap current will be developed. ELMmitigation techniques will also be investigated. The proposed research will help design efficient confinement scenarios and reduce transient heat loads on the divertor and plasma facing components. During the last two years, the principal investigator (PI) of this proposal actively participated in physics studies related to the DOE Joint Research Targets. These studies include the modeling of divertor heat load in the DIII-D, Alcator C-Mod, and NSTX tokamaks in 2010, and the modeling of H-mode pedestal structure in the DIII-D tokamak in 2011. It is proposed that this close collaboration with experimentalists from major US tokamaks continue during the next funding period. Verification and validation will be a strong component of the proposed research. During the course of the project, advances will be made in the following areas; Dynamics of the H-mode pedestal buildup and recovery after ELM crashes – The effects of neutral fueling, particle and thermal pinches will be explored; Dynamics of ELM crashes in realistic tokamak geometries – Heat loads associated with ELM crashes will be validated against experimental measurements. An improved model for ELM crashes will be developed; ELM mitigation – The effect of resonant magnetic perturbations on ELMs stability and their evolution will be investigated; Development of a new bootstrap current model – A reduced model for will be developed through careful verification of existing models for bootstrap current against first-principle kinetic neoclassical simulations. The primary focus here will be on plasma regimes where the existing reduced models deviate from each other, such as high-beta and low aspect ratio plasmas; and Results will be verified against analytical models and simulation results from other codes. Simulation results will be compared with diagnostic experimental data such as Thomson measurements (HRTC) or flux measurements if available. We have arrangements with experimentalists from DIII-D and JET to use diagnostic edge data for this validation. In summary, we propose to study the H-mode pedestal buildup, the stability and evolution of ELMs using a model that will account for dynamically evolving neoclassical sources.

42 ENGINEERING↗

Four-Dimensional Finite-Orbit-Width Fokker-Planck Code with Sources, for Neoclassical/Anomalous Transport Simulation of Ion and Electron (Final Technical Report)

Within the US Department of Energy/Office of Fusion Energy magnetic fusion research program, there is an important whole-plasma-modeling need for a radio-frequency/neutral-beam-injection (RF/NBI) transport-oriented finite-difference Fokker-Planck (FP) code with combined capabilities for 4D (2R2V) geometry near the fusion plasma periphery, and computationally less demanding 3D (1R2V) bounce-averaged capabilities for plasma in the core of fusion devices. Demonstration of proof-of-principle achievement of this goal has been carried out in research carried out under Phase I of the SBIR award. Two DOE-sponsored codes, the CQL3D bounce-average Fokker-Planck code in which CompX has specialized, and the COGENT 4D, plasma edge-oriented Fokker-Planck code which has been constructed by Lawrence Livermore National Laboratory and Lawrence Berkeley Laboratory scientists, where coupled. Coupling was achieved by using CQL3D calculated velocity distributions including an energetic tail resulting from NBI, as boundary conditions for the COGENT code over the two-dimensional velocity space on a spatial interface (flux) surface at a given radius near the plasma periphery. The finite-orbit-width fast ions from the CQL3D distributions penetrated into the peripheral plasma modeled by the COGENT code. This combined code demonstrates the feasibility of the proposed 3D/4D code. By combining these codes, the greatest computational efficiency is achieved subject to present modeling needs in toroidally symmetric magnetic fusion devices. The more efficient 3D code can be used in its regions of applicability, coupled to the more computationally demanding 4D code in higher collisionality edge plasma regions where that extended capability is necessary for accurate representation of the plasma. More efficient code leads to greater use and utility of the model. An ancillary aim of the project is to make the combined 3D/4D code user friendly. Achievement of full-coupling of these two Fokker-Planck codes will advance computational modeling of plasma devices important to the USDOE magnetic fusion energy program, in particular the DIII-D tokamak at General Atomics, San Diego, the NSTX spherical tokamak at Princeton, New Jersey, and the MST reversed-field-pinch Madison, Wisconsin. The validation studies of the code against the experiments will improve understanding of physics important for magnetic fusion, and will increase our design capabilities for achieving the goals of the International Tokamak Experimental Reactor (ITER) project in which the US is a participant and which seeks to demonstrate at least a factor of five in fusion power production divided by input power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Test for Reconditioning RA Waste with Simulated Bitumen and Concrete in a 1,2 MW Plasma Test Facility - 20092

The operation and maintenance of nuclear power plants, the non nuclear fuel cycle, etc generate low-level radioactive waste which, along with the historical radioactive waste from past nuclear activities, needs to be treated and stored, awaiting final disposal. Plasma technology offers a very effective way of treating this waste with a high volume reduction factor (VRF), free from organics, liquids and moisture, and meets without a doubt the acceptance criteria for safe storage and disposal. By means of a plasma beam of approximately 5000 deg. C, the inorganic materials are melted into a glassy slag, containing most of the radioactive isotopes while the organic material is gasified, oxidized and purified in an off-gas cleaning system. First the paper describes the new full-scale Plasma Melting Facility (PMF) at the Kozloduy Nuclear Power Plant in Bulgaria which was taken in nuclear operation In May 2018. The plant has a capacity of 250 tons per year and the maximum contact dose rates of the incoming waste is 2 mSv/h. Different mixtures of radioactive waste packed in 200 l drums were successfully treated resulting in a glassy slag free from liquids and organic material with an important volume reduction factor (VRF). The Project was co-financed through a grant by Kozloduy International Decommissioning fund (KIDSF) administrated by the EBRD through Bulgarian national funding. Plasma is a suitable technology for treatment of problematic waste or even reconditioning waste so Belgoprocess was contracted to do plasma tests with simulated conditioned waste types. One can do tests on a laboratory scale on smaller samples and torch capacities of e.g. 50 kW but Belgoprocess wanted to do more realistic and reliable tests. So Belgoprocess contracted Phoenix Solutions Co who has a full-scope test facility equipped with a 1200 kW plasma torch for full-scope treatment of simulated conditioned waste. For a first confidential contract simulated 200 l (55 gallon) bitumen drums were treated. The drums contained different pucks of compacted waste such as rags, used filters, granulates, etc. The pucks were stacked in the 200 l drums and subsequently embedded with bitumen. A total of 6 drums were treated in the plasma facility. For a second contract simulated homogeneous 200 l (55 gallon) concrete drums with on the one hand concentrates and on the other hand spent resins were selected. A total of 6 drums with concrete and spent resins were treated and melted in the plasma testing facility. The paper describe the test facility, volume reduction factor (VRF) of different waste streams and most important parameters. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

3D Gamma Source Mapping and Intervention Analysis - 19243

3D gamma source mapping and intervention analysis, through Transco Products, Inc.'s RadVision3D{sup R}, uses the N-Visage Scanner{sup TM} by Createc combining laser scanning, optical imaging, and CZT (Cadmium Zinc Telluride) technologies to capture plant geometry, images, and gamma radiation data. By capturing multiple scans in an area of interest a full 3D analysis can be performed using N-Visage Fusion{sup TM} software by Createc which develops a 3D model showing radioactive sources in the environment and allowing various interventions, such as shielding or source removal, to be modeled virtually. The output of a RadVision3D{sup R} project is a scaled virtual environment that the user can 'walk through' on a PC for training and work planning activities. Included in the viewer are dose planes for the as-scanned conditions as well as for each intervention, showing the dose rate at any point in the room. A virtual dosimeter is built into the viewer to allow users to be fully immersed in the environment in which they will be working, without being exposed to actual radiation or other hazardous conditions. This technology can help save dose at operating nuclear power plants when used for planning large work scopes in high radiation areas and is valuable as a characterization tool for decommissioning facilities. Adopters of 3D gamma source mapping and intervention analysis save time, money, and dose by evaluating options well in advance of the performance of any work, and optimizing interventions using accurate modelling techniques to plan work evolutions. (author)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

MagNetUS: a magnetized plasma research ecosystem

MagNetUS is a network of scientists and research groups that coordinates and advocates for fundamental magnetized plasma research in the USA. Its primary goal is to bring together a broad community of researchers and the experimental and numerical tools they use in order to facilitate the sharing of ideas, resources and common tasks. Discussed here are the motivation and goals for this network and details of its formation, history and structure. An overview of associated experimental facilities and numerical projects is provided, along with examples of scientific topics investigated therein. Finally, a vision for the future of the organization is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of Nonlinear Interaction and Turbulence of Alfven Waves in LAPD Experiments

The complete project had two major goals — investigate MHD turbulence generated by counterpropagating Alfven modes, and study such processes in the LAPD device. In order to study MHD turbulence in numerical simulations, two codes have been used: full MHD, and reduced MHD developed specialy for this project. Quantitative numerical results are obtained through high-resolution simulations of strong MHD turbulence, performed through the 2010 DOE INCITE allocation. We addressed the questions of the spectrum of turbulence, its universality, and the value of the so-called Kolmogorov constant (the normalization coefficient of the spectrum). In these simulations we measured with unprecedented accuracy the energy spectra of magnetic and velocity fluctuations. We also studied the so-called residual energy, that is, the difference between kinetic and magnetic energies in turbulent fluctuations. In our analytic work we explained generation of residual energy in weak MHD turbulence, in the process of random collisions of counterpropagating Alfven waves. We then generalized these results for the case of strong MHD turbulence. The developed model explained generation of residual energy is strong MHD turbulence, and verified the results in numerical simulations. We then analyzed the imbalanced case, where more Alfven waves propagate in one direction. We found that spectral properties of the residual energy are similar for both balanced and imbalanced cases. We then compared strong MHD turbulence observed in the solar wind with turbulence generated in numerical simulations. Nonlinear interaction of Alfv´en waves has been studied in the upgraded Large Plasma Device (LAPD). We have simulated the collision of the Alfven modes in the settings close to the experiment. We have created a train of wave packets with the apltitudes closed to those observed n the experiment, and allowed them to collide. We then saw the generation of the second harmonic, resembling that observed in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tunable Laser-Plasma Amplifier (Final Report)

Exploration at the laser intensity frontier has always offered new avenues for physics and reaching beyond this frontier is a grand challenge. Present-day petawatt-class lasers provide focused intensities on target of 10 22 W/cm 2 , corresponding to electric fields of 200 TV/m, while laser-plasma amplification opens a route for focused intensities well above 10 23 W/cm 2 . Intensities in this range provide the ability to test quantum electrodynamics in the unexplored low-energy, strong-field regime where signatures for new physics may arise. The behavior of matter under such extraordinary conditions is a rich and fascinating subject, not only in its own right in fundamental plasma physics, but for the many potential applications that promise to enrich the natural sciences in the future, including compact electron-beam, ion-beam particle accelerators, and ultra-bright X-ray sources. Although laser systems are now under construction internationally to access intensities of 10 23 W/cm 2 , the current technologies used appear to be fundamentally limited to these intensities. The realization of intensities beyond 10 23 W/cm 2 using parametric amplification in plasmas promises a breakthrough in high-energy density physics. Parametric amplification using Raman scattering in a plasma could provide the enabling technology for the generation of ultra-high-power laser pulses, but a more complete understanding of the nonlinear optics of plasmas is required. There is a significant gap in well-diagnosed laser-plasma instability studies of nonlinear plasma-wave phenomena, which are critical to understand for future laser-plasma devices. To achieve an efficient laser-plasma amplifier, plasma waves must be driven to large amplitude where significant energy can be rapidly transferred from the pump to the seed over the pulse duration of the seed. Simulations suggest that this nonlinear pump depletion regime can be achieved in the “pi-pulse” amplification regime. While simulations show this optimal regime with efficient amplification, it has remained elusive in experiments and there is a growing consensus within the community that thermal effects and pump beam limitations prevent laser-plasma amplifiers from progressing through the linear regime into the nonlinear pump depletion regime. Previous experiments have been significantly limited by the laser power available at the necessary wavelengths for the seed laser; therefore, the amplification is required to start in the linear regime where it is sensitive to many deleterious effects. The enabling technology (currently unique to plasma-wave amplification in the world) at the University of Rochester is the ability to provide a seed pulse with sufficient power (4 mJ/100 fs seed) to immediately drive nonlinear plasma waves into the pi-pulse regime and to tune its wavelength to optimize the efficiency of energy transfer. This in combination with the state-of-the-art OMEGA heater beams providing multiple kilojoules in a nanosecond to sufficiently heat the plasma make this system distinct from previous studies. These heater beams will provide, for the first time in Raman amplification studies, a homogeneous electron temperature high enough to prevent pump beam propagation issues that have plagued previous experiments. These systems will provide a platform for driving electron-plasma waves into the nonlinear regime where pump depletion and pulse shortening are predicted to lead to high amplification efficiencies (>30%). The Team has made significant progress through prior support from DOE Fusion Energy Sciences [DOE Office of Science Award Number DE-SC0016253 (2016-2022)]. This includes twenty-two peer-reviewed manuscripts, one patent, ten contributed talks presented at international conferences, and research that was highlighted as invited talks at fifteen international conferences. The broader impacts of this research are evident in the support of early career scientists, two Ph.D. theses, four current graduate students, a Masters Project, two undergraduate researchers, and an underrepresented minority student hired through the California Alliance for Minority Participation who now works as a Research Engineer in the group. This research met all of the funded research objectives and the highlights from primary Raman amplification thrust of this work are discussed below and form the foundation for the proposed research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collaborative research: Dynamics of electrostatic solitary waves and their effects on current layers

The project has accomplished the following achievements including the goals outlined in the original proposal. Generation and measurements of Debye-scale electron holes in laboratory: We have generated by beam injections electron solitary waves in the LAPD experiments. The measurements were made possible by the fabrication of the state-of-the-art microprobes at UCLA to measure Debye-scale electric fields [Chiang et al., 2011]. We obtained a result that challenged the state of knowledge about electron hole generation. We found that the electron holes were not due to two-stream instability, but generated by a current-driven instability that also generated whistler-mode waves [Lefebvre et al., 2011, 2010b]. Most of the grant supported a young research scientist Bertrand Lefebvre who led the dissemination of the laboratory experimental results. In addition to two publications, our work relevant to the laboratory experiments on electron holes has resulted in 7 invited talks [Chen, 2007, 2009; Pickett et al., 2009a; Lefebvre et al., 2010a; Pickett et al., 2010; Chen et al., 2011c, b] (including those given by the co-I Jolene Pickett) and 2 contributed talks [Lefebvre et al., 2009b, a]. Discovery of elecctron phase-space-hole structure in the reconnection electron layer: Our theoretical analyses and simulations under this project led to the discovery of an inversion electric field layer whose phase-space signature is an electron hole within the electron diffusion layer in 2D anti-parallel reconnection [Chen et al., 2011a]. We carried out particle tracing studies to understand the electron orbits that result in the phase-space hole structure. Most importantly, we showed that the current density in the electron layer is limited in collisionless reconnection with negligible guide field by the cyclotron turning of meandering electrons. Comparison of electrostatic solitary waves in current layers observed by Cluster and in LAPD: We compared the ESWs observed in a supersubstorm by the Cluster spacecraft and those measured in LAPD. One of the similarities in the characteristics of ESWs observed in space and in LAPD is that the time duration tends to be approximately the inverse of the electron plasma frequency [Pickett et al., 2009b]. Discovery of suprathermal electron bursts inside a series of magnetic islands: Our effort in examining the roles of ESWs in reconnection current layers resulted in the serendipitous discovery that was published in Nature Physics. In earth’s magnetosphere, we observed through the measurements from the four Cluster spacecraft, a series of magnetic islands and suprathermal electron bursts within the islands. The islands were identified to be effectively acceleration sites for electrons [Chen et al., 2008, 2009].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coulomb excitations for a short linear chain of metallic shells

A self-consistent-field theory is given for the electronic collective modes of a chain containing a finite number, N, of Coulomb-coupled spherical two-dimensional electron gases arranged with their centers along a straight line, for simulating electromagnetic response of a narrow-ribbon of metallic shells. The separation between nearest-neighbor shells is arbitrary and because of the quantization of the electron energy levels due to their confinement to the spherical surface, all angular momenta L of the Coulomb excitations, as well as their projections M on the quantization axis, are coupled. However, for incoming light with a given polarization, only one angular momentum quantum number is usually required. Therefore, the electromagnetic response of the narrow-ribbon of metallic shells is expected to be controlled externally by selecting different polarizations for incident light. We show that, when N = 3, the next-nearest-neighbor Coulomb coupling is larger than its value if they are located at opposite ends of a right-angle triangle forming the triad. Additionally, the frequencies of the plasma excitations are found to depend on the orientation of the line joining them with respect to the axis of quantization since the magnetic field generated from the induced oscillating electric dipole moment on one sphere can couple to the induced magnetic dipole moment on another. Although the transverse inter-shell electromagnetic coupling can be modeled by an effective dynamic medium, the longitudinal inter-shell Coulomb coupling, on the other hand, can still significantly modify the electromagnetic property of this effective medium between shells.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lower-hybrid drift waves and their interaction with plasmas in a 3D symmetric reconnection simulation with zero guide field

We investigate lower-hybrid drift waves (LHDW) in symmetric magnetic reconnection with zero guide field using three-dimensional particle-in-cell simulations. The long-wavelength mode with kρiρe∼1 develops in the bifurcated electron current layer around the X-line within the width of the electron meandering motion from the mid-plane, where ρi(e) is the ion (electron) gyroradius. The short-wavelength mode with kρe∼1 develops in the separatrix region downstream of the electron outflow jet, producing electron vortices in the background flow frame. Electrons follow the E × B drift with corrections from the diamagnetic drift and are heated inside the vortices with diverging electric fields. In the vortices, ions have comparable E × B and inertia drifts, which together mostly cancel the diamagnetic drift. Toward the center of diverging field vortices, ions are decelerated, leading to a decrease in the perpendicular temperature, while the loss of low-energy ions results in an increase in the parallel temperature. Parallel electric fields exist as a combination of the LHDW wave field projected to the magnetic field direction and the penetration of whistler waves that are mainly outside of the LHDW layer. The magnetic flux tube is twisted in the vortices. The twist may potentially lead to slippage reconnection, as indicated by the non-uniform parallel potential variation across field lines, while the periodic variations of the twisting directions are a limiting factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SciDAC Center for Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas

During the first year of the SciDAC gyrokinetic particle simulation (GPS) project, the GPS team (Zhihong Lin, Liu Chen, Yasutaro Nishimura, and Igor Holod) at the University of California, Irvine (UCI) studied the tokamak electron transport driven by electron temperature gradient (ETG) turbulence, and by trapped electron mode (TEM) turbulence and ion temperature gradient (ITG) turbulence with kinetic electron effects, extended our studies of ITG turbulence spreading to core-edge coupling. We have developed and optimized an elliptic solver using finite element method (FEM), which enables the implementation of advanced kinetic electron models (split-weight scheme and hybrid model) in the SciDAC GPS production code GTC. The GTC code has been ported and optimized on both scalar and vector parallel computer architectures, and is being transformed into objected-oriented style to facilitate collaborative code development. During this period, the UCI team members presented 11 invited talks at major national and international conferences, published 22 papers in peer-reviewed journals and 10 papers in conference proceedings. The UCI hosted the annual SciDAC Workshop on Plasma Turbulence sponsored by the GPS Center, 2005-2007. The workshop was attended by about fifties US and foreign researchers and financially sponsored several gradual students from MIT, Princeton University, Germany, Switzerland, and Finland. A new SciDAC postdoc, Igor Holod, has arrived at UCI to initiate global particle simulation of magnetohydrodynamics turbulence driven by energetic particle modes. The PI, Z. Lin, has been promoted to the Associate Professor with tenure at UCI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First neutral beam experiments on Wendelstein 7-X

In the previous divertor campaign, the Wendelstein 7-X (W7-X) device injected 3.6 MW of neutral beam heating power allowing for the achievement of densities approaching 2 × 10 20 m -3 , and providing the first initial assessment of fast ion confinement in a drift optimized stellarator. The neutral beam injection (NBI) system on W7-X is comprised of two beam boxes with space for four radio frequency sources each. The 3.6 MW of heating reported in this work was achieved with two sources in the NI21 beam box. The effect of combined electron-cyclotron resonance heating (ECRH) and NBI was explored through a series of discharges varying both NBI and ECRH power. Discharges without ECRH saw a linear increase in the line-integrated plasma density, and strong peaking of the core density, over the discharge duration. The presence of 1 MW of ECRH power was found to be sufficient to control a continuous density rise during NBI operation. Simulations of fast ion wall loads were found to be consistent with experimental infrared camera images during operation. In general, NBI discharges were free from the presence of fast ion induced Alfvénic activity, consistent with low beam betas. These experiments provide data for future scenario development and initial assessment of fast-ion confinement in W7-X, a key topic of the project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗