A numerical investigation of electric current conservation associated with solar wind plasma under GES-model approach

A plasma-based Gravito-Electrostatic Sheath (GES) model for theoretical description of the basic physics of the surface emission mechanism of the quasi-neutral Solar Interior Plasma (SIP) on a bounded scale, and subsequently, its transonic transition into Solar Wind Plasma (SWP) on an unbounded scale with a simplified field-free fluid model approach has already been proposed. An autonomous closed set of self-consistently coupled nonlinear eigenvalue equations for analyzing the dynamical stability of the steady GES model on both the scales is developed. The focussed aim of the present contribution is only to study the electric current profiles associated with the SIP as well as SWP, and hence to investigate conservative dynamical features thereof. Applying the developed set of coupled nonlinear dynamical evolution equations, the profiles of electric currents associated with both SIP as well as SWP have been obtained numerically under a wide range spectrum of initial values of relevant physical parameters (by nonlinear stability analysis). Interestingly, it is conjectured that the dynamical evolution of solar electric currents is found fairly to be conservative (divergence-free) in nature except some deviations near the defined solar surface boundary. Results, discussions and conclusions on the basis of the obtained numerical results are presented in brief.