Hybrid guiding-centre/full-orbit simulations in non-axisymmetric magnetic geometry exploiting general criterion for guiding-centre accuracy

To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel currents and the local shearing of field-lines. The criterion is derived for general three-dimensional magnetic equilibria including stellarator plasmas. Details are provided on how to implement it in cylindrical coordinates, and in flux coordinates that rely on the geometric toroidal angle. A means of switching between guiding-centre and full-orbit equations at first order in Larmor radius with minimal discrepancy is shown. Techniques are applied to a MAST (Mega Amp Spherical Tokamak) helical core equilibrium in which the inner kinked flux-surfaces are tightly compressed against the outer axisymmetric mantle and where the parallel current peaks at the nearly rational surface. This is put in relation with the simpler situation $\vec{B}(x,y,z) = B_0 [\sin(kx) \vec{e_y} + \cos(kx)\vec{e_z}]$, for which full orbits and lowest order drifts are obtained analytically. In the kinked equilibrium, the full orbits of NBI fast ions are solved numerically and shown to follow helical drift surfaces. This result partially explains the off-axis redistribution of NBI fast particles in the presence of MAST Long-Lived Modes (LLM).Comment: 12 pages, 12 figure

Broadband simulations for M_w 7.8 southern San Andreas earthquakes: Ground motion sensitivity to rupture speed

Using the high-performance computing resources of the Southern California Earthquake Center, we simulate broadband (0–10 Hz) ground motions for three M_w 7.8 rupture scenarios of the southern San Andreas fault. The scenarios incorporate a kinematic rupture description with the average rupture speed along the large slip portions of the fault set at 0.96, 0.89, and 0.84 times the local shear wave velocity. Consistent with previous simulations, a southern hypocenter efficiently channels energy into the Los Angeles region along the string of basins south of the San Gabriel Mountains. However, we find the basin ground motion levels are quite sensitive to the prescribed rupture speed, with peak ground velocities at some sites varying by over a factor of two for variations in average rupture speed of about 15%. These results have important implications for estimating seismic hazards in Southern California and emphasize the need for improved understanding of earthquake rupture processes

Ion cyclotron resonance heating with consistent finite orbit widths and anisotropic equilibria

Minority ion cyclotron resonance heating is studied using the self-consistent numerical model SCENIC. This model includes 3D geometries with full shaping and anisotropic pressure effects, warm contributions to the dielectric tensor and full orbit effects. It evolves the equilibrium, wave field and hot particle distribution function iteratively until a self-consistent solution is found. We will show applications to JET-like two-dimensional equilibria with minority heating scenarios. The effects due to different heating locations on the hot particle distribution function, the hot dielectric tensor and the equilibrium will be studied for symmetric wave injection. Finally, the RF-induced particle pinch is investigated using asymmetric wave injection

Drift stabilization of ballooning modes in a high-⟨β⟩\langle \beta \rangle LHD configuration

Ideal MHD yields at best inconclusive predictions about the stability of the LHD heliotron for $\langle \beta \rangle \geq 3\%$. We investigate the impact of the drift stabilization of ballooning modes for the inward shifted LHD configuration (vacuum magnetic axis $R_0 \sim 3.5m$). The background equilibrium is considered anisotropic in which the neutral beam ions contribute about $1/4$ fraction of the total diamagnetic beta, $\langle \beta_{dia} \rangle$. A drift corrected ballooning mode equation obtained from the linearized gyrokinetic equation is expanded assuming that the hot particle drifts are much larger than the mode frequency. The fast particle pressure gradients contribute weakly to both the instability drive and the diamagnetic drift stabilization (which is dominated by the thermal ion diamagnetic drifts) for $\langle \beta_{dia} \in [0,4.8] \%$. In the single fluid limit (diamagnetic drifts ignored), the thermal pressure gradients drive ballooning modes in a broad region encompassing the outer $60-90 \%$ of the plasma volume at $\langle \beta_{dia} \rangle \approx 4.8 \%$. To stabilize these modes, we find that diamagnetic drift corrections must be invoked (mainly due to the thermal ions). The energetic ion diamagnetic drifts play a role only for low wave number values, $k_{\alpha}\leq 8$. It has been verified that the fast particle drift ordering imposed by the model is amply satisfied for on-axis hot particle to thermal density $N_{h0}/N_{i0} \approx 1\%$ even at high $\langle \beta_{dia} \rangle$

Three-dimensional anisotropic pressure free boundary equilibria

