Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems
68
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Recently, a family of models that couple multifluid systems to the full Maxwell equations draw a lot of attention in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they naturally retain non-ideal effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, etc. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid kinetic processes like plasma oscillation and electron cyclotron motion, even when they are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic forms solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with $S$ species, the locally implict algorithm involves inverting a local $3S+3$ matrix only, thus is very efficient. The performance can be further elevated by using the direct update formulas to skip null calculations. Benchmarks illustrated the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and large-scale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.
rate research
Read More
Wereportonanewmultiscalemethodapproachforthestudyofsystemswith wide separation of short-range forces acting on short time scales and long-range forces acting on much slower scales. We consider the case of the Poisson-Boltzmann equation that describes the long-range forces using the Boltzmann formula (i.e. we assume the medium to be in quasi local thermal equilibrium). We developed a new approach where fields and particle information (mediated by the equations for their moments) are solved self-consistently. The new approach is implicit and numerically stable, providing exact energy conservation. We tested different implementations all leading to exact energy conservation. The new method requires the solution of a large set of non-linear equations. We considered three solution strategies: Jacobian Free Newton Krylov, an alternative, called field hiding, based on hiding part of the residual calculation and replacing them with direct solutions and a Direct Newton Schwarz solver that considers simplified single particle-based Jacobian. The field hiding strategy proves to be the most efficient approach.
We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the non linear regime. We also investigate extensions of the Schrodinger-Poisson system that employ multiple stacked cold streams, and the von Neumann-Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. We present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.
Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsim-CYL, to be used for axially symmetric problems. The main motivation for the development of the cylindrical version is to greatly improve the computational speed by utilizing cylindrical symmetry. The ECsim-CYL discretizes the field equations in two-dimensional cylindrical coordinates using the finite volume method . For the particle mover, it uses a modification of ECsims mover for cylindrical coordinates by keeping track of all three components of velocity vectors, while only keeping radial and axial coordinates of particle positions. In this paper, we describe the details of the algorithm used in the ECsim-CYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasmas inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. The ECsim-CYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.
In this work, we use the finite differences in time domain (FDTD) numerical method to compute and assess the validity of Hopf solutions, or hopfions, for the electromagnetic field equations. In these solutions, field lines form closed loops characterized by different knot topologies which are preserved during their time evolution. Hopfions have been studied extensively in the past from an analytical perspective but never, to the best of our knowledge, from a numerical approach. The implementation and validation of this technique eases the study of more complex cases of this phenomena; e.g. how these fields could interact with materials (e.g. anisotropic or non-linear), their coupling with other physical systems (e.g. plasmas), and also opens the path on their artificial generation by different means (e.g. antenna arrays or lasers).
We consider linear stability of steady states of 1(1/2) and 3D Vlasov-Maxwell systems for collisionless plasmas. The linearized systems can be written as separable Hamiltonian systems with constraints. By using a general theory for separable Hamiltonian systems, we recover the sharp linear stability criteria obtained previously by different approaches. Moreover, we obtain the exponential trichotomy estimates for the linearized Vlasov-Maxwell systems in both relativistic and nonrelativistic cases.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
