An official website of the United States government
We consider a reduced dynamics for the first four fluid moments of the one-dimensional Vlasov–Poisson equation, namely, fluid density, fluid velocity, pressure, and heat flux. This dynamics depends on an equation of state to close the system. This equation of state (closure) connects the fifth-order moment—related to the kurtosis in velocity of the Vlasov distribution—with the first four moments. By solving the Jacobi identity, we derive an equation of state, which ensures that the resulting reduced fluid model is Hamiltonian. We show that this Hamiltonian closure allows symmetric homogeneous equilibria of the reduced fluid model to be stable.
more »
« less
This content will become publicly available on February 28, 2027
A Novel Splitting Method for Vlasov-Ampère
Abstract Vlasov equations model the dynamics of plasma in the collisionless regime. A standard approach for numerically solving the Vlasov equation is to operator split the spatial and velocity derivative terms, allowing simpler time-stepping schemes to be applied to each piece separately (known as the Cheng-Knorr method). One disadvantage of such an operator split method is that the order of accuracy of fluid moments (e.g., mass, momentum, and energy) is restricted by the order of the operator splitting (second-order accuracy in the Cheng-Knorr case). In this work, we develop a novel approach that first represents the particle density function on a velocity mesh with a local fluid approximation in each discrete velocity band and then introduces an operator splitting that splits the inter-velocity band coupling terms from the dynamics within the discrete velocity band. The advantage is that the inter-velocity band coupling terms are only needed to achieve consistency of the full distribution functions, but the local fluid models within each band are sufficient to achieve high-order accuracy on global moments such as mass, momentum, and energy. The resulting scheme is verified on several standard Vlasov-Poisson test cases.
more »
« less
- Award ID(s):
- 2410538
- PAR ID:
- 10711906
- Editor(s):
- Grabe, M; Oblapenko, G; Torrilhon, M
- Publisher / Repository:
- Springer Nature Switzerland
- Date Published:
- ISBN:
- 978-3-032-00094-1
- Page Range / eLocation ID:
- 451 to 459
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Energy bands in antiferromagnets are supposed to be spin degenerate in the absence of spin–orbit coupling (SOC). Recent studies have identified formal symmetry conditions for antiferromagnetic crystals in which this degeneracy can be lifted, spin splitting,even in the vanishing SOC (i.e., non‐relativistic) limit. Materials having such symmetries could enable spin‐split antiferromagnetic spintronics without the burden of using heavy‐atom compounds. However, the symmetry conditions that involve spin and magnetic symmetry are not always effective as practical material selection filters. Furthermore, these symmetry conditions do not readily disclose trends in the magnitude and momentum dependence of the spin‐splitting energy. Here, it is shown that the formal symmetry conditions enabling spin‐split antiferromagnets can be interpreted in terms of local motif pairs, such as octahedra or tetrahedra, each carrying opposite magnetic moments. Collinear antiferromagnets with such a spin‐structure motif pair, whose components interconvert by neither translation nor spatial inversion, will show spin splitting. Such a real‐space motif‐based approach enables an easy way to identify and design materials (illustrated in real example materials) having spin splitting without the need for SOC, and offers insights into the momentum dependence and magnitude of the spin splitting.more » « less
-
Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics (arXiv:2201.10397). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.more » « less
-
The Dougherty model Fokker–Planck operator is extended to describe nonlinear full- $$f$$ ( f is the distribution function) collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the $$H$$ -theorem ( H is related to the entropy) unlike an analogous Bhatnagar–Gross–Krook operator. Formulas for both a Cartesian velocity space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.more » « less
-
We present a robust and efficient method for simulating Lagrangian solid-fluid coupling based on a new operator splitting strategy. We use variational formulations to approximate fluid properties and solid-fluid interactions, and introduce a unified two-way coupling formulation for SPH fluids and FEM solids using interior point barrier-based frictional contact. We split the resulting optimization problem into a fluid phase and a solid-coupling phase using a novel time-splitting approach with augmentedcontact proxies, and propose efficient custom linear solvers. Our technique accounts for fluids interaction with nonlinear hyperelastic objects of different geometries and codimensions, while maintaining an algorithmically guaranteed non-penetrating criterion. Comprehensive benchmarks and experiments demonstrate the efficacy of our method.more » « less
