PPPL-4161 is available in pdf format (664 KB).

Finite Element Implementation of Braginskii's Gyroviscous Stress with Application to the Gravitational Instability

Authors: N.M. Ferraro and S.C. Jardin

Date of PPPL Report: May 2006

Submitted to: Physics of Plasmas

We present a general, coordinate-independent expression for Braginskii's form of the ion gyroviscosity in the two-dimensional potential field representation, and implement this force in a full two-dimensional, two-fluid extended magnetohydrodynamic (MHD) numerical model. Our expression for the gyroviscous force requires no field to be differentiated more than twice, and thus is appropriate for finite elements with first derivatives continuous across element boundaries (C1 finite elements). We derive from our model, including the full gyroviscous stress, linear dispersion relations of a homogeneous equilibrium and of an inverted-density profile in the presence of gravity. Our treatment of the gravitational instability substantially extends previous work on the subject. Linear and nonlinear simulations of the gravitational instability are presented. The gyroviscous cancellation effect is demonstrated, and some limitations of the omega-star approximation are discussed.