Search

In this study, we find the points of transition between elliptic and hyperbolic regimes for the axisymmetric extended magnetohydrodynamic (MHD) equilibrium equations. The ellipticity condition is expressed via a single inequality but is more involved than the corresponding two-fluid ones due to the imposition of the quasineutrality condition and is also more complicated than the Hall MHD ellipticity condition, due to electron inertia. In fact, the inclusion of electron inertia is responsible for peculiar results; namely, even the static equilibrium equations can become hyperbolic.

An analytic solution to a generalized Grad-Shafranov equation with flow of arbitrary direction is obtained upon adopting the generic linearizing ansatz for the free functions related to the poloidal current, the static pressure, and the electric field. Subsequently, a D-shaped tokamak pertinent equilibrium with sheared flow is constructed using the aforementioned solution.An analytic solution to a generalized Grad-Shafranov equation with flow of arbitrary direction is obtained upon adopting the generic linearizing ansatz for the free functions related to the poloidal current, the static pressure, and the electric field. Subsequently, a D-shaped tokamak pertinent equilibrium with sheared flow is constructed using the aforementioned solution.

We employ a conformal mapping transformation to solve a generalized Grad-Shafranov equation with incompressible plasma flow of arbitrary direction and construct particular up-down asymmetric D-shaped and diverted tokamak equilibria. The proposed method can also be employed as an alternative quasi-analytic method to solving two dimensional elliptic partial differential equations., Comment: 9 pages, 4 figures, accepted in Physics of Plasmas

An alternative method based on an inverse aspect ratio (ϵ) expansion which reduces the axisymmetric equilibrium problem to a set of ODEs containing terms of arbitrary order in ϵ is employed to solve a generalized Grad-Shafranov equation with incompressible sheared flow nonparallel to the magnetic field. The method is applied to construct equilibria with either circular magnetic surfaces and reversed magnetic shear or D-shaped magnetic surfaces and normal magnetic shear. From the former equilibrium, it turns out that the electric field results in an increase of the reversed magnetic shear, thus indicating potential synergetic effects of the sheared flow and the magnetic shear in the formation of an internal transport barrier in consistent with experimental evidence.

The equilibrium of an axisymmetric magnetically confined plasma with anisotropic resistivity and incompressible flows parallel to the magnetic field is investigated within the framework of the magnetohydrodynamic (MHD) theory by keeping the convective flow term in the momentum equation. It turns out that the stationary states are determined by a second-order elliptic partial differential equation for the poloidal magnetic flux function psi along with a decoupled Bernoulli equation for the pressure identical in form with the respective ideal MHD equations; equilibrium consistent expressions for the resistivities eta(parallel to) and eta(perpendicular to) parallel and perpendicular to the magnetic field are also derived from Ohm's and Faraday's laws. Unlike in the case of stationary states with isotropic resistivity and parallel flows [G.N. Throumoulopoulos and H. Tasso, J. Plasma Phys. 64, 601 (2000)] the equilibrium is compatible with nonvanishing poloidal current densities. Also, although exactly Spitzer resistivities either eta(parallel to)(psi) or eta(perpendicular to)(psi) are not allowed, exact solutions with vanishing poloidal electric fields can be constructed with 771, and y, profiles compatible with roughly collisional resistivity profiles, i.e., profiles having a minimum close to the magnetic axis, taking very large values on the boundary and such that eta(perpendicular to)>eta(parallel to). For equilibria with vanishing flows satisfying the relation (dP/dpsi)(dI(2)/dpsi)>0, where P and I are the pressure and the poloidal current functions, the difference eta(perpendicular to)-eta(parallel to) for the reversed-field pinch scaling, B(p)approximate toB(t), is nearly two times larger than that for the tokamak scaling, B(p)approximate to0.1B(t) (B-p and B-t are the poloidal and toroidal magnetic-field components). The particular resistive equilibrium solutions obtained in the present work, inherently free of-but not inconsistent with-Pfirsch-Schluter diffusion, indicate that parallel flows might result in a reduction of the diffusion observed in magnetically confined plasmas. (C) 2003 American Institute of Physics. Physics of Plasmas

In a recent publication [H. Tasso and G. N. Throumoulopoulos, Phys. Lett. A 271, 413 (2000)] on Lyapunov stability of general mechanical systems, it is shown that “parametric excitations” can be stabilized by dissipation for positive potential energies. Specializing on the damped Mathieu equation permits one to establish its full stability chart. It is then seen that dissipation broadens the regions of stability to the extent that not only the response to parametric excitations is damped, but even “negative-energy” modes are stabilized by the combined action of the parametric excitation and the damping coefficient. The extension of this analysis to the “two-step” Hill’s equation shows that the stability regions become many times larger than those of the Mathieu equation. By analogy, these findings are a strong indication that the “resistive wall mode” could be stabilized by the joint action of a properly tailored time-dependent wall resistivity and a sufficient viscous dissipation in the plasma. Note that within this scheme neither the wall nor the plasma need to be in motion. An extension of this work to include more realistic models is in progress.

