Abstract:We consider the numerical modeling of the Farley-Buneman instability development in the earth's ionosphere plasma. The ion behavior is governed by the kinetic Landau equation in the four-dimensional phase space, and since the finite difference discretization on a tensor product grid is used, this equation becomes the most computationally challenging part of the scheme. To relax the complexity and memory consumption, an adaptive model reduction using the low-rank separation of variables, namely the Tensor Train format, is employed. The approach was verified via the prototype MATLAB implementation. Numerical experiments demonstrate the possibility of efficient separation of space and velocity variables, resulting in the solution storage reduction by a factor of order tens.

What problem does this paper attempt to address?

This paper aims to solve the numerical simulation problem of Farley - Buneman instability in the Earth's ionospheric E - region. Specifically, the paper focuses on the low - frequency plasma oscillation phenomenon that occurs when the electric field is perpendicular to the geomagnetic field under the action of the geomagnetic field, where electrons are magnetized while ions are not. This instability mainly occurs in the E - regions at the equator and the poles, and is manifested as low - frequency plasma fluctuations with wavelengths in the meter - level range. ### Core problems of the paper 1. **Challenges in numerical simulation**: The paper points out that discretizing the Landau kinetic equation in the four - dimensional phase space on a tensor - product grid using the finite - difference method is the most computationally challenging part of the entire numerical scheme. This is mainly because the solution in the four - dimensional phase space requires a large amount of memory to store and the computational complexity at each time step is very high. 2. **Need for model simplification**: In order to reduce the computational complexity and memory consumption, the paper proposes a method based on variable low - rank separation, especially using the Tensor Train format to approximate the solution in the four - dimensional phase space. This method significantly reduces the storage requirements of the solution by effectively separating the space and velocity variables. ### Solutions - **Tensor Train format**: The paper adopts the Tensor Train format to perform low - rank approximation of the solution in the four - dimensional phase space. This method can compress a four - dimensional array that originally needs to store \(n^4\) elements to only need to store \(O(dnr^2)\) elements, where \(d\) is the number of dimensions, \(n\) is the number of grid points on each dimension, and \(r\) is the Tensor Train rank. - **Adaptive model order reduction**: Through the adaptive model order reduction technique, the paper effectively separates the space and velocity variables, thereby significantly reducing the storage requirements and computational time. Experimental results show that this separation method can reduce the storage requirements by dozens of times. ### Experimental verification - **MATLAB prototype implementation**: The paper verifies the effectiveness of this method through MATLAB prototype implementation. Numerical experiments show that using the Tensor Train format can efficiently separate the space and velocity variables, thereby significantly reducing the storage requirements of the solution. - **Performance improvement**: Compared with the traditional method of directly storing four - dimensional arrays, using the Tensor Train format not only reduces memory consumption but also improves computational efficiency. Experimental results show that the storage requirements are reduced by dozens of times and the computational time is shortened by at least half. ### Conclusion The method proposed in the paper provides an effective way to solve high - dimensional problems, especially when dealing with plasma kinetic problems in the four - dimensional phase space, which can significantly improve computational efficiency and reduce memory consumption. Future work will further optimize the program code to make it suitable for parallel computing environments and explore higher - order discretization schemes to further improve performance.

Low-rank approximation in the numerical modeling of the Farley-Buneman instability in ionospheric plasma

A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations

High Fidelity simulations of the multi-species Vlasov equation in the electro-static, collisional-less limit

Numerical solution of the quantum Lenard-Balescu equation for a one-component plasma

Hermite Spectral Method for Fokker-Planck-Landau Equation Modeling Collisional Plasma

A low-rank method for two-dimensional time-dependent radiation transport calculations

Accelerating the simulation of kinetic shear Alfvén waves with a dynamical low-rank approximation

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

Higher-order moment convergent method in weakly anisotropic plasma and the NLVFP code for solution of the 0D-2V Vlasov-Fokker-Planck equation

Accelerating Kinetic Simulations of Electrostatic Plasmas with Reduced-Order Modeling

Interpolatory dynamical low-rank approximation for the 3+3d Boltzmann-BGK equation

Simulating the nonlinear interaction of relativistic electrons and tokamak plasma instabilities: Implementation and validation of a fluid model

High-order implicit solver in conservative formulation for tokamak plasma transport equations

Novel approaches for the reliable and efficient numerical evaluation of the Landau operator

Higher Dimensional Semi-Relativistic Time-Fractional Vlasov-Maxwell Code for Numerical Simulation Based on Linear Polarization and 2D Landau Damping Instability.

Plasma kinetics: Discrete Boltzmann modeling and Richtmyer–Meshkov instability

Plasma kinetics: Discrete Boltzmann modelling and Richtmyer-Meshkov instability

Low rank approximation for the numerical simulation of high dimensional Lindblad and Riccati equations

Hilbert expansion based fluid models for kinetic equations describing neutral particles in the plasma edge of a fusion device

THE RENORMALIZED STRONG TURBULENCE THEORY FOR THE LOW-FREQUENCY MAGNETIC FIELD AND THE ION ACOUSTIC WAVE EXCITED BY HIGH-FREQUENCY WAVE