Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “separatrix crossing”

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 19 records

Minimizing separatrix crossings through isoprominence

Abstract A simple property of magnetic fields that minimizes bouncing to passing type transitions of guiding center orbits is defined and discussed. This property, called isoprominence, is explored through the framework of a near-axis expansion. It is shown that isoprominent magnetic fields for a toroidal configuration exist to all orders in a formal expansion about a magnetic axis. Some key geometric features of these fields are described.

confinement↗

Chaotic jumps in the generalized first adiabatic invariant in current sheets

The present study examines how the changes in the generalized first adiabatic invariant J derived from the separatrix crossing theory can be incorporated into the drift variable approach to generating distribution functions. A method is proposed for determining distribution functions for an ensemble of particles following interaction with the tail current sheet by treating the interaction as a scattering problem characterized by changes in the invariant. Generalized drift velocities are obtained for a 1D tail configuration by using the generalized first invariant. The invariant remained constant except for the discrete changes caused by chaotic scattering as the particles cross the separatrix.

Brittnacher, M. J.↗

Topological structure of the magnetotail as a function of interplanetary magnetic field direction

Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Fedder, J. A.↗

Balmer emission measurements in JET-ILW hydrogen, deuterium, tritium and deuterium-tritium low-confinement mode plasmas

Measurements of Balmer-alpha (B-α) and Balmer-gamma (B-γ), where B = H, D, or T isotopes of hydrogen, emission across the low-field side (LFS) divertor in JET ITER-like Wall (ILW), low-confinement mode (L-mode) plasma indicate that the divertor plasma was more strongly detached in deuterium, tritium and deuterium–tritium than in hydrogen plasmas. These measurements are consistent with the lower degree of detachment in hydrogen versus tritium plasmas inferred from Langmuir probes (Groth et al., 2023). The ratio of B-γ to B-α is a suitable indicator of high density/low temperature conditions because B-α emission is dominated by electron impact excitation with a small electron–ion recombination component regardless of electron temperature and density (for conditions in tokamak plasmas) whilst for B-γ the recombination component increases and the excitation component decreases as the electron temperature decreases and the electron density increases. Within the uncertainty of the measurements, the electron density at the low-field side separatrix and the power crossing the separatrix are similar to permit these comparisons. The B-©/B-α ratio is well correlated with the onset of detachment as measured by Langmuir probes (McCracken et al., 1998; Groth et al., 2023). The ratios of B-γ and B-α varied between 0.005 and 0.035 for low-recycling, high-recycling and partially detached conditions, and were 60–70 % higher in partially detached hydrogen than in tritium plasmas. In contrast, the difference between deuterium, tritium and 40–60 % deuterium–tritium is within 10 % and thus within the uncertainties of the measurements. Profiles of the B-α and B-γ emission was measured with two filtered CCDs connected to the visible lost light leg of JET’s mirror link spectroscopy system (Meigs et al., 2010; Lomanowski et al., 2014). Spatially (radial) narrower B-α profiles were measured with increasing isotope mass, while similarly broadened B-γ profiles were measured.

Balmer↗

Ion precipitation from the inner plasma sheet due to stochastic diffusion

Plasma sheet ions do not conserve their first adiabatic invariant when the magnetic field is appreciably tail-like. They do conserve a different adiabatic invariant but only to linear, rather than exponential, accuracy in the appropriate small parameter. Thus significant stochastic diffusion can occur for particles crossing the separatrix dividing the segments of orbits on which the particles cross and do not cross the tail midplane. Such ions can escape the plasma sheet and precipitate into the atmosphere. Stochastic scattering is strongest from those field lines where the ion's Larmor period in the normal component of the neutral sheet magnetic field approximately equals its bounce period. By comparing the rates of stochastic ion loss and convection in the tail, it is possible to estimate the location and thickness of the inner edge of the ion plasma sheet created by stochastic ion loss. Ions of different masses precipitate into the atmosphere at slightly different locations. Since wave particle interactions are not needed, this precipitation will always occur and should be particularly evident during quiet geomagnetic conditions, when it is less likely to be masked by other precipitation mechanisms.

Zelenyi, L.↗

The separatrix electron density in JET, ASDEX upgrade and alcator C-Mod H-mode plasmas: A common evaluation procedure and correlation with engineering parameters

