Numerical boundary control of multi-dimensional discrete-velocity kinetic models

This paper extends our recent results on multi-dimensional discrete-velocity models to the numerical level. By adopting an operator splitting scheme and introducing a suitable discrete Lyapunov function, we derive numerical control laws that ensure the corresponding numerical solutions decay exponentially in time. To handle the stiff source term, we also use an implicit scheme for the collision part and prove the stability of the resulting schemes. The theoretical results are validated through three numerical simulations for the two-dimensional coplanar model.