Stable finite element analysis of viscous dusty plasma
Abstract
Purpose
The purpose of this paper is to provide an analysis of dust acoustic (solitary) waves including viscosity. Specifically, the authors consider a dusty unmagnetized plasma system consisting of negatively charged dust and Boltzmann electrons and ions.
Design/methodology/approach
In this paper, a Petrov–Galerkin weak form with upwinding is adopted. Nonlinearity of ion and electron number density in terms of an electrostatic potential is included. A fully implicit time integration is used (backward Euler method), which requires the first derivative of the weak form. A three-field formulation is proposed, with the dust number density, the electrostatic potential and the dust velocity being the unknown fields.
Findings
In this study, two numerical examples are introduced and results show great promise for the proposed formulation as a predictive tool in viscous dusty plasmas. Presence of solitary waves was demonstrated. Dusty plasma vortices are predicted in 2D and 3D, as mentioned in the specialized literature.
Research limitations/implications
We observed some dependence on step size, which is due to the simple time-stepping scheme. This can be solved with a higher order integration scheme, which implies an added cost to the solution.
Practical implications
Dusty plasmas are found in astrophysics (Saturn rings) and electronics industry at several scales and have high impact as a contaminant.
Originality/value
To the authors’ knowledge, this is the first paper with a simulation of dusty plasma including vortices.
Keywords
Citation
Areias, P., Sikta, J.N. and Pereira dos Santos, M. (2018), "Stable finite element analysis of viscous dusty plasma", Engineering Computations, Vol. 35 No. 3, pp. 1230-1249. https://doi.org/10.1108/EC-06-2017-0191
Publisher
:Emerald Publishing Limited
Copyright © 2018, Emerald Publishing Limited