Related papers: A fast, high-order solver for the Grad-Shafranov e…
We present an efficient and accurate algorithm for solving the Poisson equation in spherical polar coordinates with a logarithmic radial grid and open boundary conditions. The method employs a divide-and-conquer strategy, decomposing the…
Steady plasma flows have been studied almost exclusively in systems with continuous symmetry or in open domains. In the absence of continuous symmetry, the lack of a conserved quantity makes the study of flows intrinsically challenging. In…
The boundary problem about behaviour (oscillations) of the electronic plasmas with arbitrary degree of degeneration of electronic gas in half-space with diffusion boundary conditions is analytically solved. The kinetic equation of Vlasov -…
A Grad-Shafranov equation (GSE) valid for compact quasisymmetric stellarators is derived by an asymptotic expansion around a vacuum field carried to first order. We obtain an equation for the existence of flux surfaces leading up to the…
This paper presents a grid-free simulation algorithm for the fully three-dimensional Vlasov--Poisson system for collisionless electron plasmas. We employ a standard particle method for the numerical approximation of the distribution…
We consider the free boundary problem for the plasma-vacuum interface in ideal compressible magnetohydrodynamics (MHD). In the plasma region the flow is governed by the usual compressible MHD equations, while in the vacuum region we…
In three dimensional toroidal domains without symmetry, the standard magnetohydrodynamic (MHD) equilibrium model used for magnetic confinement fusion does not generally support smooth solutions. Instead, solutions have singular plasma…
For the fractional Laplacian of variable order, an efficient and accurate numerical evaluation in multi-dimension is a challenge for the nature of a singular integral. We propose a simple and easy-to-implement finite difference scheme for…
The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very…
We present the solution to an inverse problem arising in the context of finding a distribution function for a specific collisionless plasma equilibrium. The inverse problem involves the solution of two integral equations, each having the…
We present a fully pseudo-spectral scheme to solve axisymmetric hyperbolic equations of second order. With the Chebyshev polynomials as basis functions, the numerical grid is based on the Lobbato (for two spatial directions) and Radau (for…
Fast, high-order accurate algorithms for electromagnetic scattering from axisymmetric objects are of great importance when modeling physical phenomena in optics, materials science (e.g. meta-materials), and many other fields of applied…
We present a new tunably-accurate Laguerre Petrov-Galerkin spectral method for solving linear multi-term fractional initial value problems with derivative orders at most one and constant coefficients on the half line. Our method results in…
We consider the free boundary problem for relativistic plasma--vacuum interfaces in two and three spatial dimensions. The plasma flow is governed by the equations of ideal relativistic magnetohydrodynamics, while the vacuum magnetic and…
We develop an approach of Grad-Shafranov (GS) reconstruction for toroidal structures in space plasmas, based on in-situ spacecraft measurements. The underlying theory is the GS equation that describes two-dimensional magnetohydrostatic…
We are interested in developing a numerical method for capturing stationary sheaths, that a plasma forms in contact with a metallic wall. This work is based on a bi-species (ion/electron) Vlasov-Amp{\`e}re model proposed in [3]. The main…
In this paper, we consider the three dimensional Vlasov equation with an inhomogeneous, varying direction, strong magnetic field. Whenever the magnetic field has constant intensity, the oscillations generated by the stiff term are periodic.…
We perform the asymptotic analysis of parabolic equations with stiff transport terms. This kind of problem occurs, for example, in collisional gyrokinetic theory for tokamak plasmas, where the velocity diffusion of the collision mechanism…
Fast and high-order accurate algorithms for three dimensional elastic scattering are of great importance when modeling physical phenomena in mechanics, seismic imaging, and many other fields of applied science. In this paper, we develop a…
We consider the plasma confined in a general axisymmetric spatial domain with perfect conducting boundary which reflects particles specularly, and look at a certain class of equilibria, assuming axisymmetry in the problem. We prove a sharp…