اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدAn Explicitly Solvable Energy-Conserving Algorithm for Pitch-Angle Scattering in Magnetized Plasmas
264
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We develop an Explicitly Solvable Energy-Conserving (ESEC) algorithm for the Stochastic Differential Equation (SDE) describing the pitch-angle scattering process in magnetized plasmas. The Cayley transform is used to calculate both the deterministic gyromotion and stochastic scattering, affording the algorithm to be explicitly solvable and exactly energy conserving. An unusual property of the SDE for pitch-angle scattering is that its coefficients diverge at the zero velocity and do not satisfy the global Lipschitz condition. Consequently, when standard numerical methods, such as the Euler-Maruyama (EM), are applied, numerical convergence is difficult to establish. For the proposed ESEC algorithm, its energy-preserving property enables us to overcome this obstacle. We rigorously prove that the ESEC algorithm is order 1/2 strongly convergent. This result is confirmed by detailed numerical studies. For the case of pitch-angle scattering in a magnetized plasma with a constant magnetic field, the numerical solution is benchmarked against the analytical solution, and excellent agreements are found.
قيم البحث
اقرأ أيضاً
Pitch-angle scattering rates for cosmic-ray particles in magnetohydrodynamic (MHD) simulations with imbalanced turbulence are calculated for fully evolving electromagnetic turbulence. We compare with theoretical predictions derived from the quasiline
ar theory of cosmic-ray diffusion for an idealized slab spectrum and demonstrate how cross helicity affects the shape of the pitch-angle diffusion coefficient. Additional simulations in evolving magnetic fields or static field configurations provide evidence that the scattering anisotropy in imbalanced turbulence is not primarily due to coherence with propagating Alfven waves, but an effect of the spatial structure of electric fields in cross-helical MHD turbulence.
Vlasov solvers that operate on a phase-space grid are highly accurate but also numerically demanding. Coarse velocity space resolutions, which are unproblematic in particle-in-cell (PIC) simulations, lead to strong numerical heating or oscillations i
n standard continuum Vlasov methods. We present a new dual Vlasov solver which is based on an established positivity preserving advection scheme for the update of the distribution function and an energy conserving partial differential equation solver for the kinetic update of mean velocity and temperature. The solvers work together via moment fitting during which the maximum entropy part of the distribution function is replaced by the solution from the partial differential equation solver. This numerical scheme makes continuum Vlasov methods competitive with PIC methods concerning computational cost and enables us to model large scale reconnection in Earths magnetosphere with a fully kinetic continuum method. The simulation results agree well with measurements by the MMS spacecraft.
A generalized Ohms law is derived to treat strongly magnetized plasmas in which the electron gyrofrequency significantly exceeds the electron plasma frequency. The frictional drag due to Coulomb collisions between electrons and ions is found to shift
, producing an additional transverse resistivity term in the generalized Ohms law that is perpendicular to both the current ($vc{J}$) and the Hall ($vc{J} times vc{B}$) direction. In the limit of very strong magnetization, the parallel resistivity is found to increase by a factor of 3/2, and the perpendicular resistivity to scale as $ln (omega_{ce} tau_e)$, where $omega_{ce} tau_e$ is the Hall parameter. Correspondingly, the parallel conductivity coefficient is reduced by a factor of 2/3, and the perpendicular conductivity scales as $ln(omega_{ce} tau_e)/(omega_{ce} tau_e)^2$. These results suggest that strong magnetization significantly changes the magnetohydrodynamic evolution of a plasma.
The effect of a magnetic field on the characteristics of capacitively coupled radio frequency discharges is investigated and found to be substantial. A one-dimensional particle-in-cell simulation shows that geometrically symmetric discharges can be a
symmetrized by applying a spatially inhomogeneous magnetic field. This effect is similar to the recently discovered electrical asymmetry effect. Both effects act independently, they can work in the same direction or compensate each other. Also the ion energy distribution functions at the electrodes are strongly affected by the magnetic field, although only indirectly. The field influences not the dynamics of the sheath itself but rather its operating conditions, i.e., the ion flux through it and voltage drop across it. To support this interpretation, the particle-in-cell results are compared with the outcome of the recently proposed ensemble-in-spacetime algorithm. Although that scheme resolves only the sheath and neglects magnetization, it is able to reproduce the ion energy distribution functions with very good accuracy, regardless of whether the discharge is magnetized or not.
Nonlinear axisymmetric cylindrical plasma oscillations in magnetized collisionless plasmas are a model for the electron fluid collapse on the axis behind an ultrashort relativisically intense laser pulse exciting a plasma wake wave. We present an ana
lytical description of the strongly nonlinear oscillations showing that the magnetic field prevents closing of the cavity formed behind the laser pulse. This effect is demonstrated with 3D PIC simulations of the laser-plasma interaction. An analysis of the betatron oscillations of fast electrons in the presence of the magnetic field reveals a characteristic Four-Ray Star pattern which has been observed in the image of the electron bunch in experiments [T. Hosokai, et al., Phys. Rev. Lett. 97, 075004 (2006)].
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