Analytic solutions of the magnetohydrodynamic equilibrium equations for a symmetric magnetically confined plasma with sheared incompressible flows associated with electric fields similar to those observed in the transition from the low- to the high-confinement mode in tokamaks are constructed in cylindrical and toroidal geometries. In particular, an exact toroidal solution is obtained which for vanishing flows reduces to the Solovev equilibrium which has been extensively employed in tokamak confinement studies. Owing to the flow, several toroidal configurations having either one or two stagnation points are possible in addition to the usual ones with a single magnetic axis. For flows pertaining to tokamak operational regime the extremum of the electric field becomes larger as flow and its shear increase, the location of the extremum being, however, nearly independent of these variations. In addition, the flow affects the safety factor profile and the shape of the magnetic surfaces and results in an increase of the magnetic shear and a decrease of the toroidal beta. The impact of plasma elongation on the above-mentioned confinement figures of merit is also evaluated. (C) 2001 American Institute of Physics. Physics of Plasmas

A Solovev-like solution describing equilibria with field aligned incompressible flows [G. N. Throumoulopoulos and H. Tasso, Phys. Plasmas 19, 014504 (2012)] is extended to non parallel flows. The solution expressed as a superposition of Bessel functions contains an arbitrary number of free parameters which are exploited to construct a variety of configurations including ITER shaped ones. For parallel flows, application of a sufficient condition for linear stability shows that this condition is satisfied in an appreciable part of the plasma region on the high-field side mostly due to the variation of the magnetic field perpendicular to the magnetic surfaces. Also, the results indicate that depending on the shape of the Mach-function profile and the values of the free parameters the flow and flow shear may have either stabilizing or destabilizing effects.

Generic linear axisymmetric equilibria with plasma flow nonparallel to the magnetic field are obtained on the basis of a generalized Grad-Shafranov equation by employing an ansatz reducing the problem to a set of ordinary differential equations which can be solved recursively. In particular, an ITER like equilibrium with reversed magnetic shear and peaked current density is constructed and its characteristics are studied in connection with the flow. Also for parallel flows, the linear stability is examined by means of a sufficient condition. The results indicate that the flow may have a stabilizing effect.

Starting from a simple, marginally stable model considered for Lyapunov based boundary control of flexible mechanical systems, we add a term driving an instability and prove that for an appropriate control condition the system can become Lyapunov stable. A similar approximate extension is found for the general energy principle of linearized magnetohydrodynamics. The implementation of such external instantaneous actions may, however, impose challenging constraints for fusion plasmas.

Several classes of analytic solutions to a generalized Grad-Shafranov equation with incompressible plasma flow non-parallel to the magnetic field are constructed. The solutions include higher transcendental functions such as the Meijer G-function and describe D-shaped and diverted configurations with either a single or double X-points. Their characteristics are examined in particular with respect to the flow parameters associated with the electric field. It turns out that the electric field makes the safety factor flatter and increases the magnitude and shear of the toroidal velocity in qualitative agreement with experimental evidence on the formation of internal transport barriers in tokamaks, thus indicating a potential stabilizing effect of the electric field.

A general stability condition for plasma-vacuum systems with resistive walls is derived by using the Frieman Rotenberg Lagrangian stability formulation [Rev. Mod. Phys. 32, 898 ( 1960)]. It is shown that the Lyapunov stability limit for external modes does not depend upon the gyroscopic term but upon the sign of the perturbed potential energy only. In the absence of dissipation in the plasma such as viscosity, it is expected that the flow cannot stabilize the system. (C) 2011 American Institute of Physics. [doi:10.1063/1.3606470] Physics of Plasmas

A kinetic consideration of an axisymmetric equilibrium with vanishing electric field near the magnetic axis shows that del f should not vanish on axis within the framework of Vlasov theory while it can either vanish or not in the framework of both a drift kinetic and a gyrokinetic theories (f is either the pertinent particle or the guiding center distribution function). This different behavior, relating to the reduction of phase space which leads to the loss of a Vlasov constant of motion, may result in the construction of different currents in the reduced phase space than the Vlasov ones. This conclusion is indicative of some limitation on the implications of reduced kinetic theories in particular as concerns the physics of energetic particles in the central region of magnetically confined plasmas., 9 pages

Extending previous work [R. L. White and R. D. Hazeltine, Phys. Plasmas 16, 123101 (2009)] to the case of a generalized Grad-Shafranov equation (GGSE) with incompressible flow of arbitrary direction, we obtain new classes of exact solutions on the basis of Lie-point symmetries. This is done by using a previously found exact generalized Solovev solution to the GGSE. The new solutions containing five free parameters describe D-shaped toroidal configurations with plasma flow non-parallel to the magnetic field. In addition, the full symmetry group is obtained and new group-invariant solutions to the GGSE are presented.

The Hamiltonian structure of ideal translationally symmetric extended MHD (XMHD) is obtained by employing a method of Hamiltonian reduction on the three-dimensional noncanonical Poisson bracket of XMHD. The existence of the continuous spatial translation symmetry allows the introduction of Clebsch-like forms for the magnetic and velocity fields. Upon employing the chain rule for functional derivatives, the 3D Poisson bracket is reduced to its symmetric counterpart. The sets of symmetric Hall, Inertial, and extended MHD Casimir invariants are identified, and used to obtain energy-Casimir variational principles for generalized XMHD equilibrium equations with arbitrary macroscopic flows. The obtained set of generalized equations is cast into Grad-Shafranov-Bernoulli (GSB) type, and special cases are investigated: static plasmas, equilibria with longitudinal flows only, and Hall MHD equilibria, where the electron inertia is neglected. The barotropic Hall MHD equilibrium equations are derived as a limiting case of the XMHD GSB system, and a numerically computed equilibrium configuration is presented that shows the separation of ion-flow from electromagnetic surfaces.

The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton's equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.