High-resolution numerical simulation of 2D nonlinear wave structures in electromagnetic fluids with absorbing boundary conditions

Abstract

Here we show how the full set of governing equations for the dynamics of charged-particle fluids in an electromagnetic field may be solved numerically in order to model nonlinear wave structures propagating in two dimensions. We employ a source-term adaptation and two-fluid extension of the second-order high-resolution central scheme of Balbas et al. (2004) [1]. The model employed is a 2D extension of that used by Baboolal and Bharuthram (2007) [5] in studies of 1D shocks and solitons in a two-fluid plasma under 3D electromagnetic fields. Further, we outline the use of free-flow boundary conditions to obtain stable wave structures over sufficiently long modelling times. As illustrative results, we examine the formation and evolution of shock-like and soliton structures of the magnetosonic mode.