A Gas-Kinetic Scheme for Maxwell Equations

Zhigang Pu,Kun Xu

2024-12-22

Abstract:In this paper, we present a gas-kinetic scheme using discrete velocity space to solve Maxwell equations. The kinetic model recovers Maxwell equations in the zero relaxation time limit. The scheme achieves second-order spatial and temporal accuracy in structured meshes comparable to the finite-difference time-domain (FDTD) method, without requiring staggered grids or leapfrog discretization. Our kinetic scheme is inherently multidimensional due to its use of kinetic beams in multiple directions, allowing larger time steps in multidimensional computations. It demonstrates better stability than FDTD when handling discontinuities and readily extends to unstructured meshes. We validate the method through various test cases including antenna simulation, sphere scattering, and flight vehicle scattering. The results align well with Riemann-solver-based solutions. Finally, we examine charge conservation for Maxwell equations through test cases.

Numerical Analysis,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a new numerical method for solving Maxwell equations. Specifically, the author proposes a gas - kinetic scheme, using discrete velocity space to solve Maxwell equations. This method aims to overcome some limitations of existing numerical methods (such as the finite - difference time - domain method FDTD) and provide better stability and higher computational efficiency.

### Main problems and goals

1. **Improve the stability and accuracy of numerical methods**:
- The existing FDTD method faces challenges when dealing with discontinuities and complex geometries, especially when applied on unstructured grids.
- The method proposed by the author shows better stability in multi - dimensional calculations by using the gas - kinetic model and can handle larger time steps.

2. **Avoid staggered grids and leap - frog discretization**:
- The FDTD method requires staggered grids and leap - frog discretization to satisfy the elliptic divergence constraint, which increases the complexity of implementation.
- The newly proposed gas - kinetic scheme can achieve second - order space - time accuracy on structured grids without these complex discretization techniques.

3. **Expand to unstructured grids**:
- Although the FEM method is more flexible in dealing with complex geometries, it needs to solve sparse linear systems at each time step.
- The gas - kinetic scheme can be naturally extended to unstructured grids, thus better adapting to electromagnetic simulations of complex geometries.

4. **Verification and comparison**:
- The author verifies the effectiveness of this method through multiple test cases (such as antenna simulation, sphere scattering and aircraft scattering).
- The results show that this method is consistent with the results of the finite - volume method (FVM) based on the Riemann solver and shows better stability when dealing with discontinuities.

### Key formulas

The perfect hyperbolic - type equations (PHM equations) of Maxwell equations are expressed as:

\[
\begin{aligned}
& \frac{\partial E}{\partial t}-c^{2} \nabla_{x} \times B+\chi c^{2} \nabla_{x} \phi =-\frac{1}{\epsilon_{0}} J, \\
& \frac{\partial B}{\partial t}+\nabla_{x} \times E+\gamma \nabla_{x} \psi = 0, \\
& \frac{1}{\chi} \frac{\partial \phi}{\partial t}+\nabla_{x} \cdot E=\frac{q}{\epsilon_{0}}, \\
& \frac{\epsilon_{0} \mu_{0}}{\gamma} \frac{\partial \psi}{\partial t}+\nabla_{x} \cdot B = 0,
\end{aligned}
\]

where $\phi$ and $\psi$ are artificial correction potentials for dealing with divergence errors.

The discrete BGK equation of the gas - kinetic model is:

\[
\frac{\partial f_{k}}{\partial t}+v_{k} \cdot \nabla f_{k}=\frac{g_{k}-f_{k}}{\tau},
\]

where $f_{k}(x, v_{k}, t)$ is the distribution function in the discrete velocity space, $v_{k}$ is the $k$ - th discrete microscopic velocity vector, and $\tau$ is the relaxation parameter.

### Summary

The main contribution of this paper is to propose a new gas - kinetic numerical method for solving Maxwell equations. This method not only achieves second - order space - time accuracy on structured grids, but also shows significant advantages in dealing with discontinuities and expanding to unstructured grids. Through the verification of multiple practical cases, the effectiveness and superiority of this method are proved.