The separatrix electron density is an important parameter for core-edge scenario integration in tokamak devices, as it influences plasma confinement, divertor detachment and disruption avoidance. This quantity has been measured in H-mode discharges on JET, ASDEX Upgrade and Alcator C-Mod by applying the same fitting function to Thomson scattering measurements, and by employing the same analysis technique based on scrape-off layer power balance. To estimate the power crossing the separatrix, the inter-ELM time derivative of the plasma energy d W /d t has been experimentally evaluated and found to be approximately a constant fraction of the absorbed heating power. Correlations between n e,sep and engineering parameters have been investigated, revealing that n e,sep scales with the divertor neutral pressure p 0,div in a similar manner across all devices. Additionally, when n e,sep is normalized to the obtained n e,sep dependency, no clear correlation with the plasma current is found. These observations are in agreement with the 2-point model, which suggests that the upstream separatrix density is mainly set by the recycling at the divertor target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analytic model for neutral penetration and plasma fueling

Neutral atoms recycled from wall interaction interact with confined plasma, thereby refuelling it, most strongly in the region closest to the wall. This occurs near the X-point in diverted configurations, or else near the wall itself in limited configurations. A progression of analytic models is developed for neutral density in the vicinity of a planar or linear source in an ionising domain. First-principles neutral transport simulations with DEGAS2 are used throughout to test the validity and limits of the model when using equivalent sources. The model is further generalised for strong plasma gradients or the inclusion of charge exchange. An important part of the problem of neutral fuelling from recycling is thereby isolated and solved with a closed-form analytic model. A key finding is that charge exchange with the confined plasma can be significantly simplified with a reasonable sacrifice of accuracy by treating it as a loss. The several assumptions inherent to the model (and the simulations with which it is compared) can be adapted according to the particular behaviour of neutrals in the divertor and the manner in which they cross the separatrix.

fusion plasma↗

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing W sources in the all-W wall, all-RF WEST tokamak environment

In this work, experimental data, together with interpretive modeling tools, are examined to study trends in the tungsten (W) source in the all-W environment of the WEST tokamak, both from the divertor and from the main chamber. In particular, a poloidal limiter protecting an ion cyclotron resonance heating (ICRH) antenna is used as proxy for main chamber sourcing. The key study is carried out by stepping up lower hybrid current drive (LHCD) power, as the only auxiliary power source. Limiter and divertor W sources exhibit a qualitatively similar proportionality to the total power crossing the separatrix, P SEP , although the main chamber source remains substantially lower than the divertor source, for the range of P SEP accessible in the experiments. Intepretive modeling of the limiter source is carried out with a particle-in-cell (PIC) sheath model coupled to a surface sputtering model. Oxygen is used as a proxy for all light impurity species allowing for characterization of the critical W erosion regions. To get a good quantitative match to the data, it is necessary to assume that the oxygen arrives at the surface mostly at high ionization stages (4+ and above). A separate simulation with SOLEDGE-EIRENE, constrained to measured upstream scrape-off-layer plasma profiles, gives oxygen fractional abundances that are compatible with the PIC simulation result. This is understood to arise from transport processes that dominate over recombination. Substituting the LHCD by ICRH, in an equivalent experiment, the local W source exhibits a 3× enhancement. This can be matched by the simulation, by assuming local RF electric field rectification, based on ~100 eV peak-to-peak, near-antennna electric field. This work has highlighted the particular importance of understanding the ion charge state balance of light impurities as these are most likely the dominant sputtering species in fusion devices with high-Z walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predicting core transport in ITER baseline discharges with neon injections

Achieving self-consistent performance predictions for ITER requires integrated modeling of core transport and divertor power exhaust under realistic impurity conditions. We present results from a systematic power-flow and impurity-content study for the ITER 15 MA baseline scenario constrained directly by existing SOLPS-ITER neon-seeded divertor solutions. Using the OMFIT STEP workflow, stationary temperature and density profiles are predicted with TGYRO for $1.5 \unicode{x2A7D} Z_\textrm{eff} \unicode{x2A7D} 2.5$, and the corresponding power crossing the separatrix $P_\textrm{sep}$ is evaluated. We find that $P_\textrm{sep}$ varies by more than a factor of 1.7 across this scan and matches the ${\sim}100$ MW SOLPS-ITER prediction when $Z_\textrm{eff} \simeq 1.6$ or when auxiliary heating is reduced to ${\sim}75\%$ of nominal. Rotation-sensitivity studies show that plausible variations in toroidal flow magnitude modify $P_\textrm{sep}$ by $\lesssim 20\%$, while AURORA modeling confirms that charge-exchange radiation inside the separatrix is dynamically negligible under predicted ITER neutral densities. These results identify a restricted compatibility window, $Z_\textrm{eff} \approx 1.6$ –1.75 and $0.75 \lesssim f_{P_\textrm{aux}} \unicode{x2A7D} 1.0$, in which core transport predictions remain aligned with neon-seeded divertor protection targets. This self-consistent, model-constrained framework provides actionable guidance for impurity control and auxiliary-heating scheduling in early ITER operation and supports future whole-device scenario optimization.

