Related papers: A fast, high-order solver for the Grad-Shafranov e…
We present a simple and effective multigrid-based Poisson solver of second-order accuracy in both gravitational potential and forces in terms of the one, two and infinity norms. The method is especially suitable for numerical simulations…
We give an approximate algorithm of computing holonomic systems of linear differential equations for definite integrals with parameters. We show that this algorithm gives a correct answer in finite steps, but we have no general stopping…
We study smooth, spherically-symmetric solutions to the Vlasov-Poisson system and relativistic Vlasov-Poisson system in the plasma physical case. In particular, we construct solutions that initially possess arbitrarily small charge…
We show that nontrivial solutions to higher and fractional order equations with certain nonlinearity are radially symmetric and nonincreasing on geodesic balls in the hyperbolic space $\mathbb{H}^n$ as well as on the entire space…
In this study, we consider a numerical implementation of the nonlinear Rosenbluth-Trubnikov collision operator for particle simulations in plasma physics in the framework of the finite element method (FEM). The relevant particle evolution…
We analyse and compare several algorithms to compute numerically periodic solutions of high-dimensional dynamical systems and investigate their Floquet stability without building the monodromy matrix. The solution and its perturbation are…
This report addresses the boundary value problem for a second-order linear singularly perturbed FIDE. Traditional methods for solving these equations often face stability issues when dealing with small perturbation parameters. We propose an…
In this paper, a centred universal high-order finite volume method for solving hyperbolic balance laws is presented. The scheme belongs to the family of ADER methods where the Generalized Riemann Problems (GRP) is a building block. The…
An exact analytic solution is obtained for a uniformly expanding, neutral, highly conducting plasma sphere in an ambient dipole magnetic field with an arbitrary orientation of the dipole moment in the space. Based on this solution the…
We introduce a numerical method for the solution of the time-dependent Schrodinger equation with a smooth potential, based on its reformulation as a Volterra integral equation. We present versions of the method both for periodic boundary…
This paper presents a fast algorithm to solve a spectral estimation problem for two-dimensional random fields. The latter is formulated as a convex optimization problem with the Itakura-Saito pseudodistance as the objective function subject…
We introduce a high-order spline geometric approach for the initial boundary value problem for Maxwell's equations. The method is geometric in the sense that it discretizes in structure preserving fashion the two de Rham sequences of…
We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems…
Starting from the given passive particle equilibrium particle cylindrical profiles, we built self-consistent stationary conditions of the Maxwell-Vlasov equation at thermodynamic equilibrium with non-flat density profiles. The solutions to…
We present an explicit solver of the three-dimensional screened and unscreened Poisson's equation which combines accuracy, computational efficiency and versatility. The solver, based on a mixed plane-wave / interpolating scaling function…
A new highly accurate algorithm for the solution of the Falkner-Skan equation of boundary layer theory is presented. The algorithm, based on a Maclaurin series representation, finds its coefficients from recurrence. In addition,…
Outside the core of the plasma, the plasma current and pressure rapidly transition to zero in a scrape-off or edge region or plasma-vacuum interface. However, existing tools for fixed-boundary magnetohydrodynamic equilibria in 2D and 3D…
Plasma rotation in connection to both zonal and mean (equilibrium) flows can play a role in the transitions to the advanced confinement regimes in tokamaks, as the L-H transition and the formation of Internal Transport Barriers. For…
A theoretical model of the high-beta equilibrium of magnetospheric plasma was constructed by consistently connecting the (anisotropic pressure) Grad-Shafranov equation and the Vlasov equation. The Grad-Shafranov equation was used to…
We develop a fast divided-and-conquer indirect collocation method for the homogeneous Dirichlet boundary value problem of variable-order space-fractional diffusion equations. Due to the impact of the space-dependent variable order, the…