Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHybrid Vlasov-MHD models: Hamiltonian vs. non-Hamiltonian
137
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This paper investigates hybrid kinetic-MHD models, where a hot plasma (governed by a kinetic theory) interacts with a fluid bulk (governed by MHD). Different nonlinear coupling schemes are reviewed, including the pressure-coupling scheme (PCS) used in modern hybrid simulations. This latter scheme suffers from being non-Hamiltonian and to not exactly conserve total energy. Upon adopting the Vlasov description for the hot component, the non-Hamiltonian PCS and a Hamiltonian variant are compared. Special emphasis is given to the linear stability of Alfven waves, for which it is shown that a spurious instability appears at high frequency in the non-Hamiltonian version. This instability is removed in the Hamiltonian version.
rate research
Read More
Different variants of hybrid kinetic-fluid models are considered for describing the interaction of a bulk fluid plasma obeying MHD and an energetic component obeying a kinetic theory. Upon using the Vlasov kinetic theory for energetic particles, two planar Vlasov-MHD models are compared in terms of their stability properties. This is made possible by the Hamiltonian structures underlying the considered hybrid systems, whose infinite number of invariants makes the energy-Casimir method effective for determining stability. Equilibrium equations for the models are obtained from a variational principle and in particular a generalized hybrid Grad-Shafranov equation follows for one of the considered models. The stability conditions are then derived and discussed with particular emphasis on kinetic particle effects on classical MHD stability.
For the quadratic helicity $chi^{(2)}$ we present a generalization of the Arnold inequality which relates the magnetic energy to the quadratic helicity, which poses a lower bound. We then introduce the quadratic helicity density using the classical magnetic helicity density and its derivatives along magnetic field lines. For practical purposes we also compute the flow of the quadratic helicity and show that for an $alpha^2$-dynamo setting it coincides with the flow of the square of the classical helicity. We then show how the quadratic helicity can be extended to obtain an invariant even under compressible deformations. Finally, we conclude with the numerical computation of $chi^{(2)}$ which show cases the practical usage of this higher order topological invariant.
We present a new variational principle for the gyrokinetic system, similar to the Maxwell-Vlasov action presented in Ref. 1. The variational principle is in the Eulerian frame and based on constrained variations of the phase space fluid velocity and particle distribution function. Using a Legendre transform, we explicitly derive the field theoretic Hamiltonian structure of the system. This is carried out with a modified Dirac theory of constraints, which is used to construct meaningful brackets from those obtained directly from Euler-Poincar{e} theory. Possible applications of these formulations include continuum geometric integration techniques, large-eddy simulation models and Casimir type stability methods. [1] H. Cendra et. al., Journal of Mathematical Physics 39, 3138 (1998)
We introduce a new matter action principle, with a wide range of applicability, for the Vlasov equation in terms of a conjugate pair of functions. Here we apply this action principle to the study of matter in Bianchi cosmological models in general relativity. The Bianchi models are spatially-homogeneous solutions to the Einstein field equations, classified by the three-dimensional Lie algebra that describes the symmetry group of the model. The Einstein equations for these models reduce to a set of coupled ordinary differential equations. The class A Bianchi models admit a Hamiltonian formulation in which the components of the metric tensor and their time derivatives yield the canonical coordinates. The evolution of anisotropy in the vacuum Bianchi models is determined by a potential due to the curvature of the model, according to its symmetry. For illustrative purposes, we examine the evolution of anisotropy in models with Vlasov matter. The Vlasov content is further simplified by the assumption of cold, counter-streaming matter, a kind of matter that is far from thermal equilibrium and is not describable by an ordinary fluid model nor other more simplistic matter models. Qualitative differences and similarities are found in the dynamics of certain vacuum class A Bianchi models and Bianchi Type I models with cold, counter-streaming Vlasov matter potentials analogous to the curvature potentials of corresponding vacuum models.
A high-order method to evolve in time electromagnetic and velocity fields in conducting fluids with non-periodic boundaries is presented. The method has a small overhead compared with fast FFT-based pseudospectral methods in periodic domains. It uses the magnetic vector potential formulation for accurately enforcing the null divergence of the magnetic field, and allowing for different boundary conditions including perfectly conducting walls or vacuum surroundings, two cases relevant for many astrophysical, geophysical, and industrial flows. A spectral Fourier continuation method is used to accurately represent all fields and their spatial derivatives, allowing also for efficient solution of Poisson equations with different boundaries. A study of conducting flows at different Reynolds and Hartmann numbers, and with different boundary conditions, is presented to study convergence of the method and the accuracy of the solenoidal and boundary conditions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