ITER↗

Developing high performance RF heating scenarios on the WEST tokamak

High power experiments, up to 9.2 MW with LHCD and ICRH, have been carried out in the full tungsten tokamak WEST. Here, quasi non inductive discharges have been achieved allowing to extend the plasma duration to 53 s with stationary conditions in particular with respect to tungsten contamination. Transitions to H mode are observed, and H-modes lasting up to 4 s have been obtained. The increase in stored energy is weak since the power crossing the separatrix is close to the L–H threshold. Hot L mode plasmas (central temperature exceeding 3 keV) with a confinement time following the ITER L96 scaling law are routinely obtained. The weak aspect ratio dependence of this scaling law is confirmed. Tungsten accumulation is generally not an operational issue on WEST. Difficulty of burning through tungsten can prevent the discharge from accessing to a hot core plasma in the ramp-up phase, or can lead to rapid collapse of the central temperature when radiation is enhanced by a slight decrease of the temperature. Except a few pulses post-boronization, the plasma radiation is rather high (P rad /P tot ~ 50%) and is dominated by tungsten. This fraction does not vary as the RF power is ramped up and is quite similar in ICRH and/or LHCD heated plasmas. An estimate of the contribution of the RF antennas to the plasma contamination in tungsten is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced pedestal transport driven by edge collisionality on Alcator C-Mod and its role in regulating H-mode pedestal gradients

Experimental measurements of plasma and neutral profiles across the pedestal are used in conjunction with 2D edge modeling to examine pedestal stiffness in Alcator C-Mod H-mode plasmas. Enhanced D α experiments on Alcator C-Mod observed pedestal degradation and loss in confinement below a critical value of net power crossing the separatrix, P net = $P^{crit}_{net}$ ≈ 2.3 MW, in the absence of any external fueling. New analysis of ionization and particle flux profiles reveal saturation of the pedestal electron density, $n^{ped}_{e}$, despite continuous increases in ionization throughout the pedestal, inversely related to P net . A limi to the pedestal $\nabla$n e emerges as the particle flux, Γ D , continues to grow, implying increases in the effective particle diffusivity, D eff . This is well-correlated with the separatrix collisionality, $v^*_{sep}$ and a turbulence control parameter, α t , implying a possible transition in type of turbulence. The transition is well correlated with the experimentally observed value of $P^{crit}_{net}$. SOLPS-ITER modeling is performed for select discharges from the power scan, constrained with experimental electron and neutral densities, measured at the outer midpane. The modeling confirms general growth in D eff , consistent with experimental findings, and additionally suggests even larger growth in Χ e at the same $P^{crit}_{net}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observations of magnetic merging and the formation of the plasma sheet in the earth's magnetotail

Particle and magnetic field observations of the field line merging process in the earth's magnetotail are presented. By analyzing the lunar shadow pattern of electron fluxes detected by the lunar-orbiting Apollo 16 subsatellite it has been possible to identify the topology and to measure the velocity of the magnetotail field lines. The observations reported here were made as the moon crossed the separatrix between premerging and merged field lines. The measured field line velocities toward the merging region were 30-60 km/s, and the thickness of the separatrix was estimated to be about 2000 km. Most of the magnetic energy released in the merging process appears to go into the energization of particles. The length and the thickness of the merging region are inferred to be of the order of about 10 earth radii and about 4000 km, respectively. The energized particles travel away from the merging region along the separatrix. Those headed earthward may form the plasma sheet by being trapped on closed field lines. The rate of energization and the energy spectrum of those particles are consistent with those required for formation of the plasma sheet.

Lin, R. P.↗

Plasma heating by collisionless magnetic reconnection - Analysis and computation

