Search NASA⌕ Search

SEARCH · Search NASA

Results for “collisionality”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

How sheath properties change with gas pressure: modeling and simulation

Abstract Particle-in-cell simulations are used to study how neutral pressure influences plasma properties at the sheath edge. The high rate of ion–neutral collisions at pressures above several mTorr are found to cause a decrease in the ion velocity at the sheath edge (collisional Bohm criterion), a decrease in the edge-to-center density ratio ( h l factor), and an increase in the sheath width and sheath potential drop. A comparison with existing analytic models generally indicates favorable agreement, but with some distinctions. One is that models for the h l factor need to be made consistent with the collisional Bohm criterion. With this and similar corrections, a comprehensive fluid-based model of the plasma boundary transition is constructed that compares well with the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Observation of Slow Eigen-Zundel Interconversion in H + (H 2 O) 6 Clusters upon Isomer-Selective Vibrational Excitation and Buffer Gas Cooling in a Cryogenic Ion Trap

The formation of isomers when trapping floppy cluster ions in a temperature-controlled ion trap is a generally observed phenomenon. This involves collisional quenching of the ions initially formed at high temperature by buffer gas cooling until their internal energies fall below the barriers in the potential energy surface that separate them. Here we explore the kinetics at play in the case of the two isomers adopted by the H + (H 2 O) 6 cluster ion that differ in the proton accommodation motif. One of these is most like the Eigen cation with a tricoordinated hydronium motif (denoted E), and the other is most like the Zundel ion with the proton equally shared between two water molecules (denoted Z). After initial cooling to about 20 K in the radiofrequency (Paul) trap, the relative populations of these two spectroscopically distinct isomers are abruptly changed through isomer-selective photoexcitation of bands in the OH stretching region with a pulsed (~6 ns) infrared laser while the ions are in the trap. We then monitor the relaxation of the vibrationally excited clusters and reformation of the two cold isomers by recording infrared photodissociation spectra with a second IR laser as a function of delay time from the initial excitation. The latter spectra are obtained after ejecting the trapped ions into a time-of-flight photofragmentation mass spectrometer, thus enabling long (~0.1 s) delay times. Excitation of the Z isomer is observed to display long-lived vibrationally excited states that are collisionally cooled on a ms time scale, some of which quench into the E isomer. These excited E species then display spontaneous interconversion to the Z form on a ~10 ms time scale. Furthermore, these qualitative observations set the stage for a series of experimental measurements that can provide quantitative benchmarks for theoretical simulations of cluster dynamics and the potential energy surfaces that underlie them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Validation of 2D Te and ne measurements made with Helium imaging spectroscopy in the volume of the TCV divertor

Abstract Multi-spectral imaging of helium atomic emission (HeMSI) has been used to create 2D poloidal maps of T e and n e in TCV’s divertor. To achieve these measurements, TCV’s MANTIS multispectral cameras (Perek et al 2019 Rev. Sci. Instrum. 90 123514) simultaneously imaged four He I lines (two singlet and two triplet) and a He II line (468 nm) from passively present He and He + . The images, which were absolutely calibrated and covered the whole divertor region, were inverted through the assumption of toroidal symmetry to create emissivity profiles and, consequently, line-ratio profiles. A collisional-radiative model (CRM) was applied to the line-ratio profiles to produce 2D poloidal maps of T e and n e . The collisional-radiative modeling was accomplished with the Goto helium CRM code (Zholobenko et al 2018 Nucl. Fusion 58 126006, Zholobenko et al 2018 Technical Report , Goto 2003 J. Quant. Spectrosc. Radiat. Transfer 76 331–44) which accounts for electron-impact excitation (EIE) and deexcitation, and electron–ion recombination (EIR) with He + . The HeMSI T e and n e measurements were compared with co-local Thomson scattering measurements. The two sets of measurements exhibited good agreement for ionizing plasmas: ( 5 eV ⩽ T e ⩽ 60 eV , and 2 × 10 18 m − 3 ⩽ n e ⩽ 3 × 10 19 m − 3 ) in the case of majority helium plasmas, and ( 10 eV ⩽ T e ⩽ 40 eV , 2 × 10 18 m − 3 ⩽ n e ⩽ 3 × 10 19 m − 3 ) in the case of majority deuterium plasmas. However, there were instances where HeMSI measurements diverged from Thomson scattering. When T e ⩽ 10 eV in majority deuterium plasmas, HeMSI deduced inaccurately high values of T e . This disagreement cannot be rectified within the CRM’s EIE and EIR framework. Second, on sporadic occasions within the private flux region, HeMSI produced erroneously high measurements of n e . Multi-spectral imaging of Helium emission has been demonstrated to produce accurate 2D poloidal maps of T e and n e within the divertor of a tokamak for plasma conditions relevant to contemporary divertor studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Plasma heating by electron cyclotron wave and the temperature effects on lower hybrid current drive on EAST