Free boundary three-dimensional anisotropic pressure magnetohydrodynamic equilibria with nested magnetic flux surfaces are computed through the minimisation of the plasma energy functional $W={\int}_{V}{d^3}x\left[{B^2}/(2\mu_0)+p_{||}/(\Gamma-1)\right]$. The plasma–vacuum interface is varied to guarantee the continuity of the total pressure $\left[{p}_{\perp}+{B^2}/(2\mu_0)\right]$ across it and the vacuum magnetic field must satisfy the Neumann boundary condition that its component normal to this interface surface vanishes. The vacuum magnetic field corresponds to that driven by the plasma current and external coils plus the gradient of a potential function whose solution is obtained using a Green's function method. The energetic particle contributions to the pressure are evaluated analytically from the moments of the variant of a bi-Maxwellian distribution function that satisfies the constraint ${\bf B\cdot\nabla}{\cal F}_h=0$. Applications to demonstrate the versatility and reliability of the numerical method employed have concentrated on high-β off-axis energetic particle deposition with large parallel and perpendicular pressure anisotropies in a 2-field period quasiaxisymmetric stellarator reactor system. For large perpendicular pressure anisotropy, the hot particle component of the pperpendicular distribution localises in the regions where the energetic particles are deposited. For large parallel pressure anisotropy, the pressures are more uniform around the flux surfaces

An approximate single fluid 3-dimensional magnetohydrodynamic equilibrium model with toroidal flow

An approximate model for a single fluid three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium with pure isothermal toroidal flow with imposed nested magnetic flux surfaces is proposed. It recovers the rigorous toroidal rotation equilibrium description in the axisymmetric limit. The approximation is valid under conditions of nearly rigid or vanishing toroidal rotation in regions with significant 3D deformation of the equilibrium flux surfaces. Bifurcated helical core equilibrium simulations of long-lived modes in the MAST device demonstrate that the magnetic structure is only weakly affected by the flow but that the 3D pressure distortion is important. The pressure is displaced away from the major axis and therefore is not as noticeably helically deformed as the toroidal magnetic flux under the subsonic flow conditions measured in the experiment. The model invoked fails to predict any significant screening by toroidal plasma rotation of resonant magnetic perturbations in MAST free boundary computations

Western Bumble Bee: Declines in the Continental United States and Range-Wide Information Gaps

In recent decades, many bumble bee species have declined due to changes in habitat, climate, and pressures from pathogens, pesticides, and introduced species. The western bumble bee (Bombus occidentalis), once common throughout western North America, is a species of concern and will be considered for listing by the U.S. Fish and Wildlife Service (USFWS) under the Endangered Species Act (ESA). We attempt to improve alignment of data collection and research with USFWS needs to consider redundancy, resiliency, and representation in the upcoming species status assessment. We reviewed existing data and literature on B. occidentalis, highlighting information gaps and priority topics for research. Priorities include increased knowledge of trends, basic information on several life‐history stages, and improved understanding of the relative and interacting effects of stressors on population trends, especially the effects of pathogens, pesticides, climate change, and habitat loss. An understanding of how and where geographic range extent has changed for the two subspecies of B. occidentalis is also needed. We outline data that could be easily collected in other research projects that would increase their utility for understanding range‐wide trends of bumble bees. We modeled the overall trend in occupancy from 1998 to 2018 of Bombus occidentalis within the continental United States using existing data. The probability of local occupancy declined by 93% over 21 yr from 0.81 (95% CRI = 0.43, 0.98) in 1998 to 0.06 (95% CRI = 0.02, 0.16) in 2018. The decline in occupancy varied spatially by landcover and other environmental factors. Detection rates vary in both space and time, but peak detection across the continental United States occurs in mid‐July. We found considerable spatial gaps in recent sampling, with limited sampling in many regions, including most of Alaska, northwestern Canada, and the southwestern United States. We therefore propose a sampling design to address these gaps to best inform the ESA species status assessment through improved assessment of how the spatial distribution of stressors influences occupancy changes. Finally, we request involvement via data sharing, participation in occupancy sampling with repeated visits to distributed survey sites, and complementary research to address priorities outlined in this paper

System Engineering Paper

The Iowa State University team, Team LunaCY, is composed of the following sub-teams: the main student organization, the Lunabotics Club; a senior mechanical engineering design course, ME 415; a senior multidisciplinary design course, ENGR 466; and a senior design course from Wartburg College in Waverly, Iowa. Team LunaCY designed and fabricated ART-E III, Astra Robotic Tractor- Excavator the Third, for the team's third appearance in the NASA Lunabotic Mining competition. While designing ART-E III, the team had four main goals for this year's competition:to reduce the total weight of the robot, to increase the amount of regolith simulant mined, to reduce dust, and to make ART-E III autonomous. After many designs and research, a final robot design was chosen that obtained all four goals of Team LunaCY. A few changes Team LunaCY made this year was to go to the electrical, computer, and software engineering club fest at Iowa State University to recruit engineering students to accomplish the task of making ART-E III autonomous. Team LunaCY chose to use LabView to program the robot and various sensors were installed to measure the distance between the robot and the surroundings to allow ART-E III to maneuver autonomously. Team LunaCY also built a testing arena to test prototypes and ART-E III in. To best replicate the competition arena at the Kennedy Space Center, a regolith simulant was made from sand, QuickCrete, and fly ash to cover the floor of the arena. Team LunaCY also installed fans to allow ventilation in the arena and used proper safety attire when working in the arena . With the additional practice in the testing arena and innovative robot design, Team LunaCY expects to make a strong appearance at the 2012 NASA Lunabotic Mining Competition.