Analytic and numerical results on particle acceleration in 2D collisionless magnetic reconnection are presented. The particles are followed until they reach an outgoing flux surface at the same distance from the origin as the starting surface. The magnetic moment is not conserved for particles passing through the unmagnetized region around the X line at the origin. Other particles cross the separatrix without passing near the X line. The magnetic moment of the first class of outgoing particles is randomized, whereas it can be considered for the second class. The analytic model is based upon the observation of the final kinetic energy as a function of the initial conditions. Analytic results are shown to predict a Maxwellian tail for the distribution function in the perpendicular kinetic energy, with this energy much greater than the parallel kinetic energy. Numerical results showing that the predicted tail temperature agrees with the numerically computed temperature to within 10 percent over 4 orders of magnitude in the electric field are presented.

Moses, R. W.↗

Electromagnetic turbulence simulation of tokamak edge plasma dynamics and divertor heat load during thermal quench

Abstract The edge plasma turbulence and transport dynamics, as well as the divertor power loads during the thermal-quench phase of tokamak disruptions, are numerically investigated with BOUT++’s flux-driven six-field electromagnetic turbulence model. Here, transient yet intense particle and energy sources are applied at the pedestal top to mimic the plasma power drive at the edge induced by a core thermal collapse, which flattens the core temperature profile. Interesting features, such as surging of divertor heat load (up to 50 times) and broadening of heat-flux width (up to four times) on the outer-divertor target plate, are observed in the simulation, in qualitative agreement with experimental observations. The dramatic changes in divertor heat load and width are due to the enhanced plasma turbulence activities inside the separatrix. Two cross-field transport mechanisms, namely, the E × B turbulent convection and the stochastic parallel advection/conduction, are identified to play important roles in this process. First, an elevated edge pressure gradient drives instabilities and subsequent turbulence in the entire pedestal region. The enhanced turbulence not only transports particles and energy radially across the separatrix via the E × B convection, which causes the initial divertor heat-load burst, but it also induces amplified magnetic fluctuation B ˜ . Once themagnetic fluctuation is large enough to break the magnetic flux surface, magnetic flutter effect provides an additional radial transport channel. In the late stage of our simulation, | B ˜ r / B 0 | reaches to 10 −4 level that completely breaks magnetic flux surfaces such that stochastic field lines are directly connecting pedestal top plasma to the divertor target plates or first wall, further contributing to the divertor heat-flux width broadening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Some comments on the topology of the geomagnetic tail

The magnetotail topology of the noon-midnight meridian is considered, taking into account the determination of the trajectory of the X line across the tail current sheet. The X line is the line in the current sheet that separates plasma sheet field lines from interplanetary field lines. Attention is given to recent measurements of the boundary layer, the entry layer, and a crossing of the separatrix at the orbit of the moon in terms of the topological sketches.

Russell, C. T.↗

Lower-Hybrid Wave Structures and Interactions With Electrons Observed in Magnetotail Reconnection Diffusion Regions

We investigate waves close to the lower-hybrid frequency in 12 magnetotail reconnection electron diffusion region (EDR) events with guide field levels of near-zero to 30%. In about half of the events, the wave vector has a small component along the current sheet normal, consistent with known lower-hybrid drift wave properties, but the perpendicular magnetic field fluctuations can be comparable or greater than the parallel component, a feature unique to the waves inside and adjacent to EDRs. Another new wave property is that the wave vector has a significant component along the current sheet normal in some events and completely along the normal for one event. In 1/4of the events, the ∇∙𝑷(sub 𝑒) term has a significant contribution to the wave electric field, possibly a feature of lower-hybrid waves more likely to exist in the diffusion region than further away from the X-line. Electron temperature variations are correlated with the wave potential, due to wave electric field acceleration and crossings at the corrugated separatrix region with different amounts of mixing between reconnection inflowing and outflowing populations. The latter also leads to the anti-correlation between parallel and perpendicular temperature components. Using four-spacecraft measurements, the magnetic field line twisting is demonstrated by the correlated fluctuations in (∇×𝑽(sub 𝐸×𝐵))||and (∇×𝐁)||. The lower-hybrid wave in the EDR of weak guide field reconnection may be generated near separatrices and penetrate to the mid-plane or locally generated, and the latter possibility is beyond the prediction of previous reconnection simulations.

Magnetic reconnection↗