This paper presents the progress in the long—pulse operation of the electron cyclotron (EC) system and the achievements in high—electron temperature plasmas by the combined EC and lower hybrid (LH) waves heating since the EC system was built in 2015. An electron temperature of up to 12 keV with a duration over 100 s was realized by the simultaneous heating of EC and LH waves at the line-averaged density ${\bar n_{\text{e}}}$ ~1.8 × 10 19 m –3 . The plasma heating effect strongly depends on the location of EC power deposition. H-mode plasmas with solely EC wave auxiliary heating have been obtained on EAST for the first time. These H-mode discharges show an enhanced confinement factor H 98(y,2) around 1.0, which is higher than the previous H-mode using LH power alone (Xu et al 2011 Nucl. Fusion 51 072001). The total heating power is very close to the threshold value for L–H transition according to the international tokamak scaling. In addition to the temperature effects inside the separatrix, higher electron temperature produced by EC wave is found to reduce the LH power loss in the scrape-off layer due to collisional absorption, which is beneficial to further increase the LH current drive efficiency. Ray-tracing/Fokker–Planck modeling results indicate that higher electron temperature can shorten the LH wave propagation on EAST in a multiple-pass regime, thus decreasing the collisional dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Closure and transport theory for plasmas with multiple ion species

This report summarizes a DOE-funded research project on closure and transport theory for plasmas with multiple ion species. The work developed accurate moment-based closure models across the full range of collisionality, incorporating both kinetic and collisional effects. Key contributions include high-fidelity closure relations, multi-temperature models, and transport coefficients applicable to fusion, space, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine Learning Framework for Conotoxin Class and Molecular Target Prediction

Conotoxins are small and highly potent neurotoxic peptides derived from the venom of marine cone snails which have captured the interest of the scientific community due to their pharmacological potential. These toxins display significant sequence and structure diversity, which results in a wide range of specificities for several different ion channels and receptors. Despite the recognized importance of these compounds, our ability to determine their binding targets and toxicities remains a significant challenge. Predicting the target receptors of conotoxins, based solely on their amino acid sequence, remains a challenge due to the intricate relationships between structure, function, target specificity, and the significant conformational heterogeneity observed in conotoxins with the same primary sequence. We have previously demonstrated that the inclusion of post-translational modifications, collisional cross sections values, and other structural features, when added to the standard primary sequence features, improves the prediction accuracy of conotoxins against non-toxic and other toxic peptides across varied datasets and several different commonly used machine learning classifiers. Here, we present the effects of these features on conotoxin class and molecular target predictions, in particular, predicting conotoxins that bind to nicotinic acetylcholine receptors (nAChRs). We also demonstrate the use of the Synthetic Minority Oversampling Technique (SMOTE)-Tomek in balancing the datasets while simultaneously making the different classes more distinct by reducing the number of ambiguous samples which nearly overlap between the classes. In predicting the alpha, mu, and omega conotoxin classes, the SMOTE-Tomek PCA PLR model, using the combination of the SS and P feature sets establishes the best performance with an overall accuracy (OA) of 95.95%, with an average accuracy (AA) of 93.04%, and an f1 score of 0.959. Using this model, we obtained sensitivities of 98.98%, 89.66%, and 90.48% when predicting alpha, mu, and omega conotoxin classes, respectively. Similarly, in predicting conotoxins that bind to nAChRs, the SMOTE-Tomek PCA SVM model, which used the collisional cross sections (CCSs) and the P feature sets, demonstrated the highest performance with 91.3% OA, 91.32% AA, and an f1 score of 0.9131. The sensitivity when predicting conotoxins that bind to nAChRs is 91.46% with a 91.18% sensitivity when predicting conotoxins that do not bind to nAChRs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Implementation and assessment of an extended hydrogenic molecular model in UEDGE

Abstract Two sets of pure deuterium plasmas are simulated using UEDGE, one including atoms and molecules and the other including atoms only as separate fluid species. Simulation results in one and two dimensions are reported to assess the role of molecules in tokamak plasmas. It is shown that thermal coupling of the molecules to the plasma can be relevant for the onset of detachment under highly collisional conditions in simple geometries, but play a modest role at low collisionality. Ion‐molecule equipartition presents an additional ion energy loss channel, dissipating ion energy that would otherwise heat the electrons by equipartition. The resulting reduced electron temperature increases the ion‐electron recombination by an order of magnitude as the plasma density is increased, yielding earlier and deeper detachment.

Holm, A.↗
DOE OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗
DOE OSTI2020

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

SPEARS: A Database-Invariant Spectral modeling API

The Spectral Physics Environment for Advanced Remote Sensing (SPEARS) application programming interface (API) is a Python-based, line-by-line, local thermal equilibrium (LTE) spectral modeling code which is optimized for simultaneously synthesizing optical spectra from any combination of fundamental spectroscopic databases. In this article, we contribute two novel spectral modeling techniques to the scientific literature. First we describe how SPEARS integrates a physics-based collisional model for calculating pressure broadening in the absence of available broadening coefficients. With this collisional model implementation, a generalized approach to fundamental spectroscopic databases can be achieved across multiple databases. We also detail our adaptive grid mesh algorithm developed to make the code scalable for simulating large spectral bandwidths at high spectral fidelity using intuitive grid parameters. Here, we present comparisons to other modeling tools, experiments, and provide a discussion on the SPEARS user interface.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗
DOE OSTI2020

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