In this paper, a discrete Boltzmann method (DBM) for plasma kinetics is proposed and further used to investigate the non-equilibrium characteristics in Orszag–Tang (OT) vortex and Richtmyer–Meshkov instability (RMI) problems. The construction of DBM mainly considers two aspects. The first is to build a physical model with sufficient capability to capture underlying physics. The second is to devise schemes for extracting more valuable information from massive data. For the first aspect, the generated model is equivalent to a magnetohydrodynamic model, and a coarse-grained model for extracting the most relevant thermodynamic non-equilibrium (TNE) behaviors including the entropy production rate. For the second aspect, the DBM uses non-conserved kinetic moments of (f−feq) to describe the non-equilibrium states and behaviors of complex systems. It is found that (i) for OT vortex, the entropy production rate and compression difficulty first increase and then decrease with time. (ii) For RMI with interface inversion and re-shock process, the influence of magnetic field on TNE effects shows stages: before the interface inversion, the TNE strength is enhanced by delaying the interface inversion; while after the interface inversion, the TNE strength is significantly reduced. Both the global average TNE strength and entropy production rate contributed by non-organized energy flux can be used as physical criteria to identify whether or not the magnetic field is sufficient to prevent the interface inversion. In general, this paper proposes a generalized physical modeling and analysis scheme that has the potential for investigating the kinetic physics in plasma.

The effects of nonextensive electrons on nonlinear ion acoustic waves in dusty negative ion plasmas with ion–dust collisions are investigated. Analytical results show that both solitary and shock waves are supported in this system. The wave propagation is governed by a Korteweg–de Vries Burgers-type equation. The coefficients of this equation are modified by the nonextensive parameter q. Numerical calculations indicate that the amplitude of solitary wave and oscillatory shock can be obviously modified by the nonextensive electrons, but the monotonic shock is little affected.

A mesoscopic lattice Boltzmann method based on the BGK model is proposed to solve a class of two-dimensional second-order nonlinear partial differential equations by incorporating an amending function. The model provides an efficient and stable framework for simulating initial value problems of second-order nonlinear partial differential equations and is adaptable to various nonlinear systems, including strongly nonlinear cases. The numerical characteristics and evolution patterns of these nonlinear equations are systematically investigated. A D2Q4 lattice model is employed, and the kinetic moment constraints for both local equilibrium and correction distribution functions are derived in the four velocity directions. Explicit analytical expressions for these distribution functions are presented. The model is verified to recover the target macroscopic equations in the continuous limit via Chapman–Enskog analysis. Numerical experiments using exact solutions are performed to assess the model’s accuracy and stability. The results show excellent agreement with exact solutions and demonstrate the model’s robustness in capturing nonlinear dynamics.