ترغب بنشر مسار تعليمي؟ اضغط هنا

Continuum Kinetic and Multi-Fluid Simulations of Classical Sheaths

114   0   0.0 ( 0 )
 نشر من قبل Petr Cagas
 تاريخ النشر 2016
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

The kinetic study of plasma sheaths is critical, among other things, to understand the deposition of heat on walls, the effect of sputtering, and contamination of the plasma with detrimental impurities. The plasma sheath also provides a boundary condition and can often have a significant global impact on the bulk plasma. In this paper, kinetic studies of classical sheaths are performed with the continuum code, Gkeyll, that directly solves the Vlasov-Poisson/Maxwell equations. The code uses a novel version of the finite-element discontinuous Galerkin (DG) scheme that conserves energy in the continuous-time limit. The electrostatic field is computed using the Poisson equation. Ionization and scattering collisions are included, however, surface effects are neglected. The aim of this work is to introduce the continuum-kinetic method and compare its results to those obtained from an already established finite-volume multi-fluid model also implemented in Gkeyll. Novel boundary conditions on the fluids allow the sheath to form without specifying wall fluxes, so the fluids and fields adjust self-consistently at the wall. The work presented here demonstrates that the kinetic and fluid results are in agreement for the momentum flux, showing that in certain regimes, a multi-fluid model can be a useful approximation for simulating the plasma boundary. There are differences in the electrostatic potential between the fluid and kinetic results. Further, the direct solutions of the distribution function presented here highlight the non-Maxwellian distribution of electrons in the sheath, emphasizing the need for a kinetic model.



قيم البحث

اقرأ أيضاً

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma-wall interactions in thrusters and fusion devices. For mul tifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remains an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.
Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv{e}n waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of $k_{perp}^{-1.3}$. The kinetic code shows a spectral slope of $k_{perp}^{-1.5}$ for smaller simulation domain, and $k_{perp}^{-1.3}$ for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.
540 - A.O. Korotkevich , 2010
We wish to transmit messages to and from a hypersonic vehicle around which a plasma sheath has formed. For long distance transmission, the signal carrying these messages must be necessarily low frequency, typically 2 GHz, to which the plasma sheath i s opaque. The idea is to use the plasma properties to make the plasma sheath appear transparent.
Continuum kinetic simulations of plasmas, where the distribution function of the species is directly discretized in phase-space, permits fully kinetic simulations without the statistical noise of particle-in-cell methods. Recent advances in numerical algorithms have made continuum kinetic simulations computationally competitive. This work presents the first continuum kinetic description of high-fidelity wall boundary conditions that utilize the readily available particle distribution function. The boundary condition is realized through a reflection function that can capture a wide range of cases from simple specular reflection to more involved first principles models. Examples with detailed discontinuous Galerkin implementation are provided for secondary electron emission using phenomenological and first-principles quantum-mechanical models. Results presented in this work demonstrate the effect of secondary electron emission on a classical plasma sheath.
Reduced fluid models including electron inertia and ion finite Larmor radius corrections are derived asymptotically, both from fluid basic equations and from a gyrofluid model. They apply to collisionless plasmas with small ion-to-electron equilibriu m temperature ratio and low $beta_e$, where $beta_e$ indicates the ratio between the equilibrium electron pressure and the magnetic pressure exerted by a strong, constant and uniform magnetic guide field. The consistency between the fluid and gyrofluid approaches is ensured when choosing ion closure relations prescribed by the underlying ordering. A two-field reduction of the gyrofluid model valid for arbitrary equilibrium temperature ratio is also introduced, and is shown to have a noncanonical Hamiltonian structure. This model provides a convenient framework for studying kinetic Alfven wave turbulence, from MHD to sub-$d_e$ scales (where $d_e$ holds for the electron skin depth). Magnetic energy spectra are phenomenologically determined within energy and generalized helicity cascades in the perpendicular spectral plane. Arguments based on absolute statistical equilibria are used to predict the direction of the transfers, pointing out that, within the sub-ion range associated with a $k_perp^{-7/3}$ transverse magnetic spectrum, the generalized helicity could display an inverse cascade if injected at small scales, for example by reconnection processes.